Evolution and Impact of Subclonal Mutations in Papillary Thyroid Cancer

Exploration of the prognostic value of methylation regulators related to m5C in papillary thyroid carcinoma

The incidence of papillary thyroid carcinoma (PTC) has increased significantly in recent years, and for patients with metastatic and recurrent PTC, the options for treatment currently available are insufficient. To date, the exact molecular mechanism underlying PTC is still not fully understood. 5-Methylcytosine (m5C) RNA methylation is associated with the prognosis of a variety of tumors. However, the molecular mechanisms and biomarkers associated with m5C in the diagnosis, treatment, and prognosis of this disease have not been fully elucidated. Ten m5C regulators with significantly different expression levels were included in this study. Immune infiltration analysis revealed significant negative correlations between most of these regulators and regulatory T cells. TRDMT1, NSUN5, and NSUN6 had high weights and strong correlations in the protein-protein interaction network. Using gene ontology, Kyoto Encyclopedia of Genes and Genomes, and gene set enrichment analysis, 1489 differentially expressed genes were screened from The Cancer Genome Atlas messenger RNA matrix, indicating that these differentially expressed genes were significantly enriched in various pathways and functions related to cancers. Four m5C regulators, NSUN2, NSUN4, NSUN6, and DNMT3B, were screened as prognostic markers by least absolute shrinkage and selection operator regression analysis, and NSUN2 and NSUN6 were identified as risk factors for poor prognosis. We found that the prognostic prediction model constructed using the m5C regulators NSUN2, NSUN4, NSUN6, and DNMT3B showed good prognostic prediction ability and diagnostic ability. This model was applied to predict the survival probability of patients with PTC, the prediction ability of 5-year survival was the best. The multi-factor prognostic prediction model combined with the tumor node metastasis stage and risk score grouping showed better prognostic predictive power.

Mapping recurrent mosaic copy number variation in human neurons

When somatic cells acquire complex karyotypes, they often are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons are limited by relatively small sample sizes. Here, we develop an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We apply this approach to 2,125 frontal cortical neurons from a neurotypical human brain. SCOVAL identifies 226 CNV neurons, which include a subclass of 65 CNV neurons with highly aberrant karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we find that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contain fewer, but longer, genes.

Genetic alterations and allele frequency of BRAF V600E and TERT mutation in papillary thyroid carcinoma with intermediate-to-high recurrence risk: a retrospective study

The predictive value of allele frequency (AF) of BRAF V600E and TERT mutations in papillary thyroid carcinoma (PTC) remains controversial. We aimed to investigate the AF of BRAF V600E and TERT mutations in intermediate-to-high risk PTC and their association between tumor invasiveness, prognosis, and other mutations. Probe hybridization capture and high-throughput sequencing were used to quantitatively test 40 gene loci in 94 intermediate-to-high recurrence risk PTC patients, combined with clinical characteristics and follow-up for retrospective analysis. BRAF V600E mutation AF was linked to a increased risk of thyroid capsule penetration, recurrence, and concurrent mutations. Concurrent mutations could lead to a worse prognosis and increased invasiveness. TERT promoter mutation frequently accompanied other mutations and resulted in a poorer prognosis. However, there was no clear association between the TERT mutation AF and tumor invasiveness or recurrence. The sensitivity and specificity of predicting recurrence in intermediate-to-high risk PTC with BRAF V600E mutation AF > 28.2% were 60 and 80%. Although genetic alterations in PTC can differ among different ethnicities, the AF of BRAF V600E and TERT mutations may be similar. The AF of BRAF V600E has the potential to be a novel indicator in predicting PTC invasiveness and prognosis.

A Phase 2 Study of Encorafenib in Combination with Binimetinib in Patients with Metastatic BRAF-Mutated Thyroid Cancer in Japan

Background: Driver mutations at BRAF V600 are frequently identified in papillary thyroid cancer and anaplastic thyroid cancer (ATC), in which BRAF inhibitors have shown clinical effectiveness. This Japanese phase 2 study evaluated the efficacy and safety of a BRAF inhibitor, encorafenib, combined with an MEK inhibitor, binimetinib, in patients with BRAF V600-mutated thyroid cancer. Methods: This phase 2, open-label, uncontrolled study was conducted at 10 institutions targeted patients with BRAF V600-mutated locally advanced or distant metastatic thyroid cancer not amenable to curative treatment who became refractory/intolerant to ≥1 previous vascular endothelial growth factor receptor-targeted regimen(s) or were considered ineligible for those. The primary endpoint was centrally assessed objective response rate (ORR). The secondary endpoints included duration of response (DOR), progression-free survival (PFS), overall survival (OS), and safety. Results: We enrolled 22 patients with BRAFV600E-mutated thyroid cancer: 17 had differentiated thyroid cancer (DTC), and 5 had ATC. At data cutoff (October 26, 2022), the median follow-up was 11.5 (range = 3.4-19.0) months. The primary endpoint of centrally assessed ORR was 54.5% (95% confidence interval [CI] 32.2-75.6; partial response in 12 patients and stable disease in 10). The ORRs in patients with DTC and ATC were 47.1% (8 of 17) and 80.0% (4 of 5), respectively. The medians for DOR and PFS by central assessment and for OS were not reached in the overall population, the DTC subgroup, or the ATC subgroup. At 12 months, the rate of ongoing response was 90.9%, and the PFS and OS rates were 78.8% and 81.8%, respectively. All patients developed ≥1 adverse events (AEs): grade 3 AEs in 6 patients (27.3%). No patients developed grade 4-5 AEs. The most common grade 3 AE was lipase increased (4 patients [18.2%]). Those toxicities were mostly manageable with appropriate monitoring and dose adjustment. Conclusions: Treatment with encorafenib plus binimetinib met the primary endpoint criteria and demonstrated clinical benefit in patients with BRAFV600E-mutated thyroid cancer regardless of its histological type, such as DTC or ATC, with no new safety concerns identified. Encorafenib plus binimetinib could thus be a new treatment option for BRAF V600-mutated thyroid cancer. Clinical Trial Registration number: Japan Registry of Clinical Trials: jRCT2011200018.

Pathogenesis of cancers derived from thyroid follicular cells

The genomic simplicity of differentiated cancers derived from thyroid follicular cells offers unique insights into how oncogenic drivers impact tumour phenotype. Essentially, the main oncoproteins in thyroid cancer activate nodes in the receptor tyrosine kinase-RAS-BRAF pathway, which constitutively induces MAPK signalling to varying degrees consistent with their specific biochemical mechanisms of action. The magnitude of the flux through the MAPK signalling pathway determines key elements of thyroid cancer biology, including differentiation state, invasive properties and the cellular composition of the tumour microenvironment. Progression of disease results from genomic lesions that drive immortalization, disrupt chromatin accessibility and cause cell cycle checkpoint dysfunction, in conjunction with a tumour microenvironment characterized by progressive immunosuppression. This Review charts the genomic trajectories of these common endocrine tumours, while connecting them to the biological states that they confer.

Whole Exome-Wide Association Identifies Rare Variants in GALNT9 Associated with Middle Eastern Papillary Thyroid Carcinoma Risk

Papillary thyroid carcinoma (PTC) is the commonest thyroid cancer. The majority of inherited causes of PTC remain elusive. However, understanding the genetic underpinnings and origins remains a challenging endeavor. An exome-wide association study was performed to identify rare germline variants in coding regions associated with PTC risk in the Middle Eastern population. By analyzing exome-sequencing data from 249 PTC patients (cases) and 1395 individuals without any known cancer (controls), GALNT9 emerged as being strongly associated with rare inactivating variants (RIVs) (4/249 cases vs. 1/1395 controls, OR = 22.75, p = 5.09 × 10-5). Furthermore, three genes, TRIM40, ARHGAP23, and SOX4, were enriched for rare damaging variants (RDVs) at the exome-wide threshold (p < 2.5 × 10-6). An additional seven genes (VARS1, ZBED9, PRRC2A, VWA7, TRIM31, TRIM40, and COL8A2) were associated with a Middle Eastern PTC risk based on the sequence kernel association test (SKAT). This study underscores the potential of GALNT9 and other implicated genes in PTC predisposition, illuminating the need for large collaborations and innovative approaches to understand the genetic heterogeneity of PTC predisposition.

Mapping the Complex Genetic Landscape of Human Neurons

When somatic cells acquire complex karyotypes, they are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons have been limited by relatively small sample sizes. Here, we developed an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We applied this approach to 2,125 frontal cortical neurons from a neurotypical human brain. This approach identified 226 CNV neurons, as well as a class of CNV neurons with complex karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we found that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contained fewer, but longer, genes.

Tracking the evolution of esophageal squamous cell carcinoma under dynamic immune selection by multi-omics sequencing

Intratumoral heterogeneity (ITH) has been linked to decreased efficacy of clinical treatments. However, although genomic ITH has been characterized in genetic, transcriptomic and epigenetic alterations are hallmarks of esophageal squamous cell carcinoma (ESCC), the extent to which these are heterogeneous in ESCC has not been explored in a unified framework. Further, the extent to which tumor-infiltrated T lymphocytes are directed against cancer cells, but how the immune infiltration acts as a selective force to shape the clonal evolution of ESCC is unclear. In this study, we perform multi-omic sequencing on 186 samples from 36 primary ESCC patients. Through multi-omics analyses, it is discovered that genomic, epigenomic, and transcriptomic ITH are underpinned by ongoing chromosomal instability. Based on the RNA-seq data, we observe diverse levels of immune infiltrate across different tumor sites from the same tumor. We reveal genetic mechanisms of neoantigen evasion under distinct selection pressure from the diverse immune microenvironment. Overall, our work offers an avenue of dissecting the complex contribution of the multi-omics level to the ITH in ESCC and thereby enhances the development of clinical therapy.

Systematic investigation of the prognostic impact of clonal status of somatic mutations across multiple cancer types

Tumors are genetically heterogeneous and many mutations are actually present in subclonal populations. The clonal status of mutations is valuable for accurate prognosis, clinical management. The aim of this study was to identify the clonal status of somatic mutations and systematically evaluate their prognostic values across various cancer types. We totally identified 227 clonal and 432 subclonal mutations contributed to prognosis and demonstrated the importance of clonal status in improving mutation-related clinical guidance. We further developed a customized multi-step approach to identify gene-specific prognostic patterns of clonal status at pan-cancer level and found some cancer-specific prognostic patterns. The 'subclonal-dependent risk' subpattern was one of the most common subpatterns, it usually accompanied by high genomic in-stability and high extent of intra-tumor heterogeneity and could be used to improve the accuracy of prognostic analysis. Our results revealed the importance of clonal status, especially subclonal mutation in clinical survival.

APOBEC SBS13 Mutational Signature-A Novel Predictor of Radioactive Iodine Refractory Papillary Thyroid Carcinoma

Standard surgery followed by radioactive iodine (131I, RAI) therapy are not curative for 5−20% of papillary thyroid carcinoma (PTC) patients with RAI refractory disease. Early predictors indicating therapeutic response to RAI therapy in PTC are yet to be elucidated. Whole-exome sequencing was performed (at median depth 198x) on 66 RAI-refractory and 92 RAI-avid PTCs with patient-matched germline. RAI-refractory tumors were significantly associated with distinct aggressive clinicopathological features, including positive surgical margins (p = 0.016) and the presence of lymph node metastases at primary diagnosis (p = 0.012); higher nonsilent tumor mutation burden (p = 0.011); TERT promoter (TERTp) mutation (p < 0.0001); and the enrichment of the APOBEC-related single-base substitution (SBS) COSMIC mutational signatures 2 (p = 0.030) and 13 (p < 0.001). Notably, SBS13 (odds ratio [OR] 30.4, 95% confidence intervals [CI] 1.43−647.22) and TERTp mutation (OR 41.3, 95% CI 4.35−391.60) were revealed to be independent predictors of RAI refractoriness in PTC (p = 0.029 and 0.001, respectively). Although SBS13 and TERTp mutations alone highly predicted RAI refractoriness, when combined, they significantly increased the likelihood of predicting RAI refractoriness in PTC. This study highlights the APOBEC SBS13 mutational signature as a novel independent predictor of RAI refractoriness in a distinct subgroup of PTC.

Dynamic Cancer Cell Heterogeneity: Diagnostic and Therapeutic Implications

Though heterogeneity of cancers is recognized and has been much discussed in recent years, the concept often remains overlooked in different routine examinations. Indeed, in clinical or biological articles, reviews, and textbooks, cancers and cancer cells are generally presented as evolving distinct entities rather than as an independent heterogeneous cooperative cell population with its self-oriented biology. There are, therefore, conceptual gaps which can mislead the interpretations/diagnostic and therapeutic approaches. In this short review, we wish to summarize and discuss various aspects of this dynamic evolving heterogeneity and its biological, pathological, clinical, diagnostic, and therapeutic implications, using thyroid carcinoma as an illustrative example.

Intratumor Heterogeneity as a Prognostic Factor in Solid Tumors: A Systematic Review and Meta-Analysis

Background

The landscape of intratumor heterogeneity (ITH) is present from the tumor evolution. ITH is a promising clinical indicator, but the association between ITH and prognosis remains controversial. Therefore, a meta-analysis was performed to explore whether ITH can serve as a valuable prognostic indicator in solid tumors.

Methods

All included studies were from PubMed, Embase, Cochrane, and Web of Science databases up to October 10, 2020. Studies based on ITH with available prognostic information were included. Three researchers independently completed study selection and data extraction following PRISMA guidelines. The random-effect model was used for synthesis. Hazard ratio (HR) and 95% confidence intervals (CI) were used with the endpoint defined by overall survival (OS), disease-specific survival (DFS), and progression-free survival (PFS).

Results

A total of 9,804 solid tumor patients from 21 studies were included. Analysis of specific cancers in the TCGA database showed similar results based on different ITH assessment methods, which provided the logical support for data consolidation. Available evidence revealed a negative relationship between ITH and prognosis for a specific cancer (such as lung cancer). However, the OS results from 14 tumor types showed that high ITH associated with shorter survival time [HR 1.65 (95% CI, 1.42-1.91)]. PFS and DFS analyses showed similar results [HR 1.89 (95% CI, 1.41-2.54) and HR 1.87 (95% CI, 1.15-3.04)] in general. The status of tumor metastasis and sampling models were not the confounding factors.

Conclusions

High ITH is associated with worse prognosis in many solid tumors in general although this association was absent for some cancers. ITH is expected to be a promising clinical prognostic factor for the improvement of assessment, treatment, and surveillance strategy.

Gene Mutation Analysis in Papillary Thyroid Carcinoma Using a Multi-Gene Panel in China

Purpose

To detect low-frequency mutation in the 57 genes of small panels that are associated with developing thyroid cancer in papillary thyroid carcinoma (PTC) patients and provide patients with precise-targeted therapy.

Patients and Methods

This study included 144 patients diagnosed with PTC who underwent total thyroidectomy and lymph node dissection in the central area of the neck between May 2017 and October 2018. We performed ultra-deep sequencing of 57 genes from 144 patients and detected the 57 genes mutations with bioinformatics.

Results

There were 698 mutations in 45 genes from 138 PTC patients. A high frequency of mutations was detected in the RBM10 gene (44%) and TERT (43%), and some hotspot mutations, such as RBM10:p.E119D and TERT:p.P112fs, were also found.

Conclusion

Ultra-deep sequencing of small gene panels can find some low-frequency mutation genes, which can provide targeted therapy for patients.

Combined Mutational and Clonality Analyses Support the Existence of Intra-Tumor Heterogeneity in Papillary Thyroid Cancer

Despite its potential clinical impact, intra-tumor genetic heterogeneity (ITH) has been scantly investigated in papillary thyroid cancer (PTC). We studied ITH in PTC by combining, for the first time, data derived from the evaluation of the normalized allelic frequencies (NAF) of the mutation/s, using a customized MassARRAY panel, and those obtained by the HUMARA clonality assay. Among tumors with a single mutation, 80% of cases with NAF 50 ± 5% were clonal, consistent with the presence of a single mutated clone, while 20% of cases showed a polyclonal pattern, suggesting the presence of the same mutation in two or more clones. Differently, all cases with NAF < 45% were polyclonal. Among tumors with double mutation, cases with both mutations showing NAF 50 ± 5% were monoclonal, consistent with the presence of a single clone harboring both mutations. On the other hand, all cases with double mutation at NAF < 45% were polyclonal, indicating the presence of two clones with different mutations. Finally, no significant differences in the clinico-pathological characteristics were found between monoclonal and polyclonal tumors. In conclusion, the present study adds insights into the concept of ITH in PTC, which warrants attention because the occurrence of this phenomenon is likely to affect the response to targeted drugs.

Comparison of Genomic Characterization in Upper Tract Urothelial Carcinoma and Urothelial Carcinoma of the Bladder

BACKGROUND

Different genomic characterization in urothelial carcinoma (UC) by site of origin may imply contrasting therapeutic opportunities and pathogenetic mechanisms. The aim of this study was to investigate whether differences between upper tract UC (UTUC) and UC of the bladder (UCB) result from intrinsic biological diversity.

MATERIALS AND METHODS

We prospectively sequenced 118 tumors and matched blood DNA from Chinese patients with UC using next-generation sequencing techniques, including 45 UTUC and 73 UCB. Two hundred twenty-six patients with UTUC and 350 patients with UCB for The Cancer Genome Atlas were acquired from the cbioportal.

RESULTS

There were marked disparities in the mutational landscape for UC according to race and site of origin. Signature 22 for exposure to aristolochic acid was only observed in the UTUC cohort. Conversely, signature 6 for defective DNA mismatch repair only existed in the UCB cohort. Compared with UCB, UTUC had higher clonal and subclonal mutation numbers. TP53, PIK3CA, and FGFR3 mutations may be the driver genes for UTUC, whereas for UCB, the driver gene may be BRCA1. Patients with UTUC had lower PD-L1 than those with UCB. There was no significant difference in the number of DDR mutations, copy number variation counts, and tumor mutational burden between UTUC and UCB.

CONCLUSION

UTUC and UCB exhibit significant differences in the prevalence of genomic landscape and carcinogenesis. Consequently, molecular subtypes differ according to location, and these results may imply the site-specific management of patients with urothelial carcinoma. Mutational signature may be used as a screening tool to assist clinical differential diagnosis between UTUC and UCB.

IMPLICATIONS FOR PRACTICE

This study's findings lay the foundation for a deeper understanding of distinct molecular mechanisms and similar treatment opportunities between upper tract urothelial carcinoma (UTUC) and urothelial carcinoma of the bladder (UCB) and had important implications for the site-specific management of patients with urothelial carcinoma. A comprehensive understanding of the biology of UTUC and UCB is needed to identify new drug targets in order to improve clinical outcomes.

Histology-based molecular profiling improves mutation detection for advanced thyroid cancer

Advanced cancers frequently show histologic and molecular intratumoral heterogeneity. Therefore, we comprehensively characterized advanced, metastatic, radioiodine-resistant (RAIR) thyroid carcinomas at the molecular level in the context of histologic heterogeneity with the aim to identify potentially actionable mutations that may guide the use of specific tyrosine kinase inhibitor (TKI) treatment. Whole exome sequencing (WES) was applied to 29 macrodissected tissue samples of histologically heterogeneous and homogeneous areas, lymph node and lung metastases from six clinically and histologically well-characterized metastatic RAIR thyroid cancer patients with structural incomplete response to treatment. WES data were analyzed to identify potential driver mutations in oncogenic pathways, copy number alterations, microsatellite instability, mutant-allele tumor heterogeneity, and the relevance of histologic heterogeneity to molecular profiling. In addition to known driver mutations in BRAF, NRAS, EIF1AX, NCOA4-RET, and TERT, further potentially actionable drivers were identified in AKT1, ATM, E2F1, HTR2A, and MLH3. The analysis of the evolutionary history of the mutations and the reconstruction of the molecular phylogeny of the cancers show a remarkable association between histologic and molecular heterogeneity. A comprehensive molecular analysis of the primary tumor guided by histologic analysis may help to better stratify patients for precision medicine approaches. Given the association between the molecular and the histologic heterogeneity, the selection of tumor samples for molecular analysis should be based on meticulous histologic evaluation of the entire tumor.

Spatial Distribution Patterns of Clinically Relevant TERT Promoter Mutations in Follicular Thyroid Tumors of Uncertain Malignant Potential: Advantages of the Digital Droplet PCR Technique

In thyroid carcinomas, telomerase reverse transcriptase (TERT) promoter mutations C228T and C250T predict an unfavorable clinical outcome. The analysis is particularly valuable when assessing histologically equivocal follicular thyroid tumors of uncertain malignant potential (FT-UMPs). Given recent findings of TERT promoter mutational heterogeneity in thyroid cancer, we determined the frequency of this phenomenon in FT-UMPs and minimally invasive follicular thyroid carcinomas. DNA was extracted from several tissue blocks from 16 FT-UMPs as well as 10 minimally invasive follicular thyroid carcinomas, and interrogated using Sanger sequencing as well as digital droplet PCR (ddPCR). Mutational heterogeneity was observed by Sanger sequencing in four of seven (57%) FT-UMPs. In two FT-UMPs with C228T mutations, analyses of additional blocks gave wild-type results using Sanger sequencing in one or several blocks interrogated, whereas ddPCR found low-frequency C228T mutations in one of these fractions. In two additional FT-UMPs with the C228T and C250T mutation, respectively, sequencing of additional blocks revealed the opposite mutation. Moreover, in the C250T mutated area in one of these tumors, the ddPCR displayed a co-occurring C228T mutation that failed detection through Sanger sequencing. To conclude, most TERT promoter mutated FT-UMPs display mutational heterogeneity when analyzed by Sanger sequencing, thereby emphasizing the importance of the tissue sampling process. The ddPCR technique might overcome this phenomenon because of increased sensitivity and should be considered for clinical screening purposes.

The clonal origin of multifocal papillary thyroid cancer: intrathyroidal spread or independent tumors?

Multifocality is a common finding in papillary thyroid cancer but the molecular pathogenesis, prognosis and management of multifocal papillary thyroid cancer are debated. The clonal origin of multifocal papillary thyroid cancer represents a controversial aspect, as two opposite viewpoints have been proposed: independent origin or intraglandular spread. Different approaches have been used for inferring the clonality of multifocal papillary thyroid cancer, including X-chromosome inactivation, mutational analysis, determination of loss of heterozygosity and, more recently, next-generation sequencing. Next-generation sequencing, able to provide information on genetic heterogeneity and phylogenetic evolution in multifocal tumors, represents the most reliable approach. While most evidences indicated an independent origin of multifocal papillary thyroid cancer, a minority of studies suggested that multifocal papillary thyroid tumors might be monoclonally derived. This discrepancy may reflect technical limitations; nevertheless, studies based on next-generation sequencing indicated that both independent and clonal origins are possible. The co-existence of multiple tumors implies a high degree of genetic heterogeneity, which may influence the best and targeted therapeutic strategy. On the other hand, intrathyroidal dissemination may indicate metastatic potential of the dominant tumor, thereby prompting more aggressive treatments. In conclusion, data available in the literature indicated that multifocal papillary thyroid cancer may derived from both intraglandular spread and independent tumor foci. The understanding of the clonal origin of multifocal papillary thyroid tumors might represent an important issue in patient treatment.

Intratumoral Genetic Heterogeneity in Papillary Thyroid Cancer: Occurrence and Clinical Significance

Intratumoral heterogeneity (ITH) refers to a subclonal genetic diversity observed within a tumor. ITH is the consequence of genetic instability and accumulation of genetic alterations, two mechanisms involved in the progression from an early tumor stage to a more aggressive cancer. While this process is widely accepted, the ITH of early stage papillary thyroid carcinoma (PTC) is debated. By different genetic analysis, several authors reported the frequent occurrence of PTCs composed of both tumor cells with and without RET/PTC or BRAF^V600E genetic alterations. While these data, and the report of discrepancies in the genetic pattern between metastases and the primary tumor, demonstrate the existence of ITH in PTC, its extension and biological significance is debated. The ITH takes on a great significance when involves oncogenes, such as RET rearrangements and BRAF^V600E as it calls into question their role of driver genes. ITH is also predicted to play a major clinical role as it could have a significant impact on prognosis and on the response to targeted therapy. In this review, we analyzed several data indicating that ITH is not a marginal event, occurring in PTC at any step of development, and suggesting the existence of unknown genetic or epigenetic alterations that still need to be identified.