scholarly journals Integrated bioinformatics analysis and screening of hub genes in papillary thyroid carcinoma

PLoS ONE ◽  
2021 ◽  
Vol 16 (6) ◽  
pp. e0251962
Author(s):  
Rong Fan ◽  
Lijin Dong ◽  
Ping Li ◽  
Xiaoming Wang ◽  
Xuewei Chen

Background With the increasing incidence of papillary thyroid carcinoma (PTC), PTC continues to garner attention worldwide; however its pathogenesis remains to be elucidated. The purpose of this study was to explore key biomarkers and potential new therapeutic targets for, PTC. Methods GEO2R and Venn online software were used for screening of differentially expressed genes. Hub genes were screened via STRING and Cytoscape, followed by Gene Ontology and KEGG enrichment analysis. Finally, survival analysis and expression validation were performed using the UALCAN online software and immunohistochemistry. Results We identified 334 consistently differentially expressed genes (DEGs) comprising 136 upregulated and 198 downregulated genes. Gene Ontology enrichment analysis results suggested that the DEGs were mainly enriched in cancer-related pathways and functions. PPI network visualization was performed and 17 upregulated and 13 downregulated DEGs were selected. Finally, the expression verification and overall survival analysis conducted using the Gene Expression Profiling Interactive Analysis Tool (GEPIA) and UALCAN showed that LPAR5, TFPI, and ENTPD1 were associated with the development of PTC and the prognosis of PTC patients, and the expression of LPAR5, TFPI and ENTPD1 was verified using a tissue chip. Conclusions In summary, the hub genes and pathways identified in the present study not only provide information for the development of new biomarkers for PTC but will also be useful for elucidation of the pathogenesis of PTC.

Genes ◽  
2019 ◽  
Vol 10 (1) ◽  
pp. 45 ◽  
Author(s):  
Junliang Shang ◽  
Qian Ding ◽  
Shasha Yuan ◽  
Jin-Xing Liu ◽  
Feng Li ◽  
...  

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer. Identifying characteristic genes of PTC are of great importance to reveal its potential genetic mechanisms. In this paper, we proposed a framework, as well as a measure named Normalized Centrality Measure (NCM), to identify characteristic genes of PTC. The framework consisted of four steps. First, both up-regulated genes and down-regulated genes, collectively called differentially expressed genes (DEGs), were screened and integrated together from four datasets, that is, GSE3467, GSE3678, GSE33630, and GSE58545; second, an interaction network of DEGs was constructed, where each node represented a gene and each edge represented an interaction between linking nodes; third, both traditional measures and the NCM measure were used to analyze the topological properties of each node in the network. Compared with traditional measures, more genes related to PTC were identified by the NCM measure; fourth, by mining the high-density subgraphs of this network and performing Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, several meaningful results were captured, most of which were demonstrated to be associated with PTC. The experimental results proved that this network framework and the NCM measure are useful for identifying more characteristic genes of PTC.


2020 ◽  
Author(s):  
Rong Fan ◽  
Lijin Dong ◽  
Ping Li ◽  
Xiaoming Wang ◽  
Xuewei Chen

Abstract Background With the increasing incidence, papillary thyroid cancer (PTC) is receiving more and more attention, but the pathogenesis of which is still not completely elucidated. The purpose of this study was to explore key biomarkers and new therapeutic targets in PTC. Methods GEO2R and Venn online software were used for differential gene screening analysis. Hub genes were screened via STRING and Cytoscape, following Gene Ontology and KEGG enrichment analysis. Finally, survival analysis and expression validation were performed via UALCAN online software and immunohistochemistry. Results We screened 334 consistently differentially expressed genes (DEGs), composed of 136 upregulated genes and 198 downregulated genes. Gene ontology enrichment analysis suggested that DEGs mainly enriched in the cancer-related pathways and functions. PPI network visualization was performed to select 17 upregulated and 13 downregulated DEGs. Finally, the expression verification and overall survival analysis conducted in the Gene Expression Profiling Interactive Analysis Tool (GEPIA) and UALCAN showed that LPAR5, TFPI and ENTPD1 were related to the development of PTC and the prognosis of PTC patients, and the expression of LPAR5 was verified by tissue chip. Conclusions In summary, the hub genes and pathways identified in the present study not only provided new biomarkers for PTC, but also will be useful for elucidating the pathogenesis of PTC.


2020 ◽  
Vol 23 (6) ◽  
pp. 546-553
Author(s):  
Hongyuan Cui ◽  
Mingwei Zhu ◽  
Junhua Zhang ◽  
Wenqin Li ◽  
Lihui Zou ◽  
...  

Objective: Next-generation sequencing (NGS) was performed to identify genes that were differentially expressed between normal thyroid tissue and papillary thyroid carcinoma (PTC). Materials & Methods: Six candidate genes were selected and further confirmed with quantitative real-time polymerase chain reaction (qRT-PCR), and immunohistochemistry in samples from 24 fresh thyroid tumors and adjacent normal tissues. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was used to investigate signal transduction pathways of the differentially expressed genes. Results: In total, 1690 genes were differentially expressed between samples from patients with PTC and the adjacent normal tissue. Among these, SFRP4, ZNF90, and DCN were the top three upregulated genes, whereas KIRREL3, TRIM36, and GABBR2 were downregulated with the smallest p values. Several pathways were associated with the differentially expressed genes and involved in cellular proliferation, cell migration, and endocrine system tumor progression, which may contribute to the pathogenesis of PTC. Upregulation of SFRP4, ZNF90, and DCN at the mRNA level was further validated with RT-PCR, and DCN expression was further confirmed with immunostaining of PTC samples. Conclusion: These results provide new insights into the molecular mechanisms of PTC. Identification of differentially expressed genes should not only improve the tumor signature for thyroid tumors as a diagnostic biomarker but also reveal potential targets for thyroid tumor treatment.


2016 ◽  
Vol 2016 ◽  
pp. 1-9 ◽  
Author(s):  
Huairong Zhang ◽  
Bo Gao ◽  
Bingyin Shi

Aim. We aim to identify protein kinases involved in the pathophysiology of papillary thyroid carcinoma (PTC) in order to provide potential therapeutic targets for kinase inhibitors and unfold possible molecular mechanisms.Materials and Methods. The gene expression profile of GSE27155 was analyzed to identify differentially expressed genes and mapped onto human protein kinases database. Correlation of kinases with PTC was addressed by systematic literature search, GO and KEGG pathway analysis.Results. The functional enrichment analysis indicated that “mitogen-activated protein kinases pathway” expression was extremely enriched, followed by “neurotrophin signaling pathway,” “focal adhesion,” and “GnRH signaling pathway.” MAPK, SRC, PDGFRa, ErbB, and EGFR were significantly regulated to correct these pathways. Kinases investigated by the literature on carcinoma were considered to be potential novel molecular therapeutic target in PTC and application of corresponding kinase inhibitors could be possible therapeutic tool.Conclusion. SRC, MAPK, and EGFR were the most important differentially expressed kinases in PTC. Combined inhibitors may have high efficacy in PTC treatment by targeting these kinases.


Author(s):  
Zheng Zhang ◽  
Shuangshuang Zhao ◽  
Keke Wang ◽  
Mengyuan Shang ◽  
Zheming Chen ◽  
...  

Integrated analysis of accumulated data is an effective way to obtain reliable potential diagnostic molecular of cervical lymph node metastases (LNM) in papillary thyroid carcinoma (PTC). The benefits of prophylactic lymph node dissection (PLND) for these clinically node-negative (cN0) patients remained considerable controversies. Hence, elucidation of the mechanisms of LNM and exploration of potential biomarkers and prognostic indicators are essential for accurate diagnosis of LNM in PTC patients. Up to date, advanced microarray and bioinformatics analysis have advanced an understanding of the molecular mechanisms of disease occurrence and development, which are necessary to explore genetic changes and identify potential diagnostic biomarkers. In present study, we performed a comprehensive analysis of the differential expression, biological functions, and interactions of LNM-related genes. Two publicly available microarray datasets GSE60542 and GSE129562 were available from Gene Expression Omnibus (GEO) database. Differentially expressed genes between clinically node-positive (cN1) and cN0 PTC samples were screened by an integrated analysis of multiple gene expression profile after gene reannotation and batch normalization. Our results identified 48 differentially expressed genes (DEGs) genetically associated with LNM in PTC patients. Gene ontology (GO) analyses revealed the changes in the modules were mostly enriched in the regulation of MHC class II receptor activity, the immune receptor activity, and the peptide antigen binding. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis of DEGs displayed that the intestinal immune network for IgA production, staphylococcus aureus infection, and cell adhesion molecules (CAMs). To screen core genes related to LNM of PTC from the protein-protein interaction network, top 10 hub genes were identified with highest scores. Our results help us understand the exact mechanisms underlying the metastasis of cervical LNM in PTC tissues and pave an avenue for the progress of precise medicine for individual patients.


2020 ◽  
Vol 40 (8) ◽  
Author(s):  
Gang Xue ◽  
Xu Lin ◽  
Jing-Fang Wu ◽  
Da Pei ◽  
Dong-Mei Wang ◽  
...  

Abstract Background: Papillary thyroid carcinoma (PTC) is one of the fastest-growing malignant tumor types of thyroid cancer. Therefore, identifying the interaction of genes in PTC is crucial for elucidating its pathogenesis and finding more specific molecular biomarkers. Methods: Four pairs of PTC tissues and adjacent tissues were sequenced using RNA-Seq, and 3745 differentially expressed genes were screened (P<0.05, |logFC|>1). The enrichment analysis indicated that the vast majority of differentially expressed genes (DEGs) may play a positive role in the development of cancer. Then, the significant modules were analyzed using Cytoscape software in the protein–protein interaction network. Survival analysis, TNM analysis, and immune infiltration analysis of key genes were analyzed. And the expression of ADORA1, APOE, and LPAR5 genes were verified by qPCR in PTC compared with matching adjacent tissues. Results: Twenty-five genes were identified as hub genes with nodes greater than 10. The expression of 25 genes were verified by the GEPIA database, and the overall survival and disease-free survival analyses were conducted with Kaplan–Meier plotter. We found only three genes were confirmed with our validation and were statistically significant in PTC, namely ADORA1, APOE, and LPAR5. Further analysis found that the mRNA levels and methylation degree of these three genes were significantly correlated with the TNM staging of PTC. And these three genes were related to PTC immune infiltration. Verification of the expression of these three genes by RT-qPCR and Western blot further confirmed the reliability of our results. Conclusion: Our study identified three genes that may play key regulatory roles in the development, metastasis, and immune infiltration of papillary thyroid carcinoma.


2020 ◽  
Vol 15 ◽  
Author(s):  
Wei Han ◽  
Dongchen Lu ◽  
Chonggao Wang ◽  
Mengdi Cui ◽  
Kai Lu

Background: In the past decades, the incidence of thyroid cancer (TC) has been gradually increasing, owing to the widespread use of ultrasound scanning devices. However, the key mRNAs, miRNAs, and mRNA-miRNA network in papillary thyroid carcinoma (PTC) has not been fully understood. Material and Methods: In this study, multiple bioinformatics methods were employed, including differential expression analysis, gene set enrichment analysis, and miRNA-mRNA interaction network construction. Results: First, we investigated the key miRNAs that regulated significantly more differentially expressed genes based on GSEA method. Second, we searched for the key miRNAs based on the mRNA-miRNA interaction subnetwork involved in PTC. We identified hsa-mir-1275, hsa-mir-1291, hsa-mir-206 and hsa-mir-375 as the key miRNAs involved in PTC pathogenesis. Conclusion: The integrated analysis of the gene and miRNA expression data not only identified key mRNAs, miRNAs, and mRNA-miRNA network involved in papillary thyroid carcinoma, but also improved our understanding of the pathogenesis of PTC.


BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Ruoting Lin ◽  
Conor E. Fogarty ◽  
Bowei Ma ◽  
Hejie Li ◽  
Guoying Ni ◽  
...  

Abstract Background Papillary thyroid carcinoma (PTC) is the most common thyroid cancer. While many patients survive, a portion of PTC cases display high aggressiveness and even develop into refractory differentiated thyroid carcinoma. This may be alleviated by developing a novel model to predict the risk of recurrence. Ferroptosis is an iron-dependent form of regulated cell death (RCD) driven by lethal accumulation of lipid peroxides, is regulated by a set of genes and shows a variety of metabolic changes. To elucidate whether ferroptosis occurs in PTC, we analyse the gene expression profiles of the disease and established a new model for the correlation. Methods The thyroid carcinoma (THCA) datasets were downloaded from The Cancer Genome Atlas (TCGA), UCSC Xena and MisgDB, and included 502 tumour samples and 56 normal samples. A total of 60 ferroptosis related genes were summarised from MisgDB database. Gene set enrichment analysis (GSEA) and Gene set variation analysis (GSVA) were used to analyse pathways potentially involving PTC subtypes. Single sample GSEA (ssGSEA) algorithm was used to analyse the proportion of 28 types of immune cells in the tumour immune infiltration microenvironment in THCA and the hclust algorithm was used to conduct immune typing according to the proportion of immune cells. Spearman correlation analysis was performed on the ferroptosis gene expression and the correlation between immune infiltrating cells proportion. We established the WGCNA to identify genes modules that are highly correlated with the microenvironment of immune invasion. DEseq2 algorithm was further used for differential analysis of sequencing data to analyse the functions and pathways potentially involving hub genes. GO and KEGG enrichment analysis was performed using Clusterprofiler to explore the clinical efficacy of hub genes. Univariate Cox analysis was performed for hub genes combined with clinical prognostic data, and the results was included for lasso regression and constructed the risk regression model. ROC curve and survival curve were used for evaluating the model. Univariate Cox analysis and multivariate Cox analysis were performed in combination with the clinical data of THCA and the risk score value, the clinical efficacy of the model was further evaluated. Results We identify two subtypes in PTC based on the expression of ferroptosis related genes, with the proportion of cluster 1 significantly higher than cluster 2 in ferroptosis signature genes that are positively associated. The mutations of Braf and Nras are detected as the major mutations of cluster 1 and 2, respectively. Subsequent analyses of TME immune cells infiltration indicated cluster 1 is remarkably richer than cluster 2. The risk score of THCA is in good performance evaluated by ROC curve and survival curve, in conjunction with univariate Cox analysis and multivariate Cox analysis results based on the clinical data shows that the risk score of the proposed model could be used as an independent prognostic indicator to predict the prognosis of patients with papillary thyroid cancer. Conclusions Our study finds seven crucial genes, including Ac008063.2, Apoe, Bcl3, Acap3, Alox5ap, Atxn2l and B2m, and regulation of apoptosis by parathyroid hormone-related proteins significantly associated with ferroptosis and immune cells in PTC, and we construct the risk score model which can be used as an independent prognostic index to predict the prognosis of patients with PTC.


Sign in / Sign up

Export Citation Format

Share Document