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.
ObjectiveTo investigate the expression and biological function of centromere protein F (CENPF) in non-small cell lung cancer (NSCLC) and the association with prognosis.MethodsThrough retrieving and analyzing the bioinformatics data such as Oncomine database, Human Protein Atlas (HPA), Kaplan-Meier Plotter, STRING and DAVID database, the expression of CENPF in both normal tissues and cancer tissues of lung cancer patients was identified, and the protein interaction network analysis, functional annotation and pathway analysis of CENPF with its associated genes were carried out.ResultsCENPF was overexpressed in lung adenocarcinoma tissues, but not in normal tissues. The median overall survival (OS) of NSCLC patients with low expression of CENPF was significantly longer than that of patients with high expression of CENPF. Further sub-analysis showed that low expression group from lung adenocarcinoma patients had longer median disease-free survival and OS compared with high expression group patients. CENPF and its associated hub genes mainly affected the protein K11-linked ubiquitination in biological process, anaphase-promoting complex (APC) in cell composition, ATP binding in molecular function, and cell cycle in KEGG pathway.ConclusionCENPF is regulated in tumorigenesis and progression of NSCLC, and its protein expression level has the value of early diagnosis and prognosis evaluation in lung adenocarcinoma. It is suggested that CENPF gene can be a potential target for molecular targeted therapy of NSCLC.
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.
Objective To detect the expression and clinical significance of POLD1 gene in non-small cell lung cancer (NSCLC) via bioinformatics method. Methods The expression difference of POLD1 in NSCLC tissue and normal lung tissue was investigated by TIMER database. UALCAN database was used to further verify different expression of POLD1 as well as the relationship between POLD1 expression and clinicopathological characteristics of NSCLC. The correlation between POLD1 gene and prognosis of NSCLC patients was detected by GEPIA and TIMER database. cBioPortal database was used to analyze frequencies of POLD1 gene mutation. POLD1-related protein-protein interaction network was constructed by STRING database. The relationship between POLD1 and immune infiltration was based on TISIDB database. Results The expression of POLD1 gene in lung adenocarcinoma and lung squamous cell carcinoma was significantly higher than that in normal lung tissue. In lung adenocarcinoma, patients with lower POLD1 level showed better prognosis. 1.2% of lung adenocarcinoma patients and 1.8% of lung squamous cell carcinoma patients carried mutated POLD1 gene, mainly missense mutations. POLD1 may interact with POLD2, POLD3, POLD4, POLE, RPA1, PCNA, MSH6, MSH2 and FEN1. The biological processes include DNA replication, mismatch repair, etc. Besides, the expression of POLD1 in NSCLC was correlated with the number of different immune cells. Conclusions The POLD1 gene is highly expressed in NSCLC patients, and negatively related with survival prognosis in patients of lung adenocarcinoma. POLD1 gene may be a potential diagnostic target and prognostic marker in NSCLC.
To screen new tuberculosis diagnostic antigens and vaccine candidates, we predicted the epitopes of Mycobacterium tuberculosis latent infection-associated protein Rv2004c by means of bioinformatics. The homology between Rv2004c protein and human protein sequences was analyzed with BLAST method. The second structures, hydrophilicity, antigenicity, flexibility and surface probability of the protein were analyzed to predict B cell epitopes and T cell epitopes by Protean software of DNAStar software package. The Th epitopes were predicted by RANKPEP and SYFPEITHI supermotif method, the CTL epitopes were predicted by means of combination analyses of SYFPEITHI supermotif method, BIMAS quantitative motif method and NetCTL prediction method. The peptide sequences with higher scores were chosen as the candidate epitopes. Blast analysis showed that Rv2004c protein had low homology with human protein. This protein had abundant secondary structures through analysis of DNAStar software, the peptide segments with high index of hydrophilicity, antigenicity, surface probability and flexibility were widely distributed and were consistent with segments having beta turn or irregular coil. Ten candidates of B cell epitopes were predicted. The Th epitopes of Rv2004c protein were located after the 200th amino acid. Of 37 Th cell epitopes predicted, there were more epitopes of HLA-DRB1*0401 and HLA-DRB1*0701 phenotypes, and the MHC restrictive types of some Th cell epitopes exist cross overlap. Of 10 CTL epitopes predicted, there were more number and higher score of HLA-A2 restricted epitopes. Therefore Mycobacterium tuberculosis Rv2004c protein is a protein antigen with T cell and B cell epitopes, and is expected to be a new target protein candidate for tuberculosis diagnosis and vaccine.
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.
Objective To investigate the relationship between miR-3187-5p in peripheral blood and pericardial drainage after coronary artery bypass grafting (CABG) and postoperative atrial fibrillation (POAF). Methods Patients who underwent CABG in the Heart Center of Beijing Chao-Yang Hospital from March to May 2022 were enrolled. Peripheral blood and pericardial drainage were collected at 0 h after surgery (immediate time for patients to return to ICU from operating room) to detect miR-3187-5p, and perioperative confounding factors were also collected. The miR-3187-5p was measured by quantitative real-time PCR and its regulated target genes were analyzed by bioinformatics. Results A total of 15 patients were enrolled, including 9 males and 6 females with an average age of 65.6±8.2 years. The incidence rate of POAF was 40.0%. miR-3187-5p in pericardial drainage at 0 h after surgery was an independent predictor for POAF. A total of 1 642 target genes of miR-3187-5p were predicted. GO function enrichment analysis and KEGG signal pathway enrichment analysis showed that target genes of miR-3187-5p were enriched in TGF-β, MAPK, Wnt and other classical collagen metabolic signal pathways, which might activate collagen metabolism by negatively regulating SMAD6 and other inhibitors of the pathways. Conclusion This study is the first to find that miR-3187-5p in pericardial drainage at 0 h after surgery is a potential, novel, and predictive factor for POAF, which may be related to the regulation of myocardial fibrosis signal pathways like TGF-β, MAPK and Wnt pathways, promoting the early collagen metabolism imbalance after CABG, increasing the collagen deposition in the atrium, and then promoting the early structural reconstruction after CABG and leading to the occurrence of POAF. The result provides a research basis for the accurate prediction and prevention of clinical POAF.
Objective To identify potential hub genes and key pathways in the early period of septic shock via bioinformatics analysis. MethodsThe gene expression profile GSE110487 dataset was downloaded from the Gene Expression Omnibus database. Differentially expressed genes were identified by using DESeq2 package of R project. Then Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses were constructed to investigated pathways and biological processes using clusterProfiler package. Subsequently, protein-protein interaction (PPI) network was mapped using ggnetwork package and the molecular complex detection (MCODE) analysis was implemented to further investigate the interactions of differentially expressed genes using Cytoscape software. Results A total of 468 differentially expressed genes were identified in septic shock patients with different responses who accepted early supportive hemodynamic therapy, including 255 upregulated genes and 213 downregulated genes. The results of GO and the KEGG pathway enrichment analysis indicated that these up-regulated genes were highly associated with the immune-related biological processes, and the down-regulated genes are involved in biological processes related to organonitrogen compound, multicellular organismal process, ion transport. Finally, a total of 23 hub genes were identified based on PPI and the subcluster analysis through MCODE software plugin in Cytoscape, which included 19 upregulated hub genes, such as CD28, CD3D, CD8B, CD8A, CD160, CXCR6, CCR3, CCR8, CCR9, TLR3, EOMES, GZMB, PTGDR2, CXCL8, GZMA, FASLG, GPR18, PRF1, IDO1, and additional 4 downregulated hub genes, such as CNR1, GPER1, TMIGD3, GRM2. KEGG pathway enrichment analysis and GO functional annotation showed that differentially expressed genes were primarily associated with the items related to cytokine-cytokine receptor interaction, natural killer cell mediated cytotoxicity, hematopoietic cell lineage, T cell receptor signaling pathway, phospholipase D signaling pathway, cell adhesion molecules, viral protein interaction with cytokine and cytokine receptor, primary immunodeficiency, graft-versus-host disease, type 1 diabetes mellitus. Conclusions Some lymphocytes such as T cells and natural killer cells, cytokines and chemokines participate in the immune process, which plays an important role in the early treatment of septic shock, and CD160, CNR1, GPER1, and GRM2 may be considered as new biomarkers.
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 aberrantly expressed genes in hepatocellular carcinoma (HCC) and their relationship with prognosis of HCC through bioinformatics analysis. Methods Five datasets related to HCC were selected from the GeneExpression Omnibus database to explore differentially expressed genes (DEGs), followed by further gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The co-upregulated genes CNIH4 and TOMM40 were selected to explore the differences in their expressions in HCC tissues and normal tissues, and to explore the relationship between their expressions and the 5-year survival of patients by using TCGA database. Tissues and paraneoplastic tissues of eight cases of HCC who underwent surgery at the Guangdong Second Provincial General Hospital were collected to verify the expression differences of CNIH4 and TOMM40L mRNA. Results A total of 25 up-regulated genes and 21 down-regulated genes were identified in this study. The results of GO analysis and KEGG analysis indicated that DEGs were mainly related to catabolism, cell division, DNA replication and repair. The results of TCGA database analysis showed that the expression of up-regulated genes CNIH4 mRNA and TOMM40L mRNA were up-regulated in HCC tissues as compared with normal tissues (P<0.05) and that the 5-year survival of patients in the high expression group was worse than that in the low expression group (P<0.05). The results of clinical samples showed that CNIH4 mRNA and TOMM40L mRNA were up-regulated in HCC tissues as compared with paraneoplastic tissues. Conclusion CNIH4 and TOMM40L genes are up-regulated in HCC tissues, and their high expressions are associated with poor prognosis, and may be potential biomarkers and prognostic indicators for HCC.