International Journal of Genomics

International Journal of Genomics / 2019 / Article

Research Article | Open Access

Volume 2019 |Article ID 5894926 | 15 pages |

Identification and Validation of Core Genes Involved in the Development of Papillary Thyroid Carcinoma via Bioinformatics Analysis

Academic Editor: Byung-Hoon Jeong
Received25 Mar 2019
Accepted20 Jun 2019
Published08 Sep 2019


Background. Papillary thyroid carcinoma (PTC) is a common endocrine malignant neoplasm, and its incidence increases continuously worldwide in the recent years. However, efficient clinical biomarkers were still deficient; the present research is aimed at exploring significant core genes of PTC. Methods. We integrated three cohorts to identify hub genes and pathways associated with PTC by comprehensive bioinformatics analysis. Expression profiles GSE33630, GSE35570, and GSE60542, including 114 PTC tissues and 126 normal tissues, were enrolled in this research. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were utilized to search for the crucial biological behaviors and pathways involved in PTC carcinogenesis. Protein-protein interaction (PPI) network was constructed, and significant modules were deeply studied. Results. A total of 831 differentially expressed genes (DEGs) were discovered, comprising 410 upregulated and 421 downregulated genes in PTC tissues compared to normal thyroid tissues. PPI network analysis demonstrated the interactions between those DEGs, and top 10 pivotal genes (TGFB1, CXCL8, LRRK2, CD44, CCND1, JUN, DCN, BCL2, ACACB, and CXCL12) with highest degree of connectivity were extracted from the network and verified by TCGA dataset and RT-PCR experiment of PTC samples. Four of the hub genes (CXCL8, DCN, BCL2, and ACACB) were linked to the prognosis of PTC patients and considered as clinically relevant core genes via survival analysis. Conclusion. In conclusion, we propose a series of key genes associated with PTC development and these genes could serve as the diagnostic biomarkers or therapeutic targets in the future treatment for PTC.

1. Introduction

Thyroid cancer accounts for 96% malignant tumors of the head and neck, which is the most common malignant tumor in endocrine organs and the fifth most common cancer in women worldwide [1, 2]. The incidence of thyroid cancer continues to increase during the past 30 years, which has the fastest growth rate among all solid malignant tumors, generating a serious impact on the quality of life and huge social economic burden for patients [3]. The 2017 edition of American cancer statistics indicated that the incidence of thyroid cancer accounted for 3.4% of all malignant tumors. According to an analysis of thyroid cancer data from the World Health Organization during 1960 to 2012, the incidence of thyroid cancer in most countries and regions around the world presents a rapid growth in the last half century [4]. Thyroid cancer can be classified into four subtypes: papillary thyroid carcinoma (PTC), follicular thyroid carcinoma (FTC), medullary thyroid carcinoma (MTC), and anaplastic thyroid carcinoma (ATC), among which PTC is the most common type, accounting for about 80% of all thyroid cancer. Most of the patients with differentiated thyroid cancer (DTC) can obtain satisfactory therapeutic effects through surgical excision and radioactive iodine therapy and an overall 10-year survival rate exceeding 90% [5, 6]. However, a portion of thyroid cancer patients develop distant metastases and exhibit poor prognosis [7].

At present, the cause of PTC is not very clear. Several studies reported that the incidence of thyroid cancer is influenced by many factors, such as sex, age, region, and ethnic differences [8, 9]. Many studies have reported that the pathogenesis of PTC is related to BRAF mutation and RET/PTC and PAX8/PPAR rearrangement; one hypothesis speculated that the activation of serine/threonine kinase caused by mutations of BRAF at T1799A and V599E may be an alternative mechanism of oncogenic mitogen-activated protein kinase (MAPK) pathway activation [1012], and other groups found that RET/PTC rearrangement promotes the proliferation and migration of PTC cells by regulating nuclear factor-κB (NF-κB) activity and relative oncogenic protein expression [1315]. PAX8/PPAR rearrangement was also involved in thyroid tumorigenesis. PAX8 regulates thyroglobulin, thyroid peroxidase, and thyrotropin receptor gene promoters; it also regulates the differentiation and growth of thyroid cells [16]. PPARγ is one of the nuclear receptors, which participates in inflammation, cell cycle control, apoptosis, and oncogenesis. PAX/PPAR rearrangement intervenes the transcription pathway, upregulates relative gene promoter activity, and finally promotes PTC progression [1719]. However, these theories still cannot explain the pathogenesis of all PTCs. Therefore, further studies of the molecular mechanisms of PTC and seeking novel biomarkers for PTC seem to be of great significance and clinical value.

With wide application of high-throughput sequencing and public databases, the bioinformatics analysis has been broadly used as a powerful tool for cancer diagnosis, classification, and prognosis prediction, especially in obtaining gene general alterations during tumorigenesis. In this study, three original microarray datasets GSE33630, GSE35570, and GSE60542 were selected from Gene Expression Omnibus (GEO, with a total of 114 PTC samples and 126 normal thyroid tissues, which owns the largest sample size compared with a previously published study with a total of 92 PTC samples and 79 normal thyroid tissue samples [20]. With the adhibition of the GEO2R online tool, 831 DEGs (differentially expressed genes) were detected: 410 of those genes were upregulated and 421 were downregulated. Based on those DEGs, the functional roles and related pathways were jointly performed by Database for Annotation, Visualization, and Integrated Discovery (DAVID) and other online tools. In addition, we established the PPI network of the DEGs and ten hub genes (TGFB1, JUN, CXCL8, LRRK2, CD44, DCN, BCL2, ACACB, CCND1, and CXCL12) were identified and verified by RT-PCR experiments in our own thyroid samples in Huashan Hospital. These findings may provide us with a better understanding of the molecular mechanism and possible treatment targets for PTC.

2. Materials and Methods

2.1. Gene Expression Profile Data

The gene expression profile datasets GSE33630, GSE35570, and GSE60542 were downloaded from GEO database. A platform of these three GEO datasets was GPL570 (Affymetrix Human Genome U133 Plus 2.0 Array). The GSE33630 dataset included 94 samples derived from PTC (papillary thyroid cancer) patients (49 tumor tissues and 45 nontumoral thyroid tissues). In the GSE35570 dataset, there consists of 32 PTC tissues and 51 normal samples. As for the GSE60542 dataset, 63 samples were included with 33 tumor tissues and 30 noncancerous tissues. We chose these three datasets for integrated analysis in this study because all of them had a pretty large sample size compared to other datasets about PTC in GEO data repository.

2.2. DEG Identification

A GEO2R ( online program was applied to detect differentially expressed genes between PTC tumor and nontumor tissues [21]. and the adjusted value < 0.05 were set as cutoff criteria. Then, a total of 831 DEGs were identified, containing 410 upregulated genes and 421 downregulated genes, and we selected the top ten genes with highest degree of connectivity as hub genes. For the validation of the expression of these hub genes, the data from GEPIA was extracted and analyzed. The TCGA (The Cancer Genome Atlas) PTC data included 337 normal and 512 PTCs in the GEPIA. The GEPIA ( is a newly developed web-based tool to provide differential expression analysis about numerous kinds of cancer based on the data from TCGA and GTEx (Genotype-Tissue Expression) [22].

2.3. Functional Network Establishment of DEGs

To annotate the upregulated and downregulated DEGs presented from our analysis on the functional level, we discovered the gene ontology (GO) functional enrichment term including biological process (BP), cellular component (CC), molecular function (MF), and vital pathways using the online application DAVID (version 6.8,, Reactome (, and clusterProfiler package in R (, with the threshold set as [23].

2.4. PPI Network and Module Analysis

To determine the relationship between upregulated and downregulated DEGs, we uploaded all DEGs to the STRING database (version 10.5,, choosing a to construct PPI. Then, the Cytoscape software (version 3.6.0, was utilized to construct the network of PPI. The Molecular Complex Detection (MCODE) app in Cytoscape was recruited to analyze modules of the PPI network. The KEGG pathway analysis of genes in each module was performed by DAVID.

2.5. Tissue Samples

32 pairs of primary human PTC cancerous tissues and their corresponding adjacent normal thyroid tissues were obtained from Huashan Hospital of Fudan University. Samples were obtained and snap frozen at liquid nitrogen immediately and stored at -80°C until analysis. All PTC specimens were confirmed by histopathology. Ethical approval was obtained from the ethics committee of Huashan Hospital, Fudan University, and written informed consent was obtained from each patient.

2.6. RNA Extraction and RT-PCR

Total RNA was isolated from PTC cells and tissues using the TRIzol reagent (Invitrogen, USA). cDNA was synthesized using 2.0 μg of total RNA with the PrimeScript™ RT reagent kit (TaKaRa Biotechnology, Japan). The SYBR-Green Supermix kit (TaKaRa) was employed to detect the relative mRNA expression in the ABI 7900 instrument (Applied Biosystems Inc.). The relative expression levels were determined by the 2-ΔΔCt method and normalized to internal control GAPDH [24, 25]. All qPCR reactions were performed in triplicate. The primers used to explore mRNA expression of ten hub genes were shown in Table 1.

GenePrimer sequences












2.7. Statistical Analysis

The expression of related genes was presented as and analyzed using the paired -test. All data were analyzed with GraphPad Prism 7.0, and was considered as a significant difference.

3. Results

3.1. Screening of DEGs between PTC and Normal Thyroid Tissues

Three GEO datasets GSE33630, GSE35570, and GSE60542, a total of 114 PTC tissues and 126 normal thyroid tissues, were analyzed in our research. According to the cutoff criteria: and the adjusted value < 0.05, we extracted 605, 1269, and 548 upregulated DEGs and 555, 1161, and 543 downregulated DEGs from the expression profile datasets GSE33630, GSE35570, and GSE60542, respectively (Figure 1). 410 upregulated DEGs and 421 downregulated DEGs were consistently presented in all three datasets and were used for further study.

3.2. Gene Ontology Functional Analysis

In order to obtain a more in-depth understanding of the selected DEGs, GO function was performed using DAVID version 6.8 and the clusterProfiler package in R statistical software (version 3.4.2). As shown in Table 2 and Figure 2(a), the upregulated DEGs mainly enriched in biological process (BP) terms were associated with extracellular matrix organization, cell adhesion, and regulation of protein serine/threonine kinase activity. Meanwhile, the downregulated DEGs were linked to nervous system development, muscle tissue development, and regulation of phospholipase activity (Figure 2(b)). As for the cellular component (CC), the upregulated DEGs were enriched in the proteinaceous extracellular region, endoplasmic reticulum lumen, and membrane-bounded vesicle, while the downregulated DEGs were related to the proteinaceous extracellular matrix. In the molecular function (MF) group, it was concluded that the upregulated DEGs were mainly involved in receptor binding and glycosaminoglycan binding and the downregulated DEGs were closely related to sulfur compound binding and nucleoside-triphosphatase regulator activity.

ExpressionCategoryIDTermCount value

DownregulatedGOTERM_BP_FATGO:0007167Enzyme-linked receptor protein signaling pathway473.40E-08
GOTERM_BP_FATGO:0007399Nervous system development808.50E-08
GOTERM_BP_FATGO:0061061Muscle structure development341.20E-07
GOTERM_BP_FATGO:0010517Regulation of phospholipase activity114.20E-07
GOTERM_CC_FATGO:0005578Proteinaceous extracellular matrix311.50E-10
GOTERM_CC_FATGO:0031012Extracellular matrix362.20E-09
GOTERM_CC_FATGO:0044421Extracellular region part1241.80E-07
GOTERM_CC_FATGO:0005576Extracellular region1372.50E-06
GOTERM_CC_FATGO:0030018Z disc124.00E-05
GOTERM_MF_FATGO:1901681Sulfur compound binding171.90E-05
GOTERM_MF_FATGO:0005539Glycosaminoglycan binding162.00E-05
GOTERM_MF_FATGO:0008201Heparin binding139.20E-05
GOTERM_MF_FATGO:0005102Receptor binding494.30E-04
GOTERM_MF_FATGO:0019199Transmembrane receptor protein kinase activity81.30E-03

UpregulatedGOTERM_BP_FATGO:0030198Extracellular matrix organization352.10E-14
GOTERM_BP_FATGO:0007155Cell adhesion852.90E-14
GOTERM_BP_FATGO:2000026Regulation of multicellular organismal development786.70E-11
GOTERM_BP_FATGO:0051094Positive regulation of developmental process568.80E-10
GOTERM_CC_FATGO:0005576Extracellular region1602.20E-10
GOTERM_CC_FATGO:0044421Extracellular region part1412.40E-10
GOTERM_CC_FATGO:0031982Membrane-bounded vesicle1153.80E-05
GOTERM_CC_FATGO:0070062Extracellular exosome944.90E-05
GOTERM_CC_FATGO:1903561Extracellular vesicle946.00E-05
GOTERM_MF_FATGO:0005539Glycosaminoglycan binding207.60E-08
GOTERM_MF_FATGO:0005102Receptor binding605.40E-07
GOTERM_MF_FATGO:0050840Extracellular matrix binding101.30E-06
GOTERM_MF_FATGO:0042802Identical protein binding546.00E-06
GOTERM_MF_FATGO:0098772Molecular function regulator525.60E-05

Abbreviations: GO: gene ontology; BP: biological process; MF: molecular function; CC: cellular component.
3.3. KEGG Signaling Pathway Enrichment Analysis

To systematically clarify the key pathways involved in PTC pathogenesis, pathway enrichment analysis was carried out using online websites of the KEGG pathway, Reactome, and the clusterProfiler R package. The pathways enriched by upregulated DEGs were mainly related to ECM-receptor interaction, focal adhesion, and p53 signaling pathway. However, the central pathways shown in downregulated DEGs were mineral absorption, thyroid hormone synthesis, and TGF-beta signaling pathway (Table 3 and Figure 3).

ExpressionCategoryIDTermCount value

DownregulatedKEGG_PATHWAYhsa04978Mineral absorption62.90E-03
KEGG_PATHWAYhsa04918Thyroid hormone synthesis73.80E-03
KEGG_PATHWAYhsa04350TGF-beta signaling pathway63.50E-02
KEGG_PATHWAYhsa00350Tyrosine metabolism43.90E-02
KEGG_PATHWAYhsa04020Calcium signaling pathway93.90E-02

UpregulatedKEGG_PATHWAYhsa04512ECM-receptor interaction151.10E-08
KEGG_PATHWAYhsa04510Focal adhesion172.50E-05
KEGG_PATHWAYhsa04115p53 signaling pathway91.60E-04
KEGG_PATHWAYhsa05205Proteoglycans in cancer152.50E-04
KEGG_PATHWAYhsa04151PI3K-Akt signaling pathway177.20E-03

DownregulatedREACTOME_PATHWAYR-HSA-5661231Metallothioneins bind metals62.33E-05
REACTOME_PATHWAYR-HSA-5660526Response to metal ions61.03E-04
REACTOME_PATHWAYR-HSA-3560782Diseases associated with glycosaminoglycan metabolism70.001726
REACTOME_PATHWAYR-HSA-8851708Signaling by FGFR2 IIIa TM50.001887
REACTOME_PATHWAYR-HSA-2129379Molecules associated with elastic fibres60.002228

UpregulatedREACTOME_PATHWAYR-HSA-1474244Extracellular matrix organization351.39E-08
REACTOME_PATHWAYR-HSA-3000171Nonintegrin membrane-ECM interactions146.32E-08
REACTOME_PATHWAYR-HSA-6785807Interleukin-4 and 13 signaling267.33E-08
REACTOME_PATHWAYR-HSA-216083Integrin cell surface interactions161.25E-07
REACTOME_PATHWAYR-HSA-3000170Syndecan interactions91.29E-06

Abbreviations: KEGG: Kyoto Encyclopedia of Genes and Genomes; PTC: papillary thyroid carcinoma.
3.4. PPI Network Construction and Modular Analysis

The STRING online database was employed to construct the PPI network; then, the Cytoscape software was used to visualize the protein interactions. Considering as the cutoff criterion, a total of 108 DEGs (including 62 upregulated DEGs and 46 downregulated DEGs) in the 831 commonly altered DEGs were filtered to establish the PPI network, which comprises 108 nodes and 769 edges (Figure 4). The PPI network provides us with the most significant ten nodes (TGFB1, CXCL8, LRRK2, CD44, CCND1, JUN, DCN, BCL2, ACACB, and CXCL12), and they were chosen as the hub genes of PTC (Table 4). Furthermore, the whole PPI network was analyzed by MCODE based on the degree of importance. The top three modules were shown in Figure 5. Module 1 consists of 22 nodes and 71 edges, while Module 2 consists of 18 nodes and 54 edges, and there are 19 nodes and 54 edges in Module 3. Moreover, KEGG pathway enrichment analysis showed that genes in these three modules were predominantly associated with pathways in cancer, chemokine signaling pathway, and PI3K-Akt signaling pathway, respectively.

Log2FCAdj. Log2FCAdj. Log2FCAdj.


Notes: The adj. value was calculated using Student’s -test. Abbreviations: PPI, protein-protein interaction.
3.5. Validation of the DEGs by RT-PCR

To reinforce the reliability of the identified hub genes from the bioinformatics analysis, we firstly investigated the expression of these crucial genes in online website GEPIA based on TCGA and GTEx thyroid carcinoma (THCA) data. The results from TCGA and the GTEx projects containing 512 PTC samples and 337 normal samples showed that TGFB1, CXCL8, LRRK2, CD44, and CCND1 were all elevated and JUN, DCN, BCL2, ACACB, and CXCL12 were reduced to some extent in PTC tissues relative to normal tissues (Figure 6(a)). Moreover, RT-PCR experiment was conducted to reconfirm these hub genes. 32 pairs of PTC tissues and their adjacent noncancerous tissues collected in our cohort were used to measure the expression of these molecules. As indicated in Figure 6(b), the results of RT-PCR showed that the mRNA level of TGFB1, CXCL8, LRRK2, CD44, and CCND1 increased by 2.44- (), 3.64- (), 6.67- (), 5.58- (), and 3.43-fold () and the level of JUN, DCN, BCL2, ACACB, and CXCL12 decreased by 2.99- (), 4.53- (), 3.57- (), 3.89- (), and 2.37-fold (), correspondingly. Although the fold changes varied from the data in GEO datasets and TCGA database, the total tendency of these genes from our experiment results was consistent with the databases. Besides, we evaluated the relationship between these hub genes and PTC patients’ prognosis with UALCAN [26]. As shown in Figure 7, Kaplan-Meier survival curves revealed that PTC patients with high expression levels of CXCL8 (), DCN (), BCL2 (), and ACACB () were correlated with a poor overall survival rate, suggested that these four molecules may serve as the prognostic markers for PTC. In brief, our study demonstrated that these ten hub genes may strongly link to the occurrence, progression, and prognosis of PTC and serve as a solid foundation for the future genomic individualized treatment of PTC.

4. Discussion

PTC is one of the most common endocrine-related cancers in the world, and the incidence increases steadily in recent decades [4, 27]. A previous study revealed that genetic alternations such as TERT mutations, BRAF V600E mutations, and RET/PTC rearrangement all contributed to the tumor proliferation and metastasis [28]. The precise mechanism concerning the pathogenesis of PTC remains largely elusive. Therefore, it is of great importance to identify and characterize the pivotal genes associated with the development of PTC. In the current study, our results provide a comprehensive bioinformatics analysis of genes and pathways which may participated in the progression of PTC. We integrated 3 datasets including 114 PTC tissues and 126 normal tissues, which enrolled the largest number of PTC tissue samples in similar bioinformatics studies, and identified 831 DEGs (410 upregulated and 421 downregulated) at the first step [29]. The following GO analysis showed that upregulated DEGs were mainly enriched in the extracellular matrix organization, endoplasmic reticulum lumen, and glycosaminoglycan binding. In addition, downregulated DEGs were significantly involved in the organic hydroxy compound metabolic process, proteinaceous extracellular matrix, and sulfur compound binding. Subsequent KEGG pathway enrichment analysis revealed that PTC development was strongly linked to pathways, such as thyroid hormone synthesis, ECM-receptor interaction, and p53 signaling pathway. In addition, we constructed the PPI network for DEGs and discovered top ten hub genes (TGFB1, CXCL8, LRRK2, CD44, CCND1, JUN, DCN, BCL2, ACACB, and CXCL12) with highest degree of connectivity. Finally, the most significant 3 modules were filter from the PPI network. Corresponding genes among these modules were associated with pathways in cancer, chemokine signaling pathway, and PI3K-Akt signaling pathway. Furthermore, experimental validation is indispensable to confirm our results predicted by bioinformatics analysis. RT-PCR results verified the overexpression of TGFB1, CXCL8, LRRK2, CD44, and CCND1 and the downregulation of JUN, DCN, BCL2, ACACB, and CXCL12 in 32 pairs of PTC samples and their adjacent normal tissues. Among these hub genes, CXCL8, DCN, BCL2, and ACACB were identified as clinically relevant genes using the expression status from the TCGA cohort of 504 PTC samples.

CXCL8, also named IL-8, is a member of the CXC chemokine family and produced by leukocytes and other cells. As a proinflammatory/chemoattractant cytokine, CXCL8 exerts its oncogenic effect by sustaining cell proliferation, metastasis, and angiogenesis [30]. The results of our study were in consistent with several published research papers that have reported the elevated expression of CXCL8 in PTC. Basolo et al. revealed that CXCL8 was produced by TPC-1 and TT cells with activated RET/PTC1 arrangement and its expression was regulated by multiple intracellular signaling pathways [31]. Liotti et al. reported that the CXCL8 expression level was correlated with lymph node metastases and it was overexpressed in thyrospheres and critical for self-renewal ability of thyroid cancer cells [32]. Rotondi et al. discovered that CXCL8 was highly abundant in PTC-derived tumor-associated macrophages (TAMs) and promoted PTC invasion in vivo [33]. Moreover, the role of CXCL8 in thyroid cancer was recently addressed by several studies which characterized its secretion and expression; in particular, these findings were recently reviewed and CXCL8 was proposed as the chemokine playing crucial protumorigenic effects in the cancer microenvironment [34]. Due to the important position of CXCL8 in PTC, Coperchini et al. designed some experiments to test the effect of the BRAF inhibitor (PLX4720) on the basal and TNF-α-induced CXCL8 secretions in BRAFV600E-mutated (BCPAP, 8305C, and 8505C) and RET/PTC-rearranged (TPC-1) thyroid cancer cell lines and in normal human thyrocytes (NHT). The experiment results show that PLX4720 is able to inhibit the secretion of CXCL8 in BRAFV600E-mutated thyroid cancer cells, which indicated that at least some of the antitumor activities of PLX4720 could be exerted through lowering of CXCL8 in the thyroid cancer microenvironment [35, 36].

Decorin (DCN) was known as an extracellular protein and a pan-RTK inhibitor, which participates in cancer cell growth, spread, proinflammatory processes, and antifibrillogenesis [37]. DCN was found underexpressed in breast cancer, colon cancer, follicular thyroid cancer, and follicular variant of papillary thyroid carcinomas [38, 39]. Mechanically, DCN binds to EGFR/ErB2 extracellular functional domains and triggers mitogen-activated protein kinase activation and the increase of cytosolic Ca2+, which accelerates cell cycle arrest and induces intrinsic cell apoptosis, eventually leading to tumor cell growth inhibition [40]. Taken together, the above studies highlighted DCN as a candidate anticancer target for many types of solid tumors. In addition, we have reason to suspect that DCN have an important role in the carcinogenesis of papillary thyroid cancer.

Acetyl-CoA carboxylases (ACC), the rate-limiting enzymes involved in the de novo synthesis of fatty acids, encode two isoforms ACC1 and ACC2 (ACACB) [41]. AMPK directly inhibits ACC1 and ACC2 activity via phosphorylation of ACC1 at Ser117 and ACC2 at Ser222, resulting in the inhibition of fatty acid synthesis and accumulation of fatty acid β-oxidation [42, 43]. A previous study has demonstrated an upregulation of ACC1 in multiple human cancers, which may promote lipogenesis to provide the energy for cancer cell growth [44]. However, we firstly found that the expression of ACC2 was at a very low level in PTC tissues and the specific mechanism of ACC2 in PTC carcinogenesis remains unclear and needs further investigation.

BCL2 is an antiapoptotic oncogene, functioning in the regulation of cell intrinsic apoptosis. There exists numerous BCL2 family members, including the antiapoptotic proteins BCL2 and BCL-xL and proapoptotic proteins Bax, Bak, Bid, Noxa, Puma, Bad, and others [45]. Once received cellular stress, the activation of these BCL2 family members releases apoptogenic proteins from the mitochondria and initiates caspase activation that results in cell apoptosis [46]. Overexpression of prosurvival BCL2 in multiple cancers, such as prostate, ovarian, and lung cancer, is a hallmark for tumorigenesis [4749]. Nevertheless, it was significantly downregulated in colorectal cancer (CRC) and PTC [50, 51]. Wang et al. investigated the relationship between BCL2 polymorphisms and PTC susceptibility and found that BCL2 (-938 C>A) polymorphism showed a protective role in susceptibility to PTC [52]. Gupta et al. reported downregulated BCL2 in anaplastic papillary carcinoma (ATC) as compared to well-differentiated thyroid cancer [53]. Besides, loss of BCL2 was associated with loss of differentiation in thyroid cancer. Our findings suggested that aberrantly expressed BCL2 may have an unknown novel function as a potential tumor suppressor in thyroid cancer progression.

Taken together, we conducted a comprehensive bioinformatics analysis of DEGs in PTC and demonstrated that ten hub genes as the biomarkers of PTC and all these molecules were validated by our RT-PCR experiments. Our findings provide novel insights into the role of potential oncogenes and tumor suppressor genes in PTC and suggest that these findings may have a great clinical significance. Understanding of these core genes in thyroid cancer development as well as the molecular mechanisms involved in this process should be consolidated in the future investigation. Meanwhile, targeting these genes, especially the clinical relevant genes (CXCL8, DCN, BCL2, and ACACB), may provide effective therapies for the treatment of PTC. However, there is still a long way to get these genes transferred to the clinic for stratifying patients, taken as diagnostic biomarkers and even in depth for immunotherapy or for the development for oncovaccine.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors report no conflicts of interest in this work.

Authors’ Contributions

Xiaoyan Li and Jing He contributed equally to this work.


This study was supported by the Science and Technology Commission of Shanghai Municipality (No. 15411952503) and National Nature Science Foundation of China (No. 81602326).


  1. R. L. Siegel, K. D. Miller, and A. Jemal, “Cancer statistics, 2017,” CA: a Cancer Journal for Clinicians, vol. 67, no. 1, pp. 7–30, 2017. View at: Publisher Site | Google Scholar
  2. W. Chen, R. Zheng, P. D. Baade et al., “Cancer statistics in China, 2015,” CA: a Cancer Journal for Clinicians, vol. 66, no. 2, pp. 115–132, 2016. View at: Publisher Site | Google Scholar
  3. C. M. Kitahara and J. A. Sosa, “The changing incidence of thyroid cancer,” Nature Reviews Endocrinology, vol. 12, no. 11, pp. 646–653, 2016. View at: Publisher Site | Google Scholar
  4. C. La Vecchia, M. Malvezzi, C. Bosetti et al., “Thyroid cancer mortality and incidence: a global overview,” International Journal of Cancer, vol. 136, no. 9, pp. 2187–2195, 2015. View at: Publisher Site | Google Scholar
  5. S. Leboulleux, C. Rubino, E. Baudin et al., “Prognostic factors for persistent or recurrent disease of papillary thyroid carcinoma with neck lymph node metastases and/or tumor extension beyond the thyroid capsule at initial diagnosis,” The Journal of Clinical Endocrinology and Metabolism, vol. 90, no. 10, pp. 5723–5729, 2005. View at: Publisher Site | Google Scholar
  6. M. Ansari, A. A. Babaei, B. Shafiei et al., “Pathological evaluation of differentiated thyroid cancer in patients with positive serum thyroglobulin and negative iodine scan,” European Review for Medical and Pharmacological Sciences, vol. 18, no. 13, pp. 1925–1929, 2014. View at: Google Scholar
  7. M. E. Cabanillas, D. G. McFadden, and C. Durante, “Thyroid cancer,” The Lancet, vol. 388, no. 10061, pp. 2783–2795, 2016. View at: Publisher Site | Google Scholar
  8. J. F. Carrillo, M. Friasmendivil, F. J. Ochoa-Carrillo, and M. Ibarra, “Accuracy of fine-needle aspiration biopsy of the thyroid combined with an evaluation of clinical and radiologic factors,” Otolaryngology - Head and Neck Surgery, vol. 122, no. 6, pp. 917–921, 2000. View at: Publisher Site | Google Scholar
  9. C. Are and A. R. Shaha, “Anaplastic thyroid carcinoma: biology, pathogenesis, prognostic factors, and treatment approaches,” Annals of Surgical Oncology, vol. 13, no. 4, pp. 453–464, 2006. View at: Publisher Site | Google Scholar
  10. Y. Cohen, M. Xing, E. Mambo et al., “BRAF mutation in papillary thyroid carcinoma,” Journal of the National Cancer Institute, vol. 95, no. 8, pp. 625–627, 2003. View at: Publisher Site | Google Scholar
  11. P. Soares, V. Trovisco, A. S. Rocha et al., “BRAF mutations and RET/PTC rearrangements are alternative events in the etiopathogenesis of PTC,” Oncogene, vol. 22, no. 29, pp. 4578–4580, 2003. View at: Publisher Site | Google Scholar
  12. T. Fukushima, S. Suzuki, M. Mashiko et al., “BRAF mutations in papillary carcinomas of the thyroid,” Oncogene, vol. 22, no. 41, pp. 6455–6457, 2003. View at: Publisher Site | Google Scholar
  13. D. Zhou, Z. Li, and X. Bai, “BRAFV600E and RET/PTC promote proliferation and migration of papillary thyroid carcinoma cells in vitro by regulating nuclear factor-κB,” Medical Science Monitor, vol. 23, pp. 5321–5329, 2017. View at: Publisher Site | Google Scholar
  14. J. A. Knauf, H. Kuroda, S. Basu, and J. A. Fagin, “RET/PTC-induced dedifferentiation of thyroid cells is mediated through Y1062 signaling through SHC-RAS-MAP kinase,” Oncogene, vol. 22, no. 28, pp. 4406–4412, 2003. View at: Publisher Site | Google Scholar
  15. R. J. Neely, M. S. Brose, C. M. Gray et al., “The RET/PTC3 oncogene activates classical NF-κB by stabilizing NIK,” Oncogene, vol. 30, no. 1, pp. 87–96, 2011. View at: Publisher Site | Google Scholar
  16. A. Y. M. Au, C. McBride, K. G. Wilhelm Jr. et al., “PAX8-peroxisome proliferator-activated receptor γ (PPARγ) disrupts normal PAX8 or PPARγ transcriptional function and stimulates follicular thyroid cell growth,” Endocrinology, vol. 147, no. 1, pp. 367–376, 2006. View at: Publisher Site | Google Scholar
  17. L. Lacroix, V. Lazar, S. Michiels et al., “Follicular thyroid tumors with the PAX8-PPARgamma1 rearrangement display characteristic genetic alterations,” The American Journal of Pathology, vol. 167, no. 1, pp. 223–231, 2005. View at: Publisher Site | Google Scholar
  18. W. O. Lui, T. Foukakis, J. Lidén et al., “Expression profiling reveals a distinct transcription signature in follicular thyroid carcinomas with a PAX8-PPARγ fusion oncogene,” Oncogene, vol. 24, no. 8, pp. 1467–1476, 2005. View at: Publisher Site | Google Scholar
  19. J. Di Cristofaro, M. Silvy, A. Lanteaume, M. Marcy, P. Carayon, and C. De Micco, “Expression of tpo mRNA in thyroid tumors: quantitative PCR analysis and correlation with alterations of ret, Braf , ras and pax8 genes,” Endocrine-Related Cancer, vol. 13, no. 2, pp. 485–495, 2006. View at: Publisher Site | Google Scholar
  20. W. Liang and F. Sun, “Identification of key genes of papillary thyroid cancer using integrated bioinformatics analysis,” Journal of Endocrinological Investigation, vol. 41, no. 10, pp. 1237–1245, 2018. View at: Publisher Site | Google Scholar
  21. S. Davis and P. S. Meltzer, “GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor,” Bioinformatics, vol. 23, no. 14, pp. 1846-1847, 2007. View at: Publisher Site | Google Scholar
  22. Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic Acids Research, vol. 45, no. W1, pp. W98–W102, 2017. View at: Publisher Site | Google Scholar
  23. G. Yu, L. G. Wang, Y. Han, and Q. Y. He, “clusterProfiler: an R package for comparing biological themes among gene clusters,” OMICS, vol. 16, no. 5, pp. 284–287, 2012. View at: Publisher Site | Google Scholar
  24. F. He, P. Lv, X. Zhao et al., “Predictive value of circulating miR-328 and miR-134 for acute myocardial infarction,” Molecular and Cellular Biochemistry, vol. 394, no. 1-2, pp. 137–144, 2014. View at: Publisher Site | Google Scholar
  25. F. C. He, W. W. Meng, Y. H. Qu et al., “Expression of circulating microRNA-20a and let-7a in esophageal squamous cell carcinoma,” World Journal of Gastroenterology, vol. 21, no. 15, pp. 4660–4665, 2015. View at: Publisher Site | Google Scholar
  26. D. S. Chandrashekar, B. Bashel, S. A. H. Balasubramanya et al., “UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses,” Neoplasia, vol. 19, no. 8, pp. 649–658, 2017. View at: Publisher Site | Google Scholar
  27. J. He, Y. Jin, M. Zhou et al., “Solute carrier family 35 member F2 is indispensable for papillary thyroid carcinoma progression through activation of transforming growth factor-β type I receptor/apoptosis signal-regulating kinase 1/mitogen-activated protein kinase signaling axis,” Cancer Science, vol. 109, no. 3, pp. 642–655, 2018. View at: Publisher Site | Google Scholar
  28. M. Xing, “Molecular pathogenesis and mechanisms of thyroid cancer,” Nature Reviews Cancer, vol. 13, no. 3, pp. 184–199, 2013. View at: Publisher Site | Google Scholar
  29. Y. Zhao, X. Liu, L. Zhong et al., “The combined use of miRNAs and mRNAs as biomarkers for the diagnosis of papillary thyroid carcinoma,” International Journal of Molecular Medicine, vol. 36, no. 4, pp. 1097–1103, 2015. View at: Publisher Site | Google Scholar
  30. M. Rotondi, F. Coperchini, and L. Chiovato, “CXCL8 in thyroid disease: from basic notions to potential applications in clinical practice,” Cytokine & Growth Factor Reviews, vol. 24, no. 6, pp. 539–546, 2013. View at: Publisher Site | Google Scholar
  31. F. Basolo, R. Giannini, A. Toniolo et al., “Establishment of a non-tumorigenic papillary thyroid cell line (FB-2) carrying the RET/PTC1 rearrangement,” International Journal of Cancer, vol. 97, no. 5, pp. 608–614, 2002. View at: Publisher Site | Google Scholar
  32. F. Liotti, F. Collina, E. Pone et al., “Interleukin-8, but not the related chemokine CXCL1, sustains an autocrine circuit necessary for the properties and functions of thyroid cancer stem cells,” Stem Cells, vol. 35, no. 1, pp. 135–146, 2017. View at: Publisher Site | Google Scholar
  33. M. Rotondi, F. Coperchini, F. Latrofa, and L. Chiovato, “Role of chemokines in thyroid cancer microenvironment: is CXCL8 the main player?” Frontiers in Endocrinology, vol. 9, p. 314, 2018. View at: Publisher Site | Google Scholar
  34. F. Coperchini, L. Croce, M. Denegri et al., “The BRAF-inhibitor PLX4720 inhibits CXCL8 secretion in BRAFV600E mutated and normal thyroid cells: a further anti-cancer effect of BRAF-inhibitors,” Scientific Reports, vol. 9, no. 1, article 4390, 2019. View at: Publisher Site | Google Scholar
  35. O. Awwad, F. Coperchini, P. Pignatti et al., “The AMPK-activator AICAR in thyroid cancer: effects on CXCL8 secretion and on CXCL8-induced neoplastic cell migration,” Journal of Endocrinological Investigation, vol. 41, no. 11, pp. 1275–1282, 2018. View at: Publisher Site | Google Scholar
  36. W. Fang, L. Ye, L. Shen et al., “Tumor-associated macrophages promote the metastatic potential of thyroid papillary cancer by releasing CXCL8,” Carcinogenesis, vol. 35, no. 8, pp. 1780–1787, 2014. View at: Publisher Site | Google Scholar
  37. S. Goldoni and R. V. Iozzo, “Tumor microenvironment: modulation by decorin and related molecules harboring leucine-rich tandem motifs,” International Journal of Cancer, vol. 123, no. 11, pp. 2473–2479, 2008. View at: Publisher Site | Google Scholar
  38. G. Li, M. Li, X. Liang et al., “Identifying DCN and HSPD1 as potential biomarkers in colon cancer using 2D-LC-MS/MS combined with iTRAQ technology,” Journal of Cancer, vol. 8, no. 3, pp. 479–489, 2017. View at: Publisher Site | Google Scholar
  39. L. A. T. Arnaldi, R. C. Borra, R. M. B. Maciel, and J. M. Cerutti, “Gene expression profiles reveal that DCN, DIO1, and DIO2 are underexpressed in benign and malignant thyroid tumors,” Thyroid, vol. 15, no. 3, pp. 210–221, 2005. View at: Publisher Site | Google Scholar
  40. W. Zhang, Y. Ge, Q. Cheng, Q. Zhang, L. Fang, and J. Zheng, “Decorin is a pivotal effector in the extracellular matrix and tumour microenvironment,” Oncotarget, vol. 9, no. 4, pp. 5480–5491, 2018. View at: Publisher Site | Google Scholar
  41. C. Wang, S. Rajput, K. Watabe, D. F. Liao, and D. Cao, “Acetyl-CoA carboxylase-a as a novel target for cancer therapy,” Frontiers in Bioscience, vol. 2, pp. 515–526, 2010. View at: Google Scholar
  42. R. U. Svensson, S. J. Parker, L. J. Eichner et al., “Inhibition of acetyl-CoA carboxylase suppresses fatty acid synthesis and tumor growth of non-small-cell lung cancer in preclinical models,” Nature Medicine, vol. 22, no. 10, pp. 1108–1119, 2016. View at: Publisher Site | Google Scholar
  43. G. R. Steinberg and B. E. Kemp, “AMPK in health and disease,” Physiological Reviews, vol. 89, no. 3, pp. 1025–1078, 2009. View at: Publisher Site | Google Scholar
  44. J. V. Swinnen, K. Brusselmans, and G. Verhoeven, “Increased lipogenesis in cancer cells: new players, novel targets,” Current Opinion in Clinical Nutrition and Metabolic Care, vol. 9, no. 4, pp. 358–365, 2006. View at: Publisher Site | Google Scholar
  45. A. Ruefli-Brasse and J. C. Reed, “Therapeutics targeting Bcl-2 in hematological malignancies,” The Biochemical Journal, vol. 474, no. 21, pp. 3643–3657, 2017. View at: Publisher Site | Google Scholar
  46. Z. Jin and W. S. El-Deiry, “Overview of cell death signaling pathways,” Cancer Biology & Therapy, vol. 4, no. 2, pp. 139–163, 2005. View at: Publisher Site | Google Scholar
  47. S. D. Catz and J. L. Johnson, “BCL-2 in prostate cancer: a minireview,” Apoptosis, vol. 8, no. 1, pp. 29–37, 2003. View at: Publisher Site | Google Scholar
  48. N. S. Anderson, L. Turner, S. Livingston, R. Chen, S. V. Nicosia, and P. A. Kruk, “Bcl-2 expression is altered with ovarian tumor progression: an immunohistochemical evaluation,” J Ovarian Res., vol. 2, no. 1, p. 16, 2009. View at: Publisher Site | Google Scholar
  49. V. K. Anagnostou, F. J. Lowery, V. Zolota et al., “High expression of BCL-2 predicts favorable outcome in non-small cell lung cancer patients with non squamous histology,” BMC Cancer, vol. 10, no. 1, p. 186, 2010. View at: Publisher Site | Google Scholar
  50. S. Wu, F. Wu, and Z. Jiang, “Identification of hub genes, key miRNAs and potential molecular mechanisms of colorectal cancer,” Oncology Reports, vol. 38, no. 4, pp. 2043–2050, 2017. View at: Publisher Site | Google Scholar
  51. T. G. Radu, L. Mogoanta, C. J. Busuioc, C. Stanescu, and F. Grosu, “Histological and immunohistochemical aspects of papillary thyroid cancer,” Romanian Journal of Morphology and Embryology, vol. 56, 2 Supplement, pp. 789–795, 2015. View at: Google Scholar
  52. Y. X. Wang, L. Zhao, X. Y. Wang, C. M. Liu, and S. G. Yu, “Role of Caspase 8, Caspase 9 and Bcl-2 polymorphisms in papillary thyroid carcinoma risk in Han Chinese population,” Medical Oncology, vol. 29, no. 4, pp. 2445–2451, 2012. View at: Publisher Site | Google Scholar
  53. A. Gupta, S. Jain, N. Khurana, and A. K. Kakar, “Expression of p63 and Bcl-2 in malignant thyroid tumors and their correlation with other diagnostic immunocytochemical markers,” Journal of Clinical and Diagnostic Research, vol. 10, no. 7, pp. EC04–EC08, 2016. View at: Publisher Site | Google Scholar

Copyright © 2019 Xiaoyan Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

878 Views | 384 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.