Knauf, Jeffrey A.

Abstract Purpose: The determinants of response or resistance to radioiodine (RAI) are unknown. We aimed to identify genomic and transcriptomic factors associated with structural responses to RAI treatment of metastatic thyroid cancer, which occur infrequently, and to test whether high MAPK pathway output was associated with RAI refractoriness. Experimental Design: Exceptional response to RAI was defined as reduction of tumor volume based on RECIST v1.1. We performed a retrospective case–control study of genomic and transcriptomic characteristics of exceptional responders (ER; n = 8) versus nonresponders (NR; n = 16) matched by histologic type and stage at presentation on a 1:2 ratio. Results: ER are enriched for mutations that activate MAPK through RAF dimerization (RAS, class 2 BRAF, RTK fusions), whereas NR are associated with BRAFV600E, which signals as a monomer and is unresponsive to negative feedback. ER have a lower MAPK transcriptional output and a higher thyroid differentiation score (TDS) than NR (P < 0.05). NR are enriched for 1q-gain (P < 0.05) and mutations of genes regulating mRNA splicing and the PI3K pathway. BRAFV600E tumors with 1q-gain have a lower TDS than BRAFV600E/1q-quiet tumors and transcriptomic signatures associated with metastatic propensity. Conclusions: ER tumors have a lower MAPK output and higher TDS than NR, whereas NR have a high frequency of BRAFV600E and 1q-gain. Molecular profiling of thyroid cancers and further functional validation of the key findings discriminating ER from NR may help predict response to RAI therapy.

Abstract Purpose: Thyroid cancer cell lines are valuable models but have been neglected in pancancer genomic studies. Moreover, their misidentification has been a significant problem. We aim to provide a validated dataset for thyroid cancer researchers. Experimental Design: We performed next-generation sequencing (NGS) and analyzed the transcriptome of 60 authenticated thyroid cell lines and compared our findings with the known genomic defects in human thyroid cancers. Results: Unsupervised transcriptomic analysis showed that 94% of thyroid cell lines clustered distinctly from other lineages. Thyroid cancer cell line mutations recapitulate those found in primary tumors (e.g., BRAF, RAS, or gene fusions). Mutations in the TERT promoter (83%) and TP53 (71%) were highly prevalent. There were frequent alterations in PTEN, PIK3CA, and of members of the SWI/SNF chromatin remodeling complex, mismatch repair, cell-cycle checkpoint, and histone methyl- and acetyltransferase functional groups. Copy number alterations (CNA) were more prevalent in cell lines derived from advanced versus differentiated cancers, as reported in primary tumors, although the precise CNAs were only partially recapitulated. Transcriptomic analysis showed that all cell lines were profoundly dedifferentiated, regardless of their derivation, making them good models for advanced disease. However, they maintained the BRAFV600E versus RAS-dependent consequences on MAPK transcriptional output, which correlated with differential sensitivity to MEK inhibitors. Paired primary tumor-cell line samples showed high concordance of mutations. Complete loss of p53 function in TP53 heterozygous tumors was the most prominent event selected during in vitro immortalization. Conclusions: This cell line resource will help inform future preclinical studies exploring tumor-specific dependencies.