ObjectiveTo explore the differentially expressed genes (DEGs) in venous leg ulcer (VLU) by bioinformatics, and further explore the molecular mechanism of the disease, predict early diagnostic markers and treatment targets.MethodsThe expression profiles of VLU were downloaded from the gene expression omnibus (GEO) database, the DEGs of VLU and inflammatory phase of normal skin healing were identified by R software and used to perform gene ontology (GO) and KEGG pathway enrichment analysis, obtaining the key genes of the pathway. We analyzed the proteins of protein interaction (PPI) network by STRING database and Cytoscape 3.2.1 software to obtain hub genes.ResultsA total of 409 DEGs were obtained, including 173 upregulted genes and 236 downregulted genes. The GO analysis showed that the upregulated DEGs mainly distributed in collagen-containing extracellular matrix (ECM), cornified envelope and collagen trimer, involved in biological processes such as skin development, keratinocyte differentiation and cornification, which mediated molecular functions such as ECM structural constituent, ECM structural constituent conferring tensile strength and integrin binding. The downregulated DEGs mainly distributed in tertiary granule, secretory granule membrane and tertiary granule membrane cornification, involved in biological processes such as response to chemokine, leukocyte migration and neutrophil chemotaxis, which mediated molecular functions such as chemokine activity, chemokine receptor binding and cytokine activity. KEGG pathway enrichment analysis results showed that the upregulated DEGs were mainly enriched in ECM-receptor interaction and protein digestion and absorption pathways, collagen type Ⅰ alpha1 chain (COL1A1), collagen type Ⅰ alpha2 chain (COL1A2), and collagen type Ⅵ alpha 6 chain (COL6A6) were the key genes of pathway; the downregulated DEGs were mainly enriched in Staphylococcus aureus infection, Toll-like receptor signaling pathway and leukocyte transendothelial migration pathways, interleukin (IL)-1β, C-X-C motif chemokine ligand 8 (CXCL8), IL-10, matrix metalloproteinase (MMP)1, and MMP9 were the key genes of pathway. The hub core genes of the PPI network were formyl peptide receptor (FPR)1, FPR2, IL-1β, IL-10, and CXCL8.ConclusionsThe results of this study indicate that the genes and signaling pathways involved in COL1A1, COL1A2, COL6A6, IL-1β, CXCL8, IL-10, MMP1, and MMP9 affect the healing of VLU. FPR1, FPR2, IL-1β, IL-10, and CXCL8 can be used as potential therapeutic targets.
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.
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.
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 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.
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 explore the functional heterogeneity of T lymphocytes in various organs after SARS-CoV-2 infection. Methods Using the public database GEO data (GSE171668, GSE159812, GSE159556, GSE167747) and the analysis method of single-cell technology, the functional differences of T lymphocytes in various organs of patients after infection with SARS-CoV-2 were analyzed. Results Through single-cell data extraction of 16 livers, 19 hearts,2 spleens, 6 brains, 58 lungs, 21 kidneys and 5 pancreases from SARS-CoV-2 infected patients, invasion genes were relatively highly expressed in T lymphocytes of the lung and pancreas. The lung had a special ability to express the interferon signaling pathway, while the expression of other organs was relatively low; at the same time, the T lymphocytes of the lung also highly expressed fatty acid binding sites. Conclusion After SARS-CoV-2 infection, compared with other organs, the lung has a special interferon-activated signaling pathway and fatty acid binding site.
ObjectiveTo explore the clinical significance and possible potential mechanism of hepatocellular carcinoma through the screening of key genes in hepatocellular carcinoma.MethodsHepatocellular carcinoma gene chip was obtained from GEO database, differentially expressed genes (DEGs) were screened by GEO2R online tools and Venn map, GO analysis and KEGG pathway analysis were performed in DAVID database, core genes were screened by STRING and Cytscape software, core genes were analyzed in Kaplan-Meier Plotter for survival analysis, and expression was analyzed by GEPIA database. The core genes related to prognosis and highly expressed in hepatocellular carcinoma were analyzed by Metascape online tool for function and pathway enrichment analysis. Finally, the key genes were verified in hepatocellular carcinoma and paracancerous tissues.ResultsA total of 94 DEGs were screened from three gene chips GSE14520, GSE60502, and GSE102079, obtained from GEO. After the selected DEGs was analyzed by GO function analysis, KEGG pathway enrichment analysis, STRING and Cytscape software by DAVID, 19 core DEGs were screened. After 19 core DEGs were analyzed by Kaplan-Meier Plotter website, 9 genes [ribonucleotide reductase M2 (RRM2), polycomb repressive complex 1 (PRC1), topoisomerase Ⅱ alpha (TOP2A), aurora kinase A (AURKA), nucleolar spindle-associated protein 1 (NUSAP1), Rac-GTPase activating protein 1 (RACGAP1), abnormal spindle-like microcephaly-associated (ASPM), cyclin dependent kinase 1 (CDK1) and GINS complex subunit 1 (GINS1)] were found to be associated with the prognosis of hepatocellular carcinoma. The expressions of these 9 genes were analyzed by GEPIA, and the results showed that all 9 genes were highly expressed in hepatocellular carcinoma tissues. The functions and pathways of 9 highly expressed genes were analyzed by metascape website. Finally, RRM2 was selected for verification in hepatocellular carcinoma tissues and adjacent tissues, and it was found that the staining score of RRM2 in hepatocellular carcinoma tissues was (10.9±1.5) points, which was significantly higher than its staining score in adjacent tissues [(4.5±1.2) points], P<0.001.ConclusionThe nine genes identified by bioinformatics analysis may be the key genes in the occurrence and development of hepatocellular carcinoma, which can provide reference for further study on the pathogenesis, diagnosis and treatment of hepatocellular carcinoma.
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.
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.