Different Molecular Phenotypes of Progression in BRAF- and RAS-Like Papillary Thyroid Carcinoma

Abstract

Background

Papillary thyroid carcinoma (PTC) can be classified into two distinct molecular subtypes, BRAF-like (BL) and RAS-like (RL). However, the molecular characteristics of each subtype according to clinicopathological factors have not yet been determined. We aimed to investigate the gene signatures and tumor microenvironment according to clinicopathological factors, and to identify the mechanism of progression in BL-PTCs and RL-PTCs.

Methods

We analyzed RNA sequencing data and corresponding clinicopathological information of 503 patients with PTC from The Cancer Genome Atlas database. We performed differentially expressed gene (DEG), Gene Ontology, and molecular pathway enrichment analyses according to clinicopathological factors in each molecular subtype. EcoTyper and CIBERSORTx were used to deconvolve the tumor cell types and their surrounding microenvironment.

Results

Even for the same clinicopathological factors, overlapping DEGs between the two molecular subtypes were uncommon, indicating that BL-PTCs and RL-PTCs have different progression mechanisms. Genes related to the extracellular matrix were commonly upregulated in BL-PTCs with aggressive clinicopathological factors, such as old age (≥55 years), presence of extrathyroidal extension, lymph node metastasis, advanced tumor-node-metastasis (TNM) stage, and high metastasis-age-completeness of resection-invasion-size (MACIS) scores (≥6). Furthermore, in the deconvolution analysis of tumor microenvironment, cancer-associated fibroblasts were significantly enriched. In contrast, in RL-PTCs, downregulation of immune response and immunoglobulin-related genes was significantly associated with aggressive characteristics, even after adjusting for thyroiditis status.

Conclusion

The molecular phenotypes of cancer progression differed between BL-PTC and RL-PTC. In particular, extracellular matrix and cancer-associated fibroblasts, which constitute the tumor microenvironment, would play an important role in the progression of BL-PTC that accounts for the majority of advanced PTCs.