Abstract
Papillary thyroid carcinoma (PTC) is heterogeneous and its molecular characteristics remain elusive. We integrated transcriptomic sequencing, genomic analysis and clinicopathologic information from 582 tissue samples of 216 PTC and 75 benign thyroid nodule (BTN) patients. We discovered four subtypes of PTC including Immune-enriched Subtype, BRAF-enriched Subtype, Stromal Subtype and CNV-enriched Subtype. Molecular subtypes were validated in an external cohort of 497 PTC cases from the TCGA. Tumors in the Immune-enriched Subtype showed higher immune infiltration and overexpression of immune checkpoints, whilst BRAF-enriched Subtype showed a higher tendency for extrathyroidal extension and more advanced TNM stage. Key oncogenes including LRRK2, SLC34A2, MUC1, FOXQ1 and KRT19 were overexpressed and enriched in oncogenic MAPK and PI3K/AKT signaling pathways in BRAF-enriched subtype. Further analysis of BRAF-enriched Subtype identified three subclasses with different degrees of malignancies. We also uncovered the molecular link of the initiation and progression from BTN to subtypes of PTC using trajectory analysis. Moreover, a 20-gene expression signature was generated for differential diagnosis of PTC from BTN patients. Together, our work identified previously unreported molecular subtypes of PTC, offering opportunities to stratify patients into optimal treatment plans based on molecular subtyping.