The Role of Chemokines in Thyroid Carcinoma

Transcriptomic Profiling in Low-Risk Thyroid Cancer Induced by Microwave Ablation

Background

Peripheral blood mononuclear cells (PBMCs) serve as the immune system's primary transportation hub outside of the affected ablated tissue. This study aims to explore the transcriptomic profiling of the immune response in PBMCs induced by microwave ablation (MWA) in low-risk thyroid cancer.

Methods

For eight patients diagnosed with low-risk thyroid cancer, 10 ml of peripheral venous blood was collected before MWA as well as one day and one month after MWA. mRNA was extracted from PBMCs for transcriptome next-generation gene sequencing and qRT-PCR analyses. The plasma samples were used for chemokine detection purposes.

Results

One day and one month after MWA, there were significant changes in GSEA, particularly in the NF-kappa B-TNFα pathway, inflammatory response, and early and late estrogen response. Common changes in differently expressed genes resulted in a significant downregulation of tumor-promoting genes (BCL3, NR6A1, and PFKFB3). One day after low-risk thyroid cancer MWA, GO enrichment analysis mainly revealed processes related to oxygen transport and other pathways. One month after MWA, GO enrichment analysis mainly revealed regulation of toll-like receptor signaling and other pathways. Furthermore, inflammation-related cytokines and regulatory genes, as well as tumor-promoting cytokines and regulatory genes, were downregulated after MWA.

Conclusions

This study presents a comprehensive profile of the systemic immune response induced by thermal ablation for treating low-risk thyroid cancer. More significantly, this study provides valuable insight into potential references for systemic antitumor immunity of ablation against low-risk thyroid cancer. This trial is registered with ChiCTR1900024544.

Prognostic Roles of Inflammatory Biomarkers in Radioiodine-Refractory Thyroid Cancer Treated with Lenvatinib

BACKGRUOUND

Inflammatory biomarkers, such as the neutrophil-to-lymphocyte ratio (NLR), lymphocyte-to-monocyte ratio (LMR), and platelet-to-lymphocyte ratio (PLR), serve as valuable prognostic indicators in various cancers. This multicenter, retrospective cohort study assessed the treatment outcomes of lenvatinib in 71 patients with radioactive iodine (RAI)-refractory thyroid cancer, considering the baseline inflammatory biomarkers.

METHODS

This study retrospectively included patients from five tertiary hospitals in Korea whose complete blood counts were available before lenvatinib treatment. Progression-free survival (PFS) and overall survival (OS) were evaluated based on the median value of inflammatory biomarkers.

RESULTS

No significant differences in baseline characteristics were observed among patients grouped according to the inflammatory biomarkers, except for older patients with a higher-than-median NLR (≥2) compared to their counterparts with a lower NLR (P= 0.01). Patients with a higher-than-median NLR had significantly shorter PFS (P=0.02) and OS (P=0.017) than those with a lower NLR. In multivariate analysis, a higher-than-median NLR was significantly associated with poor OS (hazard ratio, 3.0; 95% confidence interval, 1.24 to 7.29; P=0.015). However, neither the LMR nor the PLR was associated with PFS. A higher-than-median LMR (≥3.9) was significantly associated with prolonged OS compared to a lower LMR (P=0.036). In contrast, a higher-than-median PLR (≥142.1) was associated with shorter OS compared to a lower PLR (P=0.039).

CONCLUSION

Baseline inflammatory biomarkers can serve as predictive indicators of PFS and OS in patients with RAI-refractory thyroid cancer treated with lenvatinib.

Increased expression of GPX4 promotes the tumorigenesis of thyroid cancer by inhibiting ferroptosis and predicts poor clinical outcomes

BACKGROUND

Ferroptosis plays a critical role in suppressing cancer progression, and its essential regulator is glutathione peroxidase 4 (GPX4). High GPX4 expression can inhibit accumulation of iron, thus suppressing ferroptosis. However, its function in thyroid cancer has not been fully illuminated. Here, we explore the effect of GPX4 on thyroid cancer tumorigenesis and prognosis.

METHODS

Based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, GPX4 expression was investigated in cancer tissues and adjacent tissues. We determined the biological functions of GPX4-associated differentially expressed genes (DEGs) by using the "clusterProfiler" R package. In addition, the predictive value of GPX4 in thyroid cancer was assessed by using Cox regression analysis and nomograms. Finally, we conducted several in vitro experiments to determine the influence of GPX4 expression on proliferation and ferroptosis in thyroid cancer cells.

RESULTS

GPX4 expression was obviously elevated in thyroid cancer tissues compared with normal tissues. Biological function analysis indicated enrichment in muscle contraction, contractile fiber, metal ion transmembrane transporter activity, and complement and coagulation cascades. GPX4 overexpression was associated with stage T3-T4 and pathologic stage III-IV in thyroid cancer patients. Cox regression analysis indicated that GPX4 may be a risk factor for the overall survival of thyroid cancer patients. In vitro research showed that knockdown of GPX4 suppressed proliferation and induced ferroptosis in thyroid cancer cells.

CONCLUSIONS

GPX4 overexpression in thyroid cancer might play an essential role in tumorigenesis and may have prognostic value for thyroid cancer patients.

Preoperative serum inflammation-based scores in medullary thyroid cancer

INTRODUCTION

Neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio (PLR) and systemic immune-inflammation index (SII) are prognostic factors in several tumours, though little is known in medullary thyroid cancer (MTC).

OBJECTIVE

To evaluate the association between preoperative NLR, PLR and SII with MTC clinicopathological and molecular features, and their predictive value for lymph node and distant metastasis.

METHODS

We retrospectively analysed 75 patients with MTC who underwent surgery at our institution. The familial form of MTC was found in 12% of patients.

RESULTS

In our cohort, 56% were females, the median age at diagnosis was 57 years (44-69), the median tumour diameter was 25mm (15-50); 21.3% were multifocal and 34.7% had extrathyroidal extension. Lymph node and distant metastasis were observed in 36 (48.0%) and 8 (10.7%) patients, respectively. Higher NLR was associated with preoperative calcitonin, angioinvasion, extrathyroidal extension, moderate/severe fibrosis; higher PLR was associated with extrathyroidal extension and advanced T stages; lower SII and NLR were associated with biochemical cure after surgery. Increased PLR, NLR and SII were associated with advanced MTC stages. In the univariate analysis, only NLR was associated with lymph node metastasis (odds ratio (OR)=2.69, 95% confidence interval (CI): 1.50-5.84; p=0.004); however, in the multivariate model, NLR was no longer a predictive factor for lymph node metastasis. None of these serum inflammatory markers predicted the occurrence of distant metastasis.

CONCLUSION

In conclusion, NLR, PLR and SII are associated with aggressive MTC, but do not predict lymph node or distant metastasis.

CCL22 and CCL26 are potential biomarkers for predicting distant metastasis in thyroid carcinoma

OBJECTIVES

Chemokines have been suggested to play significant roles in the progression of malignant cancers. This study aimed to identify the chemokines related to malignant progression in thyroid carcinoma.

METHODS

The mRNA expression levels of 52 chemokines were compared between differentiated thyroid cancer (DTC) samples and normal thyroid tissues from The Cancer Genome Atlas database; survival analysis was then performed on the basis of differentially expressed chemokines. A retrospective study was conducted on the level of differentially expressed chemokines in 76 DTC patients. Functional pathway analysis was performed to explore chemokine-related regulatory mechanisms.

RESULTS

We identified 20 chemokines with differentially expressed mRNA levels through publicly available data. High levels of CCL22 and CCL26 were found to be related with metastasis in clinical DTC samples. High levels of CCL22 were found to be significantly related to poor prognosis in DTC patients. Pathway analyses revealed that cytokines might affect cancer progression through cytokine-cytokine receptor and cytokine-interleukin interactions.

CONCLUSIONS

CCL22 and CCL26 could serve as prognostic biomarkers in thyroid carcinoma.

Comprehensive analysis of lncRNA-mediated ceRNA regulatory networks and key genes associated with papillary thyroid cancer coexistent with Hashimoto's thyroiditis

OBJECTIVE

The incidence of papillary thyroid cancer (PTC) concomitant with Hashimoto's thyroiditis (HT) is gradually increasing over the past decades. This study aims to identify differentially expressed lncRNAs between tumor tissues of PTC with or without HT and further to confer a better understanding of lncRNA-based competing endogenous RNA (ceRNA) network in PTC with HT.

METHODS

GSE138198 containing tissue mRNA data and GSE192560 containing lncRNA data were utilized to perform differentially expression analysis. The ceRNA network was constructed based on miRNA-mRNA interactions merging with lncRNA-microRNA interactions. Functional enrichment analysis and protein-protein interaction (PPI) analysis were performed. The mRNA levels of core genes in the PPI analysis in tumor tissues collected from 112 PTC patients including 35 cases coexistent with HT were determined by quantitative real-time polymerase chain reaction (qRT-PCR).

RESULTS

A total of 57 genes and 40 lncRNAs, with value of |log2 fold change (FC)|≥ 1 and the adjusted P-value < 0.05, were deemed as differentially expressed genes and lncRNAs between PTC with and without HT. The pathways most significantly enriched by differentially expressed genes between PTC with and without HT were viral protein interaction with cytokine and cytokine receptor and cytokine-cytokine receptor interaction. CXCL10, CXCL9, CCL5, FCGR3A, and CCR2 owned degree values not less than 10 were deemed as core genes differentially expressed between PTC with and without HT. A total of 76 pairs of lncRNA-miRNA-mRNA ceRNA were obtained. Results of qRT-PCR partially demonstrated the bioinformatics results that the mRNA levels of CXCL10, CXCL9, CCL5, and CCR2 were remarkably elevated in tumor tissues collected from PTC patients coexistent with HT than those without HT (P < 0.001).

CONCLUSION

Our study offers a better understanding of the lncRNA-related ceRNA network involved in PTC with HT, providing novel key genes associated with PTC coexistent with HT.

Single-Cell Transcriptome Analysis Reveals Inter-Tumor Heterogeneity in Bilateral Papillary Thyroid Carcinoma

Background

The tumor microenvironment (TME) plays a pivotal role in cancer progression in papillary thyroid carcinoma (PTC), yet the composition and the phenotype of cells within the TME in bilateral PTC are poorly understood.

Methods

We performed unbiased transcriptome-wide single-cell RNA sequencing (scRNA-seq) analysis on 29,561 cells from 3 pairs of bilateral PTC and 1 non-tumor thyroid sample. The results of the analysis were validated by a large-scale bulk transcriptomic dataset deposited in The Cancer Genome Atlas (TCGA) database.

Results

Our integrative analysis of thyroid follicular cells revealed 42 signaling pathways enriched in malignant follicular cells, including cytokine-cytokine receptor interaction, PI3K/Akt signaling pathway, mitogen-activated protein kinase (MAPK) signaling pathway, and tumor necrosis factor (TNF) signaling pathway. A 6-gene signature (CXCL3, CXCL1, IL1A, CCL5, TNFRSF12A, and IL18) in the cytokine-cytokine receptor interaction pathway was constructed to predict the prognosis of patients with PTC, with high risk scores being associated with decreased overall survival [hazard ratio (HR) = 3.863, 95% CI = 2.233-6.682, p < 0.001]. Gene set variation analysis (GSVA) indicated that the pathways enriched in bilateral PTC were significantly different, indicating great heterogeneity in bilateral PTC, even with the same BRAF V600E mutation. Comprehensive analysis of T cells revealed that the proportion of CD8+ tissue-resident memory T cells expressing IFNG decreased in tumor samples with advanced N stage. Within the myeloid compartment, the ratio of suppressive M2-like to pro-inflammatory M1-like macrophages increased with advanced disease stage, which was confirmed in the bulk dataset using transcriptomic profiles. In addition, we also identified numerous biologically critical interactions among myeloid cells, T cells, and follicular cells, which were related to T-cell recruitment, M2-like macrophage polarization, malignant follicular cell progression, and T-cell inhibitory signaling.

Conclusion

Our integrative analyses revealed great inter-tumor heterogeneity within the TME in bilateral PTC, which will offer assistance for precise diagnosis and treatment.

Transcriptome Analysis of Liver Cancer Cell Huh-7 Treated With Metformin

Metformin is a kind of widely used antidiabetic drug that regulates glucose homeostasis by inhibiting liver glucose production and increasing muscle glucose uptake. Recently, some studies showed that metformin exhibits anticancer properties in a variety of cancers. Although several antitumor mechanisms have been proposed for metformin action, its mode of action in human liver cancer remains not elucidated. In our study, we investigated the underlying molecular mechanisms of metformin's antitumor effect on Huh-7 cells of hepatocellular carcinoma (HCC) in vitro. RNA sequencing was performed to explore the effect of metformin on the transcriptome of Huh-7 cells. The results revealed that 4,518 genes (with log2 fold change > 1 or < −1, adjusted p-value < 0.05) were differentially expressed in Huh-7 cells with treatment of 25-mM metformin compared with 0-mM metformin, including 1,812 upregulated and 2,706 downregulated genes. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses identified 54 classical pathways that were significantly enriched, and 16 pathways are closely associated with cancer, such as cell cycle, DNA replication, extracellular matrix–receptor interaction, and so on. We selected 11 differentially expressed genes, which are closely associated with HCC, to validate their differential expressions through a quantitative real-time reverse transcription-polymerase chain reaction. The result exhibited that the genes of fatty acid synthase, mini-chromosome maintenance complex components 6 and 5, myristoylated alanine-rich C-kinase substrate, fatty acid desaturase 2, C-X-C motif chemokine ligand 1, bone morphogenetic protein 4, S-phase kinase-associated protein 2, kininogen 1, and proliferating cell nuclear antigen were downregulated, and Dual-specificity phosphatase-1 is significantly upregulated in Huh-7 cells with treatment of 25-mM metformin. These differentially expressed genes and pathways might play a crucial part in the antitumor effect of metformin and might be potential targets of metformin treating HCC. Further investigations are required to evaluate the metformin mechanisms of anticancer action in vivo.

Roles and new Insights of Macrophages in the Tumor Microenvironment of Thyroid Cancer

Although most thyroid cancers have a good and predictable prognosis, the anaplastic, medullary, and refractory thyroid cancers still prone to recurrence and metastasis, resulting in poor prognosis. Although a number of newly developed targeted therapies have begun to be indicated for the above types of thyroid cancer in recent years, their ability to improve overall survival remain hindered by low efficacy. As the largest component of immune cells in tumor microenvironment, tumor-associated macrophages play a key role in the invasion and metastasis of thyroid cancer. There is much evidence that the immune system, tumor microenvironment and cancer stem cell interactions may revolutionize traditional therapeutic directions. Tumor-associated macrophages have been extensively studied in a variety of tumors, however, research on the relationship between thyroid cancer and macrophages is still insufficient. In this review, we summarize the functions of tumor-associated macrophages in different types of thyroid cancer, their cytokines or chemokines effect on thyroid cancer and the mechanisms that promote tumor proliferation and migration. In addition, we discuss the mechanisms by which tumor-associated macrophages maintain the stemness of thyroid cancer and potential strategies for targeting tumor-associated macrophages to treat thyroid cancer.

CCL18 Knockdown Suppresses Cell Growth and Migration in Thyroid Cancer

BACKGROUND

Chemokine (C-C motif) ligand 18 (CCL18) is a chemokine that plays a key role in immune and inflammatory responses. In recent years, CCL18 participates in the development and progression of various cancers, but its expression and role in thyroid cancer (TC) remain unclear.

METHODS

RT-qPCR assay and Western blot assay were used to explore the expression level of CCL18 in TC tissues and cells. Cell proliferation was measured by MTT assay. Transwell assay was adopted to detect cell migration in TC cells. Dual luciferase reporter assay was performed to assess the relationship between CCL18 and miR-149-5p.

RESULTS

There was an uptrend of CCL18 in TC tissues and cells. Our findings indicated that CCL18 overexpression facilitated lymph node metastasis in patients with TC. CCL18 silencing was found to inhibit cell migration, proliferation, and EMT progression in TC cells. CCL18 was proved to be a target gene of miR-149-5p. Additionally, miR-149-5p weakened the effect of CCL18 in the progression of TC.

CONCLUSION

Therefore, our results indicated that CCL18 knockdown restrained TC progression and suggested that CCL18 might be a potential therapeutic target for TC.

Identification of Transcriptional Pattern Related to Immune Cell Infiltration With Gene Co-Expression Network in Papillary Thyroid Cancer

BACKGROUND

A growing body of evidence suggests that immune cell infiltration in cancer is closely related to clinical outcomes. However, there is still a lack of research on papillary thyroid cancer (PTC).

METHODS

Based on single-sample gene set enrichment analysis (SSGSEA) algorithm and weighted gene co-expression network analysis (WGCNA) tool, the infiltration level of immune cell and key modules and genes associated with the level of immune cell infiltration were identified using PTC gene expression data from The Cancer Genome Atlas (TCGA) database. In addition, the co-expression network and protein-protein interactions network analysis were used to identify the hub genes. Moreover, the immunological and clinical characteristics of these hub genes were verified in TCGA and GSE35570 datasets and quantitative real-time polymerase chain reaction (PCR). Finally, receiver operating characteristic (ROC) curve analysis was used to evaluate the diagnostic value of hub genes.

RESULTS

Activated B cell, activated dendritic cell, CD56bright natural killer cell, CD56dim natural killer cell, Eosinophil, Gamma delta T cell, Immature dendritic cell, Macrophage, Mast cell, Monocyte, Natural killer cell, Neutrophil and Type 17 T helper cell were significantly changed between PTC and adjacent normal groups. WGCNA results showed that the black model had the highest correlation with the infiltration level of activated dendritic cells. We found 14 hub genes whose expression correlated to the infiltration level of activated dendritic cells in both TCGA and GSE35570 datasets. After validation in the TCGA dataset, the expression level of only 5 genes (C1QA, HCK, HLA-DRA, ITGB2 and TYROBP) in 14 hub genes were differentially expressed between PTC and adjacent normal groups. Meanwhile, the expression levels of these 5 hub genes were successfully validated in GSE35570 dataset. Quantitative real-time PCR results showed the expression of these 4 hub genes (except C1QA) was consistent with the results in TCGA and GSE35570 dataset. Finally, these 4 hub genes had diagnostic value to distinguish PTC and adjacent normal controls.

CONCLUSIONS

HCK, HLA-DRA, ITGB2 and TYROBP may be key diagnostic biomarkers and immunotherapy targets in PTC.

Fish and the Thyroid: A Janus Bifrons Relationship Caused by Pollutants and the Omega-3 Polyunsaturated Fatty Acids

Benefits of the omega-3 polyunsaturated fatty acids (PUFA) on a number of clinical disorders, including autoimmune diseases, are widely reported in the literature. One major dietary source of PUFA are fish, particularly the small oily fish, like anchovy, sardine, mackerel and others. Unfortunately, fish (particularly the large, top-predator fish like swordfish) are also a source of pollutants, including the heavy metals. One relevant heavy metal is mercury, a known environmental trigger of autoimmunity that is measurable inside the thyroid. There are a number of interactions between the omega-3 PUFA and thyroid hormones, even at the level of the thyroid hormone transport proteins. Concerning the mechanisms behind the protection from/amelioration of autoimmune diseases, including thyroiditis, that are caused by the omega-3 PUFA, one can be the decreased production of chemokines, a decrease that was reported in the literature for other nutraceuticals. Recent studies point also to the involvement of resolvins. The intracellular increase in resolvins is associated with the tissue protection from inflammation that was observed in experimental animals after coadministration of omega-3 PUFA and thyroid hormone. After having presented data on fish consumption at the beginning, we conclude our review by presenting data on the market of the dietary supplements/nutraceuticals. The global omega-3 products market was valued at USD 2.10 billion in 2020, and was projected to go up at a compound annual growth rate of 7.8% from 2020 to 2028. Among supplements, fish oils, which are derived mainly from anchovies, are considered the best and generally safest source of omega-3. Taking into account (i) the anti-autoimmunity and anti-cancer properties of the omega-3 PUFA, (ii) the increasing incidence of both autoimmune thyroiditis and thyroid cancer worldwide, (iii) the predisposing role for thyroid cancer exerted by autoimmune thyroiditis, and (iv) the risk for developing metabolic and cardiovascular disorders conferred by both elevated/trendwise elevated serum TSH levels and thyroid autoimmunity, then there is enough rationale for the omega-3 PUFA as measures to contrast the appearance and/or duration of Hashimoto's thyroiditis as well as to correct the slightly elevated serum TSH levels of subclinical hypothyroidism.

Serum Inflammation-based Scores in Endocrine Tumors

CONTEXT

Serum inflammation-based scores reflect systemic inflammatory response and/or patients' nutritional status, and may predict clinical outcomes in cancer. While these are well-described and increasingly used in different cancers, their clinical usefulness in the management of patients with endocrine tumors is less known.

EVIDENCE ACQUISITION

A comprehensive PubMed search was performed using the terms "endocrine tumor," "inflammation," "serum inflammation-based score," "inflammatory-based score," "inflammatory response-related scoring," "systemic inflammatory response markers," "neutrophil-to-lymphocyte ratio," "neutrophil-to-platelet ratio," "lymphocyte-to-monocyte ratio," "Glasgow prognostic score," "neutrophil-platelet score," "Systemic Immune-Inflammation Index," and "Prognostic Nutrition Index" in clinical studies.

EVIDENCE SYNTHESIS

The neutrophil-to-lymphocyte ratio and the platelet-to-lymphocyte ratio are the ones most extensively investigated in patients with endocrine tumors. Other scores have also been considered in some studies. Several studies focused in finding whether serum inflammatory biomarkers may stratify the endocrine tumor patients' risk and detect those at risk for developing more aggressive and/or refractory disease, particularly after endocrine surgery.

CONCLUSIONS

In this review, we summarize the current knowledge on the different serum inflammation-based scores and their usefulness in predicting the phenotype, clinical aggressiveness, and disease outcomes and prognosis in patients with endocrine tumors. The value of such serum inflammation-based scores in the management of patients with endocrine tumors has been emerging over the last decade. However, further research is necessary to establish useful markers and their cut-offs for routine clinical practice for individual diseases.

Cytokines: Application in recurrence appraisal for differentiated thyroid carcinoma and their relation with radioiodine ablation

The limitations in discriminating preablation disease-active status of differentiated thyroid carcinoma (DTC) still represent a major challenge to radioiodine dose management. Cytokines, the small protein signaling molecules that constitute the thyroid tumor microenvironment, play significant roles in the facilitation of intercellular communication and the control of tumorigenesis. Also, more attention should be paid to the molecular events within the innate and adaptive immune systems that occur after the organism being exposed to ionizing radiation. Therefore, we implemented a study of 260 patients with DTC in thyroid hormone withdrawal status who were treated with total thyroidectomy to explore the relationship between cytokines and recurrence/active disease status. Besides, we made a cross-sectional study to analyze pre- and post-ablation serum concentration of cytokines of 86 patients with DTC. There was a relationship between clinicohistopathological characteristics of patients with DTC and the presence of cytokines. It is noteworthy that patients with recurrence/active disease were at a higher serum interleukin-2 receptor (IL-2R) level than the disease-free patients (213.59 ± 75.43 pg/ml vs. 186.80 ± 77.40 pg/ml, P = 0.005). Positive correlation was observed between serum IL-2R and thyroglobulin (Tg) (P = 0.003). We also found significant changes in the cytokine profile after radioiodine ablation, including the decrease of tumor necrosis factor-α and IL-8 (P < 0.001, P < 0.001, respectively), and increase of IL-2R (P < 0.001). Thus, we suggest that serum IL-2R may assist in evaluating the disease status during the post-thyroidectomy follow-up and radioiodine therapy has an immunoregulatory effect on serum cytokines.

miR-181a, delivered by hypoxic PTC-secreted exosomes, inhibits DACT2 by downregulating MLL3, leading to YAP-VEGF-mediated angiogenesis

Papillary thyroid cancer (PTC) is the most common type of thyroid cancer, and angiogenesis plays critical roles in its recurrence and metastasis. In this study, we investigated the effects of hypoxia-induced exosomal microRNA-181 (miR-181a) from PTC on tumor growth and angiogenesis. Thyroid-cancer-related differentially expressed miR-181a was identified by microarray-based analysis in the Gene Expression Omnibus (GEO) database. We validated that miR-181a was highly expressed in PTC cells and even more so in cells cultured under hypoxic conditions, which also augmented exosome secretion from PTC cells. Exosomes extracted from PTC cells with manipulated miR-181a and mixed-lineage leukemia 3 (MLL3) were subjected to normoxic or hypoxic conditions. Human umbilical vein endothelial cells (HUVECs) were transfected with miR-181a inhibitor/mimic or small interfering RNA (siRNA)-MLL3 or treated with exosomes from hypoxic PTC cells. Hypoxic exosomal miR-181a delivery promoted proliferation and capillary-like network formation in HUVECs. Mechanistically, miR-181a targeted and inhibited MLL3. Furthermore, miR-181a downregulated DACT2 and upregulated YAP and vascular endothelial growth factor (VEGF). Further, hypoxic exosomal miR-181a induced angiogenesis and tumor growth in vivo, which was reversed by hypoxic exosomal miR-181a inhibitor. In conclusion, exosomal miR-181a from hypoxic PTC cells promotes tumor angiogenesis and growth through MLL3 and DACT2 downregulation, as well as VEGF upregulation.

Long Non-coding RNA RUNDC3A-AS1 Promotes Lung Metastasis of Thyroid Cancer via Targeting the miR-182-5p/ADAM9

Long non-coding RNAs (lncRNAs) have been identified as influential indicators in variety of malignancies. Among which, LncRNA RUNDC3A-AS1 is reported to upregulate in thyroid cancer. However, the expression pattern and the pathological function of lncRNA RUNDC3A-AS1 in thyroid cancer is unclear. In this study, we examined the expression levels of lncRNA RUNDC3A-AS1 in the thyroid cancer tissues and cell lines via RT-qPCR analysis. The effects of RUNDC3A-AS1 on thyroid cancer cell metastasis were detected by transwell chamber assay, scratch assay in vitro and lung metastasis model in vivo. The results indicated that RUNDC3A-AS1 was highly expressed in the thyroid cancer tissues and cell lines. Functionally, knockdown of RUNDC3A-AS1 could repress the migration and invasion of thyroid cancer cells in vitro, and inhibit thyroid cancer metastasis to lung in vivo. Mechanistically, RUNDC3A-AS1 served as an inhibitor of miR-182-5p in tumor tissues and cell lines. RUNDC3A-AS1 inhibited the expression of miR-182-5p to increase the expression level of ADAM9, thus further aggravating the malignancy of thyroid cancer. Therefore, the RUNDC3A-AS1/miR-182-5p/ADAM9 axis may be a potential therapeutic target for the treatment of thyroid cancer metastasis.

Expression level of long non-coding RNA colon adenocarcinoma hypermethylated serves as a novel prognostic biomarker in patients with thyroid carcinoma

The present study attempts to identify the prognostic value and potential mechanism of action of colorectal adenocarcinoma hypermethylated (CAHM) in thyroid carcinoma (THCA) by using the RNA sequencing (RNA-seq) dataset from The Cancer Genome Atlas (TCGA). The functional mechanism of CAHM was explored by using RNA-seq dataset and multiple functional enrichment analysis approaches. Connectivity map (CMap) online analysis tool was also used to predict CAHM targeted drugs. Survival analysis suggests that THCA patients with high CAHM expression have lower risk of death than the low CAHM expression (log-rank P=0.022, adjusted P=0.011, HR = 0.187, 95% confidence interval (CI) = 0.051–0.685). Functional enrichment of CAHM co-expression genes suggests that CAHM may play a role in the following biological processes: DNA repair, cell adhesion, DNA replication, vascular endothelial growth factor receptor, Erb-B2 receptor tyrosine kinase 2, ErbB and thyroid hormone signaling pathways. Functional enrichment of differentially expressed genes (DEGs) between low- and high-CAHM phenotype suggests that different CAHM expression levels may have the following differences in biological processes in THCA: cell adhesion, cell proliferation, extracellular signal-regulated kinase (ERK) 1 (ERK1) and ERK2 cascade, G-protein coupled receptor, chemokine and phosphatidylinositol-3-kinase-Akt signaling pathways. Connectivity map have identified five drugs (levobunolol, NU-1025, quipazine, anisomycin and sulfathiazole) for CAHM targeted therapy in THCA. Gene set enrichment analysis (GSEA) suggest that low CAHM phenotype were notably enriched in p53, nuclear factor κB, Janus kinase-signal transducer and activators of transcription, tumor necrosis factor, epidermal growth factor receptor and other signaling pathways. In the present study, we have identified that CAHM may serve as novel prognostic biomarkers for predicting overall survival (OS) in patients with THCA.

Long non-coding RNA LINC00488 facilitates thyroid cancer cell progression through miR-376a-3p/PON2

Objective: Long non-coding RNAs (lncRNAs) recently have been identified as influential indicators in a variety of malignancies. The aim of the present study was to identify a functional lncRNA LINC00488 and its effects on thyroid cancer in the view of cell proliferation and apoptosis.

Methods: In order to evaluate the effects of LINC00488 on the cellular process of thyroid cancer, we performed a series of in vitro experiments, including cell counting kit-8 (CCK-8) assay, EdU (5-ethynyl-2′-deoxyuridine) assay, flow cytometry, transwell chamber assay, Western blot and RT-qPCR. The target gene of LINC00488 was then identified by bioinformatics analysis (DIANA and TargetScan). Finally, a series of rescue experiments was conducted to validate the effect of LINC00488 and its target genes on proliferation, migration, invasion and apoptosis of thyroid cancer.

Results: Our findings revealed that LINC00488 was highly expressed in thyroid cancer cell lines (BCPAP, BHP5-16, TPC-1 and CGTH-W3) and promoted the proliferation, migration and invasion, while inhibited the apoptosis of thyroid cancer cells (BCPAP and TPC-1). The results of bioinformatics analysis and dual luciferase reporter gene assay showed that LINC00488 could directly bind to miR-376a-3p and down-regulated the expression level of miR-376a-3p. In addition, Paraoxonase-2 (PON2) was a target gene of miR-376a-3p and negatively regulated by miR-376a-3p. Rescue experiment indicated that LINC00488 might enhance PON2 expression by sponging miR-376a-3p in thyroid cancer.

Conclusion: Taken together, our study revealed that lncRNA LINC00488 acted as an oncogenic gene in the progression of thyroid cancer via regulating miR-376a-3p/PON2 axis, which indicated that LINC00488-miR-376a-3p-PON2 axis could serve as novel biomarkers or potential targets for the treatment of thyroid cancer.

Exosomal hsa-miR-129-2 and hsa-miR-889 from a 6-microRNA signature might be a potential biomarker for predicting prognosis of papillary thyroid carcinoma

BACKGROUND

Papillary thyroid carcinoma (PTC) is a subtype of thyroid cancer with increasing incidence over time.

OBJECTIVE

This study aimed to build a risk score (RS) system for PTC patients.

METHODS

PTC microRNA (miRNA) and messenger RNA (mRNA) expression data were extracted from The Cancer Genome Atlas (TCGA) database. The 491 PTC samples were randomly divided into training and validation sets. Using the limma software package, differentially expressed mRNAs (DEGs) and miRNAs (DEMs) between the tumor and control groups were screened. In order to construct an RS system, a survival package was used to select independent miRNAs related to prognosis. Enrichment analysis was performed, and a miRNA-mRNA co-expression network was constructed. High-throughput sequencing was also used to verify the prognostic miRNAs in exosomes.

RESULTS

We found 1363 DEGs and 171 DEMs between the tumor and control groups. After identifying 26 DEMs that were significantly related to prognosis, 6 independent prognosis-associated miRNAs were selected to build an RS system. The areas under the curves of the overall survival rates of the training, validation, and entire sets were 0.847, 0.772, and 0.819, respectively. By conducting pathway analysis using the miRNA-mRNA co-expression network, one overlapping factor and five overlapping pathways were obtained. In addition, high-throughput sequencing revealed that the hsa-miR-129-2, hsa-miR-548j, hsa-miR-6734, and hsa-miR-889 expression levels in TCGA tumor tissues and exosomes were consistent, and those of hsa-miR-129-2 and hsa-miR-889 between patients and controls were significantly different in exosomes.

CONCLUSION

The six-miRNA RS system in exosomes may comprise independent signatures for predicting PTC patient prognosis.

Tumor Microenvironment-Associated Pericyte Populations May Impact Therapeutic Response in Thyroid Cancer

Thyroid cancer is the most common endocrine malignancy, and aggressive radioactive iodine refractory thyroid carcinomas still lack an effective treatment. A deeper understanding of tumor heterogeneity and microenvironment will be critical to establishing new therapeutic approaches. One of the important influencing factors of tumor heterogeneity is the diversity of cells in the tumor microenvironment. Among these are pericytes, which play an important role in blood vessel stability and angiogenesis, as well as tumor growth and metastasis. Pericytes also have stem cell-like properties and are a heterogeneous cell population, and their lineage, which has been challenging to define, may impact tumor resistance at different tumor stages. Pericytes are also important stroma cell types in the angiogenic microenvironment which express tyrosine-kinase (TK) pathways (e.g., PDGFR-β). Although TK inhibitors (TKI) and BRAF^V600E inhibitors are currently used in the clinic for thyroid cancer, their efficacy is not durable and drug resistance often develops. Characterizing the range of distinct pericyte populations and distinguishing them from other perivascular cell types may enable the identification of their specific functions in the thyroid carcinoma vasculature. This remains an essential step in developing new therapeutic strategies. Also, assessing whether thyroid tumors hold immature and/or mature vasculature with pericyte populations coverage may be key to predicting tumor response to either targeted or anti-angiogenesis therapies. It is also critical to apply different markers in order to identify pericyte populations and characterize their cell lineage. This chapter provides an overview of pericyte ontogenesis and the lineages of diverse cell populations. We also discuss the role(s) and targeting of pericytes in thyroid carcinoma, as well as their potential impact on precision targeted therapies and drug resistance.

Akt isoform-specific effects on thyroid cancer development and progression in a murine thyroid cancer model

The Akt family is comprised of three unique homologous proteins with isoform-specific effects, but isoform-specific in vivo data are limited in follicular thyroid cancer (FTC), a PI3 kinase-driven tumor. Prior studies demonstrated that PI3K/Akt signaling is important in thyroid hormone receptor βPV/PV knock-in (PV) mice that develop metastatic thyroid cancer that most closely resembles FTC. To determine the roles of Akt isoforms in this model we crossed Akt1-/-, Akt2-/-, and Akt3-/- mice with PV mice. Over 12 months, thyroid size was reduced for the Akt null crosses (p < 0.001). Thyroid cancer development and local invasion were delayed in only the PVPV-Akt1 knock out (KO) mice in association with increased apoptosis with no change in proliferation. Primary-cultured PVPV-Akt1KO thyrocytes uniquely displayed a reduced cell motility. In contrast, loss of any Akt isoform reduced lung metastasis while vascular invasion was reduced with Akt1 or 3 loss. Microarray of thyroid RNA displayed incomplete overlap between the Akt KO models. The most upregulated gene was the dendritic cell (DC) marker CD209a only in PVPV-Akt1KO thyroids. Immunohistochemistry demonstrated an increase in CD209a-expressing cells in the PVPV-Akt1KO thyroids. In summary, Akt isoforms exhibit common and differential functions that regulate local and metastatic progression in this model of thyroid cancer.

The tumour microenvironment of pituitary neuroendocrine tumours

The tumour microenvironment (TME) includes a variety of non-neoplastic cells and non-cellular elements such as cytokines, growth factors and enzymes surrounding tumour cells. The TME emerged as a key modulator of tumour initiation, progression and invasion, with extensive data available in many cancers, but little is known in pituitary tumours. However, the understanding of the TME of pituitary tumours has advanced thanks to active research in this field over the last decade. Different immune and stromal cell subpopulations, and several cytokines, growth factors and matrix remodelling enzymes, have been characterised in pituitary tumours. Studying the TME in pituitary tumours may lead to a better understanding of tumourigenic mechanisms, identification of biomarkers useful to predict aggressive disease, and development of novel therapies. This review summarises the current knowledge on the different TME cellular/non-cellular elements in pituitary tumours and provides an overview of their role in tumourigenesis, biological behaviour and clinical outcomes.

RORγt may Influence the Microenvironment of Thyroid Cancer Predicting Favorable Prognosis

We aimed to investigate the role of RORγt (Retinoic acid-related orphan receptor gamma) in the tumor microenvironment of differentiated thyroid carcinoma. We retrospectively analyzed 56 patients (48 papillary and 8 follicular thyroid carcinomas). Immunohistochemical expression of RORγt was compared to other immune markers previously investigated by our group, clinical and pathological information. All patients presented cytoplasmic expression of RORγt in thyroid tumor cells. Seven (12.5%) patients presented no nuclear expression of RORγt. Positivity was few (up to 10%) in 14 patients; 10 to 50% in 5 patients (8.9%); and more than 50% in 30 patients (53.6%). Nuclear RORγt positivity was associated with absence of distant metastasis at diagnosis (p = 0.013) and the need of less cumulative doses of radioactive iodine (p = 0.039). Patients whose tumors were positive for nuclear RORγt presented higher 10-years relapse-free survival rate than those patients who were negative for RORγt (p = 0.023). We classified the patients according to the clustering of immunological immunohistochemical markers. We were able to distinguish a subset (A) of 38 patients who presented high expression of nuclear RORγt and tended to be scarce in proinflammatory immune markers. Other 16 patients integrated a second subset (B) whose tumor microenvironment accumulated proinflammatory markers and presented low expression of nuclear nuclear RORγt. Distant metastasis at diagnosis were more frequent among patients from cluster B than from cluster A (p = 0.008). Our results reinforce that the expression of RORγt together with other immune markers might help predict the prognosis of patients with thyroid cancer and help individualize clinical management.

The function of SNHG7/miR-449a/ACSL1 axis in thyroid cancer

Thyroid cancer (TC) has been characterized as the most common malignant malady of the endocrine system. Small nucleolar RNA host gene 7 (SNHG7) has been reported to serve as a key regulator in a large number of human cancer types, but its role in TC and the underlying regulatory mechanism have never been evaluated yet. The present study indicated that the expression of SNHG7 was markedly higher in TC cell lines. Knockdown of SNHG7 led to a suppression of TC cell progression and migration. Acyl-CoA synthetase long-chain family member 1 (ACSL1) has also been demonstrated as an oncogene in many cancers. Herein an inhibition of ACSL1 after SNHG7 knockdown was captured. Further, the suppressing effects of SNHG7 knockdown on TC cell processes were counteracted by ACSL1 overexpression. Data from online bioinformatics analysis, RNA immunoprecipitation, and luciferase reporter assays validated the interaction between microRNA-449a (miR-449a) and SNHG7 or ACSL1. It was also verified that SNHG7 sequestered miR-449a and therefore elevated ACSL1 expression levels. To conclude, the current study indicated that SNHG7 promoted proliferation and migration of TC cells by sponging miR-449a and therefore upregulating ACSL1. The present study may provide more explorations about the molecular regulation mechanism of long noncoding RNAs in TC progression.

Chemokines modulate the tumour microenvironment in pituitary neuroendocrine tumours

Non-tumoural cells within the tumour microenvironment (TME) influence tumour proliferation, invasiveness and angiogenesis. Little is known about TME in pituitary neuroendocrine tumours (PitNETs). We aimed to characterise the role of TME in the aggressive behaviour of PitNETs, focusing on immune cells and cytokines. The cytokine secretome of 16 clinically non-functioning PitNETs (NF-PitNETs) and 8 somatotropinomas was assessed in primary culture using an immunoassay panel with 42 cytokines. This was correlated with macrophage (CD68, HLA-DR, CD163), T-lymphocyte (CD8, CD4, FOXP3), B-lymphocyte (CD20), neutrophil (neutrophil elastase) and endothelial cells (CD31) content, compared to normal pituitaries (NPs, n = 5). In vitro tumour-macrophage interactions were assessed by conditioned medium (CM) of GH3 (pituitary tumour) and RAW264.7 (macrophage) cell lines on morphology, migration/invasion, epithelial-to-mesenchymal transition and cytokine secretion. IL-8, CCL2, CCL3, CCL4, CXCL10, CCL22 and CXCL1 are the main PitNET-derived cytokines. PitNETs with increased macrophage and neutrophil content had higher IL-8, CCL2, CCL3, CCL4 and CXCL1 levels. CD8+ T-lymphocytes were associated to higher CCL2, CCL4 and VEGF-A levels. PitNETs had more macrophages than NPs (p < 0.001), with a 3-fold increased CD163:HLA-DR macrophage ratio. PitNETs contained more CD4+ T-lymphocytes (p = 0.005), but fewer neutrophils (p = 0.047) with a 2-fold decreased CD8:CD4 ratio. NF-PitNETs secreted more cytokines and had 9 times more neutrophils than somatotropinomas (p = 0.002). PitNETs with higher Ki-67 had more FOXP3+ T cells, as well as lower CD68:FOXP3, CD8:CD4 and CD8:FOXP3 ratios. PitNETs with "deleterious immune phenotype" (CD68^hiCD4^hiFOXP3^hiCD20^hi) had a Ki-67 ≥ 3%. CD163:HLA-DR macrophage ratio was positively correlated with microvessel density (p = 0.015) and area (p < 0.001). GH3 cell-CM increased macrophage chemotaxis, while macrophage-CM changed morphology, invasion, epithelial-to-mesenchymal transition and secreted cytokines of GH3 cells. PitNETs are characterised by increased CD163:HLA-DR macrophage and reduced CD8:CD4 and CD8:FOXP3 T cell ratios. PitNET-derived chemokines facilitate macrophage, neutrophil and T cell recruitment into the tumours which can determine aggressive behaviour.

Nevirapine inhibits migration and invasion in dedifferentiated thyroid cancer cells

Background

Metastatic or recurrent thyroid cancer often behaves aggressively, and approximately two‐thirds of patients present with radioiodine resistance. Effective therapies to suppress thyroid cancer metastasis are urgently needed. Nevirapine has been proved to suppress tumor growth and induce differentiation in several tumor cells, but has not previously been evaluated in metastasis of thyroid cancer. The present study aimed to investigate the effect of nevirapine on migration and invasion in dedifferentiated thyroid cancer cells.

Methods

Human dedifferentiated thyroid cancer cell line (WRO 82‐1) was subject to real‐time quantitative PCR, western blot and transwell migration/invasion assays. The liver metastasis in tumor xenografts of nude mice was subject to hematoxylin‐eosin (HE) staining.

Results

Nevirapine significantly repressed cell migration and invasion in WRO 82‐1 cells, and surprisingly significantly decreased liver metastatic tumor in the nude mouse model of dedifferentiated thyroid cancer compared with that of the control. Moreover, nevirapine significantly decreased the expression of IL‐6 mRNA and phosphorylation of JAK2 (Y1007+Y1008) and STAT3 (Tyr 705) in WRO 82‐1 cells compared with those in control cells.

Conclusion

Our findings suggest that nevirapine significantly repressed migration and invasion/metastasis in WRO 82‐1 cells and tumor xenografts, which may be related to inhibition of IL‐6/STAT3 signaling pathway. It promises great potential as a novel therapy for thyroid cancer, especially for those patients with metastasis.

Papillary thyroid carcinoma behavior: clues in the tumor microenvironment

The prevalence of thyroid carcinoma is increasing and represents the most common endocrine malignancy, with papillary thyroid carcinoma (PTC) being the most frequent subtype. The genetic alterations identified in PTCs fail to distinguish tumors with different clinical behaviors, such as extra-thyroidal extension and lymph node metastasis. We hypothesize that the immune microenvironment may play a critical role in tumor invasion and metastasis. Computational immunogenomic analysis was performed on 568 PTC samples in The Cancer Genome Atlas using CIBERSORT, TIMER and TIDE deconvolution analytic tools for characterizing immune cell composition. Immune cell infiltrates were correlated with histologic type, mutational type, tumor pathologic T stage and lymph node N stage. Dendritic cells (DCs) are highly associated with more locally advanced tumor T stage (T3/T4, odds ratio (OR) = 2.6, CI = 1.4-4.5, P = 5.4 × 10-4). Increased dendritic cells (OR = 3.4, CI = 1.9-6.3, P = 5.5 × 10-5) and neutrophils (OR = 10.5, CI = 2.7-44, P = 8.7 × 10-4) significantly correlate with lymph node metastasis. In addition, dendritic cells positively correlate with tall cell morphology (OR = 4.5, CI = 1.6-13, P = 4.9 × 10-3) and neutrophils negatively correlate with follicular morphology (OR = 1.3 × 10-3, CI = 5.3 × 10-5-0.031, P = 4.1 × 10-5). TIDE analysis indicates an immune-exclusive phenotype that may be mediated by increased galectin-3 found in PTCs. Thus, characterization of the PTC immune microenvironment using three computational platforms shows that specific immune cells correlate with mutational type, histologic type, local tumor extent and lymph node metastasis. Immunologic evaluation of PTCs may provide a better indication of biologic behavior, resulting in the improved diagnosis and treatment of thyroid cancer.

LRG‑1 enhances the migration of thyroid carcinoma cells through promotion of the epithelial‑mesenchymal transition by activating MAPK/p38 signaling

Leucine-rich-alpha-2-glycoprotein 1 (LRG-1) has been reported to be associated with multiple malignancies. However, its participation in thyroid carcinoma progression remains unclear. In the present study, the biological function and underlying molecular mechanisms of LRG-1 in thyroid carcinoma were investigated. It was found that LRG-1 was overexpressed in thyroid carcinoma tissues, and high LRG-1 expression predicted poor patient survival and late tumor stage. As shown in the mouse xenograft study, knockdown of LRG-1 significantly attenuated thyroid cancer growth in vivo. Based on wound healing, Transwell, proliferation and apoptosis assays, it was found that the knockdown of LRG-1, using shLRG-1, inhibited cell migration and invasion, but did not affect proliferation and apoptosis in thyroid cancer cells. Furthermore, LRG-1 also induced epithelial-mesenchymal transition (EMT) in thyroid carcinoma cells. Western blot analysis revealed that this tumor-promoting bioactivity of LRG-1 was attributed to its selective activation of MAPK/p38 signaling. All of these findings indicate that LRG-1 plays a deleterious role in the progression of thyroid carcinoma. LRG-1 may serve as a promising biomarker for predicting prognosis in thyroid carcinoma patients, and LRG-1-based therapy may be developed into a novel strategy for the treatment of thyroid carcinoma.

Identifying gene modules of thyroid cancer associated with pathological stage by weighted gene co-expression network analysis

Thyroid cancer is the most common type of endocrine tumor. The TNM classification remains a standard for treatment determination and predicting prognosis in thyroid cancer. The genes modules associated with the progression of papillary thyroid carcinoma (PTC) were not clear. We applied a weighted gene co-expression network analysis (WGCNA) and differential expression analysis to systematically identified co-expressed gene modules and hub genes associated with PTC progression based on The Cancer Genome Atlas (TCGA) PTC transcriptome sequencing data. An independent validation cohort, GSE27155, was used to evaluate the preservation of gene modules. We identified two co-expressed genes modules associated with progression of PTC. Enrichment analysis indicated that the two modules were enriched in angiogenesis and extracellular matrix organization. DCN, COL1A1, COL1A2, COL5A2 and COL3A1 were hub genes in the co-expressed network. We systematically identified co-expressed gene modules and hub genes associated with PTC progression for the first time, which provided insights into the mechanisms underlying PTC progression and some potential targets for the treatment of PTC.

[Three-dimensional reconstruction of cervical CT vs ultrasound for estimating residual thyroid volume]

OBJECTIVE

To compare the accuracy of three-dimensional reconstruction of cervical CT and ultrasound for estimating residual thyroid volume.

METHODS

We performed a retrospective analysis of 17 patients with 21 residual thyroid glands undergoing thyroidectomy surgery between February, 2017 and March, 2018 in our department. We compared the residual thyroid volume in preoperative ultrasound with the intraoperative measurement and the volume measured using threedimensional CT reconstruction before surgery.

RESULTS

The maximum vertical and anterioposterior diameters of the residual thyroid measured by preoperative ultrasound differed significantly from the volume data measured intraoperatively (P < 0.05), but the difference in the maximum left-right diameters was not statistically significant (P>0.05). The maximum vertical, leftright, and anteroposterior diameters estimated by three-dimensional reconstruction of cervical CT was all similar with those measured intraoperatively (P>0.05).

CONCLUSIONS

Compared with ultrasound examination, three-dimensional reconstruction of neck CT is more accurate for estimating the residual thyroid volume and provides more reliable evidence for clinical calculation of postoperative I131 dose for thyroid cancer.

Identification of RNA Expression Profiles in Thyroid Cancer to Construct a Competing Endogenous RNA (ceRNA) Network of mRNAs, Long Noncoding RNAs (lncRNAs), and microRNAs (miRNAs)

Background

The aims of this study were to use RNA expression profile bioinformatics data from cases of thyroid cancer from the Cancer Genome Atlas (TCGA), the Kyoto Encyclopedia of Genes and Genomes (KEGG), and the Gene Ontology (GO) databases to construct a competing endogenous RNA (ceRNA) network of mRNAs, long noncoding RNAs (lncRNAs), and microRNAs (miRNAs).

Material/Methods

TCGA provided RNA profiles from 515 thyroid cancer tissues and 56 normal thyroid tissues. The DESeq R package analyzed high-throughput sequencing data on differentially expressed RNAs. GO and KEGG pathway analysis used the DAVID 6.8 and the ClusterProfile R package. Kaplan-Meier survival statistics and Cox regression analysis were performed. The thyroid cancer ceRNA network was constructed based on the miRDB, miRTarBase, and TargetScan databases.

Results

There were 1,098 mRNAs associated with thyroid cancer; 101 mRNAs were associated with overall survival (OS). Multivariate analysis developed a risk scoring system that identified seven signature mRNAs, with a discriminative value of 0.88, determined by receiver operating characteristic (ROC) curve analysis. A ceRNA network included 13 mRNAs, 31 lncRNAs, and seven miRNAs. Four out of the 31 lncRNAs and all miRNAs were down-regulated, and the remaining RNAs were upregulated. Two lncRNAs (MIR1281A2HG and OPCML-IT1) and one miRNA (miR-184) were significantly associated with OS in patients with thyroid cancer.

Conclusions

Differential RNA expression profiling in thyroid cancer was used to construct a ceRNA network of mRNAs, lncRNAs, and miRNAs that showed potential in evaluating prognosis.

Development of a prognostic index based on an immunogenomic landscape analysis of papillary thyroid cancer

BACKGROUND

Papillary thyroid cancer (PTC) is the most common subtype of thyroid cancer, and inflammation relates significantly to its initiation and prognosis. Systematic exploration of the immunogenomic landscape therein to assist in PTC prognosis is therefore urgent. The Cancer Genome Atlas (TCGA) project provides a large number of genetic PTC samples that enable a comprehensive and reliable immunogenomic study.

METHODS

We integrated the expression profiles of immune-related genes (IRGs) and progression-free intervals (PFIs) in survival in 493 PTC patients based on the TCGA dataset. Differentially-expressed and survival-associated IRGs in PTC patients were estimated a computational difference algorithm and COX regression analysis. The potential molecular mechanisms and properties of these PTC-specific IRGs were also explored with the help of computational biology. A new prognostic index based on immune-related genes was developed by using multivariable COX analysis.

RESULTS

A total of 46 differentially expressed immune-related genes were significantly correlated with clinical outcome of PTC patients. Functional enrichment analysis revealed that these genes were actively involved in a cytokine-cytokine receptor interaction KEGG pathway. A prognostic signature based on RGs (AGTR1, CTGF, FAM3B, IL11, IL17C, PTH2R and SPAG11A) performed moderately in prognostic predictions and correlated with age, tumor stage, metastasis, number of lesions, and tumor burden. Intriguingly, the prognostic index based on IRGs reflected infiltration by several types of immune cells.

CONCLUSIONS

Together, our results screened several IRGs of clinical significance, revealed drivers of the immune repertoire, and demonstrated the importance of a personalized, IRG-based immune signature in the recognition, surveillance, and prognosis of PTC.

Molecular Network-Based Identification of Competing Endogenous RNAs in Thyroid Carcinoma

RNAs may act as competing endogenous RNAs (ceRNAs), a critical mechanism in determining gene expression regulations in many cancers. However, the roles of ceRNAs in thyroid carcinoma remains elusive. In this study, we have developed a novel pipeline called Molecular Network-based Identification of ceRNA (MNIceRNA) to identify ceRNAs in thyroid carcinoma. MNIceRNA first constructs micro RNA (miRNA)-messenger RNA (mRNA)long non-coding RNA (lncRNA) networks from miRcode database and weighted correlation network analysis (WGCNA), based on which to identify key drivers of differentially expressed RNAs between normal and tumor samples. It then infers ceRNAs of the identified key drivers using the long non-coding competing endogenous database (lnCeDB). We applied the pipeline into The Cancer Genome Atlas (TCGA) thyroid carcinoma data. As a result, 598 lncRNAs, 1025 mRNAs, and 90 microRNA (miRNAs) were inferred to be differentially expressed between normal and thyroid cancer samples. We then obtained eight key driver miRNAs, among which hsa-mir-221 and hsa-mir-222 were key driver RNAs identified by both miRNA-mRNA-lncRNA and WGCNA network. In addition, hsa-mir-375 was inferred to be significant for patients' survival with 34 associated ceRNAs, among which RUNX2, DUSP6 and SEMA3D are known oncogenes regulating cellular proliferation and differentiation in thyroid cancer. These ceRNAs are critical in revealing the secrets behind thyroid cancer progression and may serve as future therapeutic biomarkers.

Role of Chemokines in Thyroid Cancer Microenvironment: Is CXCL8 the Main Player?

Tumor-related inflammation does influence the biological behavior of neoplastic cells and ultimately the patient's outcome. With specific regard to thyroid cancer, the issue of tumor-associated inflammation has been extensively studied and recently reviewed. However, the role of chemokines, which play a crucial role in determining the immuno-phenotype of tumor-related inflammation, was not addressed in previous reviews on the topic. Experimental evidence shows that thyroid cancer cells actively secrete a wide spectrum of chemokines and, at least for some of them, solid scientific data support a role for these immune-active molecules in the aggressive behavior of the tumor. Our proposal for a review article on chemokines and thyroid cancer stems from the notion that chemokines, besides having the ability to attract and maintain immune cells at the tumor site, also produce several pro-tumorigenic actions, which include proangiogenetic, cytoproliferative, and pro-metastatic effects. Studies taking into account the role of CCL15, C-X-C motif ligand 12, CXCL16, CXCL1, CCL20, and CCL2 in the context of thyroid cancer will be reviewed with particular emphasis on CXCL8. The reason for focusing on CXCL8 is that this chemokine is the most studied one in human malignancies, displaying multifaceted pro-tumorigenic effects. These include enhancement of tumor cells growth, metastatization, and angiogenesis overall contributing to the progression of several cancers including thyroid cancer. We aim at reviewing current knowledge on the (i) ability of both normal and tumor thyroid cells to secrete CXCL8; (ii) direct/indirect pro-tumorigenic effects of CXCL8 demonstrated by in vitro and in vivo studies specifically performed on thyroid cancer cells; and (iii) pharmacologic strategies proven to be effective for lowering CXCL8 secretion and/or its effects on thyroid cancer cells.