Estrogen-related receptor alpha promotes thyroid tumor cell survival via a tumor subtype-specific regulation of target gene networks.

Mortalin (encoded by HSPA9) is a mitochondrial chaperone often overexpressed in cancer through as-yet-unknown mechanisms. By searching different RNA-sequencing datasets, we found that ESRRA is a transcription factor highly correlated with HSPA9 in thyroid cancer, especially in follicular, but not C cell-originated, tumors. Consistent with this correlation, ESRRA depletion decreased mortalin expression only in follicular thyroid tumor cells. Further, ESRRA expression and activity were relatively high in thyroid tumors with oncocytic characteristics, wherein ESRRA and mortalin exhibited relatively high functional overlap. Mechanistically, ESRRA directly regulated HSPA9 transcription through a novel ESRRA-responsive element located upstream of the HSPA9 promoter. Physiologically, ESRRA depletion suppressed thyroid tumor cell survival via caspase-dependent apoptosis, which ectopic mortalin expression substantially abrogated. ESRRA depletion also effectively suppressed tumor growth and mortalin expression in the xenografts of oncocytic or ESRRA-overexpressing human thyroid tumor cells in mice. Notably, our Bioinformatics analyses of patient data revealed two ESRRA target gene clusters that contrast oncocytic-like and anaplastic features of follicular thyroid tumors. These findings suggest that ESRRA is a tumor-specific regulator of mortalin expression, the ESRRA-mortalin axis has higher significance in tumors with oncocytic characteristics, and ESRRA target gene networks can refine molecular classification of thyroid cancer.

Clinicopathological Features and Molecular Signatures of Lateral Neck Lymph Node Metastasis in Papillary Thyroid Microcarcinoma.

BACKGRUOUND: The predictive factors for lateral neck lymph node metastasis (LLNM) in papillary thyroid microcarcinoma (PTMC) remain undetermined. This study investigated the clinicopathological characteristics, transcriptomes, and tumor microenvironment in PTMC according to the LLNM status. We aimed to identify the biomarkers associated with LLNM development. METHODS: We retrospectively reviewed the medical records of patients with PTMC from two independent institutions between 2018 and 2022 (n=597 and n=467). We compared clinicopathological features between patients without lymph node metastasis (N0) and those with LLNM (N1b). Additionally, laser capture microdissection and RNA sequencing were performed on primary tumors from both groups, including metastatic lymph nodes from the N1b group (n=30; 20 primary tumors and 10 paired LLNMs). We corroborated the findings using RNA sequencing data from 16 BRAF-like PTMCs from The Cancer Genome Atlas. Transcriptomic analyses were validated by immunohistochemical staining. RESULTS: Clinicopathological characteristics, such as male sex, multifocality, extrathyroidal extension, lymphatic invasion, and central node metastasis showed associations with LLNM in PTMCs. Transcriptomic profiles between the N0 and N1b PTMC groups were similar. However, tumor microenvironment deconvolution from RNA sequencing and immunohistochemistry revealed an increased abundance of tumor-associated macrophages, particularly M2 macrophages, in the N1b group. CONCLUSION: Patients with PTMC who have a male sex, multifocality, extrathyroidal extension, lymphatic invasion, and central node metastasis exhibited an elevated risk for LLNM. Furthermore, infiltration of M2 macrophages in the tumor microenvironment potentially supports tumor progression and LLNM in PTMCs.

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

BACKGRUOUND: Papillary thyroid carcinoma (PTC) can be classified into two distinct molecular subtypes, BRAF-like (BL) and RASlike (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.