Skip to main content

Next-generation sequencing in thyroid cancer


Next-generation sequencing (NGS) in thyroid cancer allows for simultaneous high-throughput sequencing analysis of variable genetic alterations and provides a comprehensive understanding of tumor biology. In thyroid cancer, NGS offers diagnostic improvements for fine needle aspiration (FNA) cytology of thyroid with indeterminate features. It also contributes to patient management, providing risk stratification of patients based on the risk of malignancy. Furthermore, NGS has been adopted in cancer research. It is used in molecular tumor classification, and molecular prediction of recurrence and metastasis in papillary thyroid carcinoma. This review covers previous NGS analyses in variable types of thyroid cancer, where samples including FNA cytology, fresh frozen tissue, and formalin-fixed, paraffin-embedded tissues were used. This review also focuses on the clinical and research implications of using NGS to study and treat thyroid cancer.


An understanding of the molecular mechanisms of tumor formation is mandatory for accurate diagnoses and personalized treatments. Previously, single gene assays were commonly used for finding molecular alterations in tumors. Presently, NGS technology provides the simultaneous analysis of hundreds of genes of interest, using targeted sequencing panels [1]. Thus, NGS-based molecular tests for oncology research and clinical practice appear to be rapidly evolving.

Thyroid cancer is the most common malignancy of the endocrine organs; its prevalence is increasing, more than tripling during the last three decades [2]. Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer, followed by follicular carcinoma, medullary carcinoma, poorly differentiated carcinoma, and anaplastic carcinoma. NGS assays can allow improvements in diagnostic accuracy and precise personalized treatments. Thyroid cancers harbor characteristic genetic alterations, including point mutations for proto-oncogenes (BRAF, NRAS, HRAS, KRAS) and chromosomal rearrangements (RET/PTC1, RET/PTC3, PAX8/PPARG), which vary with histologic subtype [3]. This review outlines the results of NGS assays in thyroid cancer, and highlights their clinical implications.

NGS application in the diagnosis of indeterminate cytology

A majority of previous studies using NGS in thyroid cancer analyzed variable specimen sample types and histologic subtypes (Table 1). In clinical practice, NGS assays have been used in the diagnosis of indeterminate cytology of thyroid nodules. FNA is an efficient method for evaluating thyroid nodules that has high sensitivity and specificity. However, FNA has some limitations, since 20–30% of FNA samples fall into categories of indeterminate cytology. These categories include atypia of undetermined significance/follicular lesion of undetermined significance (AUS/FLUS, category III); follicular or oncocytic (Hurthle cell) neoplasm/suspicious for a follicular or oncocytic (Hurthle cell) neoplasm (FN/SFN, category IV); and suspicious for malignant cells (SUSP) [48]. The average cancer risk for these categories is 15.9% in AUS/FLUS, 26.1% in FN/SFN, and 75.2% in SUSP [9]. NGS contributes to diagnostic decision making in patients with indeterminate cytology. Findings from previous studies, which used NGS to analyze thyroid nodules with indeterminate cytology, are summarized in Table 2.

Table 1 Studies of thyroid cancer using next-generation sequencing platforms
Table 2 Application of next-generation sequencing in the thyroid nodule with indeterminate cytology

Genetic alteration of thyroid cancer and NGS

Papillary carcinoma

Several studies have applied NGS to variable subgroups of PTCs. Nikiforova et al. analyzed FFPE or frozen tissue from 27 classic PTCs and 30 FVPTCs, using the ThyroSeq NGS panel targeting 12 genes with 34 amplicons on the Ion Torrent PGM sequencer [10]. The results showed that 70% of classic PTCs harbored mutated genes: BRAF (59%) was the most frequent, followed by PIK3A (11%), TP53 (7%), and NRAS (4%). In contrast, 83% of FVPTCs had mutated genes: RAS (73%) was the most frequent, followed by BRAF (7%) and TSHR (3%) [10]. Smallridge et al. performed RNA sequencing (RNA-Seq) using the Illumina HiSeq 2000 platform on frozen tissue from 12 BRAF V600E-mutated PTCs and 8 BRAF-wild type PTCs [11]. Among the 13,085 genes interrogated, 560 were differentially expressed between BRAF V600E-mutated PTCs and BRAF-wild type PTCs [11]. Among these 560 genes, 67 were related to immune function pathways, 51 were under-expressed in BRAF V600E-mutated PTCs, and HLAG, CXCL14, TIMP1, and IL1RAP were over-expressed. In BRAF-wild type PTCs, 4 immune function genes (IL1B, CCL19, CCL21, and CXCR4) were most significantly differentially expressed, and exhibited a high degree of correlation with lymphocytic infiltration [11]. In a study by Leeman-Neill et al. fresh frozen tissue from 62 radiation-associated PTCs and 151 sporadic PTCs was analyzed using RNA-Seq on the Illumina HiSeq 2000 platform. This identified an ETV6-NTRK3 rearrangement in 14.5% of radiation-associated PTCs and 2% of sporadic PTCs [12]. The authors suggested that an ETV6-NTRK3 rearrangement may be a key mechanism of radiation-induced carcinogenesis [12].

The Cancer Genome Atlas (TCGA) Network described the genomic characterization of 496 PTCs, and generated data using whole genome sequencing. This was done on the NGS platform and by a multiplatform analysis of SNP arrays, DNA methylation, and reverse phase protein arrays [13]. In PTCs that lacked known driver mutations, alterations of EIF1AX, PPMID, and CHEK2 were discovered as potential new tumor-initiating mutations. The TCGA project identified the TERT promoter mutation, which accounts for approximately 1% of PTCs, but shows association with a high risk of recurrence. Based on the multi-level molecular data, PTCs were separated into two groups of distinct downstream signaling pathways: the BRAFV600E-like cohort and the RAS-like cohort. Genomic, epigenomic, and proteomic differences were revealed between these two groups, and most of the RAS-like PTCs were follicular variant PTCs (FVPTCs).

Regarding pediatric thyroid carcinoma, Picarsic et al. analyzed 17 pediatric PTCs (age range 8–17 years, median 13 years) from FNA, fresh frozen tissue, and FFPE samples. Mutation analysis with a 7-gene mutation panel using real-time PCR and ThyroSeq v2 on the Ion Torrent PGM sequencer showed that: (1) The detection rate of molecular alterations was increased by up to 87% by the ThyroSeq v2 NGS assay compared to an increase of 60% by the 7-gene mutation panel. (2) In pediatric thyroid carcinoma, chromosomal rearrangement (53%) was more common than point mutation (33%). (3) ETV6-NTRK3 fusion was identified in 18% of samples, and was associated with aggressive histologic features such as non-encapsulation, solid/insular/trabecular patterns, extensive glandular involvement, and thick tumor fibrosis [14]. Ballester et al. analyzed FFPE and FNA samples from 25 pediatric PTCs (age range 10–19 years, median 14 years) using the 50-gene Ion AmpliSeq Cancer Hotspot Panel v2 on the Ion Torrent PGM sequencer [15]. No additional mutations were detected using the NGS assay on pediatric PTCs that initially were negative for BRAF V600E mutation, RET/PTC1/3 fusion, and TERT promoter mutation [15].

Follicular carcinoma

Follicular carcinoma (FC) is a well-differentiated thyroid carcinoma and the second-most common thyroid cancer after PTC. It accounts for 10% of the total thyroid cancers [16]. Although it is not firmly established, FC is classified into the minimally invasive type and the widely invasive type, according to the microscopic tumor extent [17]. It can be classified into the conventional type and the oncocytic type (Hurthle cell type) on the basis of cell type [18].

A previous study retrospectively collected FFPE or fresh frozen tissue samples from 36 FCs to study 12 cancer genes and 34 amplicons using the ThyroSeq panel on the Ion Torrent PGM sequencer. The analysis detected mutations of NRAS (N = 9), KRAS (N = 2), HRAS (N = 1), TSHR (N = 4), TP53 (N = 4), and PTEN (N = 1) [10]. Interestingly, conventional type (N = 18) and oncocytic type (N = 18) samples showed distinguished genetic alterations. In the conventional type FCs, NRAS was the most frequently affected gene, followed by TSHR and KRAS. TP53 was the most commonly mutated gene in the oncocytic type FCs, followed by HRAS, KRAS, and PTEN [19].

Swierniak et al. analyzed 26 follicular adenomas, 22 FCs, and 34 paired normal thyroid tissue samples. Targeted NGS of 372 genes using the TruSeq kit on the Illumina HiSeq 1500 platform yielded the following results [20]: (1) 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). (2) Single nucleotide variants were the most common types of mutations, and large structural variants were the least frequent. (3) A novel translocation in DERL/COX6C was detected. (4) Somatic alteration affected non-coding gene regions and exhibited high penetrance. These results suggest that FC has significant molecular heterogeneity, since FC reveals far more complex somatic alterations than PTC, and each tumor harbors distinct somatic alterations.

Poorly differentiated carcinoma and anaplastic carcinoma

Poorly differentiated carcinoma (PDC) and anaplastic carcinoma (AC) are rare types of thyroid carcinoma, each with a prevalence of 10% [21, 22], and a 1–2% among all thyroid carcinomas [23]. PDC and AC have poor prognoses, with a 5-year survival rate of 51% and 0%, respectively [24]. Since these types of cancers respond poorly to conventional treatment options (including radioiodine therapy, chemotherapy, and radiotherapy), there are already clinical trials for molecular-targeted agents underway [25, 26]. Using NGS, it may be possible to identify targetable gene alterations that can improve the course of patient treatment.

In a study using the ThyroSeq panel on the Ion Torrent PGM sequencer, 12 genes with 34 amplicons were analyzed from FFPE or fresh frozen tissue from 10 PDCs and 27 ACs. The study showed that 30% of PDCs had mutations, whereas 74% of ACs had mutations [10]. The altered genes were NRAS, PIK3CA, GNAS, and BRAF in PDCs, and TP53, BRAF, RAS, PIK3CA, PTEN, and CTNNB1 in ACs.

Sykorova et al. analyzed fresh frozen tissue samples from 3 PDCs and 5 ACs using the TruSight Cancer Panel, targeting 94 cancer-related genes on the Illumina MiSeq sequencer [27]. All PDCs and ACs showed more than one genetic alteration, and TP53 mutations were identified in all but 2 cases [27]. CDH1, FANCD2, CHECK2, ADH1B, GPC3, TP53, and PTEN genes were altered in PDCs, and ATM, HNF1A, MET, NF1, TP53, PTEN, MSH2, RB1, NBN, NF1, MUTYH, TSC2, HRAS, and EGFR genes were altered in ACs [27]. However, the study could not assess mutation of larger genes or chromosomal rearrangements in a panel of 94 known cancer genes, nor could it distinguish germline mutations among the detected genetic changes.

Landa et al. performed NGS using the MSK-IMPACT cancer exome panel, and analyzed 341 genes from FFPE (N = 80) or fresh-frozen tissue (N = 37) samples from 34 PDCs and 33 ACs [28]. The analysis revealed the following: (1) ACs showed a higher mutation number than PDCs (6 ± 5 vs. 2 ± 3, median ± interquartile range), and harbored a higher frequency of mutation in TP35, TERT promoter, PI3 K/AKT/mTOR pathway effector, SWI/SNF subunit, and histone methyltransferase. (2) In PDCs, clinicopathologic features were different based on the genetic alterations: 92% of RAS mutations were found in PDCs met the Turin criteria, whereas 81% of BRAF mutations were found in PDCs met the Memorial Sloan-Kettering Cancer Center (MSKCC) criteria. PDCs harboring BRAF mutations were smaller and had frequent nodal metastasis, whereas RAS-mutant PDCs showed larger tumor sizes and a higher rate of distant metastasis. (3) The association of mutation between EIF1AX and RAS was notable in both PDCs and ACs. EIF1AX mutations have been reported in approximately 1% of PTCs, and known to be mutually exclusive for BRAF and RAS mutations. However, EIF1AX mutations were found in 11% of PDCs and 9% of ACs, and 93% were associated with RAS mutations. (4) Chromosomal rearrangements (including RET/PTC, ALK, and PAX8-PPARG fusions) were found in 14% of PDCs, but were absent in ACs.

Medullary carcinoma

Medullary carcinoma (MC) is a neuroendocrine tumor originating from C-cells, and accounts for approximately 5% of total thyroid cancers [29]. MC is composed of 75% sporadic form and 25% hereditary form, the latter resulting from a RET proto-oncogene mutation [30, 31].

Nikiforova et al. analyzed 12 genes and 34 amplicons using the ThyroSeq panel on the Ion Torrent PGM sequencer with FFPE or fresh frozen tissue samples from 15 sporadic MCs. Mutations were identified in 11 (73%) MCs, of which 7 (47%) were RET mutations, 3 (20%) were HRAS mutations, and 1 (7%) was a KRAS mutation [10]. Simbolo et al. analyzed 50 cancer-related genes using the Ion AmpliSeq Hot Spot Cancer Panel v2 on the Ion Torrent PGM sequencer. Out of 20 retrospectively collected FFPE samples of MCs, the study found that 85% of MCs harbored mutations as follows: 13 RET mutations (60%), 3 HRAS mutations (15%), 1 KRAS mutation (5%), 1 STK11 mutation (5%), and 3 samples where mutations were undetected (15%) [32]. RET status was evaluated with both the NGS and Sanger sequencing methods, and NGS showed higher sensitivity than the Sanger method. NGS identified an additional 3 RET mutations, which were undetected by the Sanger method. Although multiple mutations were found in MCs by the NGS assay, a relevant therapeutic target has not been identified, and further investigation is required to improve the treatment of MCs.

Limitation of NGS technology for thyroid cancer

One of the most profound limitations of applying NGS for thyroid cancer is a lack of sufficient evidence-based framework applicable to the clinical practice. As discussed in this review, the number of existing studies using NGS to analyze thyroid cancer stands at less than 15. Also, most of the previous studies were carried out at single institutes using specific subtypes of thyroid cancer in small sample sizes, rather than all types of thyroid cancers; therefore, resulting data would still be insufficient for making decisions on either patient diagnosis or treatment in clinical practice. To overcome such limitations, a large-scaled, global, and multicentered NGS study for thyroid cancer is required.

In a thyroid nodule with indeterminate cytology and BRAF V600E detected by NGS, surgical resection would be the most appropriate treatment option, since BRAF V600E is a highly specific mutation for PTC [13]. In contrast, RAS mutation can be detected in FVPTC and follicular neoplasm that require surgical excision, while also being present in benign adenomatous nodules [33] that do not require excision; therefore, further studies to identify the optimal treatment plan specific to mutation is needed. Apart from mutational variant, inadequate sample preparation, of both poor quality and quantity, can lead to false negative results. In samples with low tumor purity and small amount of DNA, low coverage would not be able to detect in allele with low frequency. Although DNA quality of cytology specimen is better than that of FFPE tissue, cytology specimens would contain some amount of normal tissue component. Also, evaluation of tumor purity may be an essential step before DNA preparation, particularly if the target nodule is small or an inexperienced person performs the aspiration procedure.

In addition to well-known BRAF, RAS, and RET mutations, NGS technology facilitated detection of new somatic alterations in thyroid cancer such as MITF, MDM2, JAK3, FLI1, IDH1 etc. [20], all in which the significance of thyroid cancer has not been delineated yet. Larger scales of integrated genomic and phenotypic database should be provided to interpret NGS results. Also, an appropriate reporting system for NGS results in thyroid cancer is needed. Results from NGS analysis may encompass multiple variants, and each variant may have different clinical and biological significance. An appropriate tier system, with specific level of evidence, is required for reporting NGS results. Working group with large expertise to build a consensus guideline for reporting NGS results in thyroid cancer is requested. Besides intrinsic limitation of NGS platform, such as low detection rate of large indels, annotation errors of pipeline can be present. Clinicopathologic correlation and additional knowledge-based review of NGS report are essential for result interpretations. Currently, agents specifically targeting defined mutations are available, and patients who have the targetable mutation can benefit from optimized treatment and avoid unnecessary therapy. Inclusion of potential therapeutic target genes in the gene panel of targeted NGS, as well as the accumulation of information to build up database for future investigation, would also be required.

Most of NGS panels applied in thyroid cancer studies cover hotspot mutations, and are highly sensitive for evaluation of limited regions of selected genes; however, relevant mutation could be missing if not appropriately mapped. As described in previous studies, targeted NGS with specific gene panel showed high PPV. However, NGS analysis using 7-gene panel in previous study showed that 30–35% of thyroid cancer patients were still negative for mutation, and low NPV would require diagnostic surgery in benign nodules to prevent missing cancers [34]. High sensitivity of NGS technique showed that subclones within a nodule with mutations leading to aggressive clinical behaviors might be detected with low allele frequency [35]. Clinicians would face a dilemma with such cases, regarding whether to follow up with their patients or to refer them to surgical resection. Although different platforms and variant-calling pipelines were proven to have high concordance and sensitivity, detected mutations would be different in cases with tumor heterogeneity [36]. Genetic alteration of tumor results in tumor heterogeneity, which can be divided into intertumor heterogeneity that shows different genetic alterations based on the tumor sites, and intratumor heterogeneity that contains different genetic alterations within a same tumor. Concepts of aggressive clone and tumor heterogeneity are also present in thyroid cancer [3739]. Furthermore, complicated situations may derive from NGS results, depending on whether NGS for each tissue sample was performed at relevant site and relevant time, and this can lead to repeated aspiration and/or biopsy. Also, in metastatic cancer, tissue samples should be obtained from metastatic sites; however, sites such as the brain or a specific bone are challenging for proper tissue sampling.

Future of NGS for thyroid cancer

The fast-evolving NGS technology offers a cost-effective approach for cancer genomics, as well as in thyroid cancer. In future prospects of thyroid cancer, NGS can be used to detect circulating tumor cells or cell-free plasma DNA to identify early relapse and/or residual disease. Previous studies reported presence of circulating tumor cells or cell-free plasma DNA in thyroid cancer patients [40, 41], which appeared to be future candidates for NGS application. Furthermore, NGS can detect tumor-specific genetic alterations, which are used in follow-ups for patient monitoring. In patient monitoring, genetic alterations should be present in all tumor cells, while consistently and sustainably existing from both tumor development and during tumor progression. In thyroid cancer, BRAF V600E is the most common and the earliest genetic event in PTC [42, 43], and it appears to be a good candidate gene for monitoring. Also, during radioactive iodine and/or drug treatment of thyroid cancer, new mutation variants other than primary tumor can be recognized in NGS analysis of either circulating tumor cells or cell-free plasma DNA. This concept can be applied in identifying genetic alterations related to an acquired resistance to treatment during clinical course. Currently available studies of NGS application in thyroid cancer tend to focus on evaluating genetic alterations in specific types of thyroid cancer. In future research, an improved high-throughput pipeline should be used for a more comprehensive analysis of gene expression and DNA binding activity. In addition, a systems biology approach would also help discover the interaction and casual relationship between genes and/or proteins, introducing a new ground of thyroid cancer biology.


The emergence of NGS technology has provided in-depth analysis of multiple, diverse cancers by a number of devices and gene panels, and has led to more effective options for cancer screening, prevention, diagnosis, prognosis, and targeted therapy. The use of NGS to study thyroid cancer has improved our understanding of the molecular genetics of thyroid cancer. In thyroid nodules of indeterminate cytology, such as FN/SFN and AUS/FLUS, the NGS test detected multiple genetic alterations and identified patients with a high risk of malignancy. Risk stratification using molecular signatures offers many more precise treatment options during patient management. The application of NGS for PTC, FC, MC, PDC, and AC revealed novel genetic alterations which were not detected by past sequencing methods. Newly discovered genetic alterations include genes associated with tumor recurrence and distant metastasis, which are candidates for molecular prognostic markers. However, limitations are also present with NGS, arising from variable sample types, multiple platforms and gene panels, and variable analysis programs, each of which can confound results. Standardization for quality control and the data-analytic process is needed to minimize the discrepancies between analyses. For poor-prognostic histologic types of thyroid cancer—MC, PDC, and AC—NGS studies identified several novel genetic alterations, but drug-actionable target genes have not been identified yet, and further investigation is required. Nevertheless, development of new sequencing technologies, such as NGS, enhances the cancer genome body of knowledge, and allows for more effective cancer screening, prevention, diagnosis, and monitoring. This in turn provides for better precision medicine and more curative cancer treatments.



next-generation sequencing


fine needle aspiration


papillary thyroid carcinoma


atypia of undetermined significance/follicular lesion of undetermined significance


follicular or oncocytic (Hurthle cell) neoplasm


suspicious for malignant cells


positive predictive value


negative predictive value


The Cancer Genome Atlas


follicular variant papillary thyroid carcinoma


follicular carcinoma


poorly differentiated carcinoma


anaplastic carcinoma


Memorial Sloan-Kettering Cancer Center


medullary carcinoma


  1. LeBlanc VG, Marra MA. Next-generation sequencing approaches in cancer: where have they brought us and where will they take us? Cancers (Basel). 2015;7:1925–58.

    Article  Google Scholar 

  2. Jemal A, Siegel R, Ward E, Hao Y, Xu J, Murray T, Thun MJ. Cancer statistics, 2008. CA Cancer J Clin. 2008;58:71–96.

    Article  PubMed  Google Scholar 

  3. Hunt J. Understanding the genotype of follicular thyroid tumors. Endocr Pathol. 2005;16:311–21.

    Article  CAS  PubMed  Google Scholar 

  4. Baloch ZW, LiVolsi VA, Asa SL, Rosai J, Merino MJ, Randolph G, Vielh P, DeMay RM, Sidawy MK, Frable WJ. Diagnostic terminology and morphologic criteria for cytologic diagnosis of thyroid lesions: a synopsis of the National Cancer Institute Thyroid Fine-Needle Aspiration State of the Science Conference. Diagn Cytopathol. 2008;36:425–37.

    Article  PubMed  Google Scholar 

  5. Gharib H. Changing trends in thyroid practice: understanding nodular thyroid disease. Endocr Pract. 2004;10:31–9.

    Article  PubMed  Google Scholar 

  6. Sclabas GM, Staerkel GA, Shapiro SE, Fornage BD, Sherman SI, Vassillopoulou-Sellin R, Lee JE, Evans DB. Fine-needle aspiration of the thyroid and correlation with histopathology in a contemporary series of 240 patients. Am J Surg. 2003;186:702–9 (discussion 709–710).

    Article  PubMed  Google Scholar 

  7. Yassa L, Cibas ES, Benson CB, Frates MC, Doubilet PM, Gawande AA, Moore FD Jr, Kim BW, Nose V, Marqusee E, Larsen PR, Alexander EK. Long-term assessment of a multidisciplinary approach to thyroid nodule diagnostic evaluation. Cancer. 2007;111:508–16.

    Article  PubMed  Google Scholar 

  8. Lewis CM, Chang KP, Pitman M, Faquin WC, Randolph GW. Thyroid fine-needle aspiration biopsy: variability in reporting. Thyroid. 2009;19:717–23.

    Article  CAS  PubMed  Google Scholar 

  9. Bongiovanni M, Spitale A, Faquin WC, Mazzucchelli L, Baloch ZW. The Bethesda system for reporting thyroid cytopathology: a meta-analysis. Acta Cytol. 2012;56:333–9.

    Article  PubMed  Google Scholar 

  10. Nikiforova MN, Wald AI, Roy S, Durso MB, Nikiforov YE. Targeted next-generation sequencing panel (ThyroSeq) for detection of mutations in thyroid cancer. J Clin Endocrinol Metab. 2013;98:E1852–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Smallridge RC, Chindris AM, Asmann YW, Casler JD, Serie DJ, Reddi HV, Cradic KW, Rivera M, Grebe SK, Necela BM, Eberhardt NL, Carr JM, McIver B, Copland JA, Thompson EA. RNA sequencing identifies multiple fusion transcripts, differentially expressed genes, and reduced expression of immune function genes in BRAF (V600E) mutant vs BRAF wild-type papillary thyroid carcinoma. J Clin Endocrinol Metab. 2014;99:E338–47.

    Article  CAS  PubMed  Google Scholar 

  12. Leeman-Neill RJ, Kelly LM, Liu P, Brenner AV, Little MP, Bogdanova TI, Evdokimova VN, Hatch M, Zurnadzy LY, Nikiforova MN, Yue NJ, Zhang M, Mabuchi K, Tronko MD, Nikiforov YE. ETV6-NTRK3 is a common chromosomal rearrangement in radiation-associated thyroid cancer. Cancer. 2014;120:799–807.

    Article  CAS  PubMed  Google Scholar 

  13. Agrawal N, Akbani R, Aksoy BA, Ally A, Arachchi H, Asa SL, Auman JT, Balasundaram M, Balu S, Baylin SB, Behera M, Bernard B, Beroukhim R, Bishop JA, Black AD, Bodenheimer T, Boice L, Bootwalla MS, Bowen J, Bowlby R, Bristow CA, Brookens R, Brooks D, Bryant R, Buda E, Butterfield YSN, Carling T, Carlsen R, Carter SL, Carty SE, Chan TA, Chen AY, Cherniack AD, Cheung D, Chin L, Cho J, Chu A, Chuah E, Cibulskis K, Ciriello G, Clarke A, Clayman GL, Cope L, Copland J, Covington K, Danilova L, Davidsen T, Demchok JA, DiCara D, Dhalla N, Dhir R, Dookran SS, Dresdner G, Eldridge J, Eley G, El-Naggar AK, Eng S, Fagin JA, Fennell T, Ferris RL, Fisher S, Frazer S, Frick J, Gabriel SB, Ganly I, Gao J, Garraway LA, Gastier-Foster JM, Getz G, Gehlenborg N, Ghossein R, Gibbs RA, Giordano TJ, Gomez-Hernandez K, Grimsby J, Gross B, Guin R, Hadjipanayis A, Harper HA, Hayes DN, Heiman DI, Herman JG, Hoadley KA, Hofree M, Holt RA, Hoyle AP, Huang FW, Huang M, Hutter CM, Ideker T, Iype L, Jacobsen A, Jefferys SR, Jones CD, Jones SJM, Kasaian K, Kebebew E, Khuri FR, Kim J, Kramer R, Kreisberg R, Kucherlapati R, Kwiatkowski DJ, Ladanyi M, Lai PH, Laird PW, Lander E, Lawrence MS, Lee D, Lee E, Lee S, Lee W, Leraas KM, Lichtenberg TM, Lichtenstein L, Lin P, Ling S, Liu J, Liu W, Liu Y, LiVolsi VA, Lu Y, Ma Y, Mahadeshwar HS, Marra MA, Mayo M, McFadden DG, Meng S, Meyerson M, Mieczkowski PA, Miller M, Mills G, Moore RA, Mose LE, Mungall AJ, Murray BA, Nikiforov YE, Noble MS, Ojesina AI, Owonikoko TK, Ozenberger BA, Pantazi A, Parfenov M, Park PJ, Parker JS, Paull EO, Pedamallu CS, Perou CM, Prins JF, Protopopov A, Ramalingam SS, Ramirez NC, Ramirez R, Raphael BJ, Rathmell WK, Ren X, Reynolds SM, Rheinbay E, Ringel MD, Rivera M, Roach J, Robertson AG, Rosenberg MW, Rosenthall M, Sadeghi S, Saksena G, Sander C, Santoso N, Schein JE, Schultz N, Schumacher SE, Seethala RR, Seidman J, Senbabaoglu Y, Seth S, Sharpe S, Mills Shaw KR, Shen JP, Shen R, Sherman S, Sheth M, Shi Y, Shmulevich I, Sica GL, Simons JV, Sipahimalani P, Smallridge RC, Sofia HJ, Soloway MG, Song X, Sougnez C, Stewart C, Stojanov P, Stuart JM, Tabak B, Tam A, Tan D, Tang J, Tarnuzzer R, Taylor BS, Thiessen N, Thorne L, Thorsson V, Tuttle RM, Umbricht CB, Van Den Berg DJ, Vandin F, Veluvolu U, Verhaak RGW, Vinco M, Voet D, Walter V, Wang Z, Waring S, Weinberger PM, Weinstein JN, Weisenberger DJ, Wheeler D, Wilkerson MD, Wilson J, Williams M, Winer DA, Wise L, Wu J, Xi L, Xu AW, Yang L, Yang L, Zack TI, Zeiger MA, Zeng D, Zenklusen JC, Zhao N, Zhang H, Zhang J, Zhang J, Zhang W, Zmuda E, Zou L. Integrated genomic characterization of papillary thyroid carcinoma. Cell. 2014;159:676–90

    Article  PubMed Central  Google Scholar 

  14. Picarsic JL, Buryk MA, Ozolek J, Ranganathan S, Monaco SE, Simons JP, Witchel SF, Gurtunca N, Joyce J, Zhong S, Nikiforova MN, Nikiforov YE. Molecular characterization of sporadic pediatric thyroid carcinoma with the DNA/RNA ThyroSeq v2 next-generation sequencing assay. Pediatr Dev Pathol. 2016;19:115–22.

    Article  PubMed  Google Scholar 

  15. Ballester LY, Sarabia SF, Sayeed H, Patel N, Baalwa J, Athanassaki I, Hernandez JA, Fang E, Quintanilla NM, Roy A, Lopez-Terrada DH. Integrating molecular testing in the diagnosis and management of children with thyroid lesions. Pediatr Dev Pathol. 2016;19:94–100.

    Article  PubMed  Google Scholar 

  16. Correa P, Chen VW. Endocrine gland cancer. Cancer. 1995;75:338–52.

    Article  CAS  PubMed  Google Scholar 

  17. Delellis RA, Lloyd RV, Heitz PU. Tumours of the thyroid and parathyroid, in WHO classification of tumours: pathology and genetics of tumours of endocrine organs. Lyon: IARC Press; 2004.

  18. Kushchayeva Y, Duh QY, Kebebew E, D’Avanzo A, Clark OH. Comparison of clinical characteristics at diagnosis and during follow-up in 118 patients with Hurthle cell or follicular thyroid cancer. Am J Surg. 2008;195:457–62.

    Article  PubMed  Google Scholar 

  19. Karunamurthy A, Panebianco F, Hsiao SJ, Vorhauer J, Nikiforova MN, Chiosea S, Nikiforov YE. Prevalence and phenotypic correlations of EIF1AX mutations in thyroid nodules. Endocr Relat Cancer. 2016;23:295–301.

    Article  PubMed  Google Scholar 

  20. Swierniak M, Pfeifer A, Stokowy T, Rusinek D, Chekan M, Lange D, Krajewska J, Oczko-Wojciechowska M, Czarniecka A, Jarzab M, Jarzab B, Wojtas B. Somatic mutation profiling of follicular thyroid cancer by next generation sequencing. Mol Cell Endocrinol. 2016;433:130–7.

    Article  CAS  PubMed  Google Scholar 

  21. Burman KD, Ringel MD, Wartofsky L. Unusual types of thyroid neoplasms. Endocrinol Metab Clin North Am. 1996;25:49–68.

    Article  CAS  PubMed  Google Scholar 

  22. Sywak M, Pasieka JL, Ogilvie T. A review of thyroid cancer with intermediate differentiation. J Surg Oncol. 2004;86:44–54.

    Article  PubMed  Google Scholar 

  23. Gilliland FD, Hunt WC, Morris DM, Key CR. Prognostic factors for thyroid carcinoma. A population-based study of 15,698 cases from the surveillance, epidemiology and end results (SEER) program 1973–1991. Cancer. 1997;79:564–73.

    Article  CAS  PubMed  Google Scholar 

  24. Wreesmann VB, Ghossein RA, Patel SG, Harris CP, Schnaser EA, Shaha AR, Tuttle RM, Shah JP, Rao PH, Singh B. Genome-wide appraisal of thyroid cancer progression. Am J Pathol. 2002;161:1549–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Cabanillas ME, Zafereo M, Gunn GB, Ferrarotto R. Anaplastic thyroid carcinoma: treatment in the age of molecular targeted therapy. J Oncol Pract. 2016;12:511–8.

    Article  PubMed  Google Scholar 

  26. Smith N, Nucera C. Personalized therapy in patients with anaplastic thyroid cancer: targeting genetic and epigenetic alterations. J Clin Endocrinol Metab. 2015;100:35–42.

    Article  CAS  PubMed  Google Scholar 

  27. Sykorova V, Dvorakova S, Vcelak J, Vaclavikova E, Halkova T, Kodetova D, Lastuvka P, Betka J, Vlcek P, Reboun M, Katra R, Bendlova B. Search for new genetic biomarkers in poorly differentiated and anaplastic thyroid carcinomas using next generation sequencing. Anticancer Res. 2015;35:2029–36.

    CAS  PubMed  Google Scholar 

  28. Landa I, Ibrahimpasic T, Boucai L, Sinha R, Knauf JA, Shah RH, Dogan S, Ricarte-Filho JC, Krishnamoorthy GP, Xu B, Schultz N, Berger MF, Sander C, Taylor BS, Ghossein R, Ganly I, Fagin JA. Genomic and transcriptomic hallmarks of poorly differentiated and anaplastic thyroid cancers. J Clin Invest. 2016;126:1052–66.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Siegel R, Naishadham D, Jemal A. Cancer statistics, 2012. CA Cancer J Clin. 2012;62:10–29.

    Article  PubMed  Google Scholar 

  30. de Groot JW, Links TP, Plukker JT, Lips CJ, Hofstra RM. RET as a diagnostic and therapeutic target in sporadic and hereditary endocrine tumors. Endocr Rev. 2006;27:535–60.

    Article  PubMed  Google Scholar 

  31. Moura MM, Cavaco BM, Pinto AE, Leite V. High prevalence of RAS mutations in RET-negative sporadic medullary thyroid carcinomas. J Clin Endocrinol Metab. 2011;96:E863–8.

    Article  CAS  PubMed  Google Scholar 

  32. Simbolo M, Mian C, Barollo S, Fassan M, Mafficini A, Neves D, Scardoni M, Pennelli G, Rugge M, Pelizzo MR, Cavedon E, Fugazzola L, Scarpa A. High-throughput mutation profiling improves diagnostic stratification of sporadic medullary thyroid carcinomas. Virchows Arch. 2014;465:73–8.

    Article  CAS  PubMed  Google Scholar 

  33. Medici M, Kwong N, Angell TE, Marqusee E, Kim MI, Frates MC, Benson CB, Cibas ES, Barletta JA, Krane JF, Ruan DT, Cho NL, Gawande AA, Moore FD Jr, Alexander EK. The variable phenotype and low-risk nature of RAS-positive thyroid nodules. BMC Med. 2015;13:184.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Nikiforov YE, Carty SE, Chiosea SI, Coyne C, Duvvuri U, Ferris RL, Gooding WE, Hodak SP, LeBeau SO, Ohori NP, Seethala RR, Tublin ME, Yip L, Nikiforova MN. Highly accurate diagnosis of cancer in thyroid nodules with follicular neoplasm/suspicious for a follicular neoplasm cytology by ThyroSeq v2 next-generation sequencing assay. Cancer. 2014;120:3627–34.

    Article  CAS  PubMed  Google Scholar 

  35. Le Mercier M, D’Haene N, De Neve N, Blanchard O, Degand C, Rorive S, Salmon I. Next-generation sequencing improves the diagnosis of thyroid FNA specimens with indeterminate cytology. Histopathology. 2015;66:215–24.

    Article  PubMed  Google Scholar 

  36. Frampton GM, Fichtenholtz A, Otto GA, Wang K, Downing SR, He J, Schnall-Levin M, White J, Sanford EM, An P, Sun J, Juhn F, Brennan K, Iwanik K, Maillet A, Buell J, White E, Zhao M, Balasubramanian S, Terzic S, Richards T, Banning V, Garcia L, Mahoney K, Zwirko Z, Donahue A, Beltran H, Mosquera JM, Rubin MA, Dogan S, Hedvat CV, Berger MF, Pusztai L, Lechner M, Boshoff C, Jarosz M, Vietz C, Parker A, Miller VA, Ross JS, Curran J, Cronin MT, Stephens PJ, Lipson D, Yelensky R. Development and validation of a clinical cancer genomic profiling test based on massively parallel DNA sequencing. Nat Biotech. 2013;31:1023–31.

    Article  CAS  Google Scholar 

  37. Kuhn E, Teller L, Piana S, Rosai J, Merino MJ. Different clonal origin of bilateral papillary thyroid carcinoma, with a review of the literature. Endocr Pathol. 2012;23:101–7.

    Article  PubMed  Google Scholar 

  38. Le Pennec S, Konopka T, Gacquer D, Fimereli D, Tarabichi M, Tomas G, Savagner F, Decaussin-Petrucci M, Tresallet C, Andry G, Larsimont D, Detours V, Maenhaut C. Intratumor heterogeneity and clonal evolution in an aggressive papillary thyroid cancer and matched metastases. Endocr Relat Cancer. 2015;22:205–16.

    Article  PubMed  Google Scholar 

  39. Schopper HK, Stence A, Ma D, Pagedar NA, Robinson RA. Single thyroid tumour showing multiple differentiated morphological patterns and intramorphological molecular genetic heterogeneity. J Clin Pathol 2016.

  40. Lubitz CC, Parangi S, Holm TM, Bernasconi MJ, Schalck AP, Suh H, Economopoulos KP, Gunda V, Donovan SE, Sadow PM, Wirth LJ, Sullivan RJ, Panka DJ. Detection of circulating BRAF(V600E) in patients with papillary thyroid carcinoma. J Mol Diagn. 2016;18:100–8.

    Article  CAS  PubMed  Google Scholar 

  41. Xu JY, Handy B, Michaelis CL, Waguespack SG, Hu MI, Busaidy N, Jimenez C, Cabanillas ME, Fritsche HA, Cote GJ, Sherman SI. Detection and prognostic significance of circulating tumor cells in patients with metastatic thyroid cancer. J Clin Endocrinol Metab. 2016;101:4461–7.

    Article  PubMed  Google Scholar 

  42. Xu X, Quiros RM, Gattuso P, Ain KB, Prinz RA. High prevalence of BRAF gene mutation in papillary thyroid carcinomas and thyroid tumor cell lines. Cancer Res. 2003;63:4561–7.

    CAS  PubMed  Google Scholar 

  43. Yarchoan M, LiVolsi VA, Brose MS. BRAF mutation and thyroid cancer recurrence. J Clin Oncol. 2015;33:7–8.

    Article  PubMed  Google Scholar 

  44. Nikiforov YE, Carty SE, Chiosea SI, Coyne C, Duvvuri U, Ferris RL, Gooding WE, LeBeau SO, Ohori NP, Seethala RR, Tublin ME, Yip L, Nikiforova MN. Impact of the multi-gene ThyroSeq next-generation sequencing assay on cancer diagnosis in thyroid nodules with atypia of undetermined significance/follicular lesion of undetermined significance cytology. Thyroid. 2015;25:1217–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

YJC conducted the literature search and wrote manuscript. JSK conceived the study, participated in its design and coordination, and helped to draft the manuscript. Both authors read and approved the final manuscript.


Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and material

Data sharing not applicable to this article as no datasets were generated or analysed during the current study.


This study was supported by a grant from the National R&D Program for Cancer Control, Ministry of Health & Welfare, Republic of Korea (1420080). This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (2015R1A1A1A05001209).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Ja Seung Koo.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cha, Y.J., Koo, J.S. Next-generation sequencing in thyroid cancer. J Transl Med 14, 322 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: