Integrated Bioinformatics Analysis of Master Regulators in Anaplastic Thyroid Carcinoma
Open Access
- 28 April 2019
- journal article
- research article
- Published by Wiley in BioMed Research International
- Vol. 2019, 1-13
Abstract
Anaplastic thyroid carcinoma (ATC) is one of the most aggressive and rapidly lethal tumors. However, limited advances have been made to prolong the survival and to reduce the mortality over the last decades. Therefore, identifying the master regulators underlying ATC progression is desperately needed. In our present study, three datasets including GSE33630, GSE29265, and GSE65144 were retrieved from Gene Expression Omnibus with a total of 32 ATC samples and 78 normal thyroid tissues. A total of 1804 consistently changed differentially expressed genes (DEGs) were identified from three datasets. KEGG pathways enrichment suggested that upregulated DEGs were mainly enriched in ECM-receptor interaction, cell cycle, PI3K-Akt signaling pathway, focal adhesion, and p53 signaling pathway. Furthermore, key gene modules in PPI network were identified by Cytoscape plugin MCODE and they were mainly associated with DNA replication, cell cycle process, collagen fibril organization, and regulation of leukocyte migration. Additionally, TOP2A, CDK1, CCNB1, VEGFA, BIRC5, MAPK1, CCNA2, MAD2L1, CDC20, and BUB1 were identified as hub genes of the PPI network. Interestingly, module analysis showed that 8 out of 10 hub genes participated in Module 1 network and more than 70% genes of Module 2 consisted of collagen family members. Notably, transcription factors (TFs) regulatory network analysis indicated that E2F7, FOXM1, and NFYB were master regulators of Module 1, while CREB3L1 was the master regulator of Module 2. Experimental validation showed that CREB3L1, E2F7, and FOXM1 were significantly upregulated in ATC tissue and cell line when compared with normal thyroid group. In conclusion, the TFs regulatory network provided a more detail molecular mechanism underlying ATC occurrence and progression. TFs including E2F7, FOXM1, CREB3L1, and NFYB were likely to be master regulators of ATC progression, suggesting their potential role as molecular therapeutic targets in ATC treatment.Keywords
Funding Information
- National Natural Science Foundation of China (81703578, 81802673, LQ18H160017, 2017KY023, 2019KY037)
- National Natural Science Foundation of China (81703578, 81802673, LQ18H160017, 2017KY023, 2019KY037)
- Natural Science Foundation of Zhejiang Province (81703578, 81802673, LQ18H160017, 2017KY023, 2019KY037)
- Zhejiang Provincial Projects of Medical and Health Technology Program (81703578, 81802673, LQ18H160017, 2017KY023, 2019KY037)
- Zhejiang Provincial Projects of Medical and Health Technology Program (81703578, 81802673, LQ18H160017, 2017KY023, 2019KY037)
This publication has 40 references indexed in Scilit:
- OASIS/CREB3L1 Is Induced by Endoplasmic Reticulum Stress in Human Glioma Cell Lines and Contributes to the Unfolded Protein Response, Extracellular Matrix Production and Cell MigrationPLOS ONE, 2013
- mRNA Expression in Papillary and Anaplastic Thyroid Carcinoma: Molecular Anatomy of a Killing SwitchPLOS ONE, 2012
- FOXM1 is a molecular determinant of the mitogenic and invasive phenotype of anaplastic thyroid carcinomaEndocrine-Related Cancer, 2012
- A general method to derive robust organ-specific gene expression-based differentiation indices: application to thyroid cancer diagnosticOncogene, 2012
- Hallmarks of Cancer: The Next GenerationCell, 2011
- Cell Surface and Secreted Protein Profiles of Human Thyroid Cancer Cell Lines Reveal Distinct Glycoprotein PatternsJournal of Proteome Research, 2009
- ClueGO: a Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networksBioinformatics, 2009
- Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene listsNucleic Acids Research, 2008
- Comprehensive gene expression profiling of anaplastic thyroid cancers with cDNA microarray of 25 344 genesEndocrine-Related Cancer, 2004
- Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT MethodMethods, 2001