Identification of aberrantly methylated differentially expressed genes in melanoma

Enrichment Analysis ◽

Regulation Of Transcription ◽

Biological Processes ◽

Hub Genes ◽

Positive Regulation ◽

Upregulated Genes

Abstract Background: Skin Cutaneous Melanoma (SKCM) is known as an aggressive malignant cancer, which could be directly derived from melanocytic nevi. However, the molecular mechanisms underlying malignant transformation of melanocytes and melanoma tumor progression still remain unclear. Increasing researches showed significant roles of epigenetic modifications, especially DNA methylation, in melanoma. This study focused on identification and analysis of methylation-regulated differentially expressed genes (MeDEGs) between melanocytic nevus and malignant melanoma in genome-wide profiles. Methods: The gene expression profiling datasets (GSE3189 and GSE114445) and gene methylation profiling datasets (GSE86355 and GSE120878) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) and differentially methylated genes (DMGs) were identified via GEO2R. MeDEGs were obtained by integrating the DEGs and DMGs. Then, functional enrichment analysis of MeDEGs were performed. STRING and Cytoscape were used to describe protein-protein interaction(PPI) network. Furthermore, survival analysis was implemented to select the prognostic hub genes. Finally, we conducted gene set enrichment analysis (GSEA) of hub genes. Results: We identified 237 hypomethylated, upregulated genes and 182 hypermethylated, downregulated genes. Hypomethylation-upregulated genes were enriched in biological processes of the oxidation-reduction process, cell proliferation, cell division, phosphorylation, extracellular matrix disassembly and protein sumoylation. Pathway enrichment showed selenocompound metabolism, small cell lung cancer and lysosome. Hypermethylation-downregulated genes were enriched in biological processes of positive regulation of transcription from RNA polymerase II promoter, cell adhesion, cell proliferation, positive regulation of transcription, DNA-templated and angiogenesis. The most significantly enriched pathways involved the transcriptional misregulation in cancer, circadian rhythm, tight junction, protein digestion and absorption and Hippo signaling pathway. After PPI establishment and survival analysis, seven prognostic hub genes were CKS2, DTL, KIF2C, KPNA2, MYBL2, TPX2 and FBL. Moreover, the most involved hallmarks obtained by GSEA were E2F targets, G2M checkpoint and mitotic spindle. Conclusions: Our study identified potential aberrantly methylated-differentially expressed genes participating in the process of malignant transformation from nevus to melanoma tissues based on comprehensive genomic profiles. Transcription profiles of CKS2, DTL, KIF2C, KPNA2, MYBL2, TPX2 and FBL provided clues of aberrantly methylation-based biomarkers, which might improve the development of precise medicine.

Bioinformatics-Based Study to Investigate Potential Differentially Expressed Genes and miRNAs in Hashimoto’s Thyroiditis

10.21203/rs.3.rs-1023446/v1 ◽

2021 ◽

Author(s):

Li Guoquan ◽

Du Junwei ◽

He Qi ◽

Fu Xinghao ◽

Ji Feihong ◽

...

Keyword(s):

Gene Expression ◽

Inflammatory Responses ◽

Expression Profiles ◽

Treatment Options ◽

Enrichment Analysis ◽

Pathway Enrichment Analysis ◽

Chronic Lymphocytic Thyroiditis ◽

Hub Genes

Abstract BackgroundHashimoto's thyroiditis (HT), also known as chronic lymphocytic thyroiditis, is a common autoimmune disease, which mainly occurs in women. The early manifestation was hyperthyroidism, however, hypothyroidism may occur if HT was not controlled for a long time. Numerous studies have shown that multiple factors, including genetic, environmental, and autoimmune factors, were involved in the pathogenesis of the disease, but the exact mechanisms were not yet clear. The aim of this study was to identify differentially expressed genes (DEGs) by comprehensive analysis and to provide specific insights into HT. MethodsTwo gene expression profiles (GSE6339, GSE138198) about HT were downloaded from the Gene Expression Omnibus (GEO) database. The DEGs were assessed between the HT and normal groups using the GEO2R. The DEGs were then sent to the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. The hub genes were discovered using Cytoscape and CytoHubba. Finally, NetworkAnalyst was utilized to create the hub genes' targeted microRNAs (miRNAs). ResultsA total of 62 DEGs were discovered, including 60 up-regulated and 2 down-regulated DEGs. The signaling pathways were mainly engaged in cytokine interaction and cytotoxicity, and the DEGs were mostly enriched in immunological and inflammatory responses. IL2RA, CXCL9, IL10RA, CCL3, CCL4, CCL2, STAT1, CD4, CSF1R, and ITGAX were chosen as hub genes based on the results of the protein-protein interaction (PPI) network and CytoHubba. Five miRNAs, including mir-24-3p, mir-223-3p, mir-155-5p, mir-34a-5p, mir-26b-5p, and mir-6499-3p, were suggested as likely important miRNAs in HT. ConclusionsThese hub genes, pathways and miRNAs contribute to a better understanding of the pathophysiology of HT and offer potential treatment options for HT.

Investigation of Aberrantly Methylated Differentially Expressed Genes in Pancreatic Cancer by Integrated Bioinformatics Analysis

10.21203/rs.3.rs-132250/v1 ◽

2021 ◽

Author(s):

Cailin xue ◽

Peng gao ◽

Xudong zhang ◽

Xiaohan cui ◽

Lei jin ◽

...

Keyword(s):

Gene Expression ◽

Pancreatic Cancer ◽

Signaling Pathways ◽

Dna Sequences ◽

Receptor Interaction ◽

Biological Processes ◽

Ppi Network ◽

Hub Genes

Abstract Background: Abnormal methylation of DNA sequences plays an important role in the development and progression of pancreatic cancer (PC). The purpose of this study was to identify abnormal methylation genes and related signaling pathways in PC by comprehensive bioinformatic analysis of three datasets in the Gene Expression Omnibus (GEO). Methods: Datasets of gene expression microarrays (GSE91035, GSE15471) and gene methylation microarrays (GSE37480) were downloaded from the GEO database. Aberrantly methylated-differentially expressed genes (DEGs) were analysis by GEO2R software. GO and KEGG enrichment analyses of selected genes were performed using DAVID database. A protein–protein interaction (PPI) network was constructed by STRING and visualized in Cytoscape. Core module analysis was performed by Mcode in Cytoscape. Hub genes were obtained by CytoHubba app. in Cytoscape software. Results: A total of 267 hypomethylation-high expression genes, which were enriched in biological processes of cell adhesion, biological adhesion and regulation of signaling were obtained. KEGG pathway enrichment showed ECM-receptor interaction, Focal adhesion and PI3K-Akt signaling pathway. The top 5 hub genes of PPI network were EZH2, CCNA2, CDC20, KIF11, UBE2C. As for hypermethylation-low expression genes, 202 genes were identified, which were enriched in biological processes of cellular amino acid biosynthesis process and positive regulation of PI3K activity, etc. The pathways enriched were the pancreatic secretion and biosynthesis of amino acids pathways, etc. The five significant hub genes were DLG3, GPT2, PLCB1, CXCL12 and GNG7. In addition, five genes, including CCNA2, KIF11, UBE2C, PLCB1 and GNG7, significantly associated with patient's prognosis were also identified. Conclusion: Novel genes with abnormal expression were identified, which will help us further understand the molecular mechanism and related signaling pathways of PC, and these aberrant genes could possibly serve as biomarkers for precise diagnosis and treatment of PC.

Bioinformatics Analysis and Identification of Potential Biomarkers in Gastric Cancer

10.21203/rs.3.rs-111288/v1 ◽

2020 ◽

Author(s):

Jingdi Yang ◽

Bo Peng ◽

Xianzheng Qin ◽

Tian Zhou

Keyword(s):

Gene Expression ◽

Gastric Cancer ◽

Bioinformatics Analysis ◽

Interaction Network ◽

Enrichment Analysis ◽

Protein Protein Interaction Network

Hub Genes ◽

Protein Protein Interaction ◽

Abstract Background: Although the morbidity and mortality of gastric cancer are declining, gastric cancer is still one of the most common causes of death. Early detection of gastric cancer is of great help to improve the survival rate, but the existing biomarkers are not sensitive to diagnose early gastric cancer. The aim of this study is to identify the novel biomarkers for gastric cancer.Methods: Three gene expression profiles (GSE27342, GSE63089, GSE33335) were downloaded from Gene Expression Omnibus database to select differentially expressed genes. Then, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis were performed to explore the biological functions of differentially expressed genes. Cytoscape was utilized to construct protein-protein interaction network and hub genes were analyzed by plugin cytoHubba of Cytoscape. Furthermore, Gene Expression Profiling Interactive Analysis and Kaplan-Meier plotter were used to verify the identified hub genes.Results: 35 overlapping differentially expressed genes were screened from gene expression datasets, which consisted of 11 up-regulated genes and 24 down-regulated genes. Gene Ontology functional enrichment analysis revealed that differentially expressed genes were significantly enriched in digestion, regulation of biological quality, response to hormone and steroid hormone, and homeostatic process. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed differentially expressed genes were enriched in the secretion of gastric acid and collecting duct acid, leukocyte transendothelial migration and ECM-receptor interaction. According to protein-protein interaction network, 10 hub genes were identified by Maximal Clique Centrality method.Conclusion: By using bioinformatics analysis, COL1A1, BGN, THY1, TFF2 and SST were identified as the potential biomarkers for early detection of gastric cancer.

Screening and identification of MSX1 for predicting progestin resistance in endometrial cancer

10.21203/rs.3.rs-17361/v1 ◽

2020 ◽

Author(s):

Linlin Yang ◽

Yunxia Cui ◽

Ting Huang ◽

Xiao Sun ◽

Yudong Wang

Keyword(s):

Gene Expression ◽

Overall Survival ◽

Endometrial Cancer ◽

Endometrial Carcinoma ◽

Enrichment Analysis ◽

Functional Enrichment ◽

Pathway Enrichment Analysis ◽

Hub Genes

Abstract Background: Progestin resistance is a critical obstacle for endometrial conservative therapy. Therefore, the studies to acquire a more comprehensive understanding of the mechanisms and specific biomarkers to predict progestin resistance are very important. However, the pivotal roles of essential molecules of progestin resistance are still unexplored. Methods: We downloaded GSE121367 with gene expression profiles of medroxyprogesterone acetate (MPA) resistant and sensitive cell lines from the GEO database. The “limma” R language package was applied to identify differentially expressed genes (DEGs). Gene ontology and pathway enrichment analysis was performed through the database of DAVID. Meanwhile, we conducted GSEA analysis to identify pathway enrichments. Protein–protein interaction construction of top genes was conducted to screen hub genes by STRING and visualized in Cytoscape. A high connectivity degree of hub genes were picked out to perform the differential expression, methylation validation and overall survival analysis in the Gene Expression Profiling Interactive Analysis database, Human Protein Atlas database and Kaplan–Meier plotter online tool, respectively. In addition, microRNAs and upstream transcription factors of hub genes were predicted by miRTarBase and Network Analyst database. Results: A total number of 3282 differentially expressed genes were identified. Functional enrichment analysis demonstrated that these genes were mostly enriched in negative regulation of DNA binding, chronic inflammatory response and cell adhesion molecules pathway. We screened out ten hub genes including CDH1, JAG1, PTGES, EPCAM, CNTNAP2, TBX1, MSX1, KRT19, OAS1 and DAB2 among different groups. The genomic alteration rates of hub genes were low based on the current uterine corpus endometrial carcinoma sample sets. Their relevant microRNA and transcription factor were detected and has-miR-335-5p, has-miR-124-3p, MAZ and TFDP1 were the most prominent. The methylation status of CDH1, JAG1, EPCAM and MSX1 were decreased, corresponding to their high protein expression in endometrial cancers, which also indicated better overall survival. The homeobox protein of MSX1 showed significantly tissue specificity. Conclusions: Our study identified ten hub genes associated with progestin resistance of endometrial cancer and screened out the gene of MSX1 which promised to be the specific indicator. This would shed new light on the underlying biological marker to overcome the progestin resistance of endometrial cancer. Keywords : Bioinformatic analysis, Progestin resistance, Endometrial carcinoma, MSX1

Study on potential differentially expressed genes in stroke by bioinformatics analysis

Neurological Sciences ◽

10.1007/s10072-021-05470-1 ◽

2021 ◽

Author(s):

Xitong Yang ◽

Pengyu Wang ◽

Shanquan Yan ◽

Guangming Wang

Keyword(s):

Gene Expression ◽

Normal Control ◽

Enrichment Analysis ◽

Normal Control Group ◽

Control Group ◽

Ppi Network ◽

Hub Genes ◽

Pathway Enrichment ◽

Key Genes

AbstractStroke is a sudden cerebrovascular circulatory disorder with high morbidity, disability, mortality, and recurrence rate, but its pathogenesis and key genes are still unclear. In this study, bioinformatics was used to deeply analyze the pathogenesis of stroke and related key genes, so as to study the potential pathogenesis of stroke and provide guidance for clinical treatment. Gene Expression profiles of GSE58294 and GSE16561 were obtained from Gene Expression Omnibus (GEO), the differentially expressed genes (DEGs) were identified between IS and normal control group. The different expression genes (DEGs) between IS and normal control group were screened with the GEO2R online tool. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the DEGs were performed. Using the Database for Annotation, Visualization and Integrated Discovery (DAVID) and gene set enrichment analysis (GSEA), the function and pathway enrichment analysis of DEGS were performed. Then, a protein–protein interaction (PPI) network was constructed via the Search Tool for the Retrieval of Interacting Genes (STRING) database. Cytoscape with CytoHubba were used to identify the hub genes. Finally, NetworkAnalyst was used to construct the targeted microRNAs (miRNAs) of the hub genes. A total of 85 DEGs were screened out in this study, including 65 upward genes and 20 downward genes. In addition, 3 KEGG pathways, cytokine − cytokine receptor interaction, hematopoietic cell lineage, B cell receptor signaling pathway, were significantly enriched using a database for labeling, visualization, and synthetic discovery. In combination with the results of the PPI network and CytoHubba, 10 hub genes including CEACAM8, CD19, MMP9, ARG1, CKAP4, CCR7, MGAM, CD79A, CD79B, and CLEC4D were selected. Combined with DEG-miRNAs visualization, 5 miRNAs, including hsa-mir-146a-5p, hsa-mir-7-5p, hsa-mir-335-5p, and hsa-mir-27a- 3p, were predicted as possibly the key miRNAs. Our findings will contribute to identification of potential biomarkers and novel strategies for the treatment of ischemic stroke, and provide a new strategy for clinical therapy.

AB0005 IDENTIFICATION OF KEY GENES AND PATHWAYS FOR PSORIASIS BASED ON GEO DATABASES BY BIOINFORMATICS ANALYSIS

Annals of the Rheumatic Diseases ◽

10.1136/annrheumdis-2021-eular.1773 ◽

2021 ◽

Vol 80 (Suppl 1) ◽

pp. 1037.2-1038

Author(s):

X. Sun ◽

S. X. Zhang ◽

S. Song ◽

T. Kong ◽

C. Zheng ◽

...

Keyword(s):

Gene Expression ◽

Signaling Pathways ◽

Enrichment Analysis ◽

Shanxi Province ◽

Receptor Interaction ◽

Term Therapy ◽

Pathway Enrichment Analysis ◽

Key Genes

Background:Psoriasis is an immune-mediated, genetic disease manifesting in the skin or joints or both, and also has a strong genetic predisposition and autoimmune pathogenic traits1. The hallmark of psoriasis is sustained inflammation that leads to uncontrolled keratinocyte proliferation and dysfunctional differentiation. And it’s also a chronic relapsing disease, which often necessitates a long-term therapy2.Objectives:To investigate the molecular mechanisms of psoriasis and find the potential gene targets for diagnosis and treating psoriasis.Methods:Total 334 gene expression data of patients with psoriasis research (GSE13355 GSE14905 and GSE30999) were obtained from the Gene Expression Omnibus database. After data preprocessing and screening of differentially expressed genes (DEGs) by R software. Online toll Metascape3 was used to analyze Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs. Interactions of proteins encoded by DEGs were discovered by Protein-protein interaction network (PPI) using STRING online software. Cytoscape software was utilized to visualize PPI and the degree of each DEGs was obtained by analyzing the topological structure of the PPI network.Results:A total of 611 DEGs were found to be differentially expressed in psoriasis. GO analysis revealed that up-regulated DEGs were mostly associated with defense and response to external stimulus while down-regulated DEGs were mostly associated with metabolism and synthesis of lipids. KEGG enrichment analysis suggested they were mainly enriched in IL-17 signaling, Toll-like receptor signaling and PPAR signaling pathways, Cytokine-cytokine receptor interaction and lipid metabolism. In addition, top 9 key genes (CXCL10, OASL, IFIT1, IFIT3, RSAD2, MX1, OAS1, IFI44 and OAS2) were identified through Cytoscape.Conclusion:DEGs of psoriasis may play an essential role in disease development and may be potential pathogeneses of psoriasis.References:[1]Boehncke WH, Schon MP. Psoriasis. Lancet 2015;386(9997):983-94. doi: 10.1016/S0140-6736(14)61909-7 [published Online First: 2015/05/31].[2]Zhang YJ, Sun YZ, Gao XH, et al. Integrated bioinformatic analysis of differentially expressed genes and signaling pathways in plaque psoriasis. Mol Med Rep 2019;20(1):225-35. doi: 10.3892/mmr.2019.10241 [published Online First: 2019/05/23].[3]Zhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun 2019;10(1):1523. doi: 10.1038/s41467-019-09234-6 [published Online First: 2019/04/05].Acknowledgements:This project was supported by National Science Foundation of China (82001740), Open Fund from the Key Laboratory of Cellular Physiology (Shanxi Medical University) (KLCP2019) and Innovation Plan for Postgraduate Education in Shanxi Province (2020BY078).Disclosure of Interests:None declared

Analysis of Transcriptional Changes in Different Brassica napus Synthetic Allopolyploids

Genes ◽

10.3390/genes12010082 ◽

2021 ◽

Vol 12 (1) ◽

pp. 82

Author(s):

Yunxiao Wei ◽

Guoliang Li ◽

Shujiang Zhang ◽

Shifan Zhang ◽

Hui Zhang ◽

...

Keyword(s):

Gene Expression ◽

Brassica Napus ◽

Expression Patterns ◽

Female Parent ◽

Enrichment Analysis ◽

Transcriptional Changes ◽

Aa Genome ◽

High Degree

Allopolyploidy is an evolutionary and mechanistically intriguing process involving the reconciliation of two or more sets of diverged genomes and regulatory interactions, resulting in new phenotypes. In this study, we explored the gene expression patterns of eight F2 synthetic Brassica napus using RNA sequencing. We found that B. napus allopolyploid formation was accompanied by extensive changes in gene expression. A comparison between F2 and the parent shows a certain proportion of differentially expressed genes (DEG) and activation\silent gene, and the two genomes (female parent (AA)\male parent (CC) genomes) showed significant differences in response to whole-genome duplication (WGD); non-additively expressed genes represented a small portion, while Gene Ontology (GO) enrichment analysis showed that it played an important role in responding to WGD. Besides, genome-wide expression level dominance (ELD) was biased toward the AA genome, and the parental expression pattern of most genes showed a high degree of conservation. Moreover, gene expression showed differences among eight individuals and was consistent with the results of a cluster analysis of traits. Furthermore, the differential expression of waxy synthetic pathways and flowering pathway genes could explain the performance of traits. Collectively, gene expression of the newly formed allopolyploid changed dramatically, and this was different among the selfing offspring, which could be a prominent cause of the trait separation. Our data provide novel insights into the relationship between the expression of differentially expressed genes and trait segregation and provide clues into the evolution of allopolyploids.

Bioinformatics analysis of differentially expressed genes and pathways in the development of cervical cancer

BMC Cancer ◽

10.1186/s12885-021-08412-4 ◽

2021 ◽

Vol 21 (1) ◽

Author(s):

Baojie Wu ◽

Shuyi Xi

Keyword(s):

Gene Expression ◽

Cervical Cancer ◽

Expression Profiles ◽

Gene Expression Profiles ◽

Molecular Targets ◽

Enrichment Analysis ◽

Ppi Network ◽

Hub Genes ◽

Kegg Pathways

Abstract Background This study aimed to explore and identify key genes and signaling pathways that contribute to the progression of cervical cancer to improve prognosis. Methods Three gene expression profiles (GSE63514, GSE64217 and GSE138080) were screened and downloaded from the Gene Expression Omnibus database (GEO). Differentially expressed genes (DEGs) were screened using the GEO2R and Venn diagram tools. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed. Gene set enrichment analysis (GSEA) was performed to analyze the three gene expression profiles. Moreover, a protein–protein interaction (PPI) network of the DEGs was constructed, and functional enrichment analysis was performed. On this basis, hub genes from critical PPI subnetworks were explored with Cytoscape software. The expression of these genes in tumors was verified, and survival analysis of potential prognostic genes from critical subnetworks was conducted. Functional annotation, multiple gene comparison and dimensionality reduction in candidate genes indicated the clinical significance of potential targets. Results A total of 476 DEGs were screened: 253 upregulated genes and 223 downregulated genes. DEGs were enriched in 22 biological processes, 16 cellular components and 9 molecular functions in precancerous lesions and cervical cancer. DEGs were mainly enriched in 10 KEGG pathways. Through intersection analysis and data mining, 3 key KEGG pathways and related core genes were revealed by GSEA. Moreover, a PPI network of 476 DEGs was constructed, hub genes from 12 critical subnetworks were explored, and a total of 14 potential molecular targets were obtained. Conclusions These findings promote the understanding of the molecular mechanism of and clinically related molecular targets for cervical cancer.

Diagnostic model of combined ceRNA and DNA methylation related genes in esophageal carcinoma

PeerJ ◽

10.7717/peerj.8831 ◽

2020 ◽

Vol 8 ◽

pp. e8831 ◽

Cited By ~ 1

Author(s):

Xiaojiao Guan ◽

Yao Yao ◽

Guangyao Bao ◽

Yue Wang ◽

Aimeng Zhang ◽

...

Keyword(s):

Gene Expression ◽

Esophageal Cancer ◽

Esophageal Carcinoma ◽

Enrichment Analysis ◽

Functional Enrichment ◽

Diagnostic Model ◽

Link Type ◽

Key Genes

Esophageal cancer is a common malignant tumor in the world, and the aim of this study was to screen key genes related to the development of esophageal cancer using a variety of bioinformatics analysis tools and analyze their biological functions. The data of esophageal squamous cell carcinoma from the Gene Expression Omnibus (GEO) were selected as the research object, processed and analyzed to screen differentially expressed microRNAs (miRNAs) and differential methylation genes. The competing endogenous RNAs (ceRNAs) interaction network of differentially expressed genes was constructed by bioinformatics tools DAVID, String, and Cytoscape. Biofunctional enrichment analysis was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). The expression of the screened genes and the survival of the patients were verified. By analyzing GSE59973 and GSE114110, we found three down-regulated and nine up-regulated miRNAs. The gene expression matrix of GSE120356 was calculated by Pearson correlation coefficient, and the 11696 pairs of ceRNA relation were determined. In the ceRNA network, 643 lncRNAs and 147 mRNAs showed methylation difference. Functional enrichment analysis showed that these differentially expressed genes were mainly concentrated in the FoxO signaling pathway and were involved in the corresponding cascade of calcineurin. By analyzing the clinical data in The Cancer Genome Atlas (TCGA) database, it was found that four lncRNAs had an important impact on the survival and prognosis of esophageal carcinoma patients. QRT-PCR was also conducted to identify the expression of the key lncRNAs (RNF217-AS1, HCP5, ZFPM2-AS1 and HCG22) in ESCC samples. The selected key genes can provide theoretical guidance for further research on the molecular mechanism of esophageal carcinoma and the screening of molecular markers.

Screening the Hub Genes and Analyzing the Mechanisms in Discharged COVID-19 Patients Retesting Positive through Bioinformatics Analysis

10.21203/rs.3.rs-991635/v1 ◽

2021 ◽

Author(s):

Ke-Ying Fang ◽

Gui-Ning Liang ◽

Zhuo-Qing Zhuang ◽

Yong-Xin Fang ◽

Yu-Qian Dong ◽

...

Keyword(s):

Social Order ◽

Expression System ◽

Enrichment Analysis ◽

Virus Reactivation ◽

Hub Genes ◽

Messenger Ribonucleic Acid ◽

Ppi Networks ◽

Significant Enrichment

Abstract Background: With the worldwide spread of COVID-19, people’s health and social order have been exposed to enormous risks. After encountering patients who test positive again after discharge, our study analyzed the pathogenesis to further assess the risk and possibility of virus reactivation.Methods: A separate microarray was acquired from the Integrated Gene Expression System (GEO), and its samples were divided into two groups: a “convalescent-RTP” group consisting of recovery and “retesting-positive” (RTP) patients (group CR) and a “health-RTP” group consisting of healthy control and RTP patients (group HR). The enrichment analysis was performed with R software, obtaining the gene ontology (GO) and Kyoto pluripotent stem cells (KEGG) of the genes and genomes. Subsequently, the protein–protein interaction (PPI) networks of each group were established and the hub genes were discovered using the cytoHubba plug-in.Results: In this study, 20 differentially expressed genes were identified, and 6622 genes were identified in the group CR, consisting of 5003 up-regulated and 1619 down-regulated genes. Meanwhile, 7335 genes were screened in the group HR, including 4323 up-regulated and 3012 down-regulated ones. The GO and KEGG analysis of the two groups revealed significant enrichment of these differentially expressed genes in pathways associated with immune response and apoptosis. In the PPI network constructed, 10 hub genes in group CR were identified, including TP53BP1, SNRPD1, SNRPD2, SF3B1, SNRNP200, MRPS16, MRPS9, CALM1, PPP2R1A, YWHAZ. Similarly, TP53BP1, RPS15, EFTUD2, MRPL16, MRPL17, MRPS14, RPL35A, MRPL32, MRPS6, POLR2G were selected as hub genes.Conclusions: Using the messenger ribonucleic acid (mRNA) expression data from GSE166253, we explore the pathogenesis of retesting positive in COVID-19 from the immune mechanism and molecular level. We found TP53BP1, SNRPD1 and SNRPD2 as hub genes in RTP patients. Hence, their regulatory pathway is vital to the management and prognostic prediction of RTP patients, rendering the further study of these hub genes necessary.