Objective To explore key genes and mechanisms of depression aggravating Crohn disease. Methods In March 2023, the Public Health Genomics and Precision Health Knowledge Base and Gene Expression Omnibus database were used to identify the overlapping differentially expressed genes between Crohn disease and depression and the key genes were screened by Metascape, STRING, Cytoscape, and protein interaction network analysis. The Gene Expression Omnibus database was used to analyze the correlations between key genes and clinical pathologies such as Crohn Disease Endoscopic Index of Severity and intestinal microvilli length. Results There were 137 overlapping differentially expressed genes between Crohn disease and depression, and 25 key genes were further screened out. Among them, CREB1, FKBP5, MAPT, NTSR1, OXTR, PROK2, POMC, HTR2B, and PPARGC1A genes were significantly correlated with multiple clinical parameters. The functions of PROK2 and PROK2-related genes were mainly enriched in neutrophil and granulocyte migration, neutrophil and granulocyte chemotaxis, etc. Conclusions There are 25 key genes, especially CREB1, FKBP5, MAPT, NTSR1, OXTR, PROK2, POMC, HTR2B, and PPARGC1A, that possibly contribute to the establishment and deterioration of Crohn disease caused by depressive disorder. Among these genes, PROK2 showes the possibility of regulating immune cell (neutrophils and CD8+ T cells) infiltration.
Objective To explore the differential expression of circular RNAs (circRNAs) in polycystic ovary syndrome (PCOS) by bioinformatics, and predict the microRNAs (miRNAs) associated with them. Methods The expression profile of cumulus cells gene chip in PCOS was searched in the Gene Expression Omnibus database, and differential circRNAs were screened by GEO2R tool of the database. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes signaling pathways of different circRNA genes were analyzed using the DAVID 6.8 database. Circular RNA interactome was used to predict the potential regulated miRNAs. Cytoscape software was used to establish circRNA-miRNA network map. The potential regulatory miRNAs were predicted by the 10 circRNAs with the most significant differences in up-regulation and down-regulation. Results A total of 247 circRNAs were obtained in PCOS, and 277 miRNAs binding to up-regulated circRNA genes and 125 miRNAs binding to down-regulated circRNA genes were predicted. The top 10 miRNAs that could bind to multiple differential circRNAs were hsa-miR-557, hsa-miR-507, hsa-miR-224, hsa-miR-136, hsa-miR-127-5p, hsa-miR-579, hsa-miR-502-5p, hsa-miR-186, hsa-miR-1253, and hsa-miR-432. Conclusion The differential expression analysis of circRNAs is helpful to understand the main role of circRNAs in PCOS, and the prediction of potential regulated miRNAs can help to understand the pathogenesis of the disease.
The purpose of this paper is to present the research on the molecular biological characteristics of proto-oncogene pim-2 and to analyze the related mechanism. Proto-oncogene pim-2 was studied and analyzed by the bioinformatics method and technology. With an online server, the chromosomal localization of pim-2 gene was analyzed, and the exon, open reading frame, CpG island and miRNAs complementary fragments and the like were predicted. With bioinformatics software, the physicochemical property of transcription protein of proto-oncogene pim-2 and various modification sites of protein sequence, such as ubiquitination and glycosylation, were predicted, the antigenic index was calculated, and the spatial structural was modeled. The research findings showed that the proto-oncogene pim-2 comprised six exons, the CDS (coding sequence) transcribed a section of peptide chain including 311 amino acids, a gene promoter has a CpG island, and the 3'UTR region contains an miRNA gene. The molecular weight of the Pim-2 protein was 34, 188.47, the isoelectric point was 5.78, the instability index was 45.87, and the extinction coefficient was 279nm. A plurality of covalent modification sites, two ubiquitination sites, four glycosylation sites, an SUMO sumoylation site, a nitrosation site, two palmitoylation sites and sixteen regions with higher antigenic index were distributed in the protein sequence. This research showed that the related regions and modification sites distributed on the sequence of proto-oncogene pim-2 were closely related to the carcinogenic effect thereof.
ObjectiveTo analyze the expression and clinical significance of cyclin-dependent kinase 1 (CDK1) in lung adenocarcinoma by bioinformatics.MethodsBased on the gene expression data of lung adenocarcinoma patients in The Cancer Genome Atlas (TCGA), the differential expression of CDK1 in lung adenocarcinoma tissues and normal lung tissues was analyzed. The expression of CDK1 gene in lung adenocarcinoma was analyzed by UALCAN at different angles. Survival analysis of different levels of CDK1 gene expression in lung adenocarcinoma was performed using Kaplan-Meier Plotter. Correlation Cox analysis of CDK1 expression and overall survival was based on clinical data of lung adenocarcinoma in TCGA. Gene set enrichment analysis was performed on gene sequences related to CDK1 expression in clinical cases. The protein interaction network of CDK1 from Homo sapiens was obtained by STRING. CDK1-related gene proteins were obtained and analyzed by the web server Gene Expression Profiling Interactive Analysis (GEPIA).ResultsBased on the analysis of TCGA gene expression data, CDK1 expression in lung adenocarcinoma was higher than that in normal lung tissues. UALCAN analysis showed that high CDK1 expression may be associated with smoking. Survival analysis indicated that when CDK1 gene was highly expressed, patients with lung adenocarcinoma had a poor prognosis. Univariate and multivariate Cox regression analysis of CDK1 expression and overall survival showed that high CDK1 expression was an independent risk factor for survival of patients with lung adenocarcinoma. Gene set enrichment analysis revealed that high CDK1 expression was closely related to DNA replication, cell cycle, cancer pathway and p53 signaling pathway.ConclusionCDK1 may be a potential molecular marker for prognosis of lung adenocarcinoma. In addition, CDK1 regulation may play an important role in DNA replication, cell cycle, cancer pathway and p53 signaling pathway in lung adenocarcinoma.
ObjectiveTo explore the mechanism of DDX46 regulation of esophageal squamous cell carcinoma.MethodsPicture signals of fluorescence in gene array were scanned and differential expression of gene in two groups (a DDX46-shRNA-LV group and a control-LV group) were compared by GCOSvL.4 software. These differential expressed genes were analyzed by bioinformatics methods finally, and validated by quantitative real time polymerase chain reaction (qRT-PCR) analysis.ResultsAccording to the screening criteria of fold change ≥2 and P<0.05, 1 006 genes were differentially expressed after DDX46 knockdown, including 362 up-regulated and 644 down-regulated genes. Bioinformatics analysis and gene co-expression network building identified that these differentially expressed genes were mainly involved in cell cycle, proliferation, apoptosis, adhesion, energy metabolism, immune response, etc. Phosphatidylinositol 3-kinase (PI3K) was the key molecule in the network. The results of RT-qPCR were completely consistent with the results of gene microarra.ConclusionBioinformatics can effectively exploit the microarray data of esophageal squamous cell carcinoma after DDX46 knockdown, which provides a valuable clue for further exploration of DDX46 tumorigenesis mechanism and helps to find potential drug therapy.
ObjectiveTo explore the significance of mesenchymal epithelial transition factor (MET) as a clinical prognostic evaluation index for patients with pancreatic cancer based on bioinformatics analysis.MethodsThe GSE28735 and GSE62452 gene chips from GEO database were downloaded and the difference of MET gene expression between cancer and adjacent cancerous tissues were analyzed by bioinformatics. We downloaded pancreatic cancer gene chip from TCGA database to analyze the correlation between MET gene expression and clinicopathological features of pancreatic cancer patients and prognosis risk. Finally, the possible molecular mechanism of MET involved in pancreatic carcinogenesis was analyzed by GO and KEGG enrichment analysis.ResultsThe expression level of MET gene in pancreatic cancer tissues was significantly higher than that in adjacent cancerous tissues (P<0.001). The overall survival and disease-free survival of pancreatic cancer patients in the high MET gene expression group were lower than those in the low expression group (P<0.001). The expression level of MET gene was related to the age of pancreatic cancer patients, T stage, and histological grading of tumors (P<0.05), and high MET gene expression, age >65 years, and N1 stage were independent risk factors affecting the prognosis of pancreatic cancer patients. KEGG enrichment analysis showed that MET was mainly related to PI3K/AKT signaling pathway, FAK signaling pathway, and cancer transcription dysregulation and so on.ConclusionMET may be a valuable tumor marker for pancreatic cancer and can predict the poor prognosis of patients with pancreatic cancer.
ObjectiveTo investigate the relation between disulfidptosis-related genes (DRGs) and prognosis or immunotherapy response of patients with pancreatic cancer (PC). MethodsThe transcriptome data, somatic mutation data, and corresponding clinical information of the patients with PC in The Cancer Genome Atlas (TCGA) were downloaded. The DRGs mutated in the PC were screened out from the 15 known DRGs. The DRGs subtypes were identified by consensus clustering algorithm, and then the relation between the identified DRGs subtypes and the prognosis of patients with PC, immune cell infiltration or functional enrichment pathway was analyzed. Further, a risk score was calculated according to the DRGs gene expression level, and the patients were categorized into high-risk and low-risk groups based on the mean value of the risk score. The risk score and overall survival of the patients with high-risk and low-risk were compared. Finally, the relation between the risk score and (or) tumor mutation burden (TMB) and the prognosis of patients with PC was assessed. ResultsThe transcriptome data and corresponding clinical information of the 177 patients with PC were downloaded from TCGA, including 161 patients with somatic mutation data. A total of 10 mutated DRGs were screened out. Two DRGs subtypes were identified, namely subtype A and subtype B. The overall survival of PC patients with subtype A was better than that of patients with subtype B (χ2=8.316, P=0.003). The abundance of immune cell infiltration in the PC patients with subtype A was higher and mainly enriched in the metabolic and conduction related pathways as compaired with the patients with subtype B. The mean risk score of 177 patients with PC was 1.921, including 157 cases in the high-risk group and 20 cases in the low-risk group. The risk score of patients with subtype B was higher than that of patients with subtype A (t=14.031, P<0.001). The overall survival of the low-risk group was better than that of the high-risk group (χ2=17.058, P<0.001), and the TMB value of the PC patients with high-risk was higher than that of the PC patients with low-risk (t=5.642, P=0.014). The mean TMB of 161 patients with somatic mutation data was 2.767, including 128 cases in the high-TMB group and 33 cases in the low-TMB group. The overall survival of patients in the high-TMB group was worse than that of patients in the low-TMB group (χ2=7.425, P=0.006). ConclusionDRGs are closely related to the prognosis and immunotherapy response of patients with PC, and targeted treatment of DRGs might potentially provide a new idea for the diagnosis and treatment of PC.
ObjectiveTo identify the core genes involved in the great saphenous varicose veins (GSVVs) through bioinformatics method. MethodsThe transcriptional data of GSVVs and normal great saphenous vein tissues (control tissues) were downloaded from the gene expression omnibus database. The single sample gene set enrichment analysis (ssGSEA) was used to calculate the Hallmark score. The weighted gene co-expression network analysis (WGCNA) combined with machine learning algorithms was used to screen the key genes relevant GSVVs. The protein-protein interaction (PPI) analysis was performed using the String database, and the receiver operating characteristic (ROC) curve was used to reflect the discrimination ability of the target genes for GSVVs. ResultsCompared with the control tissues, there were 548 up-regulated genes and 706 down-regulated genes in the GSVVs tissues, the Hallmark points of KRAS signaling and apical junction were down-regulated, while which of peroxisomes, coagulation, reactive oxygen species pathways, etc. were up-regulated in the GSVVs tissues. A total of 639 differentially expressed genes relevant GSVVs were obtained and 165 interaction relations between proteins encoded by 372 genes, and the top 10 genes with the highest betweeness values, ADAM10, APP, NCBP2, SP1, ASB6, ADCY4, HP, UBE2C, QSOX1, and CXCL1, were located at the center of the interaction relation. And the core genes were mainly related to copper ion homeostasis, neutrophil degranulation G protein coupled receptor signaling, response to oxidative stress, and regulation of amide metabolism processes. The SP1 and QSOX1 were both Hub genes. The expressions of the SP1 and QSOX1 in the GSVVs tissues were significantly up-regulated as compared with the control tissues. The areas under the ROC curves of SP1 and QSOX1 in distinguishing GSVVs tissues from normal tissues were 0.972 and 1.000, respectively. ConclusionsSP1 and QSOX1 are core genes in the occurrence and development of GSVVs. Regulation of SP1 or QSOX1 gene is expected to achieve precise treatment of GSVVs.
Objective To analyze the pathways, biomarkers and diagnostic genes of systemic sclerosis associated interstitial lung disease (SSc-ILD) using bioinformatics. Methods SSc-ILD related gene data sets from April to June 2023 were downloaded from the Gene Expression Omnibus database for differential analysis and enrichment analyses including gene ontology analysis, Kyoto Encyclopedia of Genes and Genomes analysis, disease ontology analysis, and gene set enrichment analysis. Least absolute shrinkage and selection operator regression and support vector machine algorithms were applied to screen and take the intersection to get the diagnostic genes and validate the results. Disease-related data were analyzed by immune cell infiltration. Results A total of 178 differential genes were obtained, and enrichment analyses showed that they were related to 5 signaling pathways and associated with 3 diseases. The diagnostic genes screened were TNFAIP3, ID3, and NT5DC2, and immune cell infiltration showed that the diagnostic genes were associated with plasma cells, resting mast cells, activated natural killer cells, macrophage M1 and M2, resting dendritic cells, and activated dendritic cells. Conclusion The screened diagnostic genes and immune cells may be involved in the development of SSc-ILD.
ObjectiveTo investigate differentially expressed genes (DEGs) and potential molecular mechanisms between hepatitis C-related hepatocellular carcinoma (HCV-HCC) and hepatitis B-related HCC (HBV-HCC). MethodsThe data of HCV-HCC and HBV-HCC gene expressions were downloaded and integrated from the public gene expression database, and the limma package was used to investigate the DEGs between the HCV-HCC and HBV-HCC samples. The gene set enrichment analysis (GSEA) was used to explore the differences in suppressed or activated gene sets between the HCV-HCC and HBV-HCC samples, and the MCODE was used to explore the key molecular modules, and then the potential biological processes and molecular pathways of the key molecular modules were analyzed. The effect of key genes on survival of the HCC patients was analyzed by the Kaplan-Meier-Plotter database.ResultsIn this study, 119 HBV-HCC samples and 163 HCV-HCC samples were obtained, and the 199 DEGs were screened out. Compared with HBV-HCC, the activated gene sets of HCV-HCC were mainly enriched in the gene sets of inflammation, complement, up-regulation of genes in response to interferon, up-regulation of genes in response to KRAS, genes regulated by the nuclear factor- κB-tumor necrosis factor pathway, and apoptosis. However, the cell cycle-related gene sets were obviously suppressed. Eight key molecular modules enriched by DEGs were found, which included 18 key genes (IFI27, DDX60, MX1, IRF9, OAS3, OAS1, RSAD2, GBP4, HERC6, ISG15, IFIT1, CMPK2, EPSTI1, IFI44, IFI44L, HERC5, IFITM1, CXCL10). GO analysis showed that the biological process was mainly concentrated in the body response related to virus infection, the molecular component was mainly in the host cells, and the molecular function was mainly enriched in the biological combination. KEGG analysis showed that the key genes were mainly involved in the molecular signaling pathway related to virus infection. The survival analysis showed that the 9 key genes (CXCL10, HERC6, DDX60, IFITM1, IFI27, GBP4, IFI44L, IFI44, MX1) were closely related to better prognosis of patients with HCC (HR<1, P<0.05). ConclusionsThere is an essential difference between HBV-HCC and HCV-HCC. Occurrence of HCV-HCC is mainly related to virus infection and immune response induced by the virus. Therefore, for HCV infection, active antiviral treatment is necessary for avoiding hepatitis turning into chronic viral infection and preventing or blocking HCV infection converting to HCC.
