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 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 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.
ObjectiveTo explore the mechanism of paucigranulocytic asthma and to find therapeutic target for paucigranulocytic asthma.MethodsGSE143303 data and platform information were downloaded from GEO. Gene Set Enrichment Analysis were performed to construct positive and negative gene-gene interaction network correlation with paucigranulocytic asthma. Differential expression analysis, pathway commonality analysis were performed with R language.ResultsGSE143303 data set contained 47 endobronchial biopsies from adult (16 cases of paucigranulocytic asthma, 13 cases of healthy control). Compared with control group, the paucigranulocytic asthma group had 115 differential genes set (37 positive and 78 negative). The results of pathway commonality analysis showed that the crosslink existed within the negative gene-gene interaction network correlation with paucigranulocytic asthma. Among these, most of the genes belonged to the protein HLA gene family. Differential expression analysis show that HLA-DQB1, HLA-DRB5 were differential genes and TNFRSF13B was significantly downregulated genes in the intersect genes.ConclusionTNFRSF13B, HLA-DQB1, HLA-DRB5 and regulatory networks associated with them are the crucial factors contributing to paucigranulocytic asthma.
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.
ObjectiveTo construct a prognostic model of esophageal squamous cell carcinoma (ESCC) based on immune checkpoint-related genes and explore the potential relationship between these genes and the tumor microenvironment (TME). Methods The transcriptome sequencing data and clinical information of immune checkpoint genes of samples from GSE53625 in GEO database were collected. The difference of gene expression between ESCC and normal paracancerous tissues was evaluated, and the drug sensitivity of differentially expressed genes in ESCC was analyzed. We then constructed a risk model based on survival-related genes and explored the prognostic characteristics, enriched pathway, immune checkpoints, immune score, immune cell infiltration, and potentially sensitive drugs of different risk groups. ResultsA total of 358 samples from 179 patients were enrolled, including 179 ESCC samples and 179 corresponding paracancerous tissues. There were 33 males and 146 females, including 80 patients≤60 years and 99 patients>60 years. 39 immune checkpoint genes were differentially expressed in ESCC, including 14 low expression genes and 25 high expression genes. Drug sensitivity analysis of 8 highly expressed genes (TNFRSF8, CTLA4, TNFRSF4, CD276, TNFSF4, IDO1, CD80, TNFRSF18) showed that many compounds were sensitive to these immunotherapy targets. A risk model based on three prognostic genes (NRP1, ICOSLG, HHLA2) was constructed by the least absolute shrinkage and selection operator analysis. It was found that the overall survival time of the high-risk group was significantly lower than that of the low-risk group (P<0.001). Similar results were obtained in different ESCC subtypes. The risk score based on the immune checkpoint gene was identified as an independent prognostic factor for ESCC. Different risk groups had unique enriched pathways, immune cell infiltration, TME, and sensitive drugs. Conclusion A prognostic model based on immune checkpoint gene is established, which can accurately stratify ESCC and provide potential sensitive drugs for ESCC with different risks, thus providing a possibility for personalized treatment of ESCC.
Objective To explore the role of high endothelial venule (HEV) in chronic obstructive pulmonary disease (COPD) at the single cell level. Methods A total of 219257 cells from the lung tissues of 18 COPD patients and 28 healthy controls in the GEO public database (GSE136831) were used to analyze the relationship between HEV with T lymphocytes, B lymphocytes, and dendritic cells. Results Endothelial cells were extracted using single cell analysis technique, and sorting out venous endothelium, CCL14, IGFBP7, POSTN were used as marker genes for HEV endothelial cells. The ratio of HEV endothelial cells was also identified as up-regulated expression in COPD. The function of the differential genes of HEV endothelial cells was analyzed, suggesting the presence of immune regulation. By trajectory analysis, it was suggested that the differential genes of HEV endothelial cells were enriched for extracellular matrix deposition in late development. Finally, by receptor-ligand pairing, it was suggested that HEV endothelial cells was recruited through a series of ligands with T lymphocytes, B lymphocytes, and dendritic cells. Conclusions HEV endothelial cells are elevated in COPD and have an immunomodulatory role by secreting a series of ligands after recruiting T lymphocytes, B lymphocytes as well as dendritic cells for immune action. HEV may be a potential target for the study of COPD therapy.
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 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 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.