Impact of the Tumor Microenvironment on the Gene Expression Profile in Papillary Thyroid Cancer.

Transcriptome of papillary thyroid cancer (PTC) is well characterized and correlates with some prognostic and genotypic factors, but data addressing the interaction between PTC and tumor microenvironment (TME) are scarce. Therefore, in the present study, we aimed to assess the impact of TME on gene expression profile in PTC. We evaluated the gene expression profile in PTC and normal thyroid cells isolated by laser capture microdissection and in whole tissue slides corresponding to the entire tumor. We included 26 microdissected samples for gene expression analysis (HG-U133 PLUS 2.0, Affymetrix, currently Thermo Fisher Scientific USA): 15 PTC samples, 11 samples of normal thyrocytes, and 30 whole slides (15 PTC and 15 normal thyroid). Transcripts were divided into three groups: differentially expressed both in microdissected and whole slides, transcripts differently expressed in microdissected samples and not changed in whole slides, and transcripts differentially expressed in whole slides and not changed in microdissected samples. Eleven genes were selected for validation in an independent set of samples; among them, four genes differentiated only microdissected PTC and normal cells. Two genes (PTCSC and CTGF) were confirmed. One gene (FOS) was not confirmed by the validation, whereas EGR1 was also significant in whole slide analysis. The other seven genes (TFF3, FN1, MPPED2, MET, KCNJ2, TACSTD2, and GALE) showed differentiated expression in microdissected thyrocytes and in whole tumor slides. Most of identified genes were related to the tumor-microenvironment interaction and confirmed the crosstalk between TME and cancer cells.

Differences in Gene Expression Profile of Primary Tumors in Metastatic and Non-Metastatic Papillary Thyroid Carcinoma-Do They Exist?

Molecular mechanisms of distant metastases (M1) in papillary thyroid cancer (PTC) are poorly understood. We attempted to analyze the gene expression profile in PTC primary tumors to seek the genes associated with M1 status and characterize their molecular function. One hundred and twenty-three patients, including 36 M1 cases, were subjected to transcriptome oligonucleotide microarray analyses: (set A-U133, set B-HG 1.0 ST) at transcript and gene group level (limma, gene set enrichment analysis (GSEA)). An additional independent set of 63 PTCs, including 9 M1 cases, was used to validate results by qPCR. The analysis on dataset A detected eleven transcripts showing significant differences in expression between metastatic and non-metastatic PTC. These genes were validated on microarray dataset B. The differential expression was positively confirmed for only two genes: IGFBP3, (most significant) and ECM1. However, when analyzed on an independent dataset by qPCR, the IGFBP3 gene showed no differences in expression. Gene group analysis showed differences mainly among immune-related transcripts, indicating the potential influence of tumor immune infiltration or signal within the primary tumor. The differences in gene expression profile between metastatic and non-metastatic PTC, if they exist, are subtle and potentially detectable only in large datasets.

Novel TG-FGFR1 and TRIM33-NTRK1 transcript fusions in papillary thyroid carcinoma.

Papillary thyroid carcinoma (PTC) is most common among all thyroid cancers. Multiple genomic alterations occur in PTC, and gene rearrangements are one of them. Here we screened 14 tumors for novel fusion transcripts by RNA-Seq. Two samples harboring RET/PTC1 and RET/PTC3 rearrangements were positive controls whereas the remaining ones were negative regarding the common PTC alterations. We used Sanger sequencing to validate potential fusions. We detected 2 novel potentially oncogenic transcript fusions: TG-FGFR1 and TRIM33-NTRK1. We detected 4 novel fusion transcripts of unknown significance accompanying the TRIM33-NTRK1 fusion: ZSWIM5-TP53BP2, TAF4B-WDR1, ABI2-MTA3, and ARID1B-PSMA1. Apart from confirming the presence of RET/PTC1 and RET/PTC3 in positive control samples, we also detected known oncogenic fusion transcripts in remaining samples: TFG-NTRK1, ETV6-NTRK3, MKRN1-BRAF, EML4-ALK, and novel isoform of CCDC6-RET.

Heterogeneity of Thyroid Cancer.

There are 5 main histological types of thyroid cancers (TCs): papillary, follicular (also known as differentiated), poorly differentiated, anaplastic (the most aggressive form), and medullary TC, and only the latter arises from thyroid C cells. These different forms of TCs show significant variability, both among and within tumours. This great variation is particularly notable among the first 4 types, which all originate from thyroid follicular cells. Importantly, this heterogeneity is not limited to histopathological diversity only but is also manifested as variation in several genetic and/or epigenetic alterations, the numbers of interactions between the tumour and surrounding microenvironment, and interpatient differences, for example. All these factors contribute to the great complexity in the development of a tumour from cancer cells. In the present review, we summarise the knowledge accumulated about the heterogeneity of TCs. Further research in this direction should help to gain a better understanding of the underlying mechanisms contributing to the development and diversity of TCs, paving the way toward more effective treatment strategies.

Coexistence of TERT Promoter Mutations and the BRAF V600E Alteration and Its Impact on Histopathological Features of Papillary Thyroid Carcinoma in a Selected Series of Polish Patients.

TERT promoter (TERTp) mutations are important factors in papillary thyroid carcinomas (PTCs). They are associated with tumor aggressiveness, recurrence, and disease-specific mortality and their use in risk stratification of PTC patients has been proposed. In this study we investigated the prevalence of TERTp mutations in a cohort of Polish patients with PTCs and the association of these mutations with histopathological factors, particularly in coexistence with the BRAF V600E mutation. A total of 189 consecutive PTC specimens with known BRAF mutational status were evaluated. TERTp mutations were detected in 8.5% of cases (16/189) with the C228T mutation being the most frequent. In six of the PTC specimens (3.2%), four additional TERTp alterations were found, which included one known polymorphism (rs2735943) and three previously unreported alterations. The association analysis revealed that the TERTp hotspot mutations were highly correlated with the presence of the BRAF V600E mutation and their coexistence was significantly associated with gender, advanced patient age, advanced disease stage, presence of lymph node metastases, larger tumor size, and tumor-capsule infiltration. While correlations were identified, the possibility of TERTp mutations being key molecular modulators responsible for PTC aggressiveness requires further studies.

Current Advances in Thyroid Cancer Management. Are We Ready for the Epidemic Rise of Diagnoses?

A rising incidence of thyroid cancers (TCs) mainly small tumors, observed during recent years, lead to many controversies regarding treatment strategies. TCs represent a distinct molecular background and clinical outcome. Although in most cases TCs are characterized by a good prognosis, there are some aggressive forms, which do not respond to standard treatment. There are still some questions, which have to be resolved to avoid dangerous simplifications in the clinical management. In this article, we focused on the current advantages in preoperative molecular diagnostic tests and histopathological examination including noninvasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTP). We discussed the controversies regarding the extent of thyroid surgery and adjuvant radioiodine therapy, as well as new treatment modalities for radioiodine-refractory differentiated thyroid cancer (RR-DTC). Considering medullary thyroid cancer (MTC), we analyzed a clinical management based on histopathology and RET (ret proto-oncogene) mutation genotype, disease follow-up with a special attention to serum calcitonin doubling time as an important prognostic marker, and targeted therapy applied in advanced MTC. In addition, we provided some data regarding anaplastic thyroid cancer (ATC), a highly lethal neoplasm, which lead to death in nearly 100% of patients due to the lack of effective treatment options.

Gene Expression (mRNA) Markers for Differentiating between Malignant and Benign Follicular Thyroid Tumours.

Distinguishing between follicular thyroid cancer (FTC) and follicular thyroid adenoma (FTA) constitutes a long-standing diagnostic problem resulting in equivocal histopathological diagnoses. There is therefore a need for additional molecular markers. To identify molecular differences between FTC and FTA, we analyzed the gene expression microarray data of 52 follicular neoplasms. We also performed a meta-analysis involving 14 studies employing high throughput methods (365 follicular neoplasms analyzed). Based on these two analyses, we selected 18 genes differentially expressed between FTA and FTC. We validated them by quantitative real-time polymerase chain reaction (qRT-PCR) in an independent set of 71 follicular neoplasms from formaldehyde-fixed paraffin embedded (FFPE) tissue material. We confirmed differential expression for 7 genes (CPQ, PLVAP, TFF3, ACVRL1, ZFYVE21, FAM189A2, and CLEC3B). Finally, we created a classifier that distinguished between FTC and FTA with an accuracy of 78%, sensitivity of 76%, and specificity of 80%, based on the expression of 4 genes (CPQ, PLVAP, TFF3, ACVRL1). In our study, we have demonstrated that meta-analysis is a valuable method for selecting possible molecular markers. Based on our results, we conclude that there might exist a plausible limit of gene classifier accuracy of approximately 80%, when follicular tumors are discriminated based on formalin-fixed postoperative material.

Gene signature of the post-Chernobyl papillary thyroid cancer.

PURPOSE: Following the nuclear accidents in Chernobyl and later in Fukushima, the nuclear community has been faced with important issues concerning how to search for and diagnose biological consequences of low-dose internal radiation contamination. Although after the Chernobyl accident an increase in childhood papillary thyroid cancer (PTC) was observed, it is still not clear whether the molecular biology of PTCs associated with low-dose radiation exposure differs from that of sporadic PTC. METHODS: We investigated tissue samples from 65 children/young adults with PTC using DNA microarray (Affymetrix, Human Genome U133 2.0 Plus) with the aim of identifying molecular differences between radiation-induced (exposed to Chernobyl radiation, ECR) and sporadic PTC. All participants were resident in the same region so that confounding factors related to genetics or environment were minimized. RESULTS: There were small but significant differences in the gene expression profiles between ECR and non-ECR PTC (global test, p < 0.01), with 300 differently expressed probe sets (p < 0.001) corresponding to 239 genes. Multifactorial analysis of variance showed that besides radiation exposure history, the BRAF mutation exhibited independent effects on the PTC expression profile; the histological subset and patient age at diagnosis had negligible effects. Ten genes (PPME1, HDAC11, SOCS7, CIC, THRA, ERBB2, PPP1R9A, HDGF, RAD51AP1, and CDK1) from the 19 investigated with quantitative RT-PCR were confirmed as being associated with radiation exposure in an independent, validation set of samples. CONCLUSION: Significant, but subtle, differences in gene expression in the post-Chernobyl PTC are associated with previous low-dose radiation exposure.

The frequency of polymorphic variants of filaggrin gene and clinical atopic dermatitis.

INTRODUCTION: As far as pathogenesis of the atopic dermatitis (AD) is concerned, the roles of an impaired epidermal barrier and cornified cell envelope are widely emphasized. AIM: The assessment of mutations of the filaggrin gene and their connection with the clinical picture of AD as well as selected allergological and environmental indicators. MATERIAL AND METHODS: 105 patients with diagnosed AD on the basis of diagnostic criteria were included. For every patient of the examined group, quantitative determination of the total concentration of IgE and the concentration of IgE antibodies to selected allergens were examined. For all patients, studies were performed by means of analysis of two genomic gene variants of profilaggrin (FLG) - R501X and 2282del4. RESULTS: Loss-of-function mutations in the filaggrin gene were shown in 12 (11.4%) patients in the examined group. All patients in the study group who developed one of the tested loss-of-function mutations in the filaggrin gene demonstrated an extrinsic, allergic form of atopic dermatitis. A significant association (p = 0.0002) between the presence of one of the tested loss-of-function mutations in the filaggrin gene and elevated levels of total concentration of immunoglobulin E was shown. CONCLUSIONS: Patients with AD of null mutations in the filaggrin gene demonstrate a relationship with the total and specific concentration of immunoglobulin E, specifically higher concentrations of IgE against aeroallergens and alimentary allergens as well as elevated levels of total immunoglobulin E.

Impact of SNPs on methylation readouts by Illumina Infinium HumanMethylation450 BeadChip Array: implications for comparative population studies.

BACKGROUND: Infinium HumanMethylation 450 BeadChip Arrays by Illumina (Illumina HM450K) are among the most popular CpG microarray platforms widely used in biological and medical research. Several recent studies highlighted the potentially confounding impact of the genomic variation on the results of methylation studies performed using Illumina's Infinium methylation probes. However, the complexity of SNPs impact on the methylation level measurements (ß values) has not been comprehensively described. RESULTS: In our comparative study of European and Asian populations performed using Illumina HM450K, we found that the majority of Infinium probes, which differentiated two examined groups, had SNPs in their target sequence. Characteristic tri-modal or bi-modal patterns of ß values distribution among individual samples were observed for CpGs with SNPs in the first and second position, respectively. To better understand how SNPs affect methylation readouts, we investigated their impact in the context of SNP position and type, and of the Illumina probe type (Infinium I or II). CONCLUSIONS: Our study clearly demonstrates that SNP variation existing in the genome, if not accounted for, may lead to false interpretation of the methylation signal differences suggested by some of the Illumina Infinium probes. In addition, it provides important practical clues for discriminating between differences due to the methylation status and to the genomic polymorphisms, based on the inspection of methylation readouts in individual samples. This approach is of special importance when Illumina Infinium assay is used for any comparative population studies, whether related to cancer, disease, ethnicity where SNP frequencies differentiate the studied groups.

Shared and unique metabolic features of the malignant and benign thyroid lesions determined with use of 1H HR MAS NMR spectroscopy.

The purpose of this work was to investigate the distinct and common metabolic features of the malignant and benign thyroid lesions in reference to the non-transformed tissue from the contralateral gland (chronic thyroiditis and colloid goiter). 1H HR MAS NMR spectra of 38 malignant lesions, 32 benign lesions and 112 samples from the non-tumoral tissue (32 from chronic thyroiditis and 80 samples from colloid goiter) were subjected both to multivariate and univariate analysis. The increased succinate, glutamine, glutathione, serine/cysteine, ascorbate, lactate, taurine, threonine, glycine, phosphocholine/glycerophosphocholine and decreased lipids were found in both lesion types in comparison to either colloid goiter or chronic thyroiditis. The elevated glutamate and choline, and reduced citrate and glucose were additionally evident in these lesions in reference to goiter, while the increased myo-inositol-in comparison to thyroiditis. The malignant lesions were characterized by the higher alanine and lysine levels than colloid goiter and thyroiditis, while scyllo-inositol was uniquely increased in the benign lesions (not in cancer) in comparison to both non-tumoral tissue types. Moreover, the benign lesions presented with the unique increase of choline in reference to thyroiditis (not observed in the cancerous tissue). The metabolic heterogeneity of the non-tumoral tissue should be considered in the analysis of metabolic reprogramming in the thyroid lesions.

Occurrence of BRAF mutations in a Polish cohort of PTC patients - preliminary results.

INTRODUCTION: Genetic alterations involving the mitogen-activated protein kinase (MAPK) pathway are frequently demonstrated in papillary thyroid cancer (PTC). BRAF(V600E), the most frequent mutation in adult patients, is present in approximately 50% of PTC. Most clinical studies have demonstrated an association of BRAF(V600E) mutation with aggressive clinicopathological characteristics and high tumour recurrence, although the results are controversial. In this study we present the preliminary results of BRAF mutation frequence in a group of 88 Polish patients with papillary thyroid cancer (PTC) and relate it to the outcome all DTC patients operated in 2004 and 2005. BRAF (V600E) mutation was diagnosed in 38 (43%) of cases. MATERIAL AND METHODS: The presence of BRAF mutation was evaluated in 88 PTC tumours. DNA was isolated from tissue parafin blocks, and the mutation V600E was evaluated by sequence analysis with an AbiPrism 377 and 3130 xl genetic analyzer (Life Technologies). Statistical analysis was carried out with the use of SPSS 12 software. The chi² and Kaplan-Meyer survival analysis were performed. RESULTS: From all analyzed clinico-pathological factors, only older age positively correlated with BRAF mutation frequency (p = 0.0017). Lymph node/distant metastases, multifocality, and extra-thyroid extension did not correlate with BRAF status. One cancer related death and two reccurences were observed in the BRAF+ group while one relapse was diagnosed in the BRAF- group. CONCLUSIONS: Although many studies document BRAF mutation as a prognostic factor in PTC our results underline that it is too early to consider it as a routine clinical predictive factor.

European perspective on the use of molecular tests in the diagnosis and therapy of thyroid neoplasms.

Thyroid nodules are frequently observed, particularly in individuals of over 60 years of age. On the other hand, most of the detected changes are benign and they do not require surgery. Therefore, differentiation between benign and malignant lesions in preoperative diagnosis is of crucial importance. Currently, the use of fine-needle aspiration biopsy (FNAB) and cytological assessment are the gold standard in the diagnosis of thyroid nodules. This procedure significantly reduces the need for diagnostic surgical intervention. However, approximately 15-30% of cytological results are classified as indeterminate. This is mainly due to the lack of specific cytomorphologic features that would facilitate the diagnosis based on cell evaluation under microscopic assessment. For the diagnoses of atypia of undetermined significance/follicular lesion of undetermined significance (AUS/FLUS), the assessment of invasion is crucial. Such an evaluation is not possible in cytology. Recently, molecular tests have been developed. They improve cytological diagnosis, particularly in the case of indeterminate results. Commercially available tests are developed based on the North American population. It is important to assess whether such tests can be used in the evaluation of e.g., European population.

TERT Promoter Mutations and Their Impact on Gene Expression Profile in Papillary Thyroid Carcinoma.

BACKGROUND: Telomerase reverse transcriptase promoter (TERTp) mutations are related to a worse prognosis in various malignancies, including papillary thyroid carcinoma (PTC). Since mechanisms responsible for the poorer outcome of TERTp(+) patients are still unknown, searching for molecular consequences of TERTp mutations in PTC was the aim of our study. METHODS: The studied cohort consisted of 54 PTCs, among them 24 cases with distant metastases. BRAF V600E, RAS, and TERTp mutational status was evaluated in all cases. Differences in gene expression profile between TERTp(+) and TERTp(-) PTCs were examined using microarrays. The evaluation of signaling pathways and gene ontology was based on the Gene Set Enrichment Analysis. RESULTS: Fifty-nine percent (32/54) of analyzed PTCs were positive for at least one mutation: 27 were BRAF(+), among them eight were TERTp(+), and 1 NRAS(+), whereas five other samples harbored RAS mutations. Expression of four genes significantly differed in BRAF(+)TERTp(+) and BRAF(+)TERTp(-) PTCs. Deregulation of pathways involved in key cell processes was observed. CONCLUSIONS: TERTp mutations are related to higher PTC aggressiveness. CRABP2 gene was validated as associated with TERTp mutations. However, its potential use in diagnostics or risk stratification in PTC patients needs further studies.

Age at diagnosis and gender modify the risk of 9q22 and 14q13 polymorphisms for papillary thyroid carcinoma.

INTRODUCTION: Papillary thyroid cancer (PTC) shows familial occurrence, and some susceptibility single nucleotide polymorphisms (SNPs) have been identified in FOXE1 and near the NKX2-1 locus. The aim of our study was to analyse the association of PTC risk with SNPs in FOXE1 (rs965513, rs1867277, rs1443434) and near the NKX2-1 locus (rs944289) in a Polish population, and, in the second step, the interac-tion between SNPs and patient-related factors (age at diagnosis and gender). MATERIAL AND METHODS: A total of 2243 DNA samples from PTC patients and 1160 controls were included in the study. The SNP analysis was performed with the allelic discrimination technique. RESULTS: There were significant associations of all SNPs with PTC (rs965513 odds ratio [OR] = 1.72, p = 8 × 10-7; rs1867277 OR = 1.59, p = 1 × 10-6; rs1443434 OR = 1.53, p = 1 × 10-5; rs944289 OR = 1.52, p = 4 × 10-5). Logistic regression analysis revealed an increased PTC risk in the interaction of rs944289 with age at diagnosis (OR = 1.01 per year, p = 6 × 10-4) and a decreased PTC risk in the interaction of male gender with the GGT FOXE1 protective haplotype (OR = 0.69, p = 0.01). CONCLUSIONS: the association between PTC and all analysed SNPs was confirmed. It was also shown that patient-related factors modify the predisposition to PTC by increasing the risk for rs944289 per year of age, and by enhancing the protective effect of the FOXE1 GGT haplotype in men.

Differences in the transcriptome of medullary thyroid cancer regarding the status and type of RET gene mutations.

Medullary thyroid cancer (MTC) can be caused by germline mutations of the RET proto-oncogene or occurs as a sporadic form. It is well known that RET mutations affecting the cysteine-rich region of the protein (MEN2A-like mutations) are correlated with different phenotypes than those in the kinase domain (MEN2B-like mutations). Our aim was to analyse the whole-gene expression profile of MTC with regard to the type of RET gene mutation and the cancer genetic background (hereditary vs sporadic). We studied 86 MTC samples. We demonstrated that there were no distinct differences in the gene expression profiles of hereditary and sporadic MTCs. This suggests a homogeneous nature of MTC. We also noticed that the site of the RET gene mutation slightly influenced the gene expression profile of MTC. We found a significant association between the localization of RET mutations and the expression of three genes: NNAT (suggested to be a tumour suppressor gene), CDC14B (involved in cell cycle control) and NTRK3 (tyrosine receptor kinase that undergoes rearrangement in papillary thyroid cancer). This study suggests that these genes are significantly deregulated in tumours with MEN2A-like and MEN2B-like mutations; however, further investigations are necessary to demonstrate any clinical impact of these findings.

Brain Metastasis Prediction by Transcriptomic Profiling in Triple-Negative Breast Cancer.

BACKGROUND: Triple-negative breast cancer (TNBC) lacks expression of steroid hormone receptors (estrogen receptor α and progesterone) and epidermal growth factor receptor type 2. This phenotype shows high metastatic potential, with particular predilection to lungs and brain. Determination of TNBC transcriptomic profiles associated with high risk of brain metastasis (BM) might identify patients requiring alternative, more aggressive, or specific preventive and therapeutic approaches. PATIENTS AND METHODS: Using a cDNA-mediated annealing, selection, extension, and ligation assay, we investigated expression of 29,369 gene transcripts in primary TNBC tumor samples from 119 patients-71 in discovery cohort A and 48 in independent cohort B-that included best discriminating genes. Expression of mRNA was correlated with the occurrence of symptomatic BM. RESULTS: In cohort A, the difference at the noncorrected P < .005 was found for 64 transcripts (P = .23 for global test), but none showed significant difference at a preset level of false-discovery rate of < 10%. Of the 30 transcripts with the largest differences between patients with and without BM in cohort A, none was significantly associated with BM in cohort B. CONCLUSION: Analysis based on the primary tumor gene transcripts alone is unlikely to predict BM development in advanced TNBC. Despite its negative findings, the study adds to the knowledge on the biology of TNBC and paves the way for future projects using more advanced molecular assays.

[BRAF initiating mutations in the papillary thyroid carcinoma]. / Inicjujace mutacje genu BRAF w raky brodawkowatym tarczycy.

Among genetic alterations most important for the initiation of papillary thyroid carcinoma (PTC) is mutation T1799A in the BRAF gene which is the most frequent event (54.5%) in this type of thyroid cancer. It is seen in all stages, from microcarcinoma through clinically overt disease to anaplastic cancer. It has been shown that BRAF mutation is correlated with PTC histotype. It is identified most frequently in classical PTC and in tall cell variant. Moreover, BRAF mutation is described more often in older patients, whereas in young patients RET/PTC rearrangements dominate. In PTC cases with BRAF mutation V600E the prognosis is poorer, with more cancer invasiveness, metastasis and recurrence. The presence of BRAF mutation is related to the specific gene expression signature, different than in cancer cases showing RET/PTC rearrangement or no known initiating mutation.

Mouse models of papillary thyroid carcinoma - short review.

Thyroid carcinoma (TC) is the most common endocrine malignancy, and its frequency is still rising. Papillary thyroid carcinoma (PTC) accounts for 80% of all TCs and usually is related to a very good prognosis. However, the standard therapeutic approaches are not always sufficient and disease progression is sometimes observed. These data highlight the limitation of our understanding of molecular mechanisms underlying tumorigenesis and how they vary between individual patients. Over the last 19 years mouse models of thyroid cancers have been developed in order to give answers to questions about their genetic background, relations of key molecular events with pathways fundamental for cancer, and many others. Among these models genetically engineered mice were of utmost importance regarding the input of knowledge about human tumorigenesis. In the present review the most significant mouse models of PTC are described with particular emphasis on BRAFV600E-induced ones, for the sake of its frequency in PTC, relation to factors of poor prognosis, and the fact that, since its identification, it became an attractive target in novel therapies. For the presented mouse models phenotype consequences of particular genetic alterations are described as well as the limitations of the used methods.

Somatic mutation profiling of follicular thyroid cancer by next generation sequencing.

The molecular etiology of follicular thyroid tumors is largely unknown, rendering the diagnostics of these tumors challenging. The somatic alterations present in these tumors apart from RAS gene mutations and PAX8/PPARG translocations are not well described. To evaluate the profile of somatic alteration in follicular thyroid tumors, a total of 82 thyroid tissue samples derived from 48 patients were subjected to targeted Illumina HiSeq next generation sequencing of 372 cancer-related genes. New somatic alterations were identified in oncogenes (MDM2, FLI1), transcription factors and repressors (MITF, FLI1, ZNF331), epigenetic enzymes (KMT2A, NSD1, NCOA1, NCOA2), and protein kinases (JAK3, CHEK2, ALK). Single nucleotide and large structural variants were most and least frequently identified, respectively. A novel translocation in DERL/COX6C was detected. Many somatic alterations in non-coding gene regions with high penetrance were observed. Thus, follicular thyroid tumor somatic alterations exhibit complex patterns. Most tumors contained distinct somatic alterations, suggesting previously unreported heterogeneity.