Molecular classification of papillary thyroid carcinoma: distinct BRAF, RAS, and RET/PTC mutation-specific gene expression profiles discovered by DNA microarray analysis

SHMT2 promotes papillary thyroid cancer metastasis through epigenetic activation of AKT signaling

Cancer cells alter their metabolism and epigenetics to support cancer progression. However, very few modulators connecting metabolism and epigenetics have been uncovered. Here, we reveal that serine hydroxymethyltransferase-2 (SHMT2) generates S-adenosylmethionine (SAM) to epigenetically repress phosphatase and tensin homolog (PTEN), leading to papillary thyroid cancer (PTC) metastasis depending on activation of AKT signaling. SHMT2 is elevated in PTC, and is associated with poor prognosis. Overexpressed SHMT2 promotes PTC metastasis both in vitro and in vivo. Proteomic enrichment analysis shows that AKT signaling is activated, and is positively associated with SHMT2 in PTC specimens. Blocking AKT activation eliminates the effects of SHMT2 on promoting PTC metastasis. Furthermore, SHMT2 expression is negatively associated with PTEN, a negative AKT regulator, in PTC specimens. Mechanistically, SHMT2 catalyzes serine metabolism and produces activated one-carbon units that can generate SAM for the methylation of CpG islands in PTEN promoter for PTEN suppression and following AKT activation. Importantly, interference with PTEN expression affects SHMT2 function by promoting AKT signaling activation and PTC metastasis. Collectively, our research demonstrates that SHMT2 connects metabolic reprogramming and epigenetics, contributing to the poor progression of PTC.

NF-κB signaling activation and roles in thyroid cancers: implication of MAP3K14/NIK

Among follicular-derived thyroid cancers (TC), those with aggressive behavior and resistance to current treatments display poor prognosis. NF-κB signaling pathways are involved in tumor progression of various cancers. Here, we finely characterize the NF-κB pathways and their involvement in TC. By using immunoblot and gel shift assays, we demonstrated that both classical and alternative NF-κB pathways are activated in ten TC-derived cell lines, leading to activated RelA/p50 and RelB/p50 NF-κB dimers. By analyzing the RNAseq data of the large papillary thyroid carcinoma (PTC) cohort from The Cancer Genome Atlas (TCGA) project, we identified a tumor progression-related NF-κB signature in BRAFV600E mutated-PTCs. That corroborated with the role of RelA and RelB in cell migration and invasion processes that we demonstrated specifically in BRAFV600E mutated-cell lines, together with their role in the control of expression of genes implicated in invasiveness (MMP1, PLAU, LCN2 and LGALS3). We also identified NF-κB-inducing kinase (NIK) as a novel actor of the constitutive activation of the NF-κB pathways in TC-derived cell lines. Finally, its implication in invasiveness and its overexpression in PTC samples make NIK a potential therapeutic target for advanced TC treatment.

Progress in Thyroid Cancer Genomics: A 40-Year Journey

Background: Very little was known about the molecular pathogenesis of thyroid cancer until the late 1980s. As part of the Centennial celebration of the American Thyroid Association, we review the historical discoveries that contributed to our current understanding of the genetic underpinnings of thyroid cancer. Summary: The pace of discovery was heavily dependent on scientific breakthroughs in nucleic acid sequencing technology, cancer biology, thyroid development, thyroid cell signaling, and growth regulation. Accordingly, we attempt to link the primary observations on thyroid cancer molecular genetics with the methodological and scientific advances that made them possible. Conclusions: The major genetic drivers of the common forms of thyroid cancer are now quite well established and contribute to a significant extent to how we diagnose and treat the disease. However, many challenges remain. Future work will need to unravel the complexity of thyroid cancer ecosystems, which is likely to be a major determinant of their biological behavior and on how they respond to therapy.

High expression of IL4I1 is correlated with poor prognosis and immune infiltration in thyroid cancer

BACKGROUND

Thyroid cancer-related deaths mostly result from metastasis. It was reported that the immunometabolism associated enzyme interleukin-4-induced-1 (IL4I1) was related to tumor metastasis. The present study was intended to investigate the effects of IL4I1 on thyroid cancer metastasis and its relationship with the prognosis.

METHODS

Data from Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were analyzed to find out the different mRNA expressions of IL4I1 between thyroid cancer and normal tissues. And Human Protein Atlas (HPA) was used to assess IL4I1 protein expression. To further differentiate thyroid cancer from normal tissues and estimate the impact of IL4I1 on the prognosis, the receiver operating characteristic curve (ROC) and Kaplan-Meier (KM) method was performed. The protein-protein interaction (PPI) network was established using STRING, and functional enrichment analyses were conducted by "clusterProfiler" package. Then, we assayed the correlation between IL4I1 and some related molecules. The relationship between IL4I1 and immune infiltration was performed using "Gene Set Variation Analysis (GSVA)" package in TCGA and tumor-immune system interaction database (TISIDB). Finally, we did in vitro experiments in order to further prove the bioeffects of IL4I1 on metastasis.

RESULTS

The expression of IL4I1 mRNA and IL4I1 protein was significantly upregulated in thyroid cancer tissues. The increment of IL4I1 mRNA expression was related to high-grade malignancy, lymph node metastases and extrathyroidal extension. The ROC curve displayed the cutoff value of 0.782, with the sensitivity of 77.5% and the specificity of 77.8%. KM survival analysis showed that there was a worse PFS in patients with high IL4I1 expression than those with low IL4I1 expression (p = 0.013). Further study indicated that IL4I1 was associated with lactate, body fluid secretion, positive regulation of T cell differentiation, and cellular response to nutrients in Gene Ontology (GO) analysis. Moreover, IL4I1 was found correlated with immune infiltration. Finally, the in vitro experiments revealed the promotion of IL4I1 on cancer cell proliferation, migration and invasion.

CONCLUSIONS

The increased IL4I1 expression is markedly correlated with the immune imbalance in the tumor microenvironment (TME) and predicts poor survival in thyroid cancer. This study reveals the potential clinical biomarker of poor prognosis and the target of immune therapy in thyroid cancer.

Novel thrombospondin-1 transcript exhibits distinctive expression and activity in thyroid tumorigenesis

Thrombospondin 1 (TSP1) is known for its cell-specific functions in cancer progression, such as proliferation and migration. It contains 22 exons that may potentially produce several different transcripts. Here, we identified TSP1V as a novel TSP1-splicing variant produced by intron retention (IR) in human thyroid cancer cells and tissues. We observed that TSP1V functionally inhibited tumorigenesis contrary to TSP1 wild-type, as identified in vivo and in vitro. These activities of TSP1V are caused by inhibiting phospho-Smad and phospho-focal adhesion kinase. Reverse transcription polymerase chain reaction and minigene experiments revealed that some phytochemicals/non-steroidal anti-inflammatory drugs enhanced IR. We further found that RNA-binding motif protein 5 (RBM5) suppressed IR induced by sulindac sulfide treatment. Additionally, sulindac sulfide reduced phospho-RBM5 levels in a time-dependent manner. Furthermore, trans-chalcone demethylated TSP1V, thereby preventing methyl-CpG-binding protein 2 binding to TSP1V gene. In addition, TSP1V levels were significantly lower in patients with differentiated thyroid carcinoma than in those with benign thyroid nodule, indicating its potential application as a diagnostic biomarker in tumor progression.

Identification of novel characteristic biomarkers and immune infiltration profile for the anaplastic thyroid cancer via machine learning algorithms

PURPOSE

Anaplastic thyroid cancer (ATC) is a rare and lethal malignant cancer. In recent years, the application of molecular-driven targeted therapy and immunotherapy has markedly improved the prognosis of ATC. This study aimed to identify characteristic genes for ATC diagnosis and revealed the role of ATC characteristic genes in drug sensitivity and immune cell infiltration.

METHODS

We downloaded ATC RNA-sequencing data from the GEO database. Following the combination and normalization of the dataset, we first divided the combined datasets into the training cohort and the validation cohort. We identified differentially expressed genes (DEGs) in ATC by differential expression analysis in the training cohort. We used two machine learning algorithms, least absolute shrinkage and selection operator (LASSO) and support vector machine-recursive feature elimination (SVM-RFE) to identify ATC characteristic genes. The CIBERSORT algorithm was performed to calculate the abundance of various immune cells in ATC. Finally, we validated the expression of ATC characteristic genes by quantitative RT-PCR (RT-qPCR) in ATC cell lines and immunohistochemistry (IHC).

RESULTS

A total of 425 DEGs were identified in the training cohort, including 240 upregulated genes and 185 downregulated genes. Four ATC characteristic genes (ADM, PXDN, MMP1, and TFF3) were identified, and their diagnostic value was validated in the validation cohort (AUC in ROC analysis > 0.75). We established a practical gene expression-based nomogram to accurately predict the probability of ATC. We also found that ATC characteristic biomarkers are associated with the tumor immune microenvironment and drug sensitivity.

CONCLUSION

ADM, PXDN, MMP1, and TFF3 might serve as potential ATC diagnostic biomarkers and may be helpful for ATC molecular targeted therapy and immunotherapy.

Follicular Neoplasm of Thyroid Revisited: Current Differential Diagnosis and the Impact of Molecular Testing

The diagnosis of "follicular neoplasm" (FN) in thyroid cytopathology has a long history that originated not long after the practice of fine-needle aspiration (FNA) of thyroid nodules. From the outset, this interpretive category was intended to convey a set of differential diagnoses rather than a precise diagnosis, as key diagnostic features, such as capsular and vascular invasion, were not detectable on cytology preparations. Cytologic-histologic correlation studies over the past several decades have shown that FN interpretation can be applied to the spectrum of nonneoplastic tumors to carcinomas. Most tumors classified as FN include follicular adenoma, follicular carcinoma, noninvasive follicular thyroid tumor with papillary-like nuclear features, and follicular variant of papillary thyroid carcinoma. Less common entities that may be classified as FN on FNA include hyalinizing trabecular tumor (HTT), poorly differentiated thyroid carcinoma, medullary carcinoma, and nonthyroidal lesions such as parathyroid tissue, paraganglioma, and metastatic tumors. Advances in our ability to detect characteristic molecular alterations (eg, GLIS gene rearrangements for hyalinizing trabecular tumor) in FNA samples may assist in the identification of some of these entities. In this review, we summarize the pathophysiology, history, and evolution of the terminology and the current differential diagnosis according to the recently published 2022 World Health Organization classification, molecular testing, and management of nodules classified as FN.

Induction of BRAF^V600E Mutation in Thyroid Cells Leads to Frequent Hypermethylation

Objectives.

The BRAF^V600E mutation is a major driver mutation in papillary thyroid cancer. The aim of this study was to elucidate the correlation between DNA methylation and gene expression changes induced by the BRAF^V600E mutation in thyroid cells.

Methods.

We used Nthy/BRAF cell lines generated by transfection of Nthy/ori cells with the wild-type BRAF gene (Nthy/WT cells) and the V600E mutant-type BRAF gene (Nthy/V600E cells). We performed gene expression microarray and DNA methylation array analyses for Nthy/WT and Nthy/V600E cells. Two types of array data were integrated to identify inverse correlations between methylation and gene expression. The results were verified in silico using data from The Cancer Genome Atlas (TCGA) and in vivo through pyrosequencing and quantitative real-time polymerase chain reaction (qRT-PCR).

Results.

In the Nthy/V600E cells, 199,821 probes were significantly hypermethylated, and 697 genes showed a “hypermethylation-downregulation” pattern in Nthy/V600E. Tumor suppressor genes and apoptosis-related genes were included. In total, 66,446 probes were significantly hypomethylated, and 227 genes showed a “hypomethylation-upregulation” pattern in Nthy/V600E cells. Protooncogenes and developmental protein-coding genes were included. In the TCGA analysis, 491/697 (70.44%) genes showed a hypermethylation-downregulation pattern, and 153/227 (67.40%) genes showed a hypomethylation-upregulation pattern. Ten selected genes showed a similar methylation-gene expression pattern in pyrosequencing and qRT-PCR.

Conclusion.

Induction of the BRAF^V600E mutation in thyroid cells led to frequent hypermethylation. Anticancer genes, such as those involved in tumor suppression or apoptosis, were downregulated by upstream hypermethylation, whereas carcinogenic genes, such as protooncogenes, were upregulated by hypomethylation. Our results suggest that the BRAF^V600E mutation in thyroid cells modulates DNA methylation and results in cancer-related gene expression.

Targeting non-canonical pathways as a strategy to modulate the sodium iodide symporter

The sodium iodide symporter (NIS) functions to transport iodide and is critical for successful radioiodide ablation of cancer cells. Approaches to bolster NIS function and diminish recurrence post-radioiodide therapy are impeded by oncogenic pathways that suppress NIS, as well as the inherent complexity of NIS regulation. Here, we utilize NIS in high-throughput drug screening and undertake rigorous evaluation of lead compounds to identify and target key processes underpinning NIS function. We find that multiple proteostasis pathways, including proteasomal degradation and autophagy, are central to the cellular processing of NIS. Utilizing inhibitors targeting distinct molecular processes, we pinpoint combinatorial drug strategies giving robust >5-fold increases in radioiodide uptake. We also reveal significant dysregulation of core proteostasis genes in human tumors, identifying a 13-gene risk score classifier as an independent predictor of recurrence in radioiodide-treated patients. We thus propose and discuss a model for targetable steps of intracellular processing of NIS function.

Identification of a 3-Gene Prognostic Index for Papillary Thyroid Carcinoma

The accurate determination of the risk of cancer recurrence is a critical unmet need in managing thyroid cancer (TC). Although numerous studies have successfully demonstrated the use of high throughput molecular diagnostics in TC prediction, it has not been successfully applied in routine clinical use, particularly in Chinese patients. In our study, we objective to screen for characteristic genes specific to PTC and establish an accurate model for diagnosis and prognostic evaluation of PTC. We screen the differentially expressed genes by Python 3.6 in The Cancer Genome Atlas (TCGA) database. We discovered a three-gene signature Gap junction protein beta 4 (GJB4), Ripply transcriptional repressor 3 (RIPPLY3), and Adrenoceptor alpha 1B (ADRA1B) that had a statistically significant difference. Then we used Gene Expression Omnibus (GEO) database to establish a diagnostic and prognostic model to verify the three-gene signature. For experimental validation, immunohistochemistry in tissue microarrays showed that thyroid samples’ proteins expressed by this three-gene are differentially expressed. Our protocol discovered a robust three-gene signature that can distinguish prognosis, which will have daily clinical application.

Identification of Dysregulated Expression of G Protein Coupled Receptors in Endocrine Tumors by Bioinformatics Analysis: Potential Drug Targets?

Background: Many studies link G protein-coupled receptors (GPCRs) to cancer. Some endocrine tumors are unresponsive to standard treatment and/or require long-term and poorly tolerated treatment. This study explored, by bioinformatics analysis, the tumoral profiling of the GPCR transcriptome to identify potential targets in these tumors aiming at drug repurposing. Methods: We explored the GPCR differentially expressed genes (DEGs) from public datasets (Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA)). The GEO datasets were available for two medullary thyroid cancers (MTCs), eighty-seven pheochromocytomas (PHEOs), sixty-one paragangliomas (PGLs), forty-seven pituitary adenomas and one-hundred-fifty adrenocortical cancers (ACCs). The TCGA dataset covered 92 ACCs. We identified GPCRs targeted by approved drugs from pharmacological databases (ChEMBL and DrugBank). Results: The profiling of dysregulated GPCRs was tumor specific. In MTC, we found 14 GPCR DEGs, including an upregulation of the dopamine receptor (DRD2) and adenosine receptor (ADORA2B), which were the target of many drugs. In PGL, seven GPCR genes were downregulated, including vasopressin receptor (AVPR1A) and PTH receptor (PTH1R), which were targeted by approved drugs. In ACC, PTH1R was also downregulated in both the GEO and TCGA datasets and was the target of osteoporosis drugs. Conclusions: We highlight specific GPCR signatures across the major endocrine tumors. These data could help to identify new opportunities for drug repurposing.

Characterization of the CpG island methylator phenotype subclass in papillary thyroid carcinoma

CpG island methylator phenotype (CIMP), characterized by the concurrent and widespread hypermethylation of a cluster of CpGs, has been reported to play an important role in carcinogenesis. Limited studies have explored the role of CIMP in papillary thyroid carcinomas (PTCs). Here, in genome-wide DNA methylation analysis of 350 primary PTCs from the Cancer Genome Atlas database that were assessed using the Illumina HumanMethylation450K platform, our study helps to identify two subtypes displayed markedly distinct DNA methylation levels, termed CIMP (high levels of DNA methylation) and nCIMP subgroup (low levels of DNA methylation). Interestingly, PTCs with CIMP tend to have a higher degree of malignancy, since this subtype was tightly associated with older age, advanced pathological stage, and lymph node metastasis (all P < 0.05). Differential methylation analysis showed a broad methylation gain in CIMP and subsequent generalized gene set testing analysis based on the significantly methylated probes in CIMP showed remarkable enrichment in epithelial mesenchymal transition and angiogenesis hallmark pathways, confirming that the CIMP phenotype may promote the tumor progression from another perspective. Analysis of tumor microenvironment showed that CIMP PTCs are in an immune-depletion status, which may affect the effectiveness of immunotherapy. Genetically, the significantly higher tumor mutation burden and copy number alteration both at the genome and focal level confirmed the genomic heterogeneity and chromosomal instability of CIMP. tumor Corresponding to the above findings, PTC patients with CIMP showed remarkable poor clinical outcome as compared to nCIMP regarding overall survival and progression-free survival. More importantly, CIMP was associated with worse survival independent of known prognostic factors.

Characterization of immune landscape in papillary thyroid cancer reveals distinct tumor immunogenicity and implications for immunotherapy

Although the vast majority of patients with papillary thyroid cancer (PTC) have a favorable prognosis when conventional treatments are implemented, local recurrence and distant metastasis of advanced PTCs still hamper the survival and clinical management in certain patients. As immune checkpoint blockade (ICB) therapy achieves a great success in some advanced cancers, we aimed to investigate the immune landscape in PTC and its potential implications for prognosis and immunotherapy. In this study, different algorithms were conducted to estimate immune infiltration in PTC samples. A series of bioinformatic and machine learning approaches were performed to identify PTC-specific immune-related genes (IRGs) and distinct immune clusters. Differences in intrinsic tumor immunogenicity and potential immunotherapy response were observed between distinct immune clusters. A prognostic immune-related signature (IRS) was established to predict progression-free survival (PFS). IRS exhibited more powerful prognostic capacity and accurate survival prediction compared to conventional clinicopathological features. Furthermore, an integrated survival decision tree and a scoring nomogram were constructed to improve prognostic stratification and predictive accuracy for individual patients. In addition, altered pathways, mutational patterns, and potential applicable drugs were analyzed in different immune-related risk groups. Our study gained some insight into the immune landscape of PTC, and provided some useful clues for introducing immune-based molecular classification into risk stratification and guiding ICB decision-making.

Characterizing dedifferentiation of thyroid cancer by integrated analysis

Understanding of dedifferentiation, an indicator of poo prognosis for patients with thyroid cancer, has been hampered by imprecise and incomplete characterization of its heterogeneity and its attributes. Using single-cell RNA sequencing, we explored the landscape of thyroid cancer at single-cell resolution with 46,205 cells and delineated its dedifferentiation process and suppressive immune microenvironment. The developmental trajectory indicated that anaplastic thyroid cancer (ATC) cells were derived from a small subset of papillary thyroid cancer (PTC) cells. Moreover, a potential functional role of CREB3L1 on ATC development was revealed by integrated analyses of copy number alteration and transcriptional regulatory network. Multiple genes in differentiation-related pathways (e.g., EMT) were involved as the downstream targets of CREB3L1, increased expression of which can thus predict higher relapse risk of PTC. Collectively, our study provided insights into the heterogeneity and molecular evolution of thyroid cancer and highlighted the potential driver role of CREB3L1 in its dedifferentiation process.

Quercetin Induces Anticancer Activity by Upregulating Pro-NAG-1/GDF15 in Differentiated Thyroid Cancer Cells

Simple Summary

Thyroid cancer is one of the most common cancers worldwide, and its incidence has increased over the last few decades. It is difficult to diagnose different types of thyroid cancer. Tumor tissues from papillary thyroid cancer patient showed higher expression of mature NAG-1, whereas adjacent normal tissues showed higher expression of pro-NAG-1. Several anti-cancer compounds increased pro-NAG-1 expression in thyroid cancer cell line. Quercetin (3,3’,4’,5,7-pentahydroxyflavone) is a flavonoid that is a major component of various plants, including raspberries, grapes, and onions. Quercetin induced apoptosis by inducing only pro-NAG-1 expression, but not mature NAG-1, mediated by the transcription factor C/EBP. This study indicates that pro-NAG-1 could be used as a useful biomarker for thyroid cancer and also provides a potential therapeutic target for the treatment of thyroid cancer with quercetin.

Abstract

Although the treatment of thyroid cancer has improved, unnecessary surgeries are performed due to a lack of specific diagnostic and prognostic markers. Therefore, the identification of novel biomarkers should be considered in the diagnosis and treatment of thyroid cancer. In this study, antibody arrays were performed using tumor and adjacent normal tissues of patients with papillary thyroid cancer, and several potential biomarkers were identified. Among the candidate proteins chosen based on the antibody array data, mature NAG-1 exhibited increased expression in tumor tissues compared to adjacent normal tissues. In contrast, pro-NAG-1 expression increased in normal tissues, as assessed by western blot analysis. Furthermore, pro-NAG-1 expression was increased when the thyroid cancer cells were treated with phytochemicals and nonsteroidal anti-inflammatory drugs in a dose-dependent manner. In particular, quercetin highly induced the expression of pro-NAG-1 but not that of mature NAG-1, with enhanced anticancer activity, including apoptosis induction and cell cycle arrest. Examination of the NAG-1 promoter activity showed that p53, C/EBPα, or C/EBPδ played a role in quercetin-induced NAG-1 expression. Overall, our study indicated that NAG-1 may serve as a novel biomarker for thyroid cancer prognosis and may be used as a therapeutic target for thyroid cancers.

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.

A Novel Nanoproteomic Approach for the Identification of Molecular Targets Associated with Thyroid Tumors

A thyroid nodule is the most common presentation of thyroid cancer; thus, it is extremely important to differentiate benign from malignant nodules. Within malignant lesions, classification of a thyroid tumor is the primary step in the assessment of the prognosis and selection of treatment. Currently, fine-needle aspiration biopsy (FNAB) is the preoperative test most commonly used for the initial thyroid nodule diagnosis. However, due to some limitations of FNAB, different high-throughput “omics” approaches have emerged that could further support diagnosis based on histopathological patterns. In the present work, formalin-fixed paraffin-embedded (FFPE) tissue specimens from normal (non-neoplastic) thyroid (normal controls (NCs)), benign tumors (follicular thyroid adenomas (FTAs)), and some common types of well-differentiated thyroid carcinoma (follicular thyroid carcinomas (FTCs), conventional or classical papillary thyroid carcinomas (CV-PTCs), and the follicular variant of papillary thyroid carcinomas (FV-PTCs)) were analyzed. For the first time, FFPE thyroid samples were deparaffinized using an easy, fast, and non-toxic method. Protein extracts from thyroid tissue samples were analyzed using a nanoparticle-assisted proteomics approach combined with shotgun LC-MS/MS. The differentially regulated proteins found to be specific for the FTA, FTC, CV-PTC, and FV-PTC subtypes were analyzed with the bioinformatic tools STRING and PANTHER showing a profile of proteins implicated in the thyroid cancer metabolic reprogramming, cancer progression, and metastasis. These proteins represent a new source of potential molecular targets related to thyroid tumors.

The molecular and gene/miRNA expression profiles of radioiodine resistant papillary thyroid cancer

Background

Papillary thyroid cancer (PTC) is the most frequent endocrine tumor. Radioiodine (RAI) treatment is highly effective in these tumors, but up to 60% of metastatic cases become RAI-refractory. Scanty data are available on either the molecular pattern of radioiodine refractory papillary thyroid cancers (PTC) or the mechanisms responsible for RAI resistance.

Methods

We analyzed the molecular profile and gene/miRNA expression in primary PTCs, synchronous and RAI-refractory lymph node metastases (LNMs) in correlation to RAI avidity or refractoriness.

We classified patients as RAI+/D+ (RAI uptake/disease persistence), RAI−/D+ (absent RAI uptake/disease persistence), and RAI+/D- (RAI uptake/disease remission), and analyzed the molecular and gene/miRNA profiles, and the expression of thyroid differentiation (TD) related genes.

Results

A different molecular profile according to the RAI class was observed: BRAF^V600E cases were more frequent in RAI−/D+ (P = 0.032), and fusion genes in RAI+/D+ cases. RAI+/D- patients were less frequently pTERT mutations positive, and more frequently wild type for the tested mutations/fusions. Expression profiles clearly distinguished PTC from normal thyroid. On the other hand, in refractory cases (RAI+/D+ and RAI−/D+) no distinctive PTC expression patterns were associated with either tissue type, or RAI uptake, but with the driving lesion and BRAF−/RAS-like subtype. Primary tumors and RAI-refractory LNMs with BRAF^V600E mutation display transcriptome similarity suggesting that RAI minimally affects the expression profiles of RAI-refractory metastases. Molecular profiles associated with the expression of TPO, SLC26A4 and TD genes, that were found more downregulated in BRAF^V600E than in gene fusions tumors.

Conclusions

The present data indicate a different molecular profile in RAI-avid and RAI-refractory metastatic PTCs. Moreover, BRAF^V600E tumors displayed reduced differentiation and intrinsic RAI refractoriness, while PTCs with fusion oncogenes are RAI-avid but persistent, suggesting different oncogene-driven mechanisms leading to RAI refractoriness.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-020-01757-x.

Standardisation of thyroid cytology terminology and practice: are modifications necessary?-a narrative review

Universally accepted guidelines for diagnosis and management of any disease are desirable. Standardization of thyroid cytology reporting is aimed at guiding and improving clinical decision-making and management. However, socio-economic, and local factors and differences in disease prevalence and patterns require modification to suit local settings. 'One size fit all' approach is not possible for any disease diagnosis or management. The same concept is applicable in diagnosis and management of thyroid nodules. An additional special issue is the well-known high inter and intra-observer variability in the histological and cytological diagnosis of thyroid neoplasms. Despite this, thyroid cytology has a very significant influence in the management of thyroid diseases. An approach based on common principals with appropriate modifications that suits countries or continents is desirable and sustainable. The principals of TBSRTC have served as a framework for similar tiered classifications for reporting thyroid cytopathology. This article discusses globally available professional guidelines based on a common framework with appropriate modifications, with the universal aim of risk stratification of thyroid nodules.

Clinical and Translational Challenges in Thyroid Cancer

Thyroid cancer is the most common endocrine malignancy and it accounts for 1% of all newly diagnosed tumors. Approximately 10% of patients with differentiated thyroid carcinomas (DTC) and 30% with medullary thyroid carcinoma (MTC) could not be cured with locoregional treatment and could develop metastatic disease. In addition, one of the most aggressive solid tumors can arise from the thyroid gland, the anaplastic thyroid carcinoma, with a median overall survival of less than 6 months. Currently, only four drugs are approved for the treatment of DTC and MTC and several unmet needs are focusing the scientific discussions, including the resistant setting, the off-target side effects that may reduce the efficacy and the molecular knowledge-based combinations. In this review, we aimed to discuss the current molecular landscape and treatment of thyroid cancers, and the ongoing clinical and translational research lines focusing on new drugs and drug combinations to improve the inhibition of driver mutations, such as BRAF and RET, and how systemic therapies that improved outcomes of other cancer types, like immunotherapy and peptide receptor radionuclide therapy, may play a role in the future management of advanced thyroid cancers.

Vitamin D receptor knockdown attenuates the antiproliferative, pro‑apoptotic and anti‑invasive effect of vitamin D by activating the Wnt/β‑catenin signaling pathway in papillary thyroid cancer

Vitamin D and the vitamin D receptor (VDR) complex have been reported to inhibit the growth of several types of tumor; however, their function in papillary thyroid cancer (PCT) remains unknown. In addition, the Wnt/β-catenin signaling pathway was discovered to serve a critical role in the pathology of PCT. Therefore, the present study aimed to determine the role of the VDR and its association with Wnt/β-catenin signaling in vitamin D-treated PTC cells. VDR expression was detected in human PTC cells (including MDA-T120, MDA-T85, SNU-790 and IHH4 cells) and thyroid follicular cells (Nthy-ori 3–1 cells). SNU-790 and IHH4 cells were infected with KD-VDR or negative control (KD-NC) lentiviruses, treated with 1,25(OH)₂D3 (the active form of vitamin D), and subsequently referred to as the KD-VDR&vitD and KD-NC&vitD groups, respectively. Additionally, PTC cells infected with KD-NC and not treated with 1,25(OH)₂D3 were used as the normal control and referred to as the KD-NC group. VDR mRNA and protein expression levels were increased in MDA-T120, SNU-790 and MDA-T85 cells compared to Nthy-ori 3-1 cells, whereas in IHH4 cells, VDR mRNA and protein expression levels were similar to Nthy-ori 3-1 cells. In SNU-790 and IHH4 cells, cell proliferation and invasion were decreased in the KD-NC&vitD group compared with the KD-NC group, but increased in the KD-VDR&vitD group compared with the KD-NC&vitD group. Cell apoptosis was increased in the KD-NC&vitD group compared with the KD-NC group, and decreased in the KD-VDR&vitD group compared with the KD-NC&vitD group. Furthermore, the expression levels of Wnt family member 3 and catenin β1 were decreased in the KD-NC&vitD group compared with the KD-NC group, but increased in the KD-VDR&vitD group compared with the KD-NC&vitD group. In conclusion, the present study revealed that VDR-KD attenuated the antiproliferative, pro-apoptotic and anti-invasive effects of vitamin D in PTC by activating the Wnt/β-catenin signaling pathway.

DNA FISH Diagnostic Assay on Cytological Samples of Thyroid Follicular Neoplasms

Simple Summary

Cytopathology cannot distinguish benign from malignant follicular lesions in 20–30% of cases. These indeterminate cases includes the so-called follicular neoplasms (FNs) according to The Bethesda System for Reporting Thyroid Cytopathology. Frozen samples from 66 classic follicular adenomas (cFAs) and carcinomas (cFTCs) studied by array-comparative genomic hybridization identified three specific alterations of cFTCs (losses of 1p36.33-35.1 and 22q13.2-13.31, and gain of whole chromosome X) confirmed by fluorescent in situ hybridization (FISH) in a second independent series of 60 touch preparations from frozen samples of cFAs and cFTCs. In a third independent set of 27 cases of already stained pre-operative fine-needle aspiration cytology samples diagnosed as FNs and histologically verified, FISH analysis using these three markers identified half of cFTCs. Specificity of our assay for identifying cFTCs is higher than 98% which might be comparable with BRAF^600E testing in cases of suspicion of classic papillary thyroid carcinomas.

Abstract

Although fine-needle aspiration cytology (FNAC) is helpful in determining whether thyroid nodules are benign or malignant, this distinction remains a cytological challenge in follicular neoplasms. Identification of genomic alterations in cytological specimens with direct and routine techniques would therefore have great clinical value. A series of 153 cases consisting of 72 and 81 histopathologically confirmed classic follicular adenomas (cFAs) and classic follicular thyroid carcinomas (cFTCs), respectively, was studied by means of different molecular techniques in three different cohorts of patients (pts). In the first cohort (training set) of 66 pts, three specific alterations characterized by array comparative genomic hybridization (aCGH) were exclusively found in half of cFTCs. These structural abnormalities corresponded to losses of 1p36.33-35.1 and 22q13.2-13.31, and gain of whole chromosome X. The second independent cohort (validation set) of 60 pts confirmed these data on touch preparations of frozen follicular neoplasms by triple DNA fluorescent in situ hybridization using selected commercially available probes. The third cohort, consisting of 27 archived cytological samples from an equal number of pts that had been obtained for preoperative FNAC and morphologically classified as and histologically verified to be follicular neoplasms, confirmed our previous findings and showed the feasibility of the DNA FISH (DNA fluorescent in situ hybridization) assay. All together, these data suggest that our triple DNA FISH diagnostic assay may detect 50% of cFTCs with a specificity higher than 98% and be useful as a low-cost adjunct to cytomorphology to help further classify follicular neoplasms on already routinely stained cytological specimens.

Overexpression of novel long intergenic non‑coding RNA LINC02454 is associated with a poor prognosis in papillary thyroid cancer

It has been revealed from microarray data analysis that long intergenic non-coding RNA 02454 (LINC02454) is highly expressed in papillary thyroid cancer (PTC). The aim of the present study was to explore the potential role of LINC02454 in the tumorigenesis of PTC. The mRNA expression levels of LINC02454 were assessed using data from The Cancer Genome Atlas (TCGA) and the GSE66783 cohort in thyroid cancer, and were validated using reverse transcription-quantitative PCR in 104 patients with PTC recruited in the present study. The association between the LINC02454 mRNA expression levels and the clinicopathological features of the 104 patients with PTC were also analyzed. Functional enrichment analyses were conducted on the differentially expressed genes in the high and low LINC02454 expression groups that were identified from the TCGA cohort. RNA interference, using short interfering (si)RNA against LINC02454, was used to investigate the role of LINC02454 in the biological functions of PTC cells in vitro. The expression level of LINC02454 was significantly increased in PTC tissues (P=0.0011) and was significantly associated with a larger tumor size, T stage, an advanced TNM stage and an increased lymph node metastasis (P<0.05), which was consistent with that in the TCGA and GSE66783 cohort. High expression levels of LINC02454 were observed in patients with PTC that also had BRAF mutations (P<0.001), and were significantly associated with a poorer disease-free survival in the TCGA cohort (P<0.05). Functional enrichment analysis indicated that LINC02454-related genes were significantly enriched in Gene Ontology terms, such as ‘positive regulation of cell proliferation’, ‘positive regulation of cell division’ and ‘cell adhesion’, and the following Kyoto Encyclopedia of Genes and Genomes pathways: ‘Pathways in cancer’ ‘proteoglycans in cancer’ and ‘ECM-receptor interaction’. In vitro, the knockdown of LINC02454 markedly arrested the cells in the G₀/G₁ phase of the cell cycle, and also led to an overall increase in apoptosis, as well as to an unexpected decrease in cell proliferation. LINC02454 may thus potentially function as an oncogene, which inhibits the apoptosis and enhances proliferation of PTC cells. Thus, as suggested by the findings of the present study, LINC02454 may be used as a diagnostic and prognostic biomarker for PTC in the future.

Molecular trait of follicular-patterned thyroid neoplasms defined by MALDI-imaging

In the field of thyroid neoplasms, the most interesting recent change regards the introduction of a new terminology for follicular-patterned thyroid tumors, named Noninvasive Thyroid Neoplasm with Papillary-like Nuclear Features (NIFTP). This pre-malignant tumor is considered to be the putative precursor of invasive carcinoma. However, given that several issues are still unresolved, the application of ancillary tools, based on omics-techniques, may improve the clinical management of these challenging cases. The present paper highlights the proteomic profiles of a series of NIFTPs submitted to Fine Needle Aspirations (FNAs) and analysed by MALDI-imaging in order to confirm the heterogeneous phenotype of nodules included in the present NIFTP terminology and to underline the necessity of more accurate biomarkers that can be used for their characterization. Ethical and economic implications in terms of healthcare costs, operative risks, morbidity, as well as the potential need for lifelong hormone replacement therapy, seem to be significant reasons to approach the characterization of NIFTPs using alternative tools such as MALDI-MSI.

High Preponderance of BRAF V600E Mutation in Papillary Thyroid Carcinoma Among Filipinos: A Clinicopathologic Study

Purpose

BRAF mutation in papillary thyroid carcinoma (PTC) is associated with an aggressive phenotype, with varying incidence. We evaluated the prevalence of BRAF mutations in PTC among Filipino patients and their correlation with clinicopathologic characteristics.

Patients and Methods

Clinicopathologic data were retrieved from 64 sequential patients who underwent thyroidectomy from June 2016 to December 2016. BRAF mutation testing was performed using Sanger sequencing.

Results

Eighteen (28%) of 64 patients were diagnosed with PTC; 12 (70.59%) of 17 harbored a BRAF V600E mutation (no amplification in one patient). Demographics of patients with PTC were as follows: 13 women and five men, with median age of 46 years (range, 25 to 74 years). Fourteen patients had conventional subtype PTC; two, follicular variant; one, oncocytic variant; and one, tall-cell features. Tumor size ranged from 0.8 to 7.0 cm (median, 2.4 cm); extrathyroidal extension was present in seven (38.9%) of 18 patients, multifocality in six (33.33%) of eight, and lymph node involvement in eight (44.4%) of 18. Significant association between presence of a BRAF mutation and presence of extrathyroidal extension or lymph node involvement was not determined due to the limited sample size.

Conclusion

The high preponderance of BRAF mutation (70.59%) suggests some correlation with the previously reported lower 5-year survival among Filipinos. This warrants further investigation in a larger-cohort prospective study.

Aging-associated genes TNFRSF12A and CHI3L1 contribute to thyroid cancer: An evidence for the involvement of hypoxia as a driver

The prevalence of thyroid cancer (TC) is high in the elderly. The present study was based on the hypothesis that genes, which have increased activity with aging, may play a role in the development of TC. A large-scale literature-based data analysis was conducted to explore the genes that are implicated in both TC and aging. Subsequently, a mega-analysis of 16 RNA expression datasets (1,222 samples: 439 healthy controls, and 783 patients with TC) was conducted to test a set of genes associated with aging but not TC. To uncover a possible link between these genes and TC, a functional pathway analysis was conducted, and the results were validated by analysis of gene co-expression. A multiple linear regression (MLR) model was employed to study the possible influence of sample size, population region and study age on the gene expression levels in TC. A total of 262 and 816 genes were identified to have increased activity with aging and TC, respectively; with a significant overlap of 63 genes (P<3.82×10⁻³⁵). The mega-analysis revealed two aging-associated genes (CHI3L1 and TNFRSF12A) to be significantly associated with TC (P<2.05×10⁻⁸), and identified the association with multiple hypoxia-driven pathways through functional pathway analysis, also confirmed by the co-expression analysis. The MLR analysis identified population region as a significant factor contributing to the expression levels of CHI3L1 and TNFRSF12A in TC samples (P<3.24×10⁻⁴). The determination of genes that promote aging was warranted due to their possible involvement in TC. The present study suggests CHI3L1 and TNFRSF12A as novel common risk genes associated with both aging and TC.

BRAF and KRAS mutations in papillary thyroid carcinoma in the United Arab Emirates

Background

Papillary thyroid carcinoma (PTC) is the most common malignant thyroid neoplasm comprising 80–90% of all thyroid malignancies. Molecular changes in thyroid follicular cells are likely associated with the development of PTC. Mutations in serine/threonine-protein kinase (BRAF) and Rat sarcoma viral oncogene homolog (RAS) are commonly seen in PTC.

Methods

In total, 90 cases of PTC are randomly selected from archive paraffin blocks and 10μm sections were cut and processed for DNA extraction. BRAF^V600E mutation and 8 types of KRAS mutations were investigated using Real Time PCR.

Results

BRAF^V600E mutation was identified in 46% of PTC while KRAS mutations were seen in 11% of PTC. There was significant correlation between BRAF^V600E mutation and PTC larger than 5cm in diameter, positive surgical margin and lymph node metastasis. BRAF^V600E mutation was significantly higher in patients with less than 55-year of age than those more than 55-year of age. BRAF^V600E mutation was significantly higher in patients with family history of thyroid cancer than those without. There was no significant difference in BRAF^V600E mutation between males and females, PTC classic and follicular variants, unifocal and multifocal PTC. There was a significant higher percentage of BRAF^V600E mutation in classic PTC than papillary microcarcinoma variant. There was no significant age, gender, histologic type, tumor size, lymph node metastasis, tumor focality, and surgical margin status differences between KRAS mutated and non-mutated PTC.

Conclusion

BRAF^V600E and KRAS mutation are seen in a significant number of PTC in the UAE. BRAF mutation is significantly correlated with large tumor size, positive surgical margins and lymph node metastasis suggesting an association between BRAF^V600E mutation and tumor growth and spread.

Identification of hub genes in papillary thyroid carcinoma: robust rank aggregation and weighted gene co-expression network analysis

BACKGROUND

Papillary thyroid carcinoma (PTC), which is the most common endocrine malignancy, has been steadily increasing worldwide in incidence over the years, while mechanisms underlying the pathogenesis and diagnostic for PTC are incomplete. The purpose of this study is to identify potential biomarkers for diagnosis of PTC, and provide new insights into pathogenesis of PTC.

METHODS

Based on weighted gene co-expression network analysis, Robust Rank Aggregation, functional annotation, GSEA and DNA methylation, were employed for investigating potential biomarkers for diagnosis of PTC.

RESULTS

Black and turquoise modules were identified in the gene co-expression network constructed by 1807 DEGs that from 6 eligible gene expression profiles of Gene Expression Omnibus database based on Robust Rank Aggregation and weighted gene co-expression network analysis. Hub genes were significantly down-regulated and the expression levels of the hub genes were different in different stages in hub gene verification. ROC curves indicated all hub genes had good diagnostic value for PTC (except for ABCA6 AUC = 89.5%, the 15 genes with AUC > 90%). Methylation analysis showed that hub gene verification ABCA6, ACACB, RMDN1 and TFPI were identified as differentially methylated genes, and the decreased expression level of these genes may relate to abnormal DNA methylation. Moreover, the expression levels of 8 top hub genes were correlated with tumor purity and tumor-infiltrating immune cells. These findings, including functional annotations and GSEA provide new insights into pathogenesis of PTC.

CONCLUSIONS

The hub genes and methylation of hub genes may as potential biomarkers provide new insights for diagnosis of PTC, and all these findings may be the direction to study the mechanisms underlying of PTC in the future.

PTPN11 Knockdown Prevents Changes in the Expression of Genes Controlling Cell Cycle, Chemotherapy Resistance, and Oncogene-Induced Senescence in Human Thyroid Cells Overexpressing BRAF V600E Oncogenic Protein

The MAPK (RAS/BRAF/MEK/ERK) signaling pathway is a kinase cascade involved in the regulation of cell proliferation, differentiation, and survival in response to external stimuli. The V600E mutation in the BRAF gene has been detected in various tumors, resulting in a 500-fold increase in BRAF kinase activity. However, monotherapy with selective BRAF V600E inhibitors often leads to reactivation of MAPK signaling cascade and emergence of drug resistance. Therefore, new targets are being developed for the inhibition of components of the aberrantly activated cascade. It was recently discovered that resistance to BRAF V600E inhibitors may be associated with the activity of the tyrosine phosphatase SHP-2 encoded by the PTPN11 gene. In this paper, we analyzed transcriptional effects of PTPN11 gene knockdown and selective suppression of BRAF V600E in a model of thyroid follicular epithelium. We found that the siRNA-mediated knockdown of PTPN11 after vemurafenib treatment prevented an increase in the expression CCNA1 and NOTCH4 genes involved in the formation of drug resistance of tumors. On the other hand, downregulation of PTPN11 expression blocked the transcriptional activation of genes (p21, p15, p16, RB1, and IGFBP7) involved in cell cycle regulation and oncogene-induced senescence in response to BRAF V600E expression. Therefore, it can be assumed that SHP-2 participates not only in emergence of drug resistance in cancer cells, but also in oncogene-induced cell senescence.

Evaluation of CYFRA 21.1 as a Dedifferentiation Marker of Advanced Thyroid Cancer

PURPOSE OF THE STUDY

Well-differentiated thyroid carcinomas have good prognosis, but as it de-differentiates, the survival rates go down. Early identification of such patients needs a marker which indicates the dedifferentiation process. CYFRA 21.1 has also shown to be increased in patients with ¹³¹I refractory thyroid cancer. We tested whether CYFRA 21.1 can differentiate between ¹³¹I avid and refractory tumors.

METHODOLOGY

Well-differentiated thyroid cancer patients with known distant metastases were accrued and tested for stimulated and unstimulated thyroglobulin and CYFRA 21.1. All patients underwent ¹³¹I whole-body scan, ¹³¹I post therapy scan, and ¹⁸F-Fluorodeoxyglucose positron emission tomography-computed tomography. Those with even a single ¹³¹I nonavid lesion were considered ¹³¹I refractory disease. CYFRA 21.1 of both ¹³¹I avid and nonavid was compared, and CYFRA 21.1 levels against disease extent were analyzed.

RESULTS

CYFRA 21.1 levels were significantly elevated in ¹³¹I refractory group. A cutoff value of 2.07 ng/ml distinguished between ¹³¹I avid and refractory disease with high sensitivity and specificity (88% and 89. 7%, respectively). However, CYFRA 21.1 levels were similar in patients when analyzed based on disease sites.

CONCLUSION

CYFRA 21.1 can be utilized to differentiate between ¹³¹I avid and refractory diseases. Further long-term studies are required to use it as a predictive and prognostic marker.

Solid variant of papillary thyroid carcinoma: an under-recognized entity

Solid variant of papillary thyroid carcinoma (SVPTC) is a rare morphological variant of papillary thyroid carcinoma (PTC). SVPTC is histologically characterized by predominant solid, trabecular and insular nests of tumor cells while cytological features of PTC such as nuclear grooves and nuclear inclusions are preserved. In fine needle aspiration cytology smears, tumor cells of SVPTC may be presented in cohesive, syncytial or trabecular clusters accompanied by some discohesiveness in the absence of necrosis. Although SVPTC and poorly differentiated thyroid carcinoma (PDTC) share similar histological findings of solid nests, SVPTC can be differentiated from PDTC in the lack of tumor necrosis, severe nuclear atypia, and a higher mitotic index. Immunohistochemical expression of CK19 and HBME-1, common markers of PTC, is decreased in solid nests of SVPTC. In pediatric patients exposed to radiation after the Chernobyl nuclear accident, there was a higher prevalence of SVPTC with RET/PTC3 type rearrangement. BRAF mutations are also reported in a small number of adult patients with SVPTC without any prior radiation exposure. Patients with SVPTC may have a slightly higher incidence of metastasis and recurrence of the tumor compared to conventional PTC, although overall survival rate is comparable. In this article, the current knowledge of SVPTC will be reviewed and discussed with an emphasis on the histopathological feature.

Noninvasive Follicular Thyroid Neoplasm with Papillary-Like Nuclear Features (NIFTP): An Update

Based on evidence accumulated over the past three decades showing that noninvasive encapsulated follicular variant of papillary thyroid carcinoma has an indolent clinical behavior and a RAS-like molecular profile similar to follicular adenoma, the Endocrine Pathology Society working group in 2016 proposed to rename this entity as "noninvasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP)" in order to eliminate the term "carcinoma" from the diagnosis. It is a major evidence-based attempt initiated by an international group of endocrine pathologists to tackle the epidemic of thyroid cancer overdiagnosis and overtreatment. However, its creation and continuous existence are not without controversies. NIFTP has sparked a wave of follow up studies aiming to decipher the exact nature of this new entity. In this review, we summarize the rationale, diagnostic criteria, controversies and subsequent changes to the NIFTP concept, and their impact on patient care and pathology practice.

The Bethesda System for Reporting Thyroid Cytology (TBSRTC): From look-backs to look-ahead

The Bethesda System for Reporting Thyroid Cytopathology (TBSRTC) was formalized in October 2007 by experts in thyroidology at the National Institute of Health in Bethesda, Maryland. The first edition of the TBSRTC book was published in 2010 and the second edition in 2018. The TBSRTC is widely employed in cytology practices in the United States and has also served as a model for similar tiered classification schemes for reporting thyroid cytopathology specimens. The tremendous success of TBSRTC cannot be underscored, it has provided a diagnostic framework which is well aligned with the present and the future of thyroid nodule management.

Integrated bioinformatics analysis of the association between apolipoprotein E expression and patient prognosis in papillary thyroid carcinoma

The prognosis of most patients with papillary thyroid carcinoma (PTC) is excellent despite some cases of tumor progression or relapse. The present study was designed to reveal possible prognostic risk indicators for PTC. Differentially expressed genes (DEGs) extracted from 4 Gene Expression Omnibus (GEO) cohorts were subjected to functional enrichment analyses by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway analysis. A dataset from The Cancer Genome Atlas (TCGA) was obtained to filter and validate significant genes using cytoHubba, followed by analysis of their association with clinicopathological characteristics and prognosis. In total, 240 DEGs were identified after data preprocessing. These DEGs were enriched in ‘intracellular redox equilibrium’, ‘release of exosome’, ‘cell adhesion’, ‘regulation of extracellular matrix’, ‘collagen binding’ and ‘energy metabolism’ based on GO analysis which including cellular component, molecular function and biological process. KEGG pathway analysis revealed that the DEGs were enriched in thyroid hormone synthesis, pathways in cancer, focal adhesion, metabolic pathways, apoptosis, PPAR signaling pathway and PI3K/AKT signaling pathway. Using cytoHubba, the following hub genes were identified: Apolipoprotein E (APOE); hemoglobin subunit α1 (HBA1); angiotensin II receptor 1 (AGTR1); collagen I α1 (COL1A1); galectin 3 (LGALS3) and TIMP metallopeptidase inhibitor 1 (TIMP1). The expression of these genes was found to be consistent in TCGA datasets. Kaplan-Meier analysis revealed that APOE was significantly associated with overall survival (P=0.00067) and disease free survival (P=0.00220). Additionally, low expression of APOE was significantly associated with older age (P<0.001) and higher TNM stage (P<0.001) compared with the high expression group. Therefore, APOE may function as a predictive risk indicator for progression as well as prognosis of PTC.

Cancer Associated Fibroblasts and Senescent Thyroid Cells in the Invasive Front of Thyroid Carcinoma

Thyroid carcinoma (TC) comprises several histotypes with different aggressiveness, from well (papillary carcinoma, PTC) to less differentiated forms (poorly differentiated and anaplastic thyroid carcinoma, PDTC and ATC, respectively). Previous reports have suggested a functional role for cancer-associated fibroblasts (CAFs) or senescent TC cells in the progression of PTC. In this study, we investigated the presence of CAFs and senescent cells in proprietary human TCs including PTC, PDTC, and ATC. Screening for the driving lesions BRAFV600E and N/H/KRAS mutations, and gene fusions was also performed to correlate results with tumor genotype. In samples with unidentified drivers, transcriptomic profiles were used to establish a BRAF- or RAS-like molecular subtype based on a gene signature derived from The Cancer Genome Atlas. By using immunohistochemistry, we found co-occurrence of stromal CAFs and senescent TC cells at the tumor invasive front, where deposition of collagen (COL1A1) and expression of lysyl oxidase (LOX) enzyme were also detected, in association with features of local invasion. Concurrent high expression of CAFs and of the senescent TC cells markers, COL1A1 and LOX was confirmed in different TC histotypes in proprietary and public gene sets derived from Gene Expression Omnibus (GEO) repository, and especially in BRAF mutated or BRAF-like tumors. In this study, we show that CAFs and senescent TC cells co-occur in various histotypes of BRAF-driven thyroid tumors and localize at the tumor invasive front.

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.

Mouse Model of Thyroid Cancer Progression and Dedifferentiation Driven by STRN-ALK Expression and Loss of p53: Evidence for the Existence of Two Types of Poorly Differentiated Carcinoma

Background: Thyroid tumor progression from well-differentiated cancer to poorly differentiated thyroid carcinoma (PDTC) and anaplastic thyroid carcinoma (ATC) involves step-wise dedifferentiation associated with loss of iodine avidity and poor outcomes. ALK fusions, typically STRN-ALK, are found with higher incidence in human PDTC compared with well-differentiated cancer and, as previously shown, can drive the development of murine PDTC. The aim of this study was to evaluate thyroid cancer initiation and progression in mice with concomitant expression of STRN-ALK and inactivation of the tumor suppressor p53 (Trp53) in thyroid follicular cells. Methods: Transgenic mice with thyroid-specific expression of STRN-ALK and biallelic p53 loss were generated and aged on a regular diet or with methimazole and sodium perchlorate goitrogen treatment. Development and progression of thyroid tumors were monitored by using ultrasound imaging, followed by detailed histological and immunohistochemical evaluation. Gene expression analysis was performed on selected tumor samples by using RNA-Seq and quantitative RT-PCR. Results: In mice treated with goitrogen, the first thyroid cancers appeared at 6 months of age, reaching 86% penetrance by the age of 12 months, while a similar rate (71%) of tumor occurrence in mice on regular diet was observed by 18 months of age. Histological examination revealed well-differentiated papillary thyroid carcinomas (PTC) (n = 26), PDTC (n = 21), and ATC (n = 8) that frequently coexisted in the same thyroid gland. The tumors were frequently lethal and associated with the development of lung metastasis in 24% of cases. Histological and immunohistochemical characteristics of these cancers recapitulated tumors seen in humans. Detailed analysis of PDTC revealed two tumor types with distinct cell morphology and immunohistochemical characteristics, designated as PDTC type 1 (PDTC1) and type 2 (PDTC2). Gene expression analysis showed that PDTC1 tumors retained higher expression of thyroid differentiation genes including Tg and Slc5a5 (Nis) as compared with PDTC2 tumors. Conclusions: In this study, we generated a new mouse model of multistep thyroid cancer dedifferentiation with evidence of progression from PTC to PDTC and ATC. Further, PDTC in these mice showed two distinct histologic appearances correlated with levels of expression of thyroid differentiation and iodine metabolism genes, suggesting a possibility of existence of two PDTC types with different functional characteristics and potential implication for therapeutic approaches to these tumors.

Radioactive Iodine-Refractory Differentiated Thyroid Cancer and Redifferentiation Therapy

The retained functionality of the sodium iodide symporter (NIS) expressed in differentiated thyroid cancer (DTC) cells allows the further utilization of post-surgical radioactive iodine (RAI) therapy, which is an effective treatment for reducing the risk of recurrence, and even the mortality, of DTC. Whereas, the dedifferentiation of DTC could influence the expression of functional NIS, thereby reducing the efficacy of RAI therapy in advanced DTC. Genetic alternations (such as BRAF and the rearranged during transfection [RET]/papillary thyroid cancer [PTC] rearrangement) have been widely reported to be prominently responsible for the onset, progression, and dedifferentiation of PTC, mainly through activating the mitogen-activated protein kinase (MAPK) and phosphoinositide 3-kinase (PI3K) signaling cascades. These genetic alternations have been suggested to associate with the reduced expression of iodide-handling genes in thyroid cancer, especially the NIS gene, disabling iodine uptake and causing resistance to RAI therapy. Recently, novel and promising approaches aiming at various targets have been attempted to restore the expression of these iodine-metabolizing genes and enhance iodine uptake through in vitro studies and studies of RAI-refractory (RAIR)-DTC patients. In this review, we discuss the regulation of NIS, known mechanisms of dedifferentiation including the MAPK and PI3K pathways, and the current status of redifferentiation therapy for RAIR-DTC patients.

BRAFV600E mutation analysis in fine-needle aspiration cytology specimens for diagnosis of thyroid nodules: The influence of false-positive and false-negative results

Background

The accurate evaluation of BRAF^V600E mutation in preoperative fine‐needle aspiration cytology (FNAC) specimens is important for making management decisions in thyroid nodules (TNs). The aim of this study was to assess the false‐positive and false‐negative BRAF^V600E mutations in thyroid FNAC specimens and their influence on diagnosis of TN.

Methods

This prospective study enrolled 292 nodules in 269 patients who underwent BRAF^V600E mutation analysis using amplification refractory mutation system‐quantitative real‐time polymerase chain reaction (ARMS‐qPCR) both in FNAC specimens and formalin‐fixed, paraffin‐embedded (FFPE) tissue samples after surgery. The false‐positive and false‐negative mutations for BRAF^V600E analysis using ARMS‐qPCR in FNAC specimens were recorded, with reference to the results of BRAF^V600E mutation analysis using ARMS‐qPCR in FFPE tissue sample. Diagnostic performances of FNAC, BRAF^V600E mutation analysis in FNAC specimens, BRAF^V600E mutation analysis in FFPE tissue sample, and the combination of FNAC and BRAF^V600E mutation analysis for predicting thyroid malignancy were assessed.

Results

The false‐positive and false‐negative mutations for BRAF^V600E analysis using ARMS‐qPCR in FNAC specimens were 10.1% (19/189) and 7.1% (7/98), respectively. FNAC combined with preoperative BRAF^V600E mutation analysis significantly increased the diagnostic sensitivity from 75.7% to 92.3%, and accuracy from 78.7% to 90.6% in comparison with FNAC alone (both P < .001). No significant differences were found between the combination of FNAC and BRAF^V600E mutation analysis in FNAC specimens and the combination of FNAC and BRAF^V600E mutation analysis in FFPE tissue sample (sensitivity: 92.3% vs 91.9%; accuracy: 90.6% vs 91.3%; both P > .05).

Conclusions

FNAC combined with preoperative BRAF^V600E mutation analysis can significantly increase the diagnostic performance in comparison with FNAC alone. False‐positive and false‐negative BRAF^V600E mutation results are found in preoperative FNAC specimens, whereas it does not affect the overall auxiliary diagnosis of TNs.

miRNA-15a regulates the proliferation and apoptosis of papillary thyroid carcinoma via regulating AKT pathway

Aim

Aberrantly expressed microRNAs (miRNAs) are involved in many diseases including cancer. The expression of miR-15a was reported to be downregulated in papillary thyroid carcinoma (PTC) compared to control tissue. However, the mechanism underlying this downregulation remains unclear.

Methods

The effects of miR-15a on the proliferation and invasion of PTC cells were evaluated by CCK-8 and transwell assays, respectively. Expression levels of AKT and rearranged during transfection (RET) in cells were assessed using Western blotting. The correlation of RET and miR-15a was validated by luciferase reporter assay. Moreover, in vivo assay was performed to demonstrate the effect of miR-15a on tumor growth.

Results

We confirmed that the expression of miR-15a was significantly lower in PTC tissue than that in normal tissue. Overexpression of miR-15a notably inhibited PTC cell proliferation and invasion via promoting apoptosis. Additionally, RET was found to be a target of miR-15a and this correlation was confirmed by dual-luciferase assay and Western blot. Furthermore, in vivo study revealed that overexpression of miR-15a inhibited tumor growth via downregulating the levels of RET and phosphorylated AKT.

Conclusion

In the present study, we demonstrated that miR-15a played an antitumor role in regulating PTC via targeting RET/AKT pathway. Therefore, miR-15a may serve as a potential molecular target for the treatment of PTC.

Identification and interaction analysis of key miRNAs in medullary thyroid carcinoma by bioinformatics analysis

Medullary thyroid carcinoma (MTC) is an endocrine tumor and comprises 5–10% of all primary thyroid malignancies. However, the biomechanical contribution to the development and progression of MTC remains unclear. In this study, To discover the key microRNAs (miRNAs or miRs) and their potential roles in the tumorigenesis of MTC, the microarray datasets GSE97070, GSE40807 and GSE27155 were analyzed. The datasets were downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed miRNAs (DEMs) and genes (DEGs) were accessed by R. Targets of DEMs and predicted using starBase, and functional and pathway enrichment analyses were performed using Metascape. A protein-protein interaction (PPI) network and an analysis of modules were constructed using NetworkAnalyst. Finally, a network was constructed to show the regulatory association between transcription factors (TFs), DEMs and downstream genes. A total of 5 DEMs were found both in GSE97070 and GSE40807, including 3 upregulated DEMs and 2 downregulated DEMs. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses from Metascape revealed that the target genes of upregulated DEMs were significantly enriched in adherens junction, kinase and protein binding, while the target genes of downregulated DEMs were mainly involved in non-canonical Wnt signaling pathway and RNA transport. From the PPI network, 13 nodes were screened as hub genes. Pathway enrichment analysis revealed that the top 5 modules were mostly enriched in the neurotrophin signaling pathway, mRNA surveillance pathway and MAPK signaling pathway. In addition, the TF-DEMs-target gene and DEGs regulatory network revealed that 17 TFs regulated 2 miRNAs, including upregulated or downregulated DEMs, CREB1 regulated all upregulated DEMs, and TGFB1 was an activator of hsa-miR-199a-3p and a repressor of hsa-miR-429. Taken together, the present study identified several miRNAs and potential biological mechanisms involved in the tumorigenesis of MTC. This study identified the key DEMs and potential mechanisms underlying the development of MTC, and provided a series of biomarkers and targets for the management of MTC.

Gene expression profiling identifies potential molecular markers of papillary thyroid carcinoma

BACKGROUND

Thyroid cancer is the most common endocrine malignancy worldwide, with the predominant form papillary thyroid carcinoma (PTC) representing approximately 80% of cases.

OBJECTIVE

This study was addressed to identify potential genes and pathways involved in the pathogenesis of PTC and potential novel biomarkers for this disease.

METHODS

Gene expression profiling was carried out by DNA microarray technology. Validation of microarray data by qRT-PCR, western blot, and enzyme linked immunosorbent assay was also performed in a selected set of genes and gene products, with the potential to be used as diagnostic or prognostic biomarkers, such as those associated with cell adhesion, extracellular matrix (ECM) remodeling and immune/inflammatory response.

RESULTS

In this study we found that upregulation of extracellular activities, such as proteoglycans, ECM-receptor interaction, and cell adhesion molecules, were the most prominent feature of PTC. Significantly over-expressed genes included SDC1 (syndecan 1), SDC4 (syndecan 4), KLK7 (kallikrein-related peptidase 7), KLK10 (kallikrein-related peptidase 10), SLPI (secretory leukocyte peptidase inhibitor), GDF15 (growth/differentiation factor-15), ALOX5 (arachidonate 5-lipoxygenase), SFRP2 (secreted Frizzled-related protein 2), among others. Further, elevated KLK10 levels were detected in patients with PTC. Many of these genes belong to KEGG pathway "Proteoglycans in cancer".

CONCLUSIONS

Using DNA microarray analysis allowed the identification of genes and pathways with known important roles in malignant transformation, and also the discovery of novel genes that may be potential biomarkers for PTC.

Molecular Signature of Prospero Homeobox 1 (PROX1) in Follicular Thyroid Carcinoma Cells

The prospero homeobox 1 (PROX1) transcription factor is a product of one of the lymphangiogenesis master genes. It has also been suggested to play a role in carcinogenesis, although its precise role in tumour development and metastasis remains unclear. The aim of this study was to gain more knowledge on the PROX1 function in thyroid tumorigenesis. Follicular thyroid cancer-derived cells—CGTH-W-1—were transfected with PROX1-siRNA (small interfering RNA) and their proliferation, cell cycle, apoptosis and motility were then analysed. The transcriptional signature of PROX1 depletion was determined using RNA-Sequencing (RNA-Seq) and the expression of relevant genes was further validated using reverse transcriptase quantitative PCR (RT-qPCR), Western blot and immunocytochemistry. PROX1 depletion resulted in a decreased cell motility, with both migratory and invasive potential being significantly reduced. The cell morphology was also affected, while the other studied cancer-related cell characteristics were not significantly altered. RNA-seq analysis revealed significant changes in the expression of transcripts encoding genes involved in both motility and cytoskeleton organization. Our transcriptional analysis of PROX1-depleted follicular thyroid carcinoma cells followed by functional and phenotypical analyses provide, for the first time, evidence that PROX1 plays an important role in the metastasis of thyroid cancer cells by regulating genes involved in focal adhesion and cytoskeleton organization in tumour cells.

The Current Histologic Classification of Thyroid Cancer

Thyroid cancers of follicular cell derivation provide excellent phenotype-genotype correlations. Current morphologic classifications are complex and require simplification. Benign adenomas have follicular or papillary architecture and bland cytology. Well-differentiated thyroid carcinomas exhibit follicular architecture, expansile growth, and variable cytologic atypia and invasiveness; low-risk tumors have excellent prognosis after surgical resection whereas widely-invasive and angioinvasive tumors warrant total thyroidectomy and radioablation. Papillary carcinoma is less differentiated; indolent microcarcinomas can be managed by active surveillance, whereas clinical lesions with local or distant spread require therapy. Progression gives rise to poorly differentiated and anaplastic carcinomas that are less common but far more aggressive.

Low dose radiation regulates BRAF-induced thyroid cellular dysfunction and transformation

Background

The existence of differentiated thyroid cells is critical to respond radioactive iodide treatment strategy in thyroid cancer, and loss of the differentiated phenotype is a trademark of iodide-refractive thyroid disease. While high-dose therapy has been beneficial to several cancer patients, many studies have indicated this clinical benefit was limited to patients having BRAF mutation. BRAF-targeted paired box gene-8 (PAX8), a thyroid-specific transcription factor, generally dysregulated in BRAF-mutated thyroid cancer.

Methods

In this study, thyroid iodine-metabolizing gene levels were detected in BRAF-transformed thyroid cells after low and high dose of ionizing radiation. Also, an mRNA-targeted approach was used to figure out the underlying mechanism of low (0.01Gyx10 or 0.1Gy) and high (2Gy) radiation function on thyroid cancer cells after BRAF^V600E mutation.

Results

Low dose radiation (LDR)-induced PAX8 upregulation restores not only BRAF-suppressive sodium/iodide symporter (NIS) expression, one of the major protein necessary for iodine uptake in healthy thyroid, on plasma membrane but also regulate other thyroid metabolizing genes levels. Importantly, LDR-induced PAX8 results in decreased cellular transformation in BRAF-mutated thyroid cells.

Conclusion

The present findings provide evidence that LDR-induced PAX8 acts as an important regulator for suppression of thyroid carcinogenesis through novel STAT3/miR-330-5p pathway in thyroid cancers.

Electronic supplementary material

The online version of this article (10.1186/s12964-019-0322-x) contains supplementary material, which is available to authorized users.

Podoplanin (PDPN) affects the invasiveness of thyroid carcinoma cells by inducing ezrin, radixin and moesin (E/R/M) phosphorylation in association with matrix metalloproteinases

Background

Podoplanin (PDPN) is a mucin-type transmembrane glycoprotein specific to the lymphatic system. PDPN expression has been found in various human tumors and is considered to be a marker of cancer. We had previously shown that PDPN expression contributes to carcinogenesis in the TPC1 papillary thyroid cancer-derived cell line by enhancing cell migration and invasiveness. The aim of this study was to determine the effect of PDPN down-regulation in another thyroid cancer-derived cell line: BcPAP.

Methods

In order to determine the effects of PDPN on malignant features of BcPAP cells (harboring the BRAFV600E mutated allele) and TPC1 cells (carrying the RET/PTC1 rearrangement), we silenced PDPN in these cells using small interfering RNA (siRNA). The efficacy of PDPN silencing was confirmed by qRT-PCR and Western blotting. Then, we tested the motility and invasiveness of these cells (using scratch test and Transwell assay), their growth capacities F(cell cycle analysis, viability, clonogenic activity) and apoptosis assays), adhesion-independent colony-formation capacities, as well as the effect of PDPN silencing on MMPs expression and activity (zymography).

Results

We found that PDPN-induced cell phenotype depended on the genetic background of thyroid tumor cells. PDPN down-regulation in BcPAP cells was negatively correlated with the migration and invasion, in contrast to TPC1 cells in which PDPN depletion resulted in enhanced migration and invasiveness. Moreover, our results suggest that in BcPAP cells, PDPN may be involved in the epithelial-mesenchymal transition (EMT) through regulating the expression of the ezrin, radixin and moesin (E/R/M) proteins, MMPs 9 and MMP2, remodeling of actin cytoskeleton and cellular protrusions. We also demonstrated that PDPN expression is associated with the MAPK signaling pathway. The inhibition of the MAPK pathway resulted in a decreased PDPN expression, increased E/R/M phosphorylation and reduced cell migration. Additionally, PDPN depleted BcPAP cells treated with inhibitors of MEK1/2 kinases (U0126) or of the BRAF V600E protein (PLX4720) had reduced motility, similar to that previously observed in TPC1 cells after PDPN knock-down.

Conclusions

Altogether, our data suggest that PDPN may play an important role in the control of invasion and migration of papillary thyroid carcinoma cells in association with the E/R/M, MMPs and MAPK kinases.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-5239-z) contains supplementary material, which is available to authorized users.

Molecular characterization of tumors meeting diagnostic criteria for the non-invasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP)

"Follicular variant" papillary thyroid carcinomas (FV-PTC) that do not histologically invade have a miniscule risk of metastasis, and thus been reclassified as a tumor of low malignant potential, the non-invasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP). There are few molecular studies of this tumor type. We performed gene expression analysis, by RNA sequencing, on a series of FV-PTCs, NIFTPs, and follicular adenomas. A training set comprised tumors from The Cancer Genome Atlas (TCGA) repository (n = 46), digital slides from which were reviewed and classified as invasive or non-invasive FV-PTC. A validation set comprised in-house NIFTPs, invasive FV-PTCs, and follicular adenomas (n = 26). In the training set, unsupervised clustering separated tumors into three distinct expression subtypes, which associated with invasion and characteristic molecular alterations. Specifically, the "BRAF-like" subtype was enriched in invasive FV-PTCs and tumors with BRAF V600E mutations. The "THADA-like" subtype was enriched in non-invasive tumors and those with rearrangements involving THADA. The "RAS-family-like" subtype included many invasive and non-invasive FV-PTCs and was enriched in tumors with mutations in RAS family genes. In the validation set, nearest centroid analysis classified all invasive FV-PTCs as "BRAF-like" and all follicular adenomas as either "RAS-like" or "THADA-like." NIFTPs were the most molecularly diverse histologic type, with cases classified as "BRAF-like," "THADA-like," and "RAS-family-like." In conclusion, tumors fitting criteria for NIFTP are molecularly diverse, making it difficult to diagnose them with molecular studies, likely including matrial from cytopathology samples.