Dataset Information

Identification of potential functional genes in papillary thyroid cancer by co-expression network analysis.

ABSTRACT: Interactions between multiple genes are involved in the development of complex diseases. However, there are few analyses of gene interactions associated with papillary thyroid cancer (PTC). Weighted gene co-expression network analysis (WGCNA) is a novel and powerful method that detects gene interactions according to their co-expression similarities. In the present study, WGCNA was performed in order to identify functional genes associated with PTC using R package. First, differential gene expression analysis was conducted in order to identify the differentially expressed genes (DEGs) between PTC and normal samples. Subsequently, co-expression networks of the DEGs were constructed for the two sample groups, respectively. The two networks were compared in order to identify a poorly preserved module. Concentrating on the significant module, validation analysis was performed to confirm the identified genes and combined functional enrichment analysis was conducted in order to identify more functional associations of these genes with PTC. As a result, 1062 DEGs were identified for network construction. A brown module containing 118 highly related genes was selected as it exhibited the lowest module preservation. After validation analysis, 61 genes in the module were confirmed to be associated with PTC. Following the enrichment analysis, two PTC-related pathways were identified: Wnt signal pathway and transcriptional misregulation in cancer. LRP4, KLK7, PRICKLE1, ETV4 and ETV5 were predicted to be candidate genes regulating the pathogenesis of PTC. These results provide novel insights into the etiology of PTC and the identification of potential functional genes.

SUBMITTER: Ao ZX

PROVIDER: S-EPMC6144229 | biostudies-literature | 2018 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Identification of potential functional genes in papillary thyroid cancer by co-expression network analysis.

Ao Zeng-Xin ZX Chen Yuan-Cheng YC Lu Jun-Min JM Shen Jie J Peng Lin-Ping LP Lin Xu X Peng Cheng C Zeng Chun-Ping CP Wang Xia-Fang XF Zhou Rou R Chen Zhi Z Xiao Hong-Mei HM Deng Hong-Wen HW

Oncology letters 20180814 4

Interactions between multiple genes are involved in the development of complex diseases. However, there are few analyses of gene interactions associated with papillary thyroid cancer (PTC). Weighted gene co-expression network analysis (WGCNA) is a novel and powerful method that detects gene interactions according to their co-expression similarities. In the present study, WGCNA was performed in order to identify functional genes associated with PTC using R package. First, differential gene expres ...[more]

PMID: 30250553

Similar Datasets

Project description:Background The genes and genetic mechanisms underlying the occurrence and progression of papillary thyroid carcinoma (PTC) are still unknown. This study aimed to find candidate genes related to the pathogenesis and progression of PTC. Methods RNA sequencing (RNA-seq) data of PTC and normal tissues were downloaded from The Cancer Genome Atlas (TCGA) database with clinical stage data to form a test, validation, and clinical-stage data matrix. We used the test data set to analyze differentially expressed genes (DEGs) and weighted gene co-expression network analysis (WGCNA) to find those gene clusters highly correlated with PTC. We then verified the expression of genes in the interested modules using the validation matrix. The quantitative real-time polymerase chain reaction (qRT-PCR) was used to verify the reliability of the expression of selected genes. Five key genes (GDF15, LCN2, KCNN4, SH3BGRL3, and MMP2) were used to analyze the connection between gene expression and the American Joint Committee on Cancer (AJCC) stage. The upregulated and downregulated DEGs, along with the modules of interest, were subjected to Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). Results We used WGCNA to find two modules of interest, the yellow module, which was positively associated with PTC, and the blue module, which was negatively correlated with PTC. Four genes (GDF15, LCN2, KCNN4, and SH3BGRL3) from the yellow module were determined to be highly expressed in PTC in the test data matrix and were verified in both the validation data matrix and quantitative real-time PCR, which indicated that these four genes were highly correlated with the occurrence of the PTC. Furthermore, these four genes also had a significantly higher expression in the advanced levels of pathological T, N, and AJCC stage, meaning that they were correlated with the progression of PTC. Genes in the yellow module and upregulated DEGs were significantly enriched in three vital signaling pathways, including focal adhesion, extracellular matrix (ECM)-receptor interaction, and the PI3K-Akt signaling pathway. Conclusions Four candidate genes (GDF15, LCN2, KCNN4, and SH3BGRL3) may be potential biomarkers for the PTC’s pathogenesis and may be useful for predicting the disease stage.

Project description:BackgroundThyroid cancer is the most common endocrine malignancy, with a recent global increase of 20% in age-related incidence. Ultrasonography and ultrasonography-guided fine-needle aspiration biopsy (FNAB) are the most widely used diagnostic tests for thyroid nodules; however, it is estimated that up to 25% of thyroid biopsies are cytologically inconclusive. Molecular markers can help guide patient-oriented and targeted treatment of thyroid nodules and thyroid cancer.MethodsDatasets related to papillary thyroid cancer (PTC) or thyroid carcinoma (GSE129562, GSE3678, GSE54958, GSE138042, and GSE124653) were downloaded from the GEO database and analysed using the Limma package of R software. For functional enrichment analysis, the Kyoto Encyclopedia of Genes and Genomes pathway analysis and Gene Ontology were applied to differentially expressed genes (DEGs) using the Metascape website. A protein-protein interaction (PPI) network was built from the STRING database. Gene expression, protein expression, immunohistochemistry, and potential functional gene survival were analysed using the GEPIA website, the Human Protein Atlas website, and the UALCAN website. Potential target miRNAs were predicted using the miRDB and Starbase datasets.ResultsWe found 219 upregulated and 310 downregulated DEGs, with a cut-off of p < 0.01 and ∣log FC | >1.5. The DEGs in papillary thyroid cancer were mainly enriched in extracellular structural organisation. At the intersection of the PPI network and Metascape MCODEs, the hub genes in common were identified as FN1, APOE, CLU, and SDC2. In the targeted regulation network of miRNA-mRNA, the hsa-miR-424-5p was found to synchronously modulate two hub genes. Survival analysis showed that patients with high expression of CLU and APOE had better prognosis.ConclusionsCLU and APOE are involved in the molecular mechanism of papillary thyroid cancer. The hsa-miR-424-5p might have the potential to reverse the processes of papillary thyroid cancer by modulating the hub genes. These are potential targets for the treatment of patients with papillary thyroid cancer.

Project description:Background: Differentiated thyroid cancer (DTC) is the most common type of thyroid cancer. Many of them can relapse to dedifferentiated thyroid cancer (DDTC) and exhibit different gene expression profiles. The underlying mechanism of dedifferentiation and the involved genes or pathways remained to be investigated. Methods: A discovery cohort obtained from patients who received surgical resection in the Fudan University Shanghai Cancer Center (FUSCC) and two validation cohorts derived from Gene Expression Omnibus (GEO) database were used to screen out differentially expressed genes in the dedifferentiation process. Weighted gene co-expression network analysis (WGCNA) was constructed to identify modules highly related to differentiation. Gene Set Enrichment Analysis (GSEA) was used to identify pathways related to differentiation, and all differentially expressed genes were grouped by function based on the GSEA and literature reviewing data. Least absolute shrinkage and selection operator (LASSO) regression analysis was used to control the number of variables in each group. Next, we used logistic regression to build a gene signature in each group to indicate differentiation status, and we computed receiver operating characteristic (ROC) curve to evaluate the indicative performance of each signature. Results: A total of 307 upregulated and 313 downregulated genes in poorly differentiated thyroid cancer (PDTC) compared with papillary thyroid cancer (PTC) and normal thyroid (NT) were screened out in FUSCC cohort and validated in two GEO cohorts. WGCNA of 620 differential genes yielded the seven core genes with the highest correlation with thyroid differentiation score (TDS). Furthermore, 395 genes significantly correlated with TDS in univariate logistic regression analysis were divided into 11 groups. The areas under the ROC curve (AUCs) of the gene signature of group transcription and epigenetic modification, signal and substance transport, extracellular matrix (ECM), and metabolism in the training set [The Cancer Genome Atlas (TCGA) cohort] and validation set (combined GEO cohort) were both >0.75. The gene signature based on group transcription and epigenetic modification, cilia formation and movement, and proliferation can reflect the patient's disease recurrence state. Conclusion: The dedifferentiation of DTC is affected by a variety of mechanisms including many genes. The gene signature of group transcription and epigenetic modification, signal and substance transport, ECM, and metabolism can be used as biomarkers for DDTC.

Project description:BackgroundPapillary 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.MethodsThe 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.ResultsWe 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.ConclusionsOur 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.

Dataset Information

Identification of potential functional genes in papillary thyroid cancer by co-expression network analysis.

Publications

Identification of potential functional genes in papillary thyroid cancer by co-expression network analysis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets