Tumor Mutation Burden Predicts Relapse in Papillary Thyroid Carcinoma With Changes in Genes and Immune Microenvironment

Abstract

Background:

The risk factors of papillary thyroid carcinoma (PTC) recurrence are meaningful for patients and clinicians. Tumor mutation burden (TMB) has been a biomarker for the effectiveness of immune checkpoint inhibitor (ICI) and prognosis in cancer. However, the role of TMB and its latent significance with immune cell infiltration in PTC are still unclear. Herein, we aimed to explore the effect of TMB on PTC prognosis.

Material and Methods:

RNA-seq and DNA-seq datasets of PTC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The Gene Ontology (GO) and gene set enrichment analysis (GSEA 4.0.1) were applied further to explore potential differences in PTC patients’ biological functions. The differentially expressed genes (DEGs) and immune microenvironment between the high and low TMB groups were determined.

Results:

TMB had the highest AUC score than other clinical indicators in ROC analysis on recurrence-free survival, and a higher TMB score was related to a worse prognosis. Further, GSEA showed a higher level of oxidative phosphorylation (OXPHOS) in the high TMB group, and four genes correlated with recurrence-free survival rate were identified. The abundance of CD8⁺ T cells and M1 macrophages in the high TMB group was significantly lower than that in the low TMB group.

Conclusions:

Our study found that TMB was a better predictor variable at evaluating the risk of PTC recurrence. Moreover, TMB-related genes conferred dramatically correlated prognosis, which was worth exploring in guiding postoperative follow-up and predicting recurrence for PTC patients.

Introduction

Papillary thyroid carcinoma (PTC) is the most common pathological type of thyroid carcinoma, with a good prognosis (1). Over the last four decades the occurrence of thyroid cancer was growing for the increased incidence of PTC, and the increase in mortality of thyroid cancer was mainly due to the increase in mortality of advanced-stage PTC (2).

Although the recurrence rate of some low-risk and intermediate-risk PTC patients undergoing hemithyroidectomy is low, the current risk scoring system has no predictive effect on these patients, which needs to be followed up to detect and treat the recurrence early (3). The aggressive clinical pathological, age ≥55, male, tumor size >4 cm, tall cell variant, and positive lymph node metastasis are independent prognostic factors, and there may continue to be the risk of cancer death and early recurrence (4–7). Besides, many biomarkers for predicting the recurrence of PTC are reported, such as SMOC2 (8), Cyclin D1 and C-myc (9), serum fibrinogen levels (10), and BRAF^V600E mutation (11). However, most biomarkers have limitations or insufficient clinical applicability (12). It is imperative to find reliable prognostic markers for PTC.

Tumor mutation burden (TMB), which is a quantitative measurement of the total number of somatic non-synonymous mutations in each coding area of the tumor genome (13), was first reported to predict the activity of immune checkpoint inhibitor therapies in multiple cancers (14). The detection of TMB is based on massive next-generation sequencing (NGS) or whole-exome sequencing (15). In specific cancer types and ethnic groups, the number of driver mutations is significantly positively correlated with cancer incidence rates (16). One lung cancer research identified that specific mutations could result in tumorigenesis or recurrence, and patients with a high burden of tumor-specific gene mutations had poor recurrence-free survival (RFS) (17). Another study showed that further validation and research on the elevated TMB and the potential impact of tumor-initiating mutations in clinically unfavorable prostate cancer could improve the outcomes of patients (18). Higher TMB means more tumor neoantigens are exposed, so TMB has also become a predictor of response rate to immune checkpoint inhibitor therapies (19).

In this study, we investigated the features of TMB and explored the relationship of TMB with PTC recurrence. DEGs and immune cells were analyzed to explore the underlying mechanism of PTC recurrence.

Methods

Data Source

The RNA-seq and single nucleotide polymorphism (SNP) data were all obtained from the TCGA pan-cancer cohort of the GDC data portal (HTTPS://portal.gdc.cancer.gov/). The RNA-seq comprised 512 tumor samples and 58 normal samples. Out of 512 tumor samples, 487 had complete SNP data, including simple nucleotide variation data (VarScan) and masked somatic mutation data. Basic clinical information included age, gender, AJCC-TNM stages, pathological subtype, survival outcomes, and follow-up time. The “maftools” software package, which could convert mutation data in Variant Call Format (VCF) format into visualization results (20), was used to describe the basic mutation characteristics of PTC samples. The mutation information of genes in each sample was displayed in the waterfall chart.

TMB Calculation

TMB generally refers to the number of somatic non-synonymous mutations per base pair (Mb) in a specific gene group region based on whole-exome modeling calculations. We selected SNP data from genome sequencing to calculate TMB. In our study, the mutation frequency and the number of variants in each sample were calculated through the Perl script based on the JAVA8 platform, as described in the previous literature (21). Then the PTC patients were divided into two groups according to the cut-off point (the median of the TMB value).

Estimation of TMB and Prognostic Analysis

Then Kaplan–Meier analysis was utilized to compare the difference in RFS between the high and low TMB groups, evaluated by the log-rank test. Chi-square test was used to assess the differences of clinical characteristics between high and low TMB groups. Subsequently, the correlations of TMB levels with clinical characteristics were appraised. The Wilcoxon rank-sum test was used to compare the differences between the two groups of clinical variables. P-value <0.05 was considered statistically significant.

DEGs and Functional Enrichment Analysis

DEGs were selected based on the high and low TMB groups of PTC patients through the “limma” software package (22), with log2 |fold change| >1 and False Discovery Rate (FDR) <0.05. The “heatmap” package was utilized to draw DEGs’ heatmap plot, and the “vioplot” package drew the volcano plot. Then we employed the Gene Ontology (GO) analyses through the “clusterProfiler” package (23) to analyze the possible biological processes of the overlapping DEGs. We also applied gene set enrichment analysis (GSEA 4.0.1) to further explore the potential differences in PTC patients’ biological functions.

Estimation of Tumor-Infiltrating Immune Cells

We employed the downloaded RNA-seq data through the deconvolution algorithm in the CIBERSORT website (http://cibersort.stanford.edu/). CIBERSORT is a general calculation method that could combine support vector regression and prior knowledge of purified leukocyte subsets’ composition profile to accurately estimate the composition of 22 immune cells in tumor tissues, and the method could replace flow cytometry or immunohistochemistry which is difficult to standardize (24–26). We subsequently identified the differences of 22 immune cells between the two groups of PTC patients. “Vioplot” package was employed to show the differential immune infiltrates between the two groups. The data analysis used the Wilcoxon rank-sum test. The threshold p-value <0.05 was the criterion for calculating the significance of samples.

Survival Analysis

Eighteen differential genes based on the two groups were identified. The RFS of these differential genes was analyzed according to the survival package of R software to find the more valuable prognostic biomarker. Receiver operating characteristic (ROC) and the area under the curve (AUC) were presented by utilizing the “survival ROC” package to confirm the prognostic value of clinical features. All data were tested by using an independent t-test, and the P-value was set to <0.05.

Results

Gene Mutation Profile of PTC

The top 30 genes with a mutation frequency per sample were displayed, and different color annotations indicated the different types of mutations (Figure 1). The most common type of mutation in the summarized figure was a missense mutation (Figure 2A). SNPs occurred more often than insertions or deletions (Figure 2B). SNV was mainly based on C>T base substitution in PTC (Figure 2C). Also, we calculated the median of PTC variants per sample (Figure 2D) as six and summarized the variant classifications based on the median and range (Figure 2E). The top 10 most common mutated genes of PTC were exhibited (Figure 2F), in which BRAF mutations were observed in 59% of cases, NRAS mutations in 8%, HRAS mutations in 3%, TG mutations in 3%, and TTN mutations in 2%. Finally, we showed the coexistence and exclusive correlation between mutant genes, where green represented coexistence and brown described a mutually exclusive relationship (Figure 2G). Moreover, we could see the BRAF mutation had the exclusive association with NRAS mutation or HRAS mutation.

Figure 1

Figure 2

The Differences in Different Clinical Characteristics Between the High and Low TMB Groups

We calculated the TMB of 487 PTC samples, which was the total number of mutations per million bases, ranging from 0.026 to 2.210 Mut/Mb, with a median of 0.211 Mut/Mb. According to the median value of TMB, we divided the PTC patients into the high and low TMB groups. Then we compared the differences in different clinical characteristics between the high and low TMB groups, including age, gender, AJCC-TNM stages, pathological subtype, BRAF status, and survival outcomes (Table 1). Because some samples lack clinical data, the number of samples exhibited in Table 1 was not exactly complete. We could see that the high TMB group had a higher proportion of PTC patients with older than 55 years old, larger tumors, higher tumor stages, and higher recurrence rate.

Lower TMB Predicted Better RFS Rate in PTC

Patients in the high TMB group tended to have a lower recurrence-free survival rate than that in the low TMB group (Figure 3A). In addition, TMB levels were higher in patients older than 55 (Figure 3B) or male sex (Figure 3C). PTC in the T3–4 stages had higher TMB levels than that in the T1–2 stages (Figure 3D). And TMB levels in stages III–IV were increased significantly compared with that in stages I–II (Figure 3G). However, levels of TMB grouped by AJCC-N (Figure 3E), M stage (Figure 3F), histopathological type (Figure 3H), and BRAF status (Figure 3I) had no significant differences. Most importantly, the ROC AUC score of TMB was 0.780, which was the most reliable indicator in the survival prediction of the model compared with age (AUC = 0.658), gender (AUC = 0.421), AJCC-T (AUC = 0.668), N (AUC = 0.566), M (AUC = 0.611), tumor stage (AUC = 0.708), or BRAF status (AUC = 0.520) (Figure 3J).

Figure 3

Identification of DEGs Based on High and Low TMB Groups

Eighteen DEGs (Table 2) were identified from high and low TMB groups in the TCGA database, which were shown by heatmap (Figure 4A). Fourteen down-regulated and four up-regulated genes in high TMB group were shown in the volcano plot (Figure 4B). Further, we explored the function of DEGs by GO analysis and found that DEG enrichment was mainly animal organ regeneration, extracellular matrix, rough endoplasmic reticulum, receptor antagonist activity, receptor inhibitor activity, and steroid hormone receptor activity (Figure 4C). The GSEA result was related to oxidative phosphorylation (OXPHOS), which demonstrated that the TMB in the high TMB group was positively correlated with OXPHOS level, while that in the low TMB group was the opposite (Figure 4D).

Figure 4

IL37, SLC5A5, NR4A3, and ODF3L1 Associated With PTC Relapse

We performed survival analysis on the 18 TMB related DEGs. Interleukin 37 (IL37), solute carrier family 5 member (SLC5A5), nuclear receptor subfamily 4 group A member 3 (NR4A3) and outer dense fiber of sperm tails 3 like 1 (ODF3L1) had been determined to be correlated with RFS of PTC patients. The results showed that higher expression of IL37 in tumor tissues predicted worse prognosis outcomes (Figure 5A). Conversely, the other three genes with a higher level in tumor tissues predicted a better prognosis (Figures 5B–D). The relationship between the four DEGs and each clinicopathological factor per sample in the high and low TMB groups was shown by the heatmap (Figure 5E). We could observe that the high expression of IL37 was concentrated in clinical features with poor prognosis, while the SLC5A5, NR4A3, and ODF3L1 were concentrated in clinical features with good prognosis. We further confirmed the expression levels of the four differential expression genes according to the RNA-sequencing dataset from TCGA. The expression level of IL37 in the tumor was significantly higher than that in the normal tissues (Figure 6A), while SLC5A5, NR4A3, and ODF3L1 were highly expressed in normal tissues than those in tumors (Figures 6B–D).

Figure 5

Figure 6

Differences in Immune Cell Infiltration Between High and Low TMB Groups

Through the deconvolution algorithm on the CIBERSORT website, the relative scale of 22 immune cells was measured. PTC patients in the high TMB group had lower infiltration of CD8⁺ T cells and M1 macrophages but higher gamma delta T cells than those in the low TMB group (Figure 7).

Figure 7

Discussion

The previous study has shown that the Tumor Node Metastasis (TNM) staging guide cannot effectively predict the prognosis of patients with low-risk PTC during initial treatment (27). Herein, we defined the characteristics of TMB in PTC and compared TMB as the best recurrence predictor than other clinical features. Generally speaking, TMB is the number of somatic non-synonymous mutations per base pair (Mb) in a specific gene group region. Herein we selected SNP data which was similar to whole-exome sequencing to calculate TMB. However, the TMB values of patients with different cancer types are highly heterogeneous, and the high TMB thresholds screened out still show huge differences in different cancer types (28). Besides, since many of the original studies on TMB are based on whole-exome modeling calculations, it is usually visualized as the number of mutations in the coding region of the genome (exome) (29). In fact, in addition to point mutations that could produce neoantigens, frameshift mutations produced by insertions and deletions could also cause tumors to produce neoantigens, thereby changing the efficacy of immune checkpoint inhibitors (30). The TMB we calculated based on genome sequencing was caused by somatic non-synonymous mutations, and the changes in TMB which might be caused by other factors were not taken into account, such as insertions or deletions. Therefore, in future, TMB detection still needs to consider many other interference factors.

GO analysis revealed that DEG enrichment was mainly involved in metabolic pathway. Also, patients in high TMB group showed poorer RFS for an increase in the level of OXPHOS. CD8⁺ T cells and M1 macrophages were decreased, which might be responsible for the recurrence of PTC with high TMB. Together these findings indicated that TMB in the tumor could help in effectively taking initial measures to monitor recurrence probability and avoid secondary surgery in the clinic. TMB has been reported to play important roles in many cancers (31). Higher TMB has more favorable RFS and is also a good predictor of response to ICPIs within a high microsatellite instability population in colorectal cancer (32, 33). Conversely, our study showed that PTC patients in the high TMB group had worse RFS (Figure 3A). Higher TMB was also associated with age older than 55, male, tumor size, and tumor staging (Figures 3B–D, G). The result was consistent with previous reports on clear cell renal cell carcinoma (21, 34). TMB had a rising trend in tall cell subtypes compared to classic and follicular types, but there might be no statistical difference due to insufficient sample size. BRAF^V600E mutation is significantly related to aggressive clinicopathological features and PTC relapse (35, 36). Furthermore, we also showed BRAF mutations as the most frequent and more prevalent in PTC than other mutations based on the TCGA data (Figure 1), which was consistent with previous reports (37). The interaction analysis of mutant genes showed that coexistence was more common than mutual exclusion (Figure 2G). However, compared with BRAF status only, TMB was a better predictor of relapse (Figure 3H). In general, ROC analysis showed the sensitivity and specificity of TMB were better than other clinical features. In our study, we found a significant risk factor for PTC recurrence and aggressiveness.

OXPHOS is up-regulated in many cancers, and its inhibitors with appropriate therapeutic index have new applications in targeted cancer therapy, including thyroid cancer (38–40). GSEA results showed higher levels of OXPHOS in the tumors of patients in the high TMB group (Figure 4D), which might be a potential relapse mechanism in PTC with a higher mutation burden. In the immune microenvironment of the tumor, CD8⁺ T cells can recognize neoantigens on the surface of tumor to play an anti-tumor effect (41). Besides, high CD8⁺ T cell density was significantly associated with a good disease-free survival rate for PTC patients and with a reduced lymph node metastasis incidence (42), and M1 macrophages in the tumor microenvironment inhibit tumor growth (43–46). The low abundance of CD8⁺ T cells and M1 macrophages was also found in PTC with high TMB (Figure 7), which might also clarify the underlying reasons of high TMB causing relapse.

Furthermore, four TMB related differential genes were highly associated with prognosis (Figure 5). IL37 was a risky prognostic gene and might play a vital role in the development of PTC. The result was consistent with oral leukoplakia and oral squamous cell carcinoma (47). However, IL37 plays a protective role in most cancers, such as lung cancer (48) and colorectal cancer (49). At present, the function of IL37 in the carcinogenesis of PTC has not been explored. SLC5A5, which codifies sodium iodide symporter, plays an essential role in thyroid metabolism, mediating the active transport of iodine from the bloodstream into the follicular cells (50). Due to the decreased expression of SLC5A5, some patients with differentiated thyroid cancer were not sensitive to radioiodine therapy (51). NR4A3 is an orphan receptor and regulates the cellular function and inflammation reaction (52). Compared with thyroid follicular adenoma (FTA), NR4A3 was significantly down-regulated in follicular thyroid carcinoma (FTC), which led to reduced apoptosis factors (53). However, the function of NR4A3 in PTC is still unclear. As for ODF3L1, there is no research report yet; however, it is an excellent prognostic molecular marker in PTC. Though we detected four genes associated with PTC relapse, these genes’ protein expression levels still need to be verified.

In conclusion, our results provided new information for the function of TMB in PTC. We showed that TMB was an independent prognostic factor to predict the recurrence of PTC and explained the possible mechanism. We also screened out genes related to relapse in PTC patients, which could become a new prognostic marker. At present, studies have shown that the targeted sequencing panel combined with bioinformatics algorithm could effectively reduce the detection cost of TMB and sample size requirements (54). Nevertheless, our study was retrospective based on public databases, and the sample size was not large enough. The detailed mechanism needs to be further elucidated in future studies. In general, our findings strongly suggest that TMB is clinically meaningful for monitoring the recurrence to shorten the follow-up time of PTC.

Funding

This research was supported by Guangzhou medicine and healthcare technology projects (20211A011010) and Guangzhou Science and Technology Plan Project (202102080170).

References

• 1

  Cancer Genome Atlas Research, N. Integrated Genomic Characterization of Papillary Thyroid Carcinoma. Cell (2014) 159(3):676–90. doi: 10.1016/j.cell.2014.09.050

  • CrossRef
  • Google Scholar

• 2

  LimHDevesaSSSosaJACheckDKitaharaCM. Trends in Thyroid Cancer Incidence and Mortality in the United States, 1974-2013. JAMA (2017) 317(13):1338–48. doi: 10.1001/jama.2017.2719

  • CrossRef
  • Google Scholar

• 3

  AhnDLeeGJSohnJH. Recurrence Following Hemithyroidectomy in Patients With Low- and Intermediate-Risk Papillary Thyroid Carcinoma. Br J Surg (2020) 107(6):687–94. doi: 10.1002/bjs.11430

  • CrossRef
  • Google Scholar

• 4

  DongWHoriuchiKTokumitsuHSakamotoANoguchiEUedaYet al. Time-Varying Pattern of Mortality and Recurrence From Papillary Thyroid Cancer: Lessons From a Long-Term Follow-Up. Thyroid (2019) 29(6):802–8. doi: 10.1089/thy.2018.0128

  • CrossRef
  • Google Scholar

• 5

  DingJWuWFangJZhaoJJiangL. Male Sex is Associated With Aggressive Behaviour and Poor Prognosis in Chinese Papillary Thyroid Carcinoma. Sci Rep (2020) 10(1). doi: 10.1038/s41598-020-60199-9

  • CrossRef
  • Google Scholar

• 6

  WangXChengWLiuCLiJ. Tall Cell Variant of Papillary Thyroid Carcinoma: Current Evidence on Clinicopathologic Features and Molecular Biology. Oncotarget (2016) 7(26):40792–9. doi: 10.18632/oncotarget.8215

  • CrossRef
  • Google Scholar

• 7

  Llamas-OlierAECuéllarDIBuitragoG. Intermediate-Risk Papillary Thyroid Cancer: Risk Factors for Early Recurrence in Patients With Excellent Response to Initial Therapy. Thyroid (2018) 28(10):1311–7. doi: 10.1089/thy.2017.0578

  • CrossRef
  • Google Scholar

• 8

  KimHSChoiJHLeeJYKangJMyungJKKimWHet al. Downregulation of SMOC2 Expression in Papillary Thyroid Carcinoma and Its Prognostic Significance. Sci Rep (2020) 10(1):4853. doi: 10.1038/s41598-020-61828-z

  • CrossRef
  • Google Scholar

• 9

  SanjariMKordestaniZSafaviMMashroutehMFekriSoofiAbadiMGhaseminejadA. Enhanced Expression of Cyclin D1 and C-myc, a Prognostic Factor and Possible Mechanism for Recurrence of Papillary Thyroid Carcinoma. Sci Rep (2020) 10(1):5100. doi: 10.1038/s41598-020-61985-1

  • CrossRef
  • Google Scholar

• 10

  JianyongLZhihuiLRixiangGJingqiangZ. Using a Nomogram Based on Preoperative Serum Fibrinogen Levels to Predict Recurrence of Papillary Thyroid Carcinoma. BMC Cancer (2018) 18(1):390. doi: 10.1186/s12885-018-4296-7

  • CrossRef
  • Google Scholar

• 11

  EnumahSFingeretAParangiSDias–SantagataDSadowPMLubitzCCet al. Braf(V600e) Mutation is Associated With an Increased Risk of Papillary Thyroid Cancer Recurrence. World J Surg (2020) 44(8):2685–91. doi: 10.1007/s00268-020-05521-2

  • CrossRef
  • Google Scholar

• 12

  KempfEde BeyerJACookJHolmesJMohammedSNguyenTLet al. Overinterpretation and Misreporting of Prognostic Factor Studies in Oncology: A Systematic Review. Br J Cancer (2018) 119(10):1288–96. doi: 10.1038/s41416-018-0305-5

  • CrossRef
  • Google Scholar

• 13

  CrowJSamuelGGodwinAK. Beyond Tumor Mutational Burden: Potential and Limitations in Using Exosomes to Predict Response to Immunotherapy. Expert Rev Mol Diagn (2019) 19(12):1079–88. doi: 10.1080/14737159.2020.1688144

  • CrossRef
  • Google Scholar

• 14

  YarchoanMHopkinsAJaffeeEM. Tumor Mutational Burden and Response Rate to PD-1 Inhibition. N Engl J Med (2017) 377(25):2500–1. doi: 10.1056/NEJMc1713444

  • CrossRef
  • Google Scholar

• 15

  PleasanceEDStephensPJO’MearaSMcBrideDJMeynertDJonesAet al. A Small-Cell Lung Cancer Genome With Complex Signatures of Tobacco Exposure. Nature (2010) 463(7278):184–90. doi: 10.1038/nature08629

  • CrossRef
  • Google Scholar

• 16

  ZhangWEdwardsAFlemingtonEKZhangK. Racial Disparities in Patient Survival and Tumor Mutation Burden, and the Association Between Tumor Mutation Burden and Cancer Incidence Rate. Sci Rep (2017) 7(1):13639. doi: 10.1038/s41598-017-13091-y

  • CrossRef
  • Google Scholar

• 17

  HsuYCChangYHChangGCHoBCYuanSSLiYCet al. Tumor Mutation Burden and Recurrent Tumors in Hereditary Lung Cancer. Cancer Med (2019) 8(5):2179–87. doi: 10.1002/cam4.2120

  • CrossRef
  • Google Scholar

• 18

  JaratlerdsiriWChanEKFGongTPetersenDCKalsbeekAMFVenterPAet al. Whole-Genome Sequencing Reveals Elevated Tumor Mutational Burden and Initiating Driver Mutations in African Men With Treatment-Naive, High-Risk Prostate Cancer. Cancer Res (2018) 78(24):6736–46. doi: 10.1158/0008-5472.CAN-18-0254

  • CrossRef
  • Google Scholar

• 19

  NandakumarVMillsJR. The Now and Beyond of Tumor Mutational Burden as a Predictor of Response to Immune Checkpoint Inhibitors. Clin Chem (2019) 65(2):357. doi: 10.1373/clinchem.2018.295097

  • CrossRef
  • Google Scholar

• 20

  MayakondaALinDCAssenovYPlassCKoefflerHP. Maftools: Efficient and Comprehensive Analysis of Somatic Variants in Cancer. Genome Res (2018) 28(11):1747–56. doi: 10.1101/gr.239244.118

  • CrossRef
  • Google Scholar

• 21

  ZhangCLiZQiFHuXLuoJ. Exploration of the Relationships Between Tumor Mutation Burden With Immune Infiltrates in Clear Cell Renal Cell Carcinoma. Ann Transl Med (2019) 7(22):648. doi: 10.21037/atm.2019.10.84

  • CrossRef
  • Google Scholar

• 22

  RitchieMEPhipsonBWuDHuYLawCWShiWet al. Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res (2015) 43(7):e47. doi: 10.1093/nar/gkv007

  • CrossRef
  • Google Scholar

• 23

  YuGWangLGHanYHeQYet al. clusterProfiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118

  • CrossRef
  • Google Scholar

• 24

  ChenBKhodadoustMSLiuCLKhodadoustAMAlizadehAAet al. Profiling Tumor Infiltrating Immune Cells With CIBERSORT. Methods Mol Biol (2018) 1711. doi: 10.1007/978-1-4939-7493-1_12

  • CrossRef
  • Google Scholar

• 25

  DesmedtCSalgadoRForniliMPruneriGVan den EyndenGZoppoliGet al. Immune Infiltration in Invasive Lobular Breast Cancer. J Natl Cancer Inst (2018) 110(7):768–76. doi: 10.1093/jnci/djx268

  • CrossRef
  • Google Scholar

• 26

  GentlesAJNewmanAMLiuCLBratmanSVFengWKimDet al. The Prognostic Landscape of Genes and Infiltrating Immune Cells Across Human Cancers. Nat Med (2015) 21(8):938–45. doi: 10.1038/nm.3909

  • CrossRef
  • Google Scholar

• 27

  TuttleRMTalaHShahJLeboeufRGhosseinRGonenMet al. Estimating Risk of Recurrence in Differentiated Thyroid Cancer After Total Thyroidectomy and Radioactive Iodine Remnant Ablation: Using Response to Therapy Variables to Modify the Initial Risk Estimates Predicted by the New American Thyroid Association Staging System. Thyroid (2010) 20(12):1341–9. doi: 10.1089/thy.2010.0178

  • CrossRef
  • Google Scholar

• 28

  SamsteinRMLeeCHShoushtariANHellmannMDShenRJanjigianYYet al. Tumor Mutational Load Predicts Survival After Immunotherapy Across Multiple Cancer Types. Nat Genet (2019) 51(2):202–6. doi: 10.1038/s41588-018-0312-8

  • CrossRef
  • Google Scholar

• 29

  ChalmersZRConnellyCFFabrizioDGayLAliSMEnnisRet al. Analysis of 100,000 Human Cancer Genomes Reveals the Landscape of Tumor Mutational Burden. Genome Med (2017) 9(1):34. doi: 10.1186/s13073-017-0424-2

  • CrossRef
  • Google Scholar

• 30

  TurajlicSLitchfieldKXuHRosenthalRMcGranahanNet al. Insertion-and-Deletion-Derived Tumour-Specific Neoantigens and the Immunogenic Phenotype: A Pan-Cancer Analysis. Lancet Oncol (2017) 18(8):1009–21. doi: 10.1016/s1470-2045(17)30516-8

  • CrossRef
  • Google Scholar

• 31

  RomeroD. TMB is Linked With Prognosis. Nat Rev Clin Oncol (2019) 16(6):336. doi: 10.1038/s41571-019-0206-4

  • CrossRef
  • Google Scholar

• 32

  LeeDWHanSWBaeJMJangHHanHKimHet al. Tumor Mutation Burden and Prognosis in Patients With Colorectal Cancer Treated With Adjuvant Fluoropyrimidine and Oxaliplatin. Clin Cancer Res (2019) 25(20):6141–7. doi: 10.1158/1078-0432.CCR-19-1105

  • CrossRef
  • Google Scholar

• 33

  SchrockABOuyangCSandhuJSokolEJinDRossJSet al. Tumor Mutational Burden is Predictive of Response to Immune Checkpoint Inhibitors in MSI-high Metastatic Colorectal Cancer. Ann Oncol (2019) 30(7):1096–103. doi: 10.1093/annonc/mdz134

  • CrossRef
  • Google Scholar

• 34

  ZhangCHeHHuXLiuAoHuangDaXuYet al. Development and Validation of a Metastasis-Associated Prognostic Signature Based on Single-Cell RNA-Seq in Clear Cell Renal Cell Carcinoma. Aging (2019) 11(22):10183–202. doi: 10.18632/aging.102434

  • CrossRef
  • Google Scholar

• 35

  GengJWangHLiuYTaiJLiuYZhangJet al. Correlation Between BRAF (V600E) Mutation and Clinicopathological Features in Pediatric Papillary Thyroid Carcinoma. Sci China Life Sci (2017) 60(7):729–38. doi: 10.1007/s11427-017-9083-8

  • CrossRef
  • Google Scholar

• 36

  XingMAlzahraniASCarsonKAShongYKKimTYViolaDet al. Association Between Braf V600E Mutation and Recurrence of Papillary Thyroid Cancer. J Clin Oncol (2015) 33(1):42–50. doi: 10.1200/jco.2014.56.8253

  • CrossRef
  • Google Scholar

• 37

  LandaIIbrahimpasicTBoucaiLSinhaRKnaufJAShahRHet al. Genomic and Transcriptomic Hallmarks of Poorly Differentiated and Anaplastic Thyroid Cancers. J Clin Invest (2016) 126(3):1052–66. doi: 10.1172/JCI85271

  • CrossRef
  • Google Scholar

• 38

  AshtonTMMcKennaWGKunz-SchughartLAHigginsGSet al. Oxidative Phosphorylation as An Emerging Target in Cancer Therapy. Clin Cancer Res (2018) 24(11):2482–90. doi: 10.1158/1078-0432.CCR-17-3070

  • CrossRef
  • Google Scholar

• 39

  SicaVBravo-San PedroJMStollGKroemerGet al. Oxidative Phosphorylation as a Potential Therapeutic Target for Cancer Therapy. Int J Cancer (2019) 146(1):10–7. doi: 10.1002/ijc.32616

  • CrossRef
  • Google Scholar

• 40

  ThakurSDaleyBGaskinsKVaskoVVBoufraqechMPatelDet al. Metformin Targets Mitochondrial Glycerophosphate Dehydrogenase to Control Rate of Oxidative Phosphorylation and Growth of Thyroid Cancer In Vitro and In Vivo. Clin Cancer Res (2018) 24(16):4030–43. doi: 10.1158/1078-0432.ccr-17-3167

  • CrossRef
  • Google Scholar

• 41

  ReadingJLGalvez-CancinoFSwantonCLladserAPeggsKSQuezadaSAet al. The Function and Dysfunction of Memory CD8(+) T Cells in Tumor Immunity. Immunol Rev (2018) 283(1):194–212. doi: 10.1111/imr.12657

  • CrossRef
  • Google Scholar

• 42

  AghajaniMJYangTMcCaffertyCEGrahamSWuXNilesN. Predictive Relevance of Programmed Cell Death Protein 1 and Tumor-Infiltrating Lymphocyte Expression in Papillary Thyroid Cancer. Surgery (2018) 163(1):130–6. doi: 10.1016/j.surg.2017.04.033

  • CrossRef
  • Google Scholar

• 43

  PaulSChhatarSMishraALalG. Natural Killer T Cell Activation Increases Inos(+)CD206(-) M1 Macrophage and Controls the Growth of Solid Tumor. J Immunother Cancer (2019) 7(1):208. doi: 10.1186/s40425-019-0697-7

  • CrossRef
  • Google Scholar

• 44

  TraversMBrownSMDunworthMHolbertCEWiehagenKRBachmanKE. DFMO and 5-Azacytidine Increase M1 Macrophages in the Tumor Microenvironment of Murine Ovarian Cancer. Cancer Res (2019) 79(13):3445–54. doi: 10.1158/0008-5472.CAN-18-4018

  • CrossRef
  • Google Scholar

• 45

  CaoMYanHHanXWengLWeiQSunXet al. Ginseng-Derived Nanoparticles Alter Macrophage Polarization to Inhibit Melanoma Growth. J Immunother Cancer (2019) 7(1):326. doi: 10.1186/s40425-019-0817-4

  • CrossRef
  • Google Scholar

• 46

  RenFFanMMeiJWuYLiuCPuQet al. Interferon-Gamma and Celecoxib Inhibit Lung-Tumor Growth Through Modulating M2/M1 Macrophage Ratio in the Tumor Microenvironment. Drug Des Devel Ther (2014) 8. doi: 10.2147/DDDT.S66302

  • CrossRef
  • Google Scholar

• 47

  LinLWangJLiuDLiuSXuHJiNet al. Interleukin-37 Expression and its Potential Role in Oral Leukoplakia and Oral Squamous Cell Carcinoma. Sci Rep (2016) 6. doi: 10.1038/srep26757

  • CrossRef
  • Google Scholar

• 48

  GeGWangAYangJChenYYangJLiYet al. Interleukin-37 Suppresses Tumor Growth Through Inhibition of Angiogenesis in Non-Small Cell Lung Cancer. J Exp Clin Cancer Res (2016) 35(1). doi: 10.1186/s13046-016-0293-3

  • CrossRef
  • Google Scholar

• 49

  ZhuBLuoJJiangYYuLLiuMFuJ. Prognostic Significance of Nomograms Integrating IL-37 Expression, Neutrophil Level, and MMR Status in Patients With Colorectal Cancer. Cancer Med (2018) 7(8):3682–94. doi: 10.1002/cam4.1663

  • CrossRef
  • Google Scholar

• 50

  SpitzwegCMorrisJC. Genetics and Phenomics of Hypothyroidism and Goiter Due to NIS Mutations. Mol Cell Endocrinol (2010) 322(1-2):56–63. doi: 10.1016/j.mce.2010.02.007

  • CrossRef
  • Google Scholar

• 51

  DohánOBalochZBánréviZLivolsiVCarrascoN. Rapid Communication: Predominant Intracellular Overexpression of the Na(+)/I(-) Symporter (NIS) in a Large Sampling of Thyroid Cancer Cases. J Clin Endocrinol Metab (2001) 86(6):2697–700. doi: 10.1210/jcem.86.6.7746

  • CrossRef
  • Google Scholar

• 52

  PrinceLRProssedaSDHigginsKet al. NR4A Orphan Nuclear Receptor Family Members, NR4A2 and NR4A3, Regulate Neutrophil Number and Survival. Blood (2017) 130(8):1014–25. doi: 10.1182/blood-2017-03-770164

  • CrossRef
  • Google Scholar

• 53

  BorupRRossingMHenaoRYamamotoYKrogdahlAGodballeCet al. Molecular Signatures of Thyroid Follicular Neoplasia. Endocr Relat Cancer (2010) 17(3):691–708. doi: 10.1677/ERC-09-0288

  • CrossRef
  • Google Scholar

• 54

  ZehirABenayedRShahRHSyedAMiddhaSKimHRet al. Mutational Landscape of Metastatic Cancer Revealed From Prospective Clinical Sequencing of 10,000 Patients. Nat Med (2017) 23(6):703–13. doi: 10.1038/nm.4333

  • CrossRef
  • Google Scholar