Screening and identification of key genes between liver hepatocellular carcinoma (LIHC) and cholangiocarcinoma (CHOL) by bioinformatic analysis

Background:Liver hepatocellular carcinoma (LIHC) and cholangiocarcinoma (CHOL) are common primary liver cancers worldwide. Liver stem cells have biopotential to differentiate into either hepatocytes and cholangiocytes, the phenotypic overlap between LIHC and CHOL has been acceptable as a continuous liver cancer spectrum. However, few studies directly investigated the underlying molecular mechanisms between LIHC and CHOL.Method:To identify the candidate genes between LIHC and CHOL, three data series including {"type":"entrez-geo","attrs":{"text":"GSE31370","term_id":"31370"}}GSE31370, {"type":"entrez-geo","attrs":{"text":"GSE15765","term_id":"15765"}}GSE15765 and {"type":"entrez-geo","attrs":{"text":"GSE40367","term_id":"40367"}}GSE40367 were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified, and function enrichment analyses were performed. The protein-protein interaction network (PPI) was constructed and the module analysis was performed using STRING and Cytoscape.Results:A total of 171 DEGs were identified, consisting of 49 downregulated genes and 122 upregulated genes. Compared with CHOL, the enriched functions of the DEGs mainly included steroid metabolic process, acute inflammatory response, coagulation. Meanwhile, the pathway of KEGG enrichment analyses showed that the upregulated gene(s) were mainly enriched complement and coagulation cascades, cholesterol metabolism and PPAR signaling pathway, while the downregulated gene(s) were mainly enriched in ECM-receptor interaction, focal adhesion, bile secretion. Similarly, the most significant module was identified and biological process analysis revealed that these genes were mainly enriched in regulation of blood coagulation, acute inflammatory response, complement and coagulation cascades. Finally, two (ITIH2 and APOA2) of 10 hub genes had been screened out to help differential diagnosis.Conclusion:171 DEGs and two (ITIH2 and APOA2) of 10 hub genes identified in the present study help us understand the different molecular mechanisms between LIHC and CHOL, and provide candidate targets for differential diagnosis.     

Identification of candidate biomarkers correlated with the pathogenesis and prognosis of breast cancer via integrated bioinformatics analysis

Background:This study was carried out to identify potential key genes associated with the pathogenesis and prognosis of breast cancer (BC).Methods:Seven GEO datasets ({"type":"entrez-geo","attrs":{"text":"GSE24124","term_id":"24124"}}GSE24124, {"type":"entrez-geo","attrs":{"text":"GSE32641","term_id":"32641"}}GSE32641, {"type":"entrez-geo","attrs":{"text":"GSE36295","term_id":"36295"}}GSE36295, {"type":"entrez-geo","attrs":{"text":"GSE42568","term_id":"42568"}}GSE42568, {"type":"entrez-geo","attrs":{"text":"GSE53752","term_id":"53752"}}GSE53752, {"type":"entrez-geo","attrs":{"text":"GSE70947","term_id":"70947"}}GSE70947, {"type":"entrez-geo","attrs":{"text":"GSE109169","term_id":"109169"}}GSE109169) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between BC and normal breast tissue samples were screened by an integrated analysis of multiple gene expression profile datasets. Hub genes related to the pathogenesis and prognosis of BC were verified by employing protein–protein interaction (PPI) network.Results:Ten hub genes with high degree were identified, including CDK1, CDC20, CCNA2, CCNB1, CCNB2, BUB1, BUB1B, CDCA8, KIF11, and TOP2A. Lastly, the Kaplan–Meier plotter (KM plotter) online database demonstrated that higher expression levels of these genes were related to lower overall survival. Experimental validation showed that all 10 hub genes had the same expression trend as predicted.Conclusion:The findings of this research would provide some directive significance for further investigating the diagnostic and prognostic biomarkers to facilitate the molecular targeting therapy of BC, which could be used as a new biomarker for diagnosis and to guide the combination medicine of BC.     

ECM2 and GLT8D2 in human pulmonary artery hypertension: fruits from weighted gene co-expression network analysis

BackgroundPulmonary artery hypertension (PAH) is an incurable disease with a high mortality rate. Current medications ameliorate symptoms but cannot target adverse vascular remodeling. New therapeutic strategies for PAH need to be established.MethodsUsing the weighted gene coexpression network analysis (WGCNA) algorithm, we constructed a coexpression network of dataset {"type":"entrez-geo","attrs":{"text":"GSE117261","term_id":"117261"}}GSE117261 from the Gene Expression Omnibus (GEO) database. Key modules were identified by the relationship between module eigengenes and clinical traits. Hub genes were screened out based on gene significance (GS), module membership (MM), and mean pulmonary artery pressure (mPAP). External validations were conducted in {"type":"entrez-geo","attrs":{"text":"GSE48149","term_id":"48149"}}GSE48149 and {"type":"entrez-geo","attrs":{"text":"GSE113439","term_id":"113439"}}GSE113439. Functional enrichment and immune cell infiltration were analyzed using Metascape and CIBERSORTx.ResultsThe WGCNA analysis revealed 13 coexpression modules. The pink module had the highest correlation with PAH in terms of module eigengene (r=0.79; P=2e−18) and module significance (MS =0.43). Functional enrichment indicated genes in the pink module contributed to the immune system process and extracellular matrix (ECM). In the pink module, ECM2 (GS =0.65, MM =0.86, ρ=0.407, P=0.0019) and GLT8D2 (GS =0.63, MM =0.85, ρ=0.443, P=0.006) were identified as hub genes. For immune cells infiltration in PAH lung tissue, hub genes were positively correlated with M2 macrophages and resting mast cells, and were negatively correlated with monocytes, neutrophils, and CD4-naïve T cells.ConclusionsOur research identified 2 hub genes ECM2 and GLT8D2 related to PAH. The functions of these hub genes were involved in the immune process and ECM, indicating that they might serve as candidate therapeutic targets for PAH.     

Identification of potential core genes in esophageal carcinoma using bioinformatics analysis

Background:Esophageal squamous cell carcinoma (ESCC) is a common human malignancy worldwide. The tumorigenesis mechanism in ESCC is unclear.Materials and methods:To explore potential therapeutic targets for ESCC, we analyzed 3 microarray datasets ({"type":"entrez-geo","attrs":{"text":"GSE20347","term_id":"20347"}}GSE20347, {"type":"entrez-geo","attrs":{"text":"GSE38129","term_id":"38129"}}GSE38129, and {"type":"entrez-geo","attrs":{"text":"GSE67269","term_id":"67269"}}GSE67269) derived from the gene expression omnibus (GEO) database. Then, the GEO2R tool was used to screen out differently expressed genes (DEGs) between ESCC and normal tissue. Gene ontology function and kyoto encyclopedia of genes and genomes pathway enrichment analysis were performed using the database for annotation, visualization and integrated discovery to identify the pathways and functional annotation of DEGs. Protein–protein interaction of these DEGs was analyzed based on the search tool for the retrieval of interacting genes database and visualized by Cytoscape software. In addition, we used encyclopedia of RNA interactomes (ENCORI), gene expression profiling interactive analysis (GEPIA), and the human protein atlas to confirm the expression of hub genes in ESCC. Finally, GEPIA was used to evaluate the prognostic value of hub genes expression in ESCC patients and we estimated the associations between hub genes expression and immune cell populations (B Cell, CD8⁺ T Cell, CD4⁺ T Cell, Macrophage, Neutrophil, and Dendritic Cell) in esophageal carcinoma (ESCA) using tumor immune estimation resource (TIMER).Results:In this study, 707 DEGs (including 385 upregulated genes and 322 downregulated genes) and 6 hub genes (cyclin B1 [CCNB1], cyclin dependent kinase 1 [CDK1], aurora kinase A [AURKA], ubiquitin conjugating enzyme E2C [UBE2C], cyclin A2 [CCNA2], and cell division cycle 20 [CDC20]) were identified. All of the 6 hub genes were highly expressed in ESCC tissues. Among of them, only CCNB1 and CDC20 were associated with stage of ESCC and all of them were not associated with survival time of patients.Conclusion:DEGs and hub genes were confirmed in our study, providing a thorough, scientific and comprehensive research goals for the pathogenesis of ESCC.     

Integrated analysis of ceRNA network in hepatocellular carcinoma using bioinformatics analysis

Background:Long noncoding RNAs (lncRNAs) can work as microRNA (miRNA) sponges through a competitive endogenous RNA (ceRNA) mechanism. LncRNAs and miRNAs are important components of competitive endogenous binding, and their expression imbalance in hepatocellular carcinoma (HCC) is closely related to tumor development, diagnosis, and prognosis. This study explored the potential impact of the ceRNA regulatory network in HCC on the prognosis of HCC patients.Methods:We thoroughly researched the differential expression profiles of lncRNAs, miRNAs, and mRNAs from 2 HCC Gene Expression Omnibus datasets ({"type":"entrez-geo","attrs":{"text":"GSE98269","term_id":"98269"}}GSE98269 and {"type":"entrez-geo","attrs":{"text":"GSE60502","term_id":"60502"}}GSE60502). Then, a dysregulated ceRNA network was constructed by bioinformatics. In addition, hub genes in the ceRNA network were screened by Cytoscape, these hub genes functional analysis was performed by gene set enrichment analysis, and the expression of these hub genes in tumors and their correlation with patient prognosis were verified with Gene Expression Profiling Interactive Analysis.Results:A ceRNA network was successfully constructed in this study including 4 differentially expressed (DE) lncRNAs, 7 DEmiRNAs, and 166 DEmRNAs. Importantly, 4 core genes (CCNA2, CHEK1, FOXM1, and MCM2) that were significantly associated with HCC prognosis were identified.Conclusions:Our study provides comprehensive and meaningful insights into HCC tumorigenesis and the underlying molecular mechanisms of ceRNA. Furthermore, the specific ceRNAs can be further used as potential therapeutic targets and prognostic biomarkers for HCC.     

Identification of key genes and pathways in discoid lupus skin via bioinformatics analysis

Discoid lupus erythematosus (DLE) is the most common skin manifestation of lupus; however, the molecular mechanisms underlying DLE remain unknown. Therefore, we aimed to identify key differentially expressed genes (DEGs) in discoid lupus skin and investigate their potential pathways.To identify candidate genes involved in the occurrence and development of the disease, we downloaded the microarray datasets {"type":"entrez-geo","attrs":{"text":"GSE52471","term_id":"52471"}}GSE52471 and {"type":"entrez-geo","attrs":{"text":"GSE72535","term_id":"72535"}}GSE72535 from the Gene Expression Database (GEO). DEGs between discoid lupus skin and normal controls were selected using the GEO2R tool and Venn diagram software (http://bioinformatics.psb.ugent.be/webtools/Venn/). The Database for Annotation, Visualization, and Integrated Discovery (DAVID), Enrichr, and Cytoscape ClueGo were used to analyze the Kyoto Encyclopedia of Gene and Genome pathways and gene ontology. Protein-protein interactions (PPIs) of these DEGs were further assessed using the Search Tool for the Retrieval Interacting Genes version 10.0.Seventy three DEGs were co-expressed in both datasets. DEGs were predominantly upregulated in receptor signaling pathways of the immune response. In the PPI network, 69 upregulated genes were selected. Furthermore, 4 genes (CXCL10, ISG15, IFIH1, and IRF7) were found to be significantly upregulated in the RIG-I-like receptor signaling pathway, from analysis of Enrichr and Cytoscape ClueGo.The results of this study may provide new insights into the potential molecular mechanisms of DLE. However, further experimentation is required to confirm these findings.     

Bioinformatic analysis of the pathogenic mechanism of talaromyces marneffei infection

Background:Talaromyces marneffei (T marneffei), known as a significant pathogen in patients with AIDS in Southeast Asia, is a dimorphic fungus, which can cause deadly systematic infection in immunocompromised hosts. What is more, the dimorphic phase transition has been reported as a conspicuous process linked with virulence. Interestingly, the yeast form was found in infected individuals, representing the pathogenic phase. However, few researches were found to study the mechanism of dimorphic transition. Thus, a diverse insight into the dimorphic switch mechanism, is urgently needed and we are the first one to research the mechanism of dimorphism.Methods:Firstly, we investigated the microarray of T. marneffei in the Gene Expression Omnibus database (GEO) for differentially expressed genes (DEGs). Then Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 was employed to analyze the underlying enrichment and pathway in biological process of DEGs. Meanwhile, protein-protein interaction (PPI) network was constructed using STRING database. On the strength of the theory that similar amino acid sequences share similar structures, which play a decisive role on the function of protein, three dimensional structures of hub-genes were predicted to further investigate the likely function of hub-genes.Results:{"type":"entrez-geo","attrs":{"text":"GSE51109","term_id":"51109"}}GSE51109 was elected as the eligible series for the purpose of our research, including {"type":"entrez-geo","attrs":{"text":"GSM1238923","term_id":"1238923"}}GSM1238923 ({"type":"entrez-geo","attrs":{"text":"GSM23","term_id":"23"}}GSM23), {"type":"entrez-geo","attrs":{"text":"GSM1238924","term_id":"1238924"}}GSM1238924 ({"type":"entrez-geo","attrs":{"text":"GSM24","term_id":"24"}}GSM24), and {"type":"entrez-geo","attrs":{"text":"GSM1238925","term_id":"1238925"}}GSM1238925 ({"type":"entrez-geo","attrs":{"text":"GSM25","term_id":"25"}}GSM25). PMAA_012920, PMAA_028730, PMAA_068140, PMAA_092900, PMAA_032350 were the most remarkable genes in all of the three PPI networks, thus, were viewed as hub-genes. With regard to the three-dimensional construction, except that there was no significant prediction structure of PMAA_092900 with the criterion seq identify > 30%, GMQE: 0-1, QMEAN4: -4-0, the parallel templates for four structures were Crystal structure of Saccharomyces cerevesiae mitochondrial NADP(+)-dependent isocitrate dehydrogenase in complex with isocitrate, Organellar two-pore channels (TPCs), Yeast Isocitrate Dehydrogenase (Apo Form) and Crystal Structure Of ATP-Dependent Phosphoenolpyruvate Carboxykinase From Thermus thermophilus HB8 in order.Conclusion:The dimorphic transition of T. marneffei was viewed as a pathogenic factor and DEGs were observed. In-depth study of the function and pathway of DEGs revealed that PMAA_012920, PMAA_028730, PMAA_068140, PMAA_092900, PMAA_032350 were most likely acting as the hub-genes and were likely taking effect through regulating energy metabolism.     

Hub gene target of glioblastoma: LOX,SERPINH1 and TGFBI

Glioblastoma (GBM) is a malignant tumor. The long-term prognosis of the patients is poor. Therefore, it is of important clinical value to further explore the pathogenesis and look for molecular markers for early diagnosis and targeted treatment. Two expression profiling datasets [{"type":"entrez-geo","attrs":{"text":"GSE50161","term_id":"50161"}}GSE50161 ({"type":"entrez-geo","attrs":{"text":"GPL570","term_id":"570"}}GPL570 platform), {"type":"entrez-geo","attrs":{"text":"GSE116520","term_id":"116520"}}GSE116520 ({"type":"entrez-geo","attrs":{"text":"GPL10558","term_id":"10558"}}GPL10558 platform)] were respectively downloaded from the gene expression omnibus database. Volcano diagrams show the Differently expressed genes (DEGs) of {"type":"entrez-geo","attrs":{"text":"GSE50161","term_id":"50161"}}GSE50161 and {"type":"entrez-geo","attrs":{"text":"GSE116520","term_id":"116520"}}GSE116520. A Venn diagram revealed 467 common DEGs between the 2 datasets. Lysyl oxidase (LOX), serpin family H member 1 (SERPINH1) and transforming growth factor beta induced (TGFBI) were negatively correlated with the overall survival rate in patients with GBM. The hub genes are high in GBM tumor tissues. The relative expression levels of LOX, SERPINH1 and TGFBI were significantly higher in GBM samples, compared with the normal brain tissues groups. Bioinformatics technology could be a useful tool to predict progression of GBM and to explore the mechanism of GBM.LOX, SERPINH1 and TGFBI may be involved in the mechanism of the occurrence and development of GBM, and may be used as molecular targets for early diagnosis and specific treatment.DEGs identified using GEO2R. Functional annotation of DEGs using Kyoto Encyclopedia of Genes and Genomes and gene body pathway enrichment analysis. Construction of a protein-protein interaction network. The pathway and process enrichment analysis of the hub genes were performed by Metascape. Survival analysis was performed in gene expression profiling interactive analysis. Real-time fluorescent quantitative polymerase chain reaction assay was performed to verify. The animal model was established for western blot test analysis.     

Bioinformatics analysis of differentially expressed genes in tumor and paracancerous tissues of patients with lung adenocarcinoma

BackgroundLung adenocarcinoma is the main pathological type of non-small cell lung cancer (NSCLC). In this study, we analyzed the gene expression profile of lung adenocarcinoma tumor and paracancerous tissues by bioinformatics to assess the genes and signal pathways related to lung adenocarcinoma.MethodsThe expression data of {"type":"entrez-geo","attrs":{"text":"GSE7670","term_id":"7670"}}GSE7670, {"type":"entrez-geo","attrs":{"text":"GSE27262","term_id":"27262"}}GSE27262, and {"type":"entrez-geo","attrs":{"text":"GSE32863","term_id":"32863"}}GSE32863 were downloaded from the Gene Expression Omnibus (GEO) database. The three microarray data sets were integrated to obtain common differential expression genes of lung adenocarcinoma tumor and adjacent tissues. The STRING database was used to construct the protein-protein interaction (PPI) network of lung adenocarcinoma and mine the gene modules and core genes in the network, and the online tools, GEPIA and Kaplan-Meier plotter were used to further verify and analyze the core genes.ResultsThere were 109 pairs of lung adenocarcinoma tissues and matched paracancerous normal lung tissues in the three data sets. Eighty-three differentially expressed genes were identified, including 16 up-regulated and 67 down-regulated genes, and 60 differentially expressed genes were successfully incorporated into the PPI network complex. Eleven core genes were identified in the PPI network complex, including three up-regulated (COMP, SPP1, COL1A1) and eight down-regulated genes (CDH5, CAV1, CLDN5, LYVE1, IL6, VWF, TEK, PECAM1). These core genes were verified by the GEPIA tumor database. Survival analysis showed that expression of the core genes was significantly related to the prognosis of lung adenocarcinoma. KEGG pathway analysis of core genes showed six genes (COMP, SPP1, COL1A1, IL6, VWF, TEK) were significantly enriched in the PI3K-Akt signaling-pathway (P=1.62E-06).ConclusionsBy analyzing the differential expression genes of lung adenocarcinoma and paracancerous normal tissues with bioinformatics, 11 genes with significant differential expression and significant influence on prognosis were identified. The findings may provide new concepts for developing diagnosis and treatment targets and prognosis markers for lung adenocarcinoma.     

Comparison of critical biomarkers in 2 erectile dysfunction models based on GEO and NOS-cGMP-PDE5 pathway

Background:Erectile dysfunction is a disease commonly caused by diabetes mellitus (DMED) and cavernous nerve injury (CNIED). Bioinformatics analyses including differentially expressed genes (DEGs), enriched functions and pathways (EFPs), and protein-protein interaction (PPI) networks were carried out in DMED and CNIED rats in this study. The critical biomarkers that may intervene in nitric oxide synthase (NOS, predominantly nNOS, ancillary eNOS, and iNOS)-cyclic guanosine monophosphate (cGMP)-phosphodiesterase 5 enzyme (PDE5) pathway, an important mechanism in erectile dysfunction treatment, were then explored for potential clinical applications.Methods:{"type":"entrez-geo","attrs":{"text":"GSE2457","term_id":"2457"}}GSE2457 and {"type":"entrez-geo","attrs":{"text":"GSE31247","term_id":"31247"}}GSE31247 were downloaded. Their DEGs with a |logFC (fold change)| > 0 were screened out. Database for Annotation, Visualization and Integrated Discovery (DAVID) online database was used to analyze the EFPs in Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes networks based on down-regulated and up-regulated DEGs respectively. PPI analysis of 2 datasets was performed in Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and Cytoscape. Interactions with an average score greater than 0.9 were chosen as the cutoff for statistical significance.Results:From a total of 1710 DEGs in {"type":"entrez-geo","attrs":{"text":"GSE2457","term_id":"2457"}}GSE2457, 772 were down-regulated and 938 were up-regulated, in contrast to the 836 DEGs in {"type":"entrez-geo","attrs":{"text":"GSE31247","term_id":"31247"}}GSE31247, from which 508 were down-regulated and 328 were up-regulated. The 25 common EFPs such as aging and response to hormone were identified in both models. PPI results showed that the first 10 hub genes in DMED were all different from those in CNIED.Conclusions:The intervention of iNOS with the hub gene complement component 3 in DMED and the aging process in both DMED and CNIED deserves attention.     

Screening of sphingolipid metabolism-related genes associated with immune cells in myocardial infarction: a bioinformatics analysis

BackgroundInflammation and immune cell infiltration in infarcted myocardial tissue are critical to myocardial infarction (MI) prognosis, and alterations in sphingolipid metabolism (SM) have been shown to potentially influence the inflammatory response and induce cardioprotection, but the underlying mechanisms are unclear. We therefore performed bioinformatics analysis to screen for key genes of SM in MI immune cells.MethodsThree matrix files including {"type":"entrez-geo","attrs":{"text":"GSE61145","term_id":"61145"}}GSE61145, {"type":"entrez-geo","attrs":{"text":"GSE23294","term_id":"23294"}}GSE23294, and {"type":"entrez-geo","attrs":{"text":"GSE71906","term_id":"71906"}}GSE71906 were downloaded from the Gene Expression Omnibus (GEO) database. {"type":"entrez-geo","attrs":{"text":"GSE61145","term_id":"61145"}}GSE61145 was a human peripheral blood database, and {"type":"entrez-geo","attrs":{"text":"GSE23294","term_id":"23294"}}GSE23294 and {"type":"entrez-geo","attrs":{"text":"GSE71906","term_id":"71906"}}GSE71906 were 2 mouse myocardial tissue databases. R and annotation packages were used to screen for differentially-expressed genes (DEGs). Datasets of human and mouse cardiac tissues were downloaded from the GEO database for subsequent validation. The downloaded platform and matrix files were processed using R language and annotation packages. Key targets and enrichment pathways were identified using Gene Ontology (GO) term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The Wilcoxon test was performed on the genes involved in SM pathways in neutrophils.ResultsA total of 261 DEGs were obtained from human peripheral blood datasets, among which 101 were immune-related. GO analysis revealed that neutrophil activation, T cell activation, and T cell differentiation were significantly enriched in the immune-related DEGs. Three types of immune cells were identified in infarcted myocardial tissues. In addition, 194 DEGs were obtained from mouse myocardial tissue data, among which 6 SM-related genes (Asah1, Degs1, Neu1, Sptlc2, Sphk1, and Gba2) were significantly associated with MI. Evaluation of the relationships between these DEGs and neutrophils showed that the expression of the Sptlc2 gene was significantly upregulated in neutrophils of the MI group, while the expression levels of the Asah1 and Degs1 genes were downregulated.ConclusionsWe identified 3 SM-related genes that were highly associated with neutrophils in MI, which may advance our understanding of SM in immune cells after MI.     

Identification of endometriosis-associated genes and pathways based on bioinformatic analysis

Endometriosis is associated with dysmenorrhea, chronic pelvic pain, and infertility. The specific mechanism of endometriosis remains unclear. The aim of this study was to apply a bioinformatics approach to reveal related pathways or genes involved in the development of endometriosis.The gene expression profiles of {"type":"entrez-geo","attrs":{"text":"GSE25628","term_id":"25628"}}GSE25628, {"type":"entrez-geo","attrs":{"text":"GSE5108","term_id":"5108"}}GSE5108, and {"type":"entrez-geo","attrs":{"text":"GSE7305","term_id":"7305"}}GSE7305 were downloaded from the gene expression omnibus (GEO) database. Differentially expressed gene (DEG) analysis was performed using GEO2R. The database for annotation, visualization, and integrated discovery (DAVID) was utilized to analyze the functional enrichment, gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) pathway of the differentially expressed genes. A protein-protein interaction (PPI) network was constructed and module analysis was performed using search tool for the retrieval of interacting genes and cytoscape.A total of 119 common differentially expressed genes were extracted, consisting of 51 downregulated genes and 68 upregulated genes. The enriched functions and pathways of the DEGs and hub genes include DNA strand separation, cellular proliferation, degradation of the extracellular matrix, encoding of smooth muscle myosin as a major contractile protein, exiting the proliferative cycle and entering quiescence, growth regulation, and implication in a wide variety of biological processes.A bioinformatics approach combined with cell experiments in this study revealed that identifying DEGs and hub genes leads to better understanding of the molecular mechanisms underlying the progression of endometriosis, and efficient biomarkers underlying this pathway need to be further investigated.     

Identifying key genes and small molecule compounds for nasopharyngeal carcinoma by various bioinformatic analysis

Nasopharyngeal carcinoma (NPC) is one of the most prevalent head and neck cancer in southeast Asia. It is necessary to proceed further studies on the mechanism of occurrence and development of NPC.In this study, we employed the microarray dataset {"type":"entrez-geo","attrs":{"text":"GSE12452","term_id":"12452"}}GSE12452 and {"type":"entrez-geo","attrs":{"text":"GSE53819","term_id":"53819"}}GSE53819 including 28 normal samples and 49 nasopharyngeal carcinoma samples downloaded from the Gene Expression Omnibus(GEO) to analysis. R software, STRING, CMap, and various databases were used to screen differentially expressed genes (DEGs), construct the protein–protein interaction (PPI) network, and proceed small molecule compounds analysis, among others.Totally, 424 DEGs were selected from the dataset. DEGs were mainly enriched in extracellular matrix organization, cilium organization, PI3K-Akt signaling pathway, collagen-containing extracellular matrix, and extracellular matrix-receptor interaction, among others. Top 10 upregulated and top 10 downregulated hub genes were identified as hub DEGs. Piperlongumine, apigenin, menadione, 1,4-chrysenequinone, and chrysin were identified as potential drugs to prevent and treat NPC. Besides, the effect of genes CDK1, CDC45, RSPH4A, and ZMYND10 on survival of NPC was validated in GEPIA database.The data revealed novel aberrantly expressed genes and pathways in NPC by bioinformatics analysis, potentially providing novel insights for the molecular mechanisms governing NPC progression. Although further studies needed, the results demonstrated that the expression levels of CDK1, CDC45, RSPH4A, and ZMYND10 probably affected survival of NPC patients.     

Analysis of differentially expressed genes in individuals with noninfectious uveitis based on data in the gene expression omnibus database

Noninfectious uveitis (NIU), an intraocular inflammation caused by immune-mediated reactions to eye antigens, is associated with systemic rheumatism and several autoimmune diseases. However, the mechanisms underlying the pathogenesis of uveitis are poorly understood. Therefore, we aimed to identify differentially expressed genes (DEGs) in individuals with NIU and to explore its etiologies using bioinformatics tools.{"type":"entrez-geo","attrs":{"text":"GSE66936","term_id":"66936"}}GSE66936 and {"type":"entrez-geo","attrs":{"text":"GSE18781","term_id":"18781"}}GSE18781 datasets from the gene expression omnibus (GEO) database were merged and analyzed. Functional enrichment analysis was performed, and protein-protein interaction (PPI) networks were constructed.A total of 89 DEGs were identified. Gene ontology (GO) enrichment analysis identified 21 enriched gene sets. Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis identified four core enriched pathways: antigen processing and expression signaling, natural killer (NK) cell-mediated cytotoxicity signaling, glutathione metabolic signal transduction, and arachidonic acid metabolism pathways. PPI network analysis revealed an active component-target network with 40 nodes and 132 edges, as well as several hub genes, including CD27, LTF, NCR3, SLC4A1, CD69, KLRB1, KIR2DL3, KIR3DL1, and GZMK.The eight potential hub genes may be associated with the risk of developing NIU. NK cell-mediated cytotoxicity signaling might be the key molecular mechanism in the occurrence and development of NIU. Our study provided new insights on NIU, its genetics, molecular pathogenesis and new therapeutic targets.     

Identification of key immune-related genes and immune infiltration in atrial fibrillation with valvular heart disease based on bioinformatics analysis

BackgroundAtrial fibrillation (AF) is the most common persistent arrhythmia. Valvular heart disease (VHD) and AF frequently coexist. In our study, from performing bioinformatics analysis, we sought to identify immune-related genes (IRGs) and explore the role of immune cell infiltration in AF-VHD in depth, aiming at investigating the potential molecular mechanism and developing new therapeutic targets for AF, including AF-VHD.MethodsThe gene expression of the {"type":"entrez-geo","attrs":{"text":"GSE41177","term_id":"41177"}}GSE41177 and {"type":"entrez-geo","attrs":{"text":"GSE79768","term_id":"79768"}}GSE79768 datasets were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were analyzed via the limma package in Bioconductor with R software. Differentially expressed immune-related genes (DEIRGs) were selected via combination ImmPort database with DEGs, and the enrichment function and pathway analysis were explored. A protein-protein interaction (PPI) network was built with a Search Tool for the Retrieval of Interacting Genes/Proteins plugin in Cytoscape. The CIBERSORT algorithm was used to evaluate immune infiltration in the left atrial (LA) tissues between AF-VHD and sinus rhythm (SR) patients. Finally, a correlation analysis between key DEIRGs and infiltrating immune cells was performed.ResultsA total of 130 DEIRGs were detected. Enrichment function of DEIRGs demonstrated that they are significant in immune and inflammatory responses. The key DEIRGs assessed by the PPI network and involved in both the immune and inflammatory responses were the C-X-C motif chemokine ligand (CXCL) 1, pro-platelet basic protein (PPBP), CXCL12, and C-C motif chemokine ligand 4 (CCL4). The immune infiltration findings indicated that, compared with the LA tissues from SR patients, the tissues from AF-VHD patients contained a higher proportion of gamma delta T cells, but a lower proportion of CD8 and regulatory T cells. The results of correlation analysis demonstrated that CXCL1 was positively correlated with activated mast cells and significantly negatively correlated with resting mast cells. PPBP, CXCL12, and CCL4 were positively correlated with the infiltration of various immune cells, such as neutrophils, plasma cells, and resting dendritic cells.ConclusionsThe key immune-related genes and the differences in immune infiltration in LA tissues play an essential role in the occurrence and progression of AF-VHD.     

Exploring new targets for the treatment of hepatitis-B virus and hepatitis-B virus-associated hepatocellular carcinoma: A new perspective in bioinformatics

Background:Hepatitis B Virus (HBV) infection is a global public health problem. After infection, patients experience a natural course from chronic hepatitis to cirrhosis and even Hepatitis B associated Hepatocellular Carcinoma (HBV-HCC). With the multi-omics research, many differentially expressed genes from chronic hepatitis to HCC stages have been discovered. All these provide important clues for new biomarkers and therapeutic targets. The purpose of this study is to explore the differential gene expression of HBV and HBV-related liver cancer, and analyze their enrichments and significance of related pathways.Methods:In this study, we downloaded four microarray datasets {"type":"entrez-geo","attrs":{"text":"GSE121248","term_id":"121248"}}GSE121248, {"type":"entrez-geo","attrs":{"text":"GSE67764","term_id":"67764"}}GSE67764, {"type":"entrez-geo","attrs":{"text":"GSE55092","term_id":"55092"}}GSE55092, {"type":"entrez-geo","attrs":{"text":"GSE55092","term_id":"55092"}}GSE55092 and {"type":"entrez-geo","attrs":{"text":"GSE83148","term_id":"83148"}}GSE83148 from the Gene Expression Omnibus (GEO) database. Using these four datasets, patients with chronic hepatitis B (CHB) differentially expressed genes (CHB DEGs) and patients with HBV-related HCC differentially expressed genes (HBV-HCC DEGs) were identified. Then Protein–protein Interaction (PPI) network analysis, Gene Ontology (GO) functional analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to excavate the functional interaction of these two groups of DEGs and the common DEGs. Finally, the Kaplan website was used to analyze the role of these genes in HCC prognostic.Results:A total of 241 CHB DEGs, 276 HBV-HCC DEGs, and 4 common DEGs (cytochrome P450 family 26 subfamily A member 1 (CYP26A1), family with sequence similarity 110 member C(FAM110C), SET and MYND domain containing 3(SMYD3) and zymogen granule protein 16(ZG16)) were identified. CYP26A1, FAM110C, SMYD3 and ZG16 exist in 4 models and interact with 33 genes in the PPI network of CHB and HBV-HCC DEGs,. GO function analysis showed that: CYP26A1, FAM110C, SMYD3, ZG16, and the 33 genes in their models mainly affect the regulation of synaptic vesicle transport, tangential migration from the subventricular zone to the olfactory bulb, cellular response to manganese ion, protein localization to mitochondrion, cellular response to dopamine, negative regulation of neuron death in the biological process of CHB. In the biological process of HBV-HCC, they mainly affect tryptophan catabolic process, ethanol oxidation, drug metabolic process, tryptophan catabolic process to kynurenine, xenobiotic metabolic process, retinoic acid metabolic process, steroid metabolic process, retinoid metabolic process, steroid catabolic process, retinal metabolic process, and rogen metabolic process. The analysis of the 4 common DEGs related to the prognosis of liver cancer showed that: CYP26A1, FAM110C, SMYD3 and ZG16 are closely related to the development of liver cancer and patient survival. Besides, further investigation of the research status of the four genes showed that CYP26A1 and SMYD3 could also affect HBV replication and the prognosis of liver cancer.Conclusion:CYP26A1, FAM110C, SMYD3 and ZG16 are unique genes to differentiate HBV infection and HBV-related HCC, and expected to be novel targets for HBV-related HCC occurrence and prognostic judgement.