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 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.     

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.     

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.     

Network pharmacology-based identification of miRNA expression of Astragalus membranaceus in the treatment of diabetic nephropathy

Diabetic nephropathy (DN) is a common microvascular complication of diabetic patients, along with hypertension, hyperlipemia, proteinuria, edema, and other clinical manifestations. Astragalus membranaceus (AM) is a traditional Chinese medicine and has shown significant clinical efficacy against DN. However, the overall molecular mechanism of this therapeutic effect has not been entirely elucidated. Using network pharmacology, we aimed to identify the key active ingredients and potential pharmacological mechanisms of AM in treating DN and provide scientific evidence of its clinical efficacy.The active ingredients of AM were obtained from the traditional Chinese medicine systems pharmacology database, and the potential targets of AM were identified using the therapeutic target database. DN-related target genes were acquired from the Gene Expression Omnibus microarray dataset {"type":"entrez-geo","attrs":{"text":"GSE1009","term_id":"1009"}}GSE1009 and 3 widely used databases-DisGeNET, GeneCards, and Comparative Toxicogenomics Database. The DN–AM common target protein interaction network was established by using the STRING database. Active ingredients candidate targets proteins networks were constructed using Cytoscape software for visualization. Additionally, gene ontology (GO) and Kyoto encyclopedia of genes and genomes pathway analyses were performed using the Database for Annotation, Visualization, and Integrated Discovery database. Target-regulating microRNAs (miRNAs) of these hub genes were obtained from the therapeutic target database, which could then be used for further identification of AM-regulated key miRNAs.A total of 17 active ingredients and 214 target proteins were screened from AM. 61 candidate co-expressed genes with therapeutic effects against DN were obtained and considered as potential therapeutic targets. GO and Kyoto encyclopedia of genes and genomes enrichment analysis showed that these genes were mainly involved in inflammatory response, angiogenesis, oxidative stress reaction, HIF signaling pathway, tumor necrosis factor signaling pathway, and VEGF signaling pathway. In all, 636 differentially expressed genes were identified between the DN patients and control group by using microarray data, {"type":"entrez-geo","attrs":{"text":"GSE1009","term_id":"1009"}}GSE1009. Lastly, VEGFA, epidermal growth factor receptor, STAT1, and GJA1 were screened as hub genes. The relationships between miRNAs and hub genes were constructed, which showed that miR-302-3p, miR-372-3p, miR-373-3p, and miR-520-3p were regulated by VEGFA and epidermal growth factor receptor. Meanwhile, VEGFA also influenced miR-15-5p, miR-16-5p, miR-17-5p, miR-20-5p, miR-93-5p, miR-106-5p, miR-195-5p, miR-424-5p, miR-497-5p, and miR-519-3p. In addition, miR-1-3p and miR-206 were regulated by VEGFA and GJA1, and miR-23-3p was regulated by STAT1 and GJA1.To our knowledge, this study revealed for the first time the characteristic multiple components, multiple targets, and multiple pathways of AM that seem to be the underlying mechanisms of action of AM in the treatment of DN with respect to miRNAs.Private information from individuals will not be published. This systematic review also does not involve endangering participant rights. Ethical approval will not be required. The results may be published in a peer-reviewed journal or disseminated at relevant conferences.     

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.     

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.     

Identification of key genes and microRNAs for multiple sclerosis using bioinformatics analysis

To better understand the molecular mechanism underlying the pathogenesis of multiple sclerosis (MS), we aimed to identify the key genes and microRNAs (miRNA) associated with MS and analyze their interactions. Differentially expressed genes (DEGs) and miRNAs (DEMs) based on the gene miRNA dataset {"type":"entrez-geo","attrs":{"text":"GSE17846","term_id":"17846"}}GSE17846 and mRNA dataset {"type":"entrez-geo","attrs":{"text":"GSE21942","term_id":"21942"}}GSE21942 were determined using R software. Next, we performed functional enrichment analysis and constructed a protein–protein interaction network. Data validation was performed to ensure the reliability of hub genes. The miRNA-mRNA regulatory network was constructed. In total, 47 DEMs and 843 DEGs were identified. Protein–protein interaction network analysis identified several hub genes, including JUN, FPR2, AKT1, POLR2L, LYZ, CXCL8, HBB, CST3, CTSZ, and MMP9, especially LYZ and CXCL8. We constructed an miRNA-mRNA regulatory network and found that hsa-miR-142-3p, hsa-miR-107, hsa-miR-140-5p, and hsa-miR-613 were the most important miRNAs. This study reveals some key genes and miRNAs that may be involved in the pathogenesis of MS, providing potential targets for the diagnosis and treatment of MS.     

CD74, a novel predictor for bronchopulmonary dysplasia in preterm infants

Bronchopulmonary dysplasia (BPD) remains a major complication and accounts for high morbidity and mortality of preterm infants. The present study aimed to identify the key genes in the development of BPD and to provide some new insights into the pathogenesis of BPD. The {"type":"entrez-geo","attrs":{"text":"GSE108754","term_id":"108754"}}GSE108754 dataset was downloaded from Gene Expression Omnibus database containing 5 samples of BPD patients and 6 of non-BPD infants. The differentially expressed genes (DEGs) between BPD and non-BPD patients were identified by R software. The pathway and function enrichment analyses were performed through Database for Annotation Visualization and Integrated Discovery website. The protein-protein interaction network for DEGs was established by Cytoscape software and the most highly connected module was selected through MCODE plugin. Furthermore, the clinical sample verification among 25 BPD patients and 10 non-BPD infants was carried out in our center. Finally, based on the results above, the gene set enrichment analysis focusing on CD74 upregulated status was employed. Totally, 189 DEGs including 147 upregulated genes and 42 downregulated genes between BPD and non-BPD patients were screened out. The pathway and function enrichments revealed these DEGs were mainly enriched in asthma, intestinal immune network for IgA production, antigen processing and presentation and immune response. Thirteen DEGs (CD74, HLA-DMA, HLA-DRA, HLA-DMB, HLA-DOB, HLA-DQA1, HLA-DRB5, HLA-DPA1, HLA-DOA, HLA-DPB1, HLA-DQB2, HLA-DQA2, and HLA-DQB1) were determined as hub genes. The mRNA expression levels of the 13 hub genes were tested by quantitative real-time polymerase chain reaction among our clinical samples. Eventually, CD74 was confirmed to be the most significant highly expressed in BPD samples (P < .001) and its expression level was negatively correlated with gestational age (r = –0.653) and birth weight (r = –0.675). The gene set enrichment analysis results showed the gene sets associated with lupus erythematosus, viral myocarditis, immune network for IgA production, graft versus host disease, cell adhesion molecules and so no were differentially enriched with the phenotype of high-expression CD74. In conclusion, CD74 may serve to predict the BPD development and provide a new therapeutic target for BPD.     

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.     

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.     

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.     

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.     

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.