ObjectiveTo analyze the expression of cold-induced RNA-binding protein (CIRBP) in lung adenocarcinoma and its clinical significance based on bioinformatics, in order to provide a new direction for the study of therapeutic targets for lung adenocarcinoma.MethodsThe CIRBP gene expression data and patient clinical information data in lung adenocarcinoma tissues and adjacent tissues were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus databases. The expression of CIRBP in lung adenocarcinoma was analyzed. Furthermore, its relationship with clinicopathological features and prognosis in patients with lung adenocarcinoma was analyzed. GO and KEGG enrichment analysis were carried out for the screened genes. The CIRBP protein interaction network was constructed by STRING, and the correlation analysis was carried out using the GEPIA online website.ResultsThe expression level of CIRBP gene in lung adenocarcinoma tissues was significantly lower than that in adjacent tissues (P<0.01), and its expression level was correlated with T stage and N stage in clinicopathological features. The prognosis of patients with high CIRBP expression in lung adenocarcinoma was significantly better than that with low CIRBP expression. Univariate and multivariate Cox regression analysis showed that CIRBP was an independent prognostic factor in patients with lung adenocarcinoma. GO functional annotation showed its enrichment in organelle fission, nuclear fission, chromosome separation, and DNA replication, etc. KEGG analysis showed that it was mainly involved in cell cycle and DNA replication. Protein interaction network and GEPIA online analysis showed that the expression level of CIRBP was negatively correlated with the expression level of cyclin B2.ConclusionCIRBP gene is down-regulated in lung adenocarcinoma tissues, and its expression level is closely related to patient prognosis. CIRBP gene may be a potential therapeutic target and prognostic marker for lung adenocarcinoma.
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 relationship between the expression of carbonic anhydrase 3 (CA3) in breast cancer tissues, its prognostic potential and the number of immune cells by a variety of online databases. Methods GEPIA2.0 and TIMER databases were used to analyze the difference of CA3 mRNA expression in breast cancer tissues. Bc-GenExMinerv4.7 database was used to analyze the difference of CA3 mRNA expression in breast cancer subcategories. Kaplan-Meier plotter, Bc-GenExMinerv4.7 and PrognoScan databases were used to analyze the effect of CA3 mRNA expression levels on prognosis of patient. LinkedOmics database was used to analyze of the biological behavior involved in CA3 co-expressed genes. TIMER database was used to analyze the relationship between CA3 mRNA expression and immune cells infiltration in breast cancer tissues. Results The expression of CA3 mRNA in breast cancer tissues was lower than that in normal breast tissues (P<0.05), and the expression levels of CA3 mRNA were higher in ER negative (P<0.05), PR negative (P<0.05), HER2 negative (P<0.05) and no lymphatic metastasis (P<0.05). In addition, the expression level of CA3 in breast cancer patients with high Ki67 expression was lower (P<0.05) and closely related to SBR and NPI grade (P<0.05). Breast cancer patients with low expression of CA3 mRNA had lower overall survivall, recurrence free survival, and disease free survival ( P<0.05). Ten of the top 50 positively correlated co-expressed genes screened out had low risk ratio (P<0.05), and 11 of the top 50 negatively correlated co-expressed genes screened out had high risk ratio (P<0.05). The expression of CA3 mRNA was positively correlated with CD4+ T cells and CD8+ T cells in breast cancer tissues (rs=0.175, P<0.001; rs=0.137, P<0.001), and negatively correlated with T cell failure markers LAG3, TIM-3 and PVRL2 (rs=–0.100, P<0.01; rs=–0.143, P<0.001; rs=–0.082, P<0.05). Conclusions The low expression of CA3 mRNA in breast cancer tissues is correlated with the occurrence, development and prognosis of breast cancer. CA3 can be used as a potential independent prognostic marker for breast cancer and may be related to immune infiltration.
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.
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.
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 role and clinical significance of cell-cycle dependent kinase 1 (CDK1) and its upstream and downstream molecules in the development of malignant peripheral nerve sheath tumor (MPNST) through the analysis of clinical tissue samples. Methods A total of 56 tumor samples from MPNST patients (“Tianjin” dataset) who underwent surgical resection, confirmed by histology and pathology between September 2011 and March 2020, along with 17 normal tissue samples, were selected as the research subjects. MPNST-related hub genes were identified through transcriptome sequencing, bioinformatics analysis, immunohistochemistry staining, and survival analysis, and their expression levels and prognostic associations were analyzed. Results Transcriptome sequencing and bioinformatics analysis revealed that upregulated genes in MPNST were predominantly enriched in cell cycle-related pathways, with CDK1 occupying a central position among all differentially expressed genes. Further differential analysis demonstrated that CDK1 mRNA expression in sarcoma tissues was significantly higher than in normal tissues [based on searching the cancer genome atlas (TCGA) dataset, P<0.05]. In MPNST tissues, CDK1 mRNA expression was not only significantly higher than in normal tissues (based on Tianjin, GSE141438 datasets, P<0.05), but also significantly higher than in neurofibromatosis (NF) and plexiform neurofibromas (PNF) (based on GSE66743 and GSE145064 datasets, P<0.05). Immunohistochemical staining results indicated that the expression rate of CDK1 protein in MPNST tissues was 40.31%. Survival analysis results demonstrated that CDK1 expression was associated with poor prognosis. The survival time of MPNST patients with high CDK1 mRNA expression was significantly lower than that of the low expression group (P<0.05), and the overall survival trend of patients with positive CDK1 protein expression was worse than that of patients with negative CDK1 expression. Additionally, differential analysis of CDK family genes (CDK1-8) revealed that only CDK1 was significantly upregulated in MPNST, NF, and PNF. Conclusion Increased expression of CDK1 is associated with poor prognosis in MPNST patients. Compared to other CDK family members, CDK1 exhibits a unique expression pattern, suggesting its potential as a therapeutic target for MPNST.
Objective To explore the key genes, pathways and immune cell infiltration of bicuspid aortic valve (BAV) with ascending aortic dilation by bioinformatics analysis. Methods The data set GSE83675 was downloaded from the Gene Expression Omnibus database (up to May 12th, 2022). Differentially expressed genes (DEGs) were analyzed and gene set enrichment analysis (GSEA) was conducted using R language. STRING database and Cytoscape software were used to construct protein-protein interaction (PPI) network and identify hub genes. The proportion of immune cells infiltration was calculated by CIBERSORT deconvolution algorithm. Results There were 199 DEGs identified, including 19 up-regulated DEGs and 180 down-regulated DEGs. GSEA showed that the main enrichment pathways were cytokine-cytokine receptor interaction, pathways in cancer, regulation of actin cytoskeleton, chemokine signaling pathway and mitogen-activated protein kinase signaling pathway. Ten hub genes (EGFR, RIMS3, DLGAP2, RAPH1, CCNB3, CD3E, PIK3R5, TP73, PAK3, and AGAP2) were identified in PPI network. CIBERSORT analysis showed that activated natural killer cells were significantly higher in dilated aorta with BAV. Conclusions These identified key genes and pathways provide new insights into BAV aortopathy. Activated natural killer cells may participate in the dilation of ascending aorta with BAV.
Objective To explore depression-related biomarkers and potential therapeutic drugs in order to alleviate depression symptoms and improve patients’ quality of life. Methods From November 2022 to January 2024, gene expression profiles of depression patients and healthy volunteers were downloaded from the Gene Expression Omnibus database. Differential expression analysis was performed to identify differentially expressed genes. Enrichment analysis of these genes was conducted, followed by the construction of a protein-protein interaction network. Finally, Cytoscape software with the Cytohubba plugin was used to identify potential key genes, and drug prediction was performed. Results Through differential expression analysis, a total of 110 differentially expressed genes (74 upregulated and 36 downregulated) were identified. Protein-protein interaction network identified 10 key genes, and differential expression analysis showed that 8 of these genes (CPA3, HDC, IL3RA, ENPP3, PTGDR2, VTN, SPP1, and SERPINE1) exhibited significant differences in expression levels between healthy volunteers and patients with depression (P<0.05). Enrichment analysis revealed that the upregulated genes were significantly enriched in pathways related to circadian rhythm, niacin and nicotinamide metabolism, and pyrimidine metabolism, while the downregulated genes were primarily enriched in extracellular matrix-receptor interaction and interleukin-17 signaling pathways. Six overlapping verification genes (SALL2, AKAP12, GCSAML, CPA3, FCRL3, and MS4A3) were obtained across two datasets using the Wayn diagram. Single-cell sequencing analysis indicated that these genes were significantly expressed in astrocytes and neurons. Mendelian randomization analysis suggested that the FCRL3 gene might play a critical role in the development of depression. Drug prediction analysis revealed several potential antidepressant agents, such as cefotiam, harmol, lincomycin, and ribavirin. Conclusions Circadian rhythm, nicotinate and nicotinamide metabolism, and pyrimidine metabolism pathways may represent potential pathogenic mechanisms in depression. Harmol may be a potential therapeutic drug for the treatment of depression.
Objective To screen the key genes in childhood therapy-resistant asthma by bioinformatic method, and to verify its expression and diagnostic value in peripheral blood of children with therapy-resistant asthma. Methods The transcriptome dataset GSE27011 of peripheral blood mononuclear cells from healthy children (healthy control group), mild asthma (MA) children (MA group) and severe asthma (SA) children (SA group) was downloaded from the Gene Expression Omnibus of the National Center for Biotechnology Information of the United States. Key genes were obtained by using R software for gene differential expression analysis, weighted gene co-expression network analysis (WGCNA) and clinical phenotypic correlation analysis. The differential expression levels of key genes were verified in children with asthma and immune cell transcriptome datasets. Seventy-eight children with asthma and 30 healthy children who were diagnosed in the Department of Pediatrics of Tangshan People’s Hospital between September 2020 and September 2021 were selected and divided into control group, MA group and SA group. Peripheral blood samples from children with asthma and healthy children who underwent physical examination were collected to detect the expression levels of key genes and inflammatory factors interleukin (IL)-4 and IL-17 in peripheral blood of children. Receiver operating characteristic curve was used to evaluate the sensitivity, specificity and accuracy of key genes in predicting childhood therapy-resistant asthma. Results The key gene GNA15 was obtained by bioinformatic analysis. Analysis of asthma validation dataset showed that GNA15 was up-regulated in asthma groups, and was specifically expressed in eosinophils. Clinical results showed that the expression levels of IL-4, IL-17 and GNA15 among the three groups were significantly different (P<0.05). The expression levels of IL-4 and IL-17 in the MA group and the SA group were higher than those in the control group (P<0.05). Compared with the control group and the MA group, the expression level of GNA15 in the SA group was up-regulated (P<0.05). Neither the difference in the expression level of IL-4 or IL-17 between the MA group and the SA group, nor the difference in the expression level of GNA15 between the control group and the MA group was statistically significant (P>0.05). The specificity, sensitivity and accuracy of GNA15 in predicting SA were 92.90%, 80.00% and 86.10%, respectively. Conclusion GNA15 has a significant clinical value in predicting the childhood therapy-resistant asthma, and may become a potential diagnostic marker for predicting the severity of asthma in children.