Dataset Information

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 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.

SUBMITTER: Shen X

PROVIDER: S-EPMC8320213 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The morphologic assessment of uterine leiomyosarcoma (LMS) may be challenging, and diagnostic immunohistochemical (IHC) analysis is currently lacking. We evaluated the genomic landscape of 167 uterine LMS by targeted next-generation sequencing (NGS) to identify common genomic alterations. IHC analyses corresponding to these genomic landmarks were applied to a test cohort of 16 uterine LMS, 6 smooth muscle tumors of uncertain malignant potential (STUMP), and 6 leiomyomas with NGS data and a validation cohort of 8 uterine LMS, 12 STUMP, 21 leiomyomas and leiomyoma variants, 7 low-grade endometrial stromal sarcomas, and 2 diagnostically challenging uterine smooth muscle tumors. IHC results were individually interpreted by 3 pathologists blinded to NGS data. Overall, 94% of LMS showed ≥1 genomic alteration involving TP53, RB1, ATRX, PTEN, CDKN2A, or MDM2, with 80% showing alterations in ≥2 of these genes. In the test cohort, an initial panel of p53, Rb, PTEN, and ATRX was applied, followed by a panel of DAXX, MTAP, and MDM2 in cases without abnormalities. Abnormal p53, Rb, PTEN, and ATRX IHC expression was seen in 75%, 88%, 44%, and 38% of LMS, respectively, in the test cohort. Two or more abnormal IHC results among these markers were seen in 81% of LMS. STUMPs demonstrated only 1 IHC abnormality involving these markers. No IHC abnormalities were seen in leiomyomas. In the validation cohort, abnormal p53, Rb, and PTEN IHC results were seen in LMS, whereas rare STUMP or leiomyomas with bizarre nuclei showed IHC abnormalities involving only 1 of the markers. Abnormalities in ≥2 markers were present in both diagnostically challenging smooth muscle tumors, confirming LMS. Concordance was excellent among pathologists in the interpretation of IHC (κ = 0.97) and between IHC and NGS results (κ = 0.941). Uterine LMS exhibit genomic landmark alterations for which IHC surrogates exist, and a diagnostic algorithm involving molecular-based IHC may aid in the evaluation of unusual uterine smooth muscle tumors.

Project description:ObjectiveTo investigate the differential expression gene modules and hub genes associated with Alzheimer's disease (AD) by weighted gene co-expression network analysis (WGCNA) and annotate the biological functions of these modules.MethodsWe downloaded transcriptome sequencing data from the GEO database, and according to the correlation of the genes, a gene co-expression network was constructed with the parameter setting of β=8 and a correlation coefficient threshold of 0.85. Pearson correlation test was used to calculate the correlation between the module genes and clinical traits to screen the gene modules significantly associated with AD and identify the hub genes according to the connectivity within modules. GO functional enrichment analysis and KEGG pathway analysis were used to annotate the functions of the modules. A cell model of AD was established in SH-SY5Y cells by Aβ1-42 treatment, and the mRNA expression levels of the hub genes were compared between the Aβ1-42-treated cells and the control cells.ResultsTen gene co-expression modules were constructed based on the correlations of gene expression, in which the brown (r=0.66, P < 0.001) and turquoise modules (r=-0.68, P < 0.001) were significantly correlated with the AD group. Forty-eight genes were identified as the hub genes in the co-expression network. Function annotation revealed that the genes in both modules were mainly enriched in DNA damage and repair pathways and metabolism-related pathways. Differential expression analysis of the genes revealed that the genes DNASE1, TEKT2 and MTSS1L were highly expressed while ACP2, LANCL2 and GMPR2 were lowly expressed in AD group. The results of cell experiment confirmed the up-regulation of DNASE1, TEKT2 and MTSS1L genes and the down-regulation of ACP2, LANCL2, and GMPR2 in Aβ1-42-treated SH-SY5Y cells (P < 0.01).ConclusionThe brown and turquoise modules are closely correlated with AD. The hub genes including MTSS1L, GMPR2, ACP2, ACTG1 and LANCL2 selected from the modules may participate in AD pathogenesis by regulating DNA damage and repair.