Comprehensive gene expression analysis for exploring the association between glucose metabolism and differentiation of thyroid cancer

doi:10.21203/rs.2.10560/v5

Download PDF

Research article

Comprehensive gene expression analysis for exploring the association between glucose metabolism and differentiation of thyroid cancer

https://doi.org/10.21203/rs.2.10560/v5

This work is licensed under a CC BY 4.0 License

Journal Publication

published 30 Dec, 2019

Read the published version in BMC Cancer →

You are reading this older preprint version

Read the latest preprint version →

Background: The principle of loss of iodine uptake and increased glucose metabolism according to dedifferentiation of thyroid cancer is clinically assessed by imaging. Though these biological properties are widely applied to appropriate iodine therapy, the understanding of the genomic background of this principle is still lacking. We investigated the association between glucose metabolism and differentiation in advanced thyroid cancer as well as papillary thyroid cancer (PTC). Methods: We used RNA sequencing of 505 patients with PTC obtained from the Cancer Genome Archives and microarray data of poorly-differentiated and anaplastic thyroid cancer (PDTC/ATC). The signatures of GLUT and glycolysis were estimated to assess glucose metabolic profiles. The glucose metabolic profiles were associated with tumor differentiation score (TDS) and BRAFV600E mutation status. In addition, survival analysis of glucose metabolic profiles was performed for predicting recurrence-free survival. Results: In PTC, the glycolysis signature was positively correlated with TDS, while the GLUT signature was inversely correlated with TDS. These correlations were significantly stronger in the BRAFV600E negative group than the positive group. Meanwhile, both GLUT and glycolysis signatures were negatively correlated with TDS in advanced thyroid cancer. The high glycolysis signature was significantly associated with poor prognosis in PTC in spite of high TDS. The glucose metabolic profiles are intricately associated with tumor differentiation in PTC and PDTC/ATC. Conclusions: As glycolysis was an independent prognostic marker, we suggest that the glucose metabolism features of thyroid cancer could be another biological progression marker different from differentiation and provide clinical implications for risk stratification.

Cancer Biology

glucose metabolism

thyroid cancer

tumor differentiation

glycolysis

GLUT

Molecular imaging has been widely used for the biological characterization of thyroid cancer in the clinic. Iodine and glucose metabolism have been respectively identified by radioactive iodine imaging and 2-deoxy-2-[¹⁸F]fluoro-D-glucose (FDG) positron emission tomography (PET), which have played important role in therapeutic plan (1, 2). The fact that thyroid cancers with poor differentiation have lower I-131 avidity and higher F-18 FDG avidity have long been used to presume the aggressiveness of thyroid cancer (3). This ‘flip-flop phenomenon’ was thought to have resulted from the loss of I-131 concentration capacity with increasing demand for glucose during the dedifferentiation of tumor cells (4).

As metabolic features of thyroid cancer could be associated with genomic alteration as well as tumor differentiation, previous studies have revealed the trend of relationship by using imaging studies and pathologic profiles (3, 5, 6). Among the major cancer drivers of differentiated thyroid cancers, B-type Raf kinase (BRAF) mutation exclusively occurred in papillary thyroid cancer (PTC) and PTC-derived anaplastic thyroid cancer (ATC) is associated with low iodine avidity (3). Thyroid cancer with BRAF mutation acquires more aggressive features through BRAF mutation-associated silencing of thyroid-specific genes, differentiation markers of iodine metabolism (6). The dedifferentiation according to BRAF mutation leads to increased glucose consumption of cancer cells (1), which is supported by the finding that PTC with BRAF mutation showed a trend of more FDG avidity (7). In spite of these cross-sectional studies using imaging, there is a lack of comprehensive understanding of the association between tumor glucose metabolism, differentiation and BRAF mutation. Meanwhile, recent comprehensive analyses using next generation sequencing have revealed new integrative molecular subtypes as well as cancer drivers (8). In this regard, the integrative analysis of iodine and glucose metabolism based on the systemic gene expression data is needed for the molecular background of imaging and therapy of thyroid cancer.

The aim of this study is to identify the association between the glucose metabolism profile and the differentiation of thyroid cancer using transcriptomic signatures. As tumor glucose metabolism clinically assessed by FDG PET is associated with the expression of glucose transporters (GLUT) and glycolysis pathways (9-11), gene signatures of GLUT and glycolysis were used for defining glucose metabolic profiles of the tumors. We investigated whether the glucose metabolic profiles were correlated with tumor differentiation and BRAF mutations in advanced thyroid cancer as well as PTC. We also investigated whether these metabolic profiles were associated with clinical outcome, and eventually aimed to redefine the biological progression of thyroid cancer in terms of glucose metabolism as well as differentiation related with iodine uptake.

Subjects and data acquisition

RNA sequencing data of well differentiated papillary thyroid cancers (THCA) were obtained from The Cancer Genome Atlas (TCGA) which was publicly available. The gene expression data of 21,022 genes and clinical information from 505 patients with PTC were downloaded from the TCGA portal (http://portal.gd.cancer.gov/) using ‘TCGABiolinks’ R/Bioconductor package. Microarray data of poorly differentiated thyroid cancer (PDTC) and ATC were also obtained in NCBI’s Gene Expression Omnibus (GEO) Series accession number GSE76039. Clinical information of PDTC/ATC patients was downloaded from cBioPortal (http://www.cbioportal. GSE76039org/public-portal). Data from patients with PTC and PDTC/ATC were all collected from primary tumor, except 7 samples from recurred or metastatic tumors of PDTC/ATC (6 samples from recurred tumor, one sample from metastatic tumor).

Glucose metabolic profiles

We normalized mRNA transcripts counts to adjust GC-content effects and between-lane distributional difference. Using normalized gene counts, we calculated the enrichment scores of glucose metabolic signatures for GLUT and glycolysis. Glycolysis signatures were defined by Reactome pathway (12). To estimate the enrichment score of glycolysis, a single sample gene set enrichment analysis (ssGSEA) was performed using the curated gene sets of glycolysis signatures(13). This analysis was performed with GSVA R/Bioconductor package (14). For estimating GLUT score, the sum of log-expression values of GLUT1 and GLUT3 was used as FDG uptake is dependent on these two GLUT subtypes (15). Tumor differentiation score (TDS) were also calculated by expression values of a set of specific genes (8). Notably, the set of genes to calculate TDS did not include genes related to glucose transporters and glycolysis.

Statistical analysis

Student’s t-test was used to analyze whether the signatures of glucose metabolism were significantly different according to the genetic mutation status of BRAFV600E. Pearson’s correlation test was performed to evaluate the associations between TDS and signatures of glucose metabolism. The difference between two correlation coefficients regarding BRAFV600E was calculated by Fisher z-transformation. In addition, the difference between two correlation coefficients according to the histologic cell types was calculated by Fisher z-transformation. The effect of the signatures of glucose metabolism on patients’ recurrence-free survival was analyzed using the Cox proportional hazard model. Patients were divided into two groups; lower or higher than the median value of GLUT, and glycolysis unit. Kaplan-Meier survival curves were demonstrated to compare recurrence-free survival between two groups. The statistical significance was evaluated using the Log rank test. The differences between groups were considered statistically significant at p-value < 0.05. All statistical analyses were performed with R 3.4.1 (https://www.r-project.org/).

To evaluate the association between tumor differentiation and glucose metabolism in the thyroid cancer, we used two cohorts, PTC samples of TCGA and PDTC/ATC samples obtained by GEO (GSE76039). The clinical and pathological characteristics of PTC patients and PDTC/ATC patients are summarized in Table 1 and Table 2. Notably, there were unavailable clinical data such as race.

Firstly, GLUT and glycolysis enrichment scores were estimated by gene expression data of the two cohorts and associated with the mutation status, BRAFV600E. The enrichment scores of GLUT and glycolysis were significantly different in PTC according to BRAFV600E mutation status (Figure 1). BRAFV600E positive PTC have higher GLUT signature and lower glycolysis signature (BRAFV600E positive group 20.53±1.08 vs. BRAFV600E negative group 19.13±1.46, t = -12.09, p < 0.0001 for GLUT; BRAFV600E positive group -0.28±0.77 vs. BRAFV600E negative group 0.24±1.11, t = 6.06, p < 0.0001 for glycolysis) than BRAFV600E negative PTC (Figure 1).

To evaluate the tumor differentiation, we employed tumor differentiation score, TDS, calculated by sixteen genes related to thyroid functions. The TDS was compared with glucose metabolic profiles. In PTC group, the TDS was negatively correlated with GLUT signature, but positively correlated with glycolysis signature (r = -0.59, p < 0.0001 for GLUT; Figure 2A, r = 0.33, p < 0.0001 for glycolysis; Figure 2B). We divided PTC into two subgroups according to the BRAFV600E mutation status. PTC without BRAFV600E showed a negative correlation between TDS and GLUT (r = -0.57, p < 0.001), while it showed a positive correlation between TDS and glycolysis (r = 0.35, p < 0.001). PTC with BRAFV600E showed a trend of negative correlation between GLUT and TDS (r = -0.18, p = 0.065) and no significant correlation between glycolysis and TDS. These correlations were significantly stronger in BRAFV600E negative group than positive group (z = 5.22 for GLUT, z = -4.48 for glycolysis, and p < 0.001 for both GLUT and glycolysis; Figure 2C-F). In PDTC/ATC group, the relationship between the TDS and GLUT signature showed a significant negative correlation (r = -0.68, p < 0.0001; Figure 2G). The TDS also showed a significant negative correlation to the glycolysis signature (r = -0.38, p = 0.019; Figure 2H). We also evaluated whether different cell types of PTC were associated with glucose metabolic signatures. Classical cell type PTC have higher GLUT signature and lower glycolysis signature than follicular cell type PTC (Classical cell type 20.13±1.32 vs. Follicular cell type 18.42±1.20, t = 11.55, p < 0.001 for GLUT; Classical cell type -0.13±0.94 vs. Follicular cell type 0.52±1.12, t = -5.65, p < 0.001 for glycolysis; Additional file 1: Figure S1A, B). The negative correlation between TDS and GLUT and the positive correlation between TDS and glycolysis were found in both cell types, classical and follicular types. The strength of their correlations showed no significant difference when we compare the correlation coefficients by Fisher z-transformation (Classical cell type: r = -0.47, p < 0.001 for GLUT, r = 0.23, p < 0.001 for glycolysis; Follicular cell type: r = -0.43, p < 0.001 for GLUT, r = 0.34, p = 0.001 for glycolysis; Additional file 1: Figure S1C, D). Individual genes that constitute TDS and glucose metabolic pathway were represented by heatmaps (Figure 3).

We assessed the association between signatures of glucose metabolism and patients’ prognosis. Kaplan-Meier survival curves of both PTC and PDTC/ATC patients with low glycolysis signature showed significantly better recurrence-free survival than the other group (p = 0.045 and 0.015, respectively, estimated by log-rank tests; Figure 4A, B). The glycolysis signature of the primary tumor showed a significant correlation with the N-stage (N negative group 0.13±1.03 vs. N positive group -0.24±0.89, t = 3.86, p = 0.0001), but not with the M-stage (M negative group -0.08±1.03 vs. M positive group -0.16±0.89, t = 0.22, p = 0.823; Additional file 2: Figure S2). We then assessed the predictive value of glycolysis signature for recurrence using a cox proportional hazard model. The hazard ratios (H.R.) of the influence of the glycolysis signature on the recurrence of PTC are estimated by uni- and multi-variate analyses (Table 3). In PTC, a univariate cox hazard proportional model revealed that the high glycolysis signature alone had a significant influence on the recurrence-free survival (H.R. = 1.50; C.I. = 1.03–2.17; p = 0.033). The adjustment for age and gender maintained its influence on the recurrence-free survival (H.R. = 1.50; C.I. = 1.04–2.16; p = 0.031). Additional adjustment for T-stage and N-stage still demonstrated worse prognosis on PTC patients with high glycolysis (H.R. = 1.92; C.I. = 1.22–3.00; p = 0.005). Although TDS itself did not influence the recurrence-free survival (H.R. = 0.76; C.I. = 0.55 – 1.06; p = 0.102), the adjustment for TDS as well as age, gender, T-stager, and N-stage maintained the influence of glycolysis signature on the recurrence-free survival (H.R. = 1.98; C.I. = 1.26–3.10; p = 0.003). The clinicopathological characteristics between PTC patients with low glycolysis and PTC patients with high glycolysis showed no significant difference except for age (Table 4). GLUT score showed no significant correlation with PTC patients’ prognosis (p = 0.85; Figure 4C). On the contrary, we found that PDTC/ATC patients with low GLUT signature have significantly longer recurrence-free survival than the other group (p = 0.0063; Figure 4D).

In this study, we used transcriptome data of thyroid cancers of two cohorts, PTC and PDTC/ATC to evaluate glucose metabolic profiles and tumor differentiation. In PTC, a trend of higher GLUT and lower glycolysis was found in tumors with BRAFV600E mutation and those with relatively poor differentiation. There is a controversy regarding the correlation of BRAFV600E, GLUT1, and tumor differentiation genes according to several previous reports (16-18). Our results are consistent with other results of higher GLUT1 in less differentiated thyroid cancer (16) and negative correlation between BRAFV600E mutation and tumor differentiation genes in PTC (17, 18). The results of TDS positively correlated with glycolysis and negatively correlated with GLUT were particularly found in PTC without BRAFV600E mutation. Moreover, the high glycolysis enrichment was significantly associated with poor clinical outcome, even it was associated with well differentiation in PTC. However, this paradoxical opposite direction of correlation was not found in PDTC/ATC cohorts, which showed both GLUT and glycolysis were negatively correlated with TDS and associated with poor clinical outcome.

We have found that thyroid cancers with poor differentiation showed higher GLUT expression (Figure 2). Cancer cells demand a higher amount of energy according to the progression, which is associated with enhanced aerobic glycolysis in the advanced cancer (19). The glucose demand of cancers cause overexpression of GLUT1 and/or GLUT3 to increase glucose influx (20). Our results were compatible with the previous studies since poorly differentiated thyroid cancers would need higher glucose uptake through GLUT expression (21). On the other hand, glycolysis signatures of thyroid cancers with poor differentiation were inconsistent between PTC and PDTC/ATC (Figure 2). In PTC, relatively well-differentiated tumors, glycolysis was positively associated with TDS. In PDTC/ATC, a negative correlation was shown between TDS and glycolysis. Considering PTC of TCGA data are relatively well-differentiated tumors compared with PDTC/ATC cohort, the association between the differentiation and glycolysis might have ‘U shape’ pattern; high glycolysis signatures were shown in ATC and some types of well-differentiated PTC. Moreover, increased glycolysis was associated with poor clinical outcome in spite of high TDS. It implies that the differentiation of thyroid cancer may not be the only factor that reflects the biological progression of thyroid cancer. Instead, in addition to differentiation, glucose metabolic profiles represented by glycolysis should be further considered to infer the progression of thyroid cancer, particularly for the specific subtypes of PTC, BRAFV600E negative and/or follicular variants.

Our results demonstrated that the signatures of GLUT and glycolysis can act as prognostic factors in predicting recurrence-free survival of thyroid cancer patients (Figure 4). The multivariate analysis revealed that the two variables, N-stage and the glycolysis signature, were significantly associated with the recurrence. As N-stage has been regarded as an important conventional biomarker related to prognosis which can be confirmed by surgical exploration, our results emphasized again the importance of the lymph node status in thyroid cancer as postoperative risk stratification. In terms of another prognostic marker in our results, glycolysis signature, it is notable that glucose metabolism profiles can be noninvasively estimated by FDG PET. GLUT score, calculated by GLUT1 and GLUT3, was reported as the major deterministic factor for the FDG uptake in various studies (10, 15, 21, 22), while a recent study showed a moderate correlation between FDG uptake and GLUT regarding a complex mechanism of glucose metabolism (23). Furthermore, glycolysis activity is also associated with FDG uptake in vivo (22, 24). According to the kinetic model of FDG, glycolysis activity is associated with FDG retention, which can be visualized by dual-time FDG PET (25, 26). In general, poorly differentiated thyroid cancers are known to have a worse outcome as compared to well differentiated thyroid cancers (27). However, a subset of well-differentiated carcinoma shows relatively poorly outcome in tumors with increased glycolysis with well-differentiated type. As GLUT and glycolytic activity were differently associated with TDS, noninvasive characterization using FDG PET and radioactive iodine imaging could play a role in risk stratification when considered with other prognostic factors as well as biological characterization of the tumor. As a future work, more specifically, dual-time FDG PET could be used for estimating glycolysis activity, which might identify a subtype of the tumor with enhanced glycolysis and well-differentiation.

Although, we analyzed glucose profiles of thyroid cancer in the two different cohorts with different differentiation, there are some limitations. Firstly, the PDTC/ATC sample size was small. Though we found the expected role of GLUT and glycolysis on the prognosis of PDTC/ATC patients, further studies with a larger group are needed. Another limitation was that the transcriptome data from two cohorts were not be combined and analyzed since those data were obtained from different resources, RNA sequencing and microarray. For PDTC/ATC patients, a few samples were acquired from recurred or metastatic tumors (recurred tumors for 6 samples and one sample from metastatic tumor), which might affect our analysis though recurred and metastatic tumors have similar glucose metabolic profiles with their primary tumors. Furthermore, noninvasive imaging biomarkers using iodine scan and FDG PET integrated with transcriptomic data could clarify our results of the association of differentiation and glucose metabolic profiles. We could expect the clinical application of our results as a future study, such as the estimation of glycolysis score based on gene expression profiles in fine-needle aspiration samples. In addition, according to our results, the expression of genes related to glycolysis may be examined by tissue samples in the clinic to stratify patients’ outcome, even though further studies focusing on clinical outcomes and clinical decision according to the glycolysis are needed.

According to the integrative analysis of iodine and glucose metabolism based on the systemic gene expression data, metabolic profiles were not simply associated with tumor differentiation. Cancer cellular GLUT expression was negatively associated with tumor differentiation in both PTC and PDTC/ATC. The enrichment of glycolysis was positively associated with the differentiation in well-differentiated PTC, while it was negatively associated with the PDTC/ATC. Overall, there might be ‘U-shape’ pattern for the association of the differentiation and glycolysis. Furthermore, increased glycolysis was poor prognosis in spite of the well-differentiated tumor. We anticipate that the biological and prognostic characteristics of glucose metabolic profiles could provide insight for biomarker using FDG PET and appropriate therapeutic plan.

FDG: 2-deoxy-2-[18F]fluoro-D-glucose

PET: positron emission tomography

BRAF: B-type Raf kinase

PTC: papillary thyroid cancer

ATC: anaplastic thyroid cancer

GLUT: glucose transporter

THCA: well differentiated papillary thyroid cancer

TCGA: The Cancer Genome Atlas

PDTC: poorly differentiated thyroid cancer

GEO: Gene Expression Omnibus

ssGSEA: single sample gene set enrichment analysis

TSD: tumor differentiation score

H.R.: hazard ratio

HK: hexokinase

C.I. = Confidence Interval

Ethics approval and consent to participate

Patient data we used were acquired by a publicly available dataset that removed patient identifiers. The original datasets were approval-free to use based on public databases.

Consent for publication

Not applicable.

Availability of data and materials

The datasets used and/or analyzed during the current study are available in the TCGA Research Networks: http://cancergenome.nih.gov/ and GEO: http://www.ncbi.nlm.nih.gov/geo.

Competing interests

The authors declare that they have no competing interests.

Funding

Funding was provided by Seoul National University Hospital Research Fund (grant no. 2620180060). The funder had no role in designing the study, in collection, analysis, and interpretation of data, or in writing the manuscript.

Authors’ contributions

HC designed the study. HYS and HC acquired and analyzed data, and drafted the manuscript. JCP, GJC, JC, and KWK participated in discussion and provided critical revision. All authors read and confirmed the manuscript.

Acknowledgements

The results here are in based upon data generated by the TCGA Research Network: “http://cancergenome.nih.gov/”.

Feine U, Lietzenmayer R, Hanke JP, Held J, Wöhrle H, Müller-Schauenburg W. Fluorine-18-FDG and iodine-131-iodide uptake in thyroid cancer. Journal of Nuclear Medicine. 1996;37(9):1468-72.
Schlumberger M, Brose M, Elisei R, Leboulleux S, Luster M, Pitoia F, et al. Definition and management of radioactive iodine-refractory differentiated thyroid cancer. The lancet Diabetes & endocrinology. 2014;2(5):356-8.
Rivera M, Ghossein RA, Schoder H, Gomez D, Larson SM, Tuttle RM. Histopathologic characterization of radioactive iodine‐refractory fluorodeoxyglucose‐positron emission tomography‐positive thyroid carcinoma. Cancer. 2008;113(1):48-56.
Wang W, Macapinlac H, Larson SM, Yeh SDJ, Akhurst T, Finn RD, et al. [18F]-2-Fluoro-2-Deoxy-d-Glucose Positron Emission Tomography Localizes Residual Thyroid Cancer in Patients with Negative Diagnostic 131I Whole Body Scans and Elevated Serum Thyroglobulin Levels. The Journal of Clinical Endocrinology & Metabolism. 1999;84(7):2291-302.
Xing M. BRAF mutation in thyroid cancer. Endocrine-related cancer. 2005;12(2):245-62.
Xing M. BRAF mutation in papillary thyroid cancer: pathogenic role, molecular bases, and clinical implications. Endocrine reviews. 2007;28(7):742-62.
Nagarajah J, Ho AL, Tuttle RM, Weber WA, Grewal RK. Correlation of BRAFV600E mutation and glucose metabolism in thyroid cancer patients: an 18F-FDG PET study. Journal of nuclear medicine. 2015;56(5):662.
Agrawal N, Akbani R, Aksoy BA, Ally A, Arachchi H, Asa SL, et al. Integrated genomic characterization of papillary thyroid carcinoma. Cell. 2014;159(3):676-90.
Mamede M, Higashi T, Kitaichi M, Ishizu K, Ishimori T, Nakamoto Y, et al. [18F] FDG uptake and PCNA, Glut-1, and Hexokinase-II expressions in cancers and inflammatory lesions of the lung. Neoplasia. 2005;7(4):369-79.
Haberkorn U, Ziegler SI, Oberdorfer F, Trojan H, Haag D, Peschke P, et al. FDG uptake, tumor proliferation and expression of glycolysis associated genes in animal tumor models. Nuclear medicine and biology. 1994;21(6):827-34.
Choi H, Na KJ. Integrative analysis of imaging and transcriptomic data of the immune landscape associated with tumor metabolism in lung adenocarcinoma: Clinical and prognostic implications. Theranostics. 2018;8(7):1956.
Croft D, Mundo AF, Haw R, Milacic M, Weiser J, Wu G, et al. The Reactome pathway knowledgebase. Nucleic acids research. 2013;42(D1):D472-D7.
Barbie DA, Tamayo P, Boehm JS, Kim SY, Moody SE, Dunn IF, et al. Systematic RNA interference reveals that oncogenic KRAS-driven cancers require TBK1. Nature. 2009;462(7269):108.
Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC bioinformatics. 2013;14(1):7.
Mochizuki T, Tsukamoto E, Kuge Y, Kanegae K, Zhao S, Hikosaka K, et al. FDG uptake and glucose transporter subtype expressions in experimental tumor and inflammation models. Journal of Nuclear Medicine. 2001;42(10):1551-5.
Kim S, Chung JK, Min HS, Kang JH, Park DJ, Jeong JM, et al. Expression patterns of glucose transporter-1 gene and thyroid specific genes in human papillary thyroid carcinoma. Nuclear medicine and molecular imaging. 2014;48(2):91-7.
Durante C, Puxeddu E, Ferretti E, Morisi R, Moretti S, Bruno R, et al. BRAF mutations in papillary thyroid carcinomas inhibit genes involved in iodine metabolism. The Journal of Clinical Endocrinology & Metabolism. 2007;92(7):2840-3.
Romei C, Ciampi R, Faviana P, Agate L, Molinaro E, Bottici V, et al. BRAFV600E mutation, but not RET/PTC rearrangements, is correlated with a lower expression of both thyroperoxidase and sodium iodide symporter genes in papillary thyroid cancer. Endocrine-related cancer. 2008;15(2):511-20.
Choi H, Na KJ. Pan-cancer analysis of tumor metabolic landscape associated with genomic alterations. Molecular cancer. 2018;17(1):150.
Macheda ML, Rogers S, Best JD. Molecular and cellular regulation of glucose transporter (GLUT) proteins in cancer. Journal of cellular physiology. 2005;202(3):654-62.
Grabellus F, Nagarajah J, Bockisch A, Schmid KW, Sheu SY. Glucose transporter 1 expression, tumor proliferation, and iodine/glucose uptake in thyroid cancer with emphasis on poorly differentiated thyroid carcinoma. Clinical nuclear medicine. 2012;37(2):121-7.
de Geus-Oei L-F, van Krieken JHJ, Aliredjo RP, Krabbe PF, Frielink C, Verhagen AF, et al. Biological correlates of FDG uptake in non-small cell lung cancer. Lung cancer. 2007;55(1):79-87.
Meyer HJ, Wienke A, Surov A. Associations between GLUT expression and SUV values derived from FDG-PET in different tumors—A systematic review and meta analysis. PloS one. 2019;14(6):e0217781.
Bos R, Van der Hoeven JJ, Van der Wall E, Van der Groep P, Van Diest PJ, Comans EF, et al. Biologic correlates of 18fluorodeoxyglucose uptake in human breast cancer measured by positron emission tomography. Journal of Clinical Oncology. 2002;20(2):379-87.
Higashi T, Saga T, Nakamoto Y, Ishimori T, Mamede MH, Wada M, et al. Relationship between retention index in dual-phase 18F-FDG PET, and hexokinase-II and glucose transporter-1 expression in pancreatic cancer. Journal of Nuclear Medicine. 2002;43(2):173-80.
Zhao S, Kuge Y, Mochizuki T, Takahashi T, Nakada K, Sato M, et al. Biologic correlates of intratumoral heterogeneity in 18F-FDG distribution with regional expression of glucose transporters and hexokinase-II in experimental tumor. Journal of Nuclear Medicine. 2005;46(4):675.
HaugenBryan R, AlexanderErik K, BibleKeith C, DohertyGerard M, MandelSusan J, NikiforovYuri E, et al. 2015 American Thyroid Association management guidelines for adult patients with thyroid nodules and differentiated thyroid cancer: the American Thyroid Association guidelines task force on thyroid nodules and differentiated thyroid cancer. Thyroid. 2016.

Table 1. Clinicopathological characteristics of the 505 PTC included in the study

Variables	PTC	No. of available data
Gender (M:F)	136 : 369	505
Age (years, mean ± SD)	46.89±15.84 (15.38–89.54)	505
Race
White / Asian / Black / Etc.	334 / 51 / 27 / 1	413
Histology
Classical / Follicular / Tall cell / Other	323 / 99 / 34 / 9	465
BRAFV600E (Positive:Negative)	235 : 270	505
Stage
1 / 2 / 3 / 4	284 / 52 / 112 / 55	503
Recurrence-free survival time (months)	14.26 (0–158.71)	423

PTC = Papillary Thyroid Cancer

Data in parentheses are ranges.

Table 2. Clinicopathological characteristics of the 10 PDTC and 27 ATC included in the study

Variables	PDTC/ATC	No. of available data
Gender (M:F)	10 : 27	37
Age (years, mean ± SD)	69±13.74 (27–87)	37
Histology
Papillary / Follicular / Tall cell / Other	14 / 5 / 5 / 5	29
BRAFV600E (Positive:Negative)	NA
Stage
1 / 2 or 3 / 4	1 / 5 / 24	30
Recurrence-free survival time (months)	10.3 (0.23–223.49)	35

PDTC/ATC = Poorly-Differentiated Thyroid Cancer/Anaplastic Thyroid Cancer

Data in parentheses are ranges.

Table 3. Univariate and multivariate analysis of influence on the recurrence of PTC

	Univariate analysis		Multivariate analysis
Variables *(The first group is the reference group)*	H.R. (95% C.I.)	p-value	H.R. (95% C.I.)	p-value
Age (Continuous variable)	1.01 (0.99–1.04)	0.351	1.01 (0.99–1.04)	0.265
Gender (Female vs. Male)	1.56 (0.71–3.45)	0.272	1.98 (0.83–4.71)	0.124
T-stage (T1+T2 vs. ≥T3)	1.36 (0.63–2.95)	0.435	1.41 (0.61–3.21)	0.420
N-stage (N0 vs. ≥N1)	3.76 (1.40–10.09)	0.008	4.24 (1.50–11.96)	0.006
TDS (Continuous variable)	0.76 (0.55–1.06)	0.102	0.83 (0.58–1.18)	0.293
Glycolysis signature (Continuous variable)	1.50 (1.03–2.17)	0.033	1.98 (1.26–3.10)	0.003

PTC = Papillary Thyroid Cancer; H.R. = Hazard Ratio; C.I. = Confidence Interval; TDS = Tumor Differentiation Score

Table 4. Clinicopathological characteristics of the 423 PTC

Variables	PTCs with low glycolysis signature	PTCs with high glycolysis signature	χ^2*or t^† or H.R.^§/ p-value
Gender (M:F)	59:153	50:161	0.74^*/0.389
Age (years)	43 (15 – 89)	47 (19 – 83)	3.06^†/0.002
Stage
1 / 2 / 3 / 4	131 / 15 / 45 / 20	119 / 23 / 48 / 20	0.995^‡
Recurrence-free survival time (months)	14.59 (0 – 158.71)	14.06 (0 – 157.03)	1.50^§/0.033

^‡p-value from Analysis of Variance

PTC = Papillary Thyroid Cancer; H.R. = Hazard Ratio

Data in parentheses are ranges.

Additional file 1: Figure S1. Plots for GLUT and glycolysis signatures in PTC with different cell types. (A) Box and whiskers plot of GLUT signature in PTC according to cell type (Classical cell type 20.13±1.32 vs. Follicular cell type 18.42±1.20, t = 11.55, p < 0.001) (B) Box and whiskers plot of glycolysis signature in PTC according to cell type (Classical cell type -0.13±0.94 vs. Follicular cell type 0.52±1.12, t = -5.65, p < 0.001) (C) Scatter plot of TDS versus glucose metabolism signature in classical cell type PTC (r = -0.47, p < 0.001 for GLUT; r = 0.23, p < 0.001 for glycolysis) (D) Scatter plot of TDS versus glucose metabolism signature in follicular cell type PTC (r = 0.43, p < 0.001 for GLUT; r = 0.34, p = 0.001 for glycolysis)

Additional file 2: Figure S2. Box and whisker plot of glycolysis signatures by N-stage or M-stage in PTC. The line across each box represent the median, and the top edge, and the bottom edge represents the first quartile, and the third quartile, respectively. Student’s t-test showed significant difference of signatures of glycolysis between N positive and N negative groups (N negative group 0.13±1.03 vs. N positive group -0.24±0.89, t = 3.86, p = 0.0001). No significant difference of signatures of glycolysis were found between M positive and M negative groups (M negative group -0.08±1.03 vs. M positive group -0.16±0.89, t = 0.22, p = 0.82). (*** = p < 0.001)

Download PDF

Journal Publication

published 30 Dec, 2019

Read the published version in BMC Cancer →

Editorial decision: Minor revision
16 Dec, 2019
Submission checks completed at journal
13 Dec, 2019

You are reading this older preprint version

Read the latest preprint version →

Comprehensive gene expression analysis for exploring the association between glucose metabolism and differentiation of thyroid cancer

Status:

Journal Publication

Version 5

Abstract

Figures

Background

Methods

Results

Discussion

Conclusions

List of Abbreviations

Declarations

References

Tables

Additional Files

Supplementary Files

Status:

Journal Publication

Version 5