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.
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.
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 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.
Objective A series of bioinformatics methods were used to identify ferroptosis related biomarkers in lupus nephritis (LN). Methods We retrieved sequencing data of GSE112943 from the GEO (Gene Expression Omnibus) database and screened LN differentially expressed genes. We searched for ferroptosis-related gene (FRG) through FerrDb database, and screened LN-FRG. We conducted enrichment analysis on the LN-FRGs using David online bioinformatics database and screened the core LN-FRG using cytoHubba. We used external data sets to verify the core LN-FRGs, constructed competing endogenous RNA networks, and conducted molecular docking analysis. Results A total of 37 LN-FRGs were selected through screening. These genes are mainly enriched in inflammation, immune regulation and ferroptosis related signaling pathways. Through the cytoHubba and external dataset validation, the key core LN-FRG of ATF3 (activating transcription factor 3) was ultimately identified, and its expression was significantly increased in LN (P<0.05). Molecular docking analysis showed that ATF3 was closely bound to SLC7A11 and NRF2, and may participate in the occurrence and development of LN through the microRNA-27-ATF3 regulation axis. Conclusion The pivotal gene ATF3 may participate in the inflammation and immune injury of LN through ferroptosis.
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 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.
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.
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.