Whole-Exome Sequencing of Matched Primary and Metastatic Papillary Thyroid Cancer

Arid4a Suppresses Breast Tumor Metastasis by Enhancing MTSS1 Expression via mRNA Stability

BACKGROUND

Tumor metastasis is one of the main causes of death in cancer patients; however, the mechanism controlling metastasis is unclear. The posttranscriptional regulation of metastasis-related genes mediated by AT-rich interactive domain-containing protein 4A (Arid4a), an RNA-binding protein (RBP), has not been elucidated.

METHODS

Bioinformatic analysis, qRT-PCR, immunohistochemistry, and immunoblotting were employed to determine the expression of Arid4a in breast tumor tissues and its association with the survival of cancer patients. In vitro and in vivo cellular experiments were used to assess the function of Arid4a in breast tumor metastasis. PCR array, RNA immunoprecipitation (RIP), luciferase, mRNA stability, RIP-ChIP, and EMSA were conducted to elucidate the potential mechanism of Arid4a.

RESULTS

Reduced expression of Arid4a in breast tumor samples was detected via bioinformatic analyses and experimental methods. Low Arid4a expression was significantly correlated with poor prognosis in breast cancer patients. Gain-of-function and silencing experiments confirmed the inhibitory effect of Arid4a on tumor metastasis in vitro and in vivo. Mechanistically, Arid4a preferentially stabilizes metastasis-suppressing transcripts, including metastasis suppressor 1 (MTSS1), tissue inhibitor of metalloproteinase 2 (TIMP2), retinoblastoma 1 (Rb1), and phosphatase and tensin homolog (PTEN), through binding to a conserved structural RNA element localized in the 3' untranslated region (3'UTR). The Arid domain of Arid4a is required for its mRNA stabilization and metastasis inhibition. Notably, the expression of Arid4a and metastasis-suppressing genes was positively correlated in human breast tumor tissues.

CONCLUSIONS

Arid4a was confirmed to suppress breast tumor metastasis progression by stabilizing the transcripts of tumor metastasis-suppressing genes, suggesting that Arid4a might be a potential therapeutic target for breast cancer treatment.

Genomic Differences in Thyroid Cancers From Primary Sites Versus Distant Metastases in Individual Patients: A Clinical Perspective and Preliminary Report

BACKGROUND

Distant metastasis is a leading cause of thyroid cancer (TC)-related deaths. Genetic profiling is typically limited to one sample per patient due to cost and sampling-risk concerns. Differences between samples from thyroid and distant metastasis within individual patients are unclear.

METHODS

Patients with TC and distant metastasis were recruited for genetic analysis.

RESULTS

Using a TC-specific NGS panel, 66 specimens from 29 patients were analyzed, identifying 16 mutations and 4 fusions, including two novel fusions (FGFR2-SHTN1 and RFTN1-BRAF). Genetic alterations differed between primary and metastatic sites in nine patients (31%), predominantly in additional oncogenic alterations (89%). More genetic alterations were found at the primary site in three patients and metastatic sites in four. Distinct mutations were found in two patients. A longer time interval between specimen acquisitions was significantly associated with genetic discrepancies (p = 0.032).

CONCLUSION

Patterns of genetic discrepancies between primary and metastatic TC vary, offering valuable insights for clinical practice.

TERT upstream promoter methylation regulates TERT expression and acts as a therapeutic target in TERT promoter mutation-negative thyroid cancer

BACKGROUND

DNA hypermethylation and hotspot mutations were frequently observed in the upstream and core promoter of telomerase reverse transcriptase (TERT), respectively, and they were associated with increased TERT expression and adverse clinical outcomes in thyroid cancer. In TERT promoter mutant cancer cells, the hypomethylated TERT mutant allele was active and the hypermethylated TERT wild-type allele was silenced. However, whether and how the upstream promoter methylation regulates TERT expression in TERT mutation-negative cells were largely unknown.

METHODS

DNA demethylating agents 5-azacytidine and decitabine and a genomic locus-specific demethylation system based on dCas9-TET1 were used to assess the effects of TERT upstream promoter methylation on TERT expression, cell growth and apoptosis of thyroid cancer cells. Regulatory proteins binding to TERT promoter were identified by CRISPR affinity purification in situ of regulatory elements (CAPTURE) combined with mass spectrometry. The enrichments of selected regulatory proteins and histone modifications were evaluated by chromatin immunoprecipitation.

RESULTS

The level of DNA methylation at TERT upstream promoter and expression of TERT were significantly decreased after treatment with 5-azacytidine or decitabine in TERT promoter wild-type thyroid cancer cells. Genomic locus-specific demethylation of TERT upstream promoter induced TERT downregulation, along with cell apoptosis and growth inhibition. Consistently, demethylating agents sharply inhibited the growth of thyroid cancer cells harboring hypermethylated TERT but had little effect on cells with TERT hypomethylation. Moreover, we identified that the chromatin remodeling protein CHD4 binds to methylated TERT upstream promoter and promotes its transcription by suppressing the enrichment of H3K9me3 and H3K27me3 at TERT promoter.

CONCLUSIONS

This study uncovered the mechanism of promoter methylation mediated TERT activation in TERT promoter mutation-negative thyroid cancer cells and indicated TERT upstream promoter methylation as a therapeutic target for thyroid cancer.

Deep learning model accurately classifies metastatic tumors from primary tumors based on mutational signatures

Metastatic propagation is the leading cause of death for most cancers. Prediction and elucidation of metastatic process is crucial for the treatment of cancer. Even though somatic mutations have been linked to tumorigenesis and metastasis, it is less explored whether metastatic events can be identified through genomic mutational signatures, which are concise descriptions of the mutational processes. Here, we developed MetaWise, a Deep Neural Network (DNN) model, by applying mutational signatures as input features calculated from Whole-Exome Sequencing (WES) data of TCGA and other metastatic cohorts. This model can accurately classify metastatic tumors from primary tumors and outperform traditional machine learning (ML) models and a deep learning (DL) model, DiaDeL. Signatures of non-coding mutations also have a major impact on the model's performance. SHapley Additive exPlanations (SHAP) and Local Surrogate (LIME) analyses identify several mutational signatures which are directly correlated to metastatic spread in cancers, including APOBEC-mutagenesis, UV-induced signatures, and DNA damage response deficiency signatures.

The epigenetic factor CHD4 contributes to metastasis by regulating the EZH2/β-catenin axis and acts as a therapeutic target in ovarian cancer

BACKGROUND

The overall survival rate of patients with advanced ovarian cancer (OC) has remained static for several decades. Advanced ovarian cancer is known for its poor prognosis due to extensive metastasis. Epigenetic alterations contribute to tumour progression and therefore are of interest for potential therapeutic strategies.

METHODS

Following our previous study, we identified that CHD4, a chromatin remodelling factor, plays a strong role in ovarian cancer cell metastasis. We investigated the clinical significance of CHD4 through TCGA and GEO database analyses and explored the effect of CHD4 expression modulation and romidepsin treatment on the biological behaviour of ovarian cancer through CCK-8 and transwell assays. Bioluminescence imaging of tumours in xenografted mice was applied to determine the therapeutic effect of romidepsin. GSEA and western blotting were used to screen the regulatory mechanism of CHD4.

RESULTS

In ovarian cancer patient specimens, high CHD4 expression was associated with a poor prognosis. Loss of function of CHD4 in ovarian cancer cells induced suppression of migration and invasion. Mechanistically, CHD4 knockdown suppressed the expression of EZH2 and the nuclear accumulation of β-catenin. CHD4 also suppressed the metastasis of ovarian cancer cells and prevented disease progression in a mouse model. To inhibit the functions of CHD4 that are mediated by histone deacetylase, we evaluated the effect of the HDAC1/2 selective inhibitor romidepsin. Our findings indicated that treatment with romidepsin suppressed the progression of metastases in vitro and in vivo.

CONCLUSIONS

Collectively, our results uncovered an oncogenic function of CHD4 in ovarian cancer and provide a rationale for clinical trials of romidepsin in ovarian cancer patients.

Thyroglobulin expression, Ki-67 index, and lymph node ratio in the prognostic assessment of papillary thyroid cancer

The clinical significance of thyroglobulin (Tg) expression in papillary thyroid cancer (PTC) has not been systematically explored in relation to the Ki-67 index, lymph node ratio (LNR), or other conventional prognostic predictors. In this retrospective study of 327 patients with PTC, we investigated the immunohistochemical expression of Tg in both primary tumors and their matching lymph node metastases in relation to the Ki-67 index, LNR, and clinical data. Tumoral Tg immunoreactivity was inversely correlated to the Ki-67 index and tumor recurrence. The Ki-67 index was higher in lymph node metastases (mean 4%) than in the primary tumors (mean 3%). Reduced Tg expression, estimated as 0-25% Tg positive tumor cells, was more common in lymph node metastases compared to primary tumors. In addition to advanced metastatic burden (defined as N1b stage and LNR ≥ 21%), low Tg expression (0-25% positive tumor cells) in lymph node metastases had a significant prognostic impact with shorter recurrence-free survival. These findings support the potential value of histopathological assessment of Tg expression and Ki-67 index in lymph node metastases as complementary predictors to anticipate the prognosis of PTC patients better.

In silico Screening and Validation of Achyranthes aspera as a Potential Inhibitor of BRAF and NRAS in Controlling Thyroid Cancer

BACKGROUND

Thyroid carcinoma (THCA) is one of the most prevalent endocrine tumors, accounting for 3.4% of all cancers diagnosed annually. Single Nucleotide Polymorphisms (SNPs) are the most prevalent genetic variation associated with thyroid cancer. Understanding thyroid cancer genetics will enhance diagnosis, prognosis, and treatment.

METHODS

This TCGA-based study analyzes thyroid cancer-associated highly mutated genes through highly robust in silico techniques. Pathway, gene expression, and survival studies were performed on the top 10 highly mutated genes (BRAF, NRAS, TG, TTN, HRAS, MUC16, ZFHX3, CSMD2, EIFIAX, SPTA1). Novel natural compounds from Achyranthes aspera Linn were discovered to target two highly mutated genes. The natural compounds and synthetic drugs used to treat thyroid cancer were subjected to comparative molecular docking against BRAF and NRAS targets. The ADME characteristics of Achyranthes aspera Linn compounds were also investigated.

RESULTS

The gene expression analysis revealed that the expression of ZFHX3, MCU16, EIF1AX, HRAS, and NRAS was up-regulated in tumor cells while BRAF, TTN, TG, CSMD2, and SPTA1 were down-regulated in tumor cells. In addition, the protein-protein interaction network demonstrated that HRAS, BRAF, NRAS, SPTA1, and TG proteins have strong interactions with each other as compared to other genes. The ADMET analysis shows that seven compounds have druglike properties. These compounds were further studied for molecular docking studies. The compounds MPHY012847, IMPHY005295, and IMPHY000939 show higher binding affinity with BRAF than pimasertib. In addition, IMPHY000939, IMPHY000303, IMPHY012847, and IMPHY005295 showed a better binding affinity with NRAS than Guanosine Triphosphate.

CONCLUSION

The outcomes of docking experiments conducted on BRAF and NRAS provide insight into natural compounds with pharmacological characteristics. These findings indicate that natural compounds derived from plants as a more promising cancer treatment option. Thus, the results of docking investigations conducted on BRAF and NRAS substantiate the conclusions that the molecule possesses the most suited drug-like qualities. Compared to other compounds, natural compounds are superior, and they are also druggable. This demonstrates that natural plant compounds can be an excellent source of potential anti-cancer agents. The preclinical research will pave the road for a possible anti-cancer agent.

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.

Positive BRAFV600E mutation of primary tumor influences radioiodine avidity but not prognosis of papillary thyroid cancer with lung metastases

PURPOSE

This study investigated the relationship between BRAFV600E mutation of the primary tumor and radioiodine avidity in lung metastases (LMs) and then further evaluated the impact of BRAFV600E mutation and radioiodine avidity status on the prognosis of papillary thyroid cancer (PTC) with LMs.

METHODS

Ninety-four PTC patients with LMs after total thyroidectomy and cervical lymph node dissection between January 2012 and September 2021 were retrospectively included. All patients received BRAFV600E mutation examination of primary tumors and radioactive iodine (RAI) therapy. The therapeutic response was evaluated by Response Evaluation Criteria in Solid Tumors (RECIST) assessments (version 1.1). For patients with target lesions, the response was divided into complete response (CR), partial response (PR), stable disease (SD), and progressive disease (PD); for patients without target lesions, the response was divided into CR, non-CR/non-PD, and PD. In therapeutic response, PR and SD were classified as non-CR/non-PD for analysis. The chi-square test and logistic regression were used to analyze the impact factor on PD and mortality. Progression-free survival (PFS) and overall survival (OS) curves were constructed by the Kaplan-Meier method.

RESULTS

It was found that 21.2% (7/33) of patients with positive BRAFV600E mutation and 62.3% (38/61) of patients with negative BRAFV600E mutation had radioiodine-avid LMs (χ² = 14.484, p = 0.000). Patients with positive BRAFV600E mutation are more likely to lose radioiodine avidity; the odds ratios (ORs) were 5.323 (95% CI: 1.953-14.514, p = 0.001). Finally, 25 patients had PD, and six patients died; loss of radioiodine avidity was the independent predictor for PD, and the ORs were 10.207 (95% CI: 2.629-39.643, p = 0.001); BRAFV600E mutation status was not correlated with PD (p = 0.602), whether in the radioiodine avidity group (p = 1.000) or the non-radioiodine avidity group (p = 0.867). Similarly, BRAFV600E mutation status was not correlated with mortality; only loss of radioiodine avidity was the unfavorable factor associated with mortality in univariate analyses (p = 0.030).

CONCLUSION

Patients with LMs of PTC were more likely to lose radioiodine avidity when their primary tumor had positive BRAFV600E mutation; however, only radioiodine avidity and not BRAFV600E mutation status affected the clinical outcome of patients with lung metastatic PTC.

Identification of Prognostic Biomarkers in Papillary Thyroid Cancer and Developing Non-Invasive Diagnostic Models Through Integrated Bioinformatics Analysis

BACKGROUND

Papillary thyroid cancer (PTC) is the most frequent subtype of thyroid carcinoma, mainly detected in patients with benign thyroid nodules (BTN). Due to the invasiveness of accurate diagnostic tests, there is a need to discover applicable biomarkers for PTC. So, in this study, we aimed to identify the genes associated with prognosis in PTC. Besides, we performed a machine learning tool to develop a non-invasive diagnostic approach for PTC.

METHODS

For the study purposes, the miRNA dataset GSE130512 was downloaded from the GEO database and then analyzed to identify the common differentially expressed miRNAs in patients with non-metastatic PTC (nm-PTC)/metastatic PTC (m-PTC) compared with BTNs. The SVM was also applied to differentiate patients with PTC from those patients with BTN using the common DEMs. A protein-protein interaction network was also constructed based on the targets of the common DEMs. Next, functional analysis was performed, the hub genes were determined, and survival analysis was then executed.

RESULTS

A total of three common miRNAs were found to be differentially expressed among patients with nm-PTC/m-PTC compared with BTNs. In addition, it was established that the autophagosome maturation, ciliary basal body-plasma membrane docking, antigen processing as ubiquitination & proteasome degradation, and class I MHC mediated antigen processing & presentation are associated with the pathogenesis of PTC. Furthermore, it was illustrated that RPS6KB1, CCNT1, SP1, and CHD4 might serve as new potential biomarkers for PTC prognosis.

CONCLUSION

RPS6KB1, CCNT1, SP1, and CHD4 may be considered new potential biomarkers used for prognostic aims in PTC. However, performing validation tests is inevitable in the future.

Whole-Exome Sequencing in Papillary Microcarcinoma: Potential Early Biomarkers of Lateral Lymph Node Metastasis

BACKGROUND

Early identification of patients with high-risk papillary thyroid microcarcinoma (PTMC) that is likely to progress has become a critical challenge. We aimed to identify somatic mutations associated with lateral neck lymph node (LN) metastasis (N1b) in patients with PTMC.

METHODS

Whole-exome sequencing (WES) of 14 PTMCs with no LN metastasis (N0) and 13 N1b PTMCs was performed using primary tumors and matched normal thyroid tissues.

RESULTS

The mutational burden was comparable in N0 and N1b tumors, as the median number of mutations was 23 (range, 12 to 46) in N0 and 24 (range, 12 to 50) in N1b PTMC (P=0.918). The most frequent mutations were detected in PGS1, SLC4A8, DAAM2, and HELZ in N1b PTMCs alone, and the K158Q mutation in PGS1 (four patients, Fisher's exact test P=0.041) was significantly enriched in N1b PTMCs. Based on pathway analysis, somatic mutations belonging to the receptor tyrosine kinase-RAS and NOTCH pathways were most frequently affected in N1b PTMCs. We identified four mutations that are predicted to be pathogenic in four genes based on Clinvar and Combined Annotation-Dependent Depletion score: BRAF, USH2A, CFTR, and PHIP. A missense mutation in CFTR and a nonsense mutation in PHIP were detected in N1b PTMCs only, although in one case each. BRAF mutation was detected in both N0 and N1b PTMCs.

CONCLUSION

This first comprehensive WES analysis of the mutational landscape of N0 and N1b PTMCs identified pathogenic genes that affect biological functions associated with the aggressive phenotype of PTMC.

Targeted Therapy of Papillary Thyroid Cancer: A Comprehensive Genomic Analysis

Background

A limited number of targeted therapy options exist for papillary thyroid cancer (PTC) to date. Based on genetic alterations reported by the "The Cancer Genome Atlas (TCGA)", we explored whether PTC shows alterations that may be targetable by drugs approved by the FDA for other solid cancers.

Methods

Databases of the National Cancer Institute and MyCancerGenome were screened to identify FDA-approved drugs for targeted therapy. Target genes were identified using Drugbank. Genetic alterations were classified into conferring drug sensitivity or resistance using MyCancerGenome, CiViC, TARGET, and OncoKB. Genomic data for PTC were extracted from TCGA and mined for alterations predicting drug response.

Results

A total of 129 FDA-approved drugs with 128 targetable genes were identified. One hundred ninety-six (70%) of 282 classic, 21 (25%) of 84 follicular, and all 30 tall-cell variant PTCs harbored druggable alterations: 259 occurred in 29, 39 in 19, and 31 in 2 targetable genes, respectively. The BRAF V600 mutation was seen in 68% of classic, 16% of follicular variant, and 93% of tall-cell variant PTCs. The RET gene fusion was seen in 8% of classic PTCs, NTRK1 and 3 gene fusions in 3%, and other alterations in <2% of classic variant PTCs. Ninety-nine of 128 (77%) FDA-approved targetable genes did not show any genetic alteration in PTC. Beside selective and non-selective BRAF-inhibitors, no other FDA-approved drug showed any frequent predicted drug sensitivity (<10%).

Conclusion

Treatment strategies need to focus on resistance mechanisms to BRAF inhibition and on genetic alteration-independent alternatives rather than on current targeted drugs.

Identification and validation of potential novel biomarkers to predict distant metastasis in differentiated thyroid cancer

Background

Distant metastasis (DM) is not common in differentiated thyroid cancer (DTC). However, it is associated with a significantly poor prognosis. Early detection of high-risk DTC patients is difficult, and the molecular mechanism is still unclear. Therefore, the present study aims to establish a novel predictive model based on clinicopathological parameters and DM-related gene signatures to provide guidelines for clinicians in decision making.

Methods

Weighted gene co-expression network analysis (WGCNA) was performed to discover co-expressed gene modules and hub genes associated with DM. Univariate and multivariate analyses were carried out to identify independent clinicopathological risk factors based on The Cancer Genome Atlas (TCGA) database. An integrated nomogram prediction model was established. Finally, real hub genes were validated using the GSE60542 database and various thyroid cell lines.

Results

The midnightblue module was most significantly positively correlated with DM (R=0.56, P=9e-06) by as per WGCNA. DLX5 (AUC: 0.769), COX6B2 (AUC: 0.764), and LYPD1 (AUC: 0.760) were determined to be the real hub genes that play a crucial role in predicting DM. Meanwhile, univariate and multivariate analyses demonstrated that T-stage (OR, 15.03; 95% CI, 1.75-319.40; and P=0.024), histologic subtype (OR, 0.17; 95% CI, 0.03-0.92; and P=0.042) were the independent predictors of DM. Subsequently, a nomogram model was constructed based on gene signatures and independent clinical risk factors exhibited good performance. Additionally, the mRNA expressions of real hub genes in the GSE60542 dataset were consistent with TCGA.

Conclusions

The present study has provided a reliable model to predict DM in patients with DTC. This model is likely to serve as an individual risk assessment tool in therapeutic decision-making.

Pre-Therapeutic Measurements of Iodine Avidity in Papillary and Poorly Differentiated Thyroid Cancer Reveal Associations with Thyroglobulin Expression, Histological Variants and Ki-67 Index

Papillary thyroid cancer (PTC) and poorly differentiated thyroid cancer (PDTC) are treated with radioiodine to reduce recurrence and to treat the spread of disease. Adequate iodine accumulation in cancer tissue, iodine avidity, is important for treatment effect. This study investigated which clinical and histological tumour characteristics correlate with avidity. To quantify avidity in cancer tissue, tracer amounts of iodine-131 were given to 45 patients with cytologically confirmed thyroid cancer. At pathology grossing, representative samples of tumour and lymph nodes were taken and subjected to radioactivity quantification ex vivo to determine avidity. Afterwards, samples underwent extended pathology work-up and analysis. We found that tumoural Tg expression and Ki-67 index were correlated with avidity, whereas tumour size and pT stage were not. The histological variant of thyroid cancer was also correlated with iodine avidity. Variants associated with worse clinical prognoses displayed lower avidity than variants with better prognoses. This work provides new information on which tumours have low iodine avidity. Lower avidity in aggressive histological PTC variants may explain their overall poorer prognoses. Our findings also suggest that radioiodine dosage could be adapted to Tg expression, Ki-67 index or histological variant instead of pT stage, potentially improving the efficacy of radioiodine therapy.

Afirma Genomic Sequencing Classifier and Xpression Atlas Molecular Findings in Consecutive Bethesda III-VI Thyroid Nodules

CONTEXT

Broad genomic analyses among thyroid histologies have been described from relatively small cohorts.

OBJECTIVE

Investigate the molecular findings across a large, real-world cohort of thyroid fine-needle aspiration (FNA) samples.

DESIGN

Retrospective analysis of RNA sequencing data files.

SETTING

Clinical Laboratory Improvement Amendments laboratory performing Afirma Genomic Sequencing Classifier (GSC) and Xpression Atlas (XA) testing.

PARTICIPANTS

A total of 50 644 consecutive Bethesda III-VI nodules.

INTERVENTION

None.

MAIN OUTCOME MEASURES

Molecular test results.

RESULTS

Of 48 952 Bethesda III/IV FNAs studied, 66% were benign by Afirma GSC. The prevalence of BRAF V600E was 2% among all Bethesda III/IV FNAs and 76% among Bethesda VI FNAs. Fusions involving NTRK, RET, BRAF, and ALK were most prevalent in Bethesda V (10%), and 130 different gene partners were identified. Among small consecutive Bethesda III/IV sample cohorts with one of these fusions and available surgical pathology excision data, the positive predictive value of an NTRK or RET fusion for carcinoma or noninvasive follicular thyroid neoplasm with papillary-like nuclear features was >95%, whereas for BRAF and ALK fusions it was 81% and 67%, respectively. At least 1 genomic alteration was identified by the expanded Afirma XA panel in 70% of medullary thyroid carcinoma classifier-positive FNAs, 44% of Bethesda III or IV Afirma GSC suspicious FNAs, 64% of Bethesda V FNAs, and 87% of Bethesda VI FNAs.

CONCLUSIONS

This large study demonstrates that almost one-half of Bethesda III/IV Afirma GSC suspicious and most Bethesda V/VI nodules had at least 1 genomic variant or fusion identified, which may optimize personalized treatment decisions.

Longitudinal Genomic Evolution of Conventional Papillary Thyroid Cancer With Brain Metastasis

BACKGROUND

Brain metastasis is extremely rare but predicts dismal prognosis in papillary thyroid cancer (PTC). Dynamic evaluation of stepwise metastatic lesions was barely conducted to identify the longitudinal genomic evolution of brain metastasis in PTC.

METHOD

Chronologically resected specimen was analyzed by whole exome sequencing, including four metastatic lymph nodes (lyn 1-4) and brain metastasis lesion (BM). Phylogenetic tree was reconstructed to infer the metastatic pattern and the potential functional mutations.

RESULTS

Contrasting with lyn1, ipsilateral metastatic lesions (lyn2-4 and BM) with shared biallelic mutations of TSC2 indicated different genetic originations from multifocal tumors. Lyn 3/4, particularly lyn4 exhibited high genetic similarity with BM. Besides the similar mutational compositions and signatures, shared functional mutations (CDK4 R24C , TP53R342*) were observed in lyn3/4 and BM. Frequencies of these mutations gradually increase along with the metastasis progression. Consistently, TP53 knockout and CDK4 R24C introduction in PTC cells significantly decreased radioiodine uptake and increased metastatic ability.

CONCLUSION

Genomic mutations in CDK4 and TP53 during the tumor evolution may contribute to the lymph node and brain metastasis of PTC.

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.

Pan-genomic characterization of high-risk pediatric papillary thyroid carcinoma

Pediatric papillary thyroid carcinomas (pPTCs) are often indolent tumors with excellent long-term outcome, although subsets of cases are clinically troublesome and recur. Although it is generally thought to exhibit similar molecular aberrancies as their counterpart tumors in adults, the pan-genomic landscape of clinically aggressive pPTCs has not been previously described. In this study, five pairs of primary and synchronously metastatic pPTC from patients with high-risk phenotypes were characterized using parallel whole-genome and -transcriptome sequencing. Primary tumors and their metastatic components displayed an exceedingly low number of coding somatic mutations and gross chromosomal alterations overall, with surprisingly few shared mutational events. Two cases exhibited one established gene fusion event each (SQSTM1-NTRK3 and NCOA4-RET) in both primary and metastatic tissues, and one case each was positive for a BRAF V600E mutation and a germline truncating CHEK2 mutation, respectively. One single case was without apparent driver events and was considered as a genetic orphan. Non-coding mutations in cancer-associated regions were generally not present. By expressional analyses, fusion-driven primary and metastatic pPTC clustered separately from the mutation-driven cases and the sole genetic orphan. We conclude that pPTCs are genetically indolent tumors with exceedingly stable genomes. Several mutations found exclusively in the metastatic samples which may represent novel genetic events that drive the metastatic behavior, and the differences in mutational compositions suggest early clonal divergence between primary tumors and metastases. Moreover, an overrepresentation of mutational and expressional dysregulation of immune regulatory pathways was noted among fusion-positive pPTC metastases, suggesting that these tumors might facilitate spread through immune evasive mechanisms.

Intratumoral Heterogeneity in Differentiated Thyroid Tumors: An Intriguing Reappraisal in the Era of Personalized Medicine

Differentiated thyroid tumors (DTTs) are characterized by significant molecular variability in both spatial and temporal intra-tumoral heterogeneity (ITH), that could influence the therapeutic management. ITH phenomenon appears to have a relevant role in tumor growth, aggressive behavior and drug resistance. Accordingly, characteristics and consequences of ITH in DTTs should be better analyzed and understood in order to guide clinical practice, improving survival. Consequently, in the present review, we investigated morphological and molecular ITH of DTTs in benign, borderline neoplasms and in malignant entities, summarizing the most significant data. Molecular testing in DTTs documents a high risk for recurrence of cancer associated with BRAF^V600E, RET/PTC 1/3, ALK and NTRK fusions, while the intermediate risk may be related to BRAF^K601E, H/K/N RAS and PAX8/PPARγ. In addition, it may be suggested that tumor genotype is associated with peculiar phenotype.

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.

Whole-Exome Sequencing Reveals Novel Variations in Patients with Familial Von Hippel-Lindau Syndrome

OBJECTIVE

Von Hippel-Lindau (VHL) syndrome is a rare disease that occurs in an autosomal-dominant genetic pattern. Due to the high genetic variability of VHL diseases, current studies have limited clinical value. Moreover, casual genetic variations in patients with VHL syndrome are still unclear.

METHODS

Here, we performed whole-exome sequencing of 25 individuals to identify reliable disease-related variations. Systemic computational analysis was performed for variant detection, and Sanger sequencing was used to validate detected mutations.

RESULTS

Most of the known mutations in the VHL gene were observed in the studied population. In addition, a large fragment deletion in VHL exon 2 in the immediate family members of the last family was detected. This had not been reported earlier. Moreover, we identified 3 novel mutation sites in the MAP2K3 gene that may be involved in the occurrence and development of the VHL disease.

CONCLUSIONS

These results demonstrated that the heterogeneous nature of VHL syndrome and novel mutational signatures may help to improve the diagnostic ability of VHL syndrome.

Do Molecular Profiles of Primary Versus Metastatic Radioiodine Refractory Differentiated Thyroid Cancer Differ?

Management of metastatic radioiodine refractory differentiated thyroid cancer (DTC) can be a therapeutic challenge. Generally, little is known about the paired molecular profile of the primary tumor and the metastases and whether they harbor the same genetic abnormalities. The present study compared the molecular profile of paired tumor specimens (primary tumor/metastatic sites) from patients with radioiodine refractory DTC in order to gain insight into a possible basis for resistance to radioiodine. Twelve patients with radioiodine refractory metastases were studied; median age at diagnosis of 61 years (range, 25-82). Nine patients had papillary TC (PTC), one had follicular TC (FTC), and two had Hürthle cell TC (HTC). Distant metastases were present in the lungs (n = 10), bones (n = 4), and liver (n = 1). The molecular profiling of paired tumors was performed with a panel of 592 genes for Next Generation Sequencing, RNA-sequencing, and immunohistochemistry. Digital microfluidic PCR was used to investigate TERT promoter mutations. The genetic landscape of all paired sites comprised BRAF, NRAS, HRAS, TP53, ATM, MUTYH, POLE, and NTRK genes, including BRAF and NTRK fusions. BRAF V600E was the most common point mutation in the paired specimens (5/12). TERT promoter mutation C228T was detected in one case. PD-L1 expression at metastatic sites was highly positive (95%) for one patient with HTC. All specimens were stable for microsatellite instability testing, and the tumor mutation burden was low to intermediate. Therefore, the molecular profile of DTC primary and metastatic lesions can show heterogeneity, which may help explain some altered responses to therapeutic intervention.

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.

Progression and dormancy in metastatic thyroid cancer: concepts and clinical implications

Distant metastasis classically has been defined as a late-stage event in cancer progression. However, it has become clear that metastases also may occur early in the "lifetime" of a cancer and that they may remain stable at distant sites. This stability of metastatic cancer deposits has been termed "metastatic dormancy" or, as we term it, "metastatic progression dormancy" as the progression either may reflect growth of already existing metastases or new cancer spread. Biologically, dormancy is the presence of nongrowing, static metastatic cells that survive over time. Clinically, dormancy is defined by stability in tumor markers, imaging, and clinical course. Metastatic well-differentiated thyroid cancer offers an excellent tumor type to understand these processes for several reasons: (1) primary therapy often includes removal of the entire gland with ablation of residual normal tissue thereby removing one source for new metastases; (2) the presence of a sensitive biochemical and radiographic monitoring tests enabling monitoring of metastasis throughout the progression process; and (3) its tendency toward prolonged clinical dormancy that can last for years or decades be followed by progression. This latter factor provides opportunities to define therapeutic targets and/or markers of progression. In this review, we will discuss concepts of metastatic progression dormancy and the factors that drive both long-term stability and loss of dormancy with a focus on thyroid cancer.

Differences in Gene Expression Profile of Primary Tumors in Metastatic and Non-Metastatic Papillary Thyroid Carcinoma-Do They Exist?

Molecular mechanisms of distant metastases (M1) in papillary thyroid cancer (PTC) are poorly understood. We attempted to analyze the gene expression profile in PTC primary tumors to seek the genes associated with M1 status and characterize their molecular function. One hundred and twenty-three patients, including 36 M1 cases, were subjected to transcriptome oligonucleotide microarray analyses: (set A-U133, set B-HG 1.0 ST) at transcript and gene group level (limma, gene set enrichment analysis (GSEA)). An additional independent set of 63 PTCs, including 9 M1 cases, was used to validate results by qPCR. The analysis on dataset A detected eleven transcripts showing significant differences in expression between metastatic and non-metastatic PTC. These genes were validated on microarray dataset B. The differential expression was positively confirmed for only two genes: IGFBP3, (most significant) and ECM1. However, when analyzed on an independent dataset by qPCR, the IGFBP3 gene showed no differences in expression. Gene group analysis showed differences mainly among immune-related transcripts, indicating the potential influence of tumor immune infiltration or signal within the primary tumor. The differences in gene expression profile between metastatic and non-metastatic PTC, if they exist, are subtle and potentially detectable only in large datasets.

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.