Identification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and CIBERSORT algorithm
Abstract Background While large-scale genomic analyses symbolize a precious attempt to decipher the molecular foundation of uterine leiomyosarcoma (ULMS), bioinformatics results associated with the occurrence of ULMS based totally on WGCNA and CIBERSORT have not yet been reported. This study aimed t...
Main Authors: | , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
BMC
2021-07-01
|
Series: | World Journal of Surgical Oncology |
Subjects: | |
Online Access: | https://doi.org/10.1186/s12957-021-02333-z |
id |
doaj-11db03835387436a915e5620276a5487 |
---|---|
record_format |
Article |
spelling |
doaj-11db03835387436a915e5620276a54872021-08-01T11:36:31ZengBMCWorld Journal of Surgical Oncology1477-78192021-07-0119111210.1186/s12957-021-02333-zIdentification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and CIBERSORT algorithmXiaoqing Shen0Zhujuan Yang1Songwei Feng2Yi Li3Department of Gynecology and Obstetrics, The Second Affiliated Hospital of Soochow UniversityDepartment of Gynecology and Obstetrics, The Second Affiliated Hospital of Soochow UniversityDepartment of Gynecology and Obstetrics, The Second Affiliated Hospital of Soochow UniversityDepartment of Gynecology, The Affiliated Jiangsu Shengze Hospital of Nanjing Medical University and Jiangsu Shengze HospitalAbstract Background While large-scale genomic analyses symbolize a precious attempt to decipher the molecular foundation of uterine leiomyosarcoma (ULMS), bioinformatics results associated with the occurrence of ULMS based totally on WGCNA and CIBERSORT have not yet been reported. This study aimed to screen the hub genes and the immune cell infiltration pattern in ULMS by bioinformatics methods. Methods Firstly, the GSE67463 dataset, including 25 ULMS tissues and 29 normal myometrium (NL) tissues, was downloaded from the public database. The differentially expressed genes (DEGs) were screened by the ‘limma’ package and hub modules were identified by weighted gene co-expression network analysis (WGCNA). Subsequently, gene function annotations were performed to investigate the biological role of the genes from the intersection of two groups (hub module and DEGs). The above genes were calculated in the protein–protein interaction (PPI) network to select the hub genes further. The hub genes were validated using external data (GSE764 and GSE68295). In addition, the differential immune cell infiltration between UL and ULMS tissues was investigated using the CIBERSORT algorithm. Finally, we used western blot to preliminarily detect the hub genes in cell lines. Results WGCNA analysis revealed a green-yellow module possessed the highest correlation with ULMS, including 1063 genes. A total of 172 DEGs were selected by thresholds set in the ‘limma’ package. The above two groups of genes were intersected to obtain 72 genes for functional annotation analysis. Interestingly, it indicated that 72 genes were mainly involved in immune processes and the Neddylation pathway. We found a higher infiltration of five types of cells (memory B cells, M0-type macrophages, mast cells activated, M1-type macrophages, and T cells follicular helper) in ULMS tissues than NL tissues, while the infiltration of two types of cells (NK cells activated and mast cells resting) was lower than in NL tissues. In addition, a total of five genes (KDR, CCL21, SELP, DPT, and DCN) were identified as the hub genes. Internal and external validation demonstrated that the five genes were over-expressed in NL tissues compared with USML tissues. Finally, the correlation analysis results indicate that NK cells activated and mast cells activated positively correlated with the hub genes. However, M1-type macrophages had a negative correlation with the hub genes. Moreover, only the DCN may be associated with the Neddylation pathway. Conclusion A series of evidence confirm that the five hub genes and the infiltration of seven types of immune cells are related to USML occurrence. These hub genes may affect the occurrence of USML through immune-related and Neddylation pathways, providing molecular evidence for the treatment of USML in the future.https://doi.org/10.1186/s12957-021-02333-zUterine leiomyosarcomaCIBERSORTWeighted co-expression networkNeddylation |
collection |
DOAJ |
language |
English |
format |
Article |
sources |
DOAJ |
author |
Xiaoqing Shen Zhujuan Yang Songwei Feng Yi Li |
spellingShingle |
Xiaoqing Shen Zhujuan Yang Songwei Feng Yi Li Identification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and CIBERSORT algorithm World Journal of Surgical Oncology Uterine leiomyosarcoma CIBERSORT Weighted co-expression network Neddylation |
author_facet |
Xiaoqing Shen Zhujuan Yang Songwei Feng Yi Li |
author_sort |
Xiaoqing Shen |
title |
Identification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and CIBERSORT algorithm |
title_short |
Identification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and CIBERSORT algorithm |
title_full |
Identification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and CIBERSORT algorithm |
title_fullStr |
Identification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and CIBERSORT algorithm |
title_full_unstemmed |
Identification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and CIBERSORT algorithm |
title_sort |
identification of uterine leiomyosarcoma-associated hub genes and immune cell infiltration pattern using weighted co-expression network analysis and cibersort algorithm |
publisher |
BMC |
series |
World Journal of Surgical Oncology |
issn |
1477-7819 |
publishDate |
2021-07-01 |
description |
Abstract Background While large-scale genomic analyses symbolize a precious attempt to decipher the molecular foundation of uterine leiomyosarcoma (ULMS), bioinformatics results associated with the occurrence of ULMS based totally on WGCNA and CIBERSORT have not yet been reported. This study aimed to screen the hub genes and the immune cell infiltration pattern in ULMS by bioinformatics methods. Methods Firstly, the GSE67463 dataset, including 25 ULMS tissues and 29 normal myometrium (NL) tissues, was downloaded from the public database. The differentially expressed genes (DEGs) were screened by the ‘limma’ package and hub modules were identified by weighted gene co-expression network analysis (WGCNA). Subsequently, gene function annotations were performed to investigate the biological role of the genes from the intersection of two groups (hub module and DEGs). The above genes were calculated in the protein–protein interaction (PPI) network to select the hub genes further. The hub genes were validated using external data (GSE764 and GSE68295). In addition, the differential immune cell infiltration between UL and ULMS tissues was investigated using the CIBERSORT algorithm. Finally, we used western blot to preliminarily detect the hub genes in cell lines. Results WGCNA analysis revealed a green-yellow module possessed the highest correlation with ULMS, including 1063 genes. A total of 172 DEGs were selected by thresholds set in the ‘limma’ package. The above two groups of genes were intersected to obtain 72 genes for functional annotation analysis. Interestingly, it indicated that 72 genes were mainly involved in immune processes and the Neddylation pathway. We found a higher infiltration of five types of cells (memory B cells, M0-type macrophages, mast cells activated, M1-type macrophages, and T cells follicular helper) in ULMS tissues than NL tissues, while the infiltration of two types of cells (NK cells activated and mast cells resting) was lower than in NL tissues. In addition, a total of five genes (KDR, CCL21, SELP, DPT, and DCN) were identified as the hub genes. Internal and external validation demonstrated that the five genes were over-expressed in NL tissues compared with USML tissues. Finally, the correlation analysis results indicate that NK cells activated and mast cells activated positively correlated with the hub genes. However, M1-type macrophages had a negative correlation with the hub genes. Moreover, only the DCN may be associated with the Neddylation pathway. Conclusion A series of evidence confirm that the five hub genes and the infiltration of seven types of immune cells are related to USML occurrence. These hub genes may affect the occurrence of USML through immune-related and Neddylation pathways, providing molecular evidence for the treatment of USML in the future. |
topic |
Uterine leiomyosarcoma CIBERSORT Weighted co-expression network Neddylation |
url |
https://doi.org/10.1186/s12957-021-02333-z |
work_keys_str_mv |
AT xiaoqingshen identificationofuterineleiomyosarcomaassociatedhubgenesandimmunecellinfiltrationpatternusingweightedcoexpressionnetworkanalysisandcibersortalgorithm AT zhujuanyang identificationofuterineleiomyosarcomaassociatedhubgenesandimmunecellinfiltrationpatternusingweightedcoexpressionnetworkanalysisandcibersortalgorithm AT songweifeng identificationofuterineleiomyosarcomaassociatedhubgenesandimmunecellinfiltrationpatternusingweightedcoexpressionnetworkanalysisandcibersortalgorithm AT yili identificationofuterineleiomyosarcomaassociatedhubgenesandimmunecellinfiltrationpatternusingweightedcoexpressionnetworkanalysisandcibersortalgorithm |
_version_ |
1721245746281316352 |