Key modules and hub genes identified by coexpression network analysis for revealing novel biomarkers for larynx squamous cell carcinoma

Larynx squamous cell carcinoma (LSCC) is the second most aggressive head and neck squamous cell carcinoma. Numerous genes have been identified to be aberrantly expressed during the development of LSCC. However, currently, researchers focus more on the individual molecule and downstream genes, leaving the coexpression among genes and key upstream disease driver genes unexploited. In this study, we applied weighted gene coexpression analysis (WGCNA) to decipher potential hub genes driving the development of LSCC. After downloading of LSCC microarray profile from gene expression omnibus, different expression analysis was performed, which was used to conduct functional enrichment analysis. Then, we applied WGCNA to highlight the hub genes which were relevant to the carcinogenesis and progression. A total of 2858 differentially expressed genes were identified in LSCC samples compared with adjacent non-neoplastic tissues. WGCNA revealed three LSCC set-specific modules having significant Kyoto Encyclopedia of Genes and Genomes enrichment effect, including pink, cyan, and black module. Nine hub genes were identified to be crucial in LSCC onset and progression, which may assist clinical decisions and serve as potential targets for LSCC treatment.     

Identification of biomarkers correlated with hypertrophic cardiomyopathy with co-expression analysis

Hypertrophic cardiomyopathy (HCM) is reported to be the most common genetic heart disease. To identify key module and candidate biomarkers correlated with clinical prognosis of patients with HCM, we carried out this study with co-expression analysis. To construct a co-expression network of hub genes correlated with HCM, the Weighted Gene Co-expression Network Analysis (WGCNA) was performed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed by Database for Annotation, Visualization and Integrated Discovery (DAVID). The protein-protein interaction network analysis of central genes was performed to recognize the interactions of central genes. Gene set enrichment analyses were carried out to discover the possible mechanisms involved in the pathways promoted by hub genes. To validate the hub genes, quantitative real-time polymerase chain reaction (RT-PCR) was performed. Based on the results of topological overlap measure based clustering, 2,351 differentially expressed genes (DEGs) were identified. Those genes were included in six different modules. Of these modules, the yellow and the blue modules showed a pivotal correlation with HCM. DEGs were enriched in immune system procedure associated GO terms and KEGG pathways. We identified nine hub genes (TYROBP, STAT3, CSF1R, ITGAM, SYK, ITGB2, LILRB2, LYN, and HCK) affected the immune system significantly. Among the genes we validated with RT-PCR, TYROBP, CSF1R, and SYK showed significant increasing expression levels in model HCM rats. In conclusion, we identified two modules and nine hub genes, which were prominently associated with HCM. We found that immune system may play a crucial role in the HCM. Accordingly, those genes and pathways might become therapeutic targets with clinical usefulness in the future.     

Exploring the mechanism of ellagic acid against gastric cancer based on bioinformatics analysis and network pharmacology

Ellagic acid (EA) is a natural polyphenolic compound. Recent studies have shown that EA has potential anticancer properties against gastric cancer (GC). This study aims to reveal the potential targets and mechanisms of EA against GC. This study adopted methods of bioinformatics analysis and network pharmacology, including the weighted gene co-expression network analysis (WGCNA), construction of protein–protein interaction (PPI) network, receiver operating characteristic (ROC) and Kaplan–Meier (KM) survival curve analysis, Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, molecular docking and molecular dynamics simulations (MDS). A total of 540 EA targets were obtained. Through WGCNA, we obtained a total of 2914 GC clinical module genes, combined with the disease database for screening, a total of 606 GC-related targets and 79 intersection targets of EA and GC were obtained by constructing Venn diagram. PPI network was constructed to identify 14 core candidate targets; TP53, JUN, CASP3, HSP90AA1, VEGFA, HRAS, CDH1, MAPK3, CDKN1A, SRC, CYCS, BCL2L1 and CDK4 were identified as the key targets of EA regulation of GC by ROC and KM curve analysis. The enrichment analysis of GO and KEGG pathways of key targets was performed, and they were mainly enriched in p53 signalling pathway, PI3K-Akt signalling pathway. The results of molecular docking and MDS showed that EA could effectively bind to 13 key targets to form stable protein–ligand complexes. This study revealed the key targets and molecular mechanisms of EA against GC and provided a theoretical basis for further study of the pharmacological mechanism of EA against GC.     

Differential coexpression networks in bronchiolitis and emphysema phenotypes reveal heterogeneous mechanisms of chronic obstructive pulmonary disease

Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease with multiple molecular mechanisms. To investigate and contrast the molecular processes differing between bronchiolitis and emphysema phenotypes of COPD, we downloaded the GSE69818 microarray data set from the Gene Expression Omnibus (GEO), which based on lung tissues from 38 patients with emphysema and 32 patients with bronchiolitis. Then, weighted gene coexpression network analysis (WGCNA) and differential coexpression (DiffCoEx) analysis were performed, followed by gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes enrichment analysis (KEGG) analysis. Modules and hub genes for bronchiolitis and emphysema were identified, and we found that genes in modules linked to neutrophil degranulation, Rho protein signal transduction and B cell receptor signalling were coexpressed in emphysema. DiffCoEx analysis showed that four hub genes (IFT88, CCDC103, MMP10 and Bik) were consistently expressed in emphysema patients; these hub genes were enriched, respectively, for functions of cilium assembly and movement, proteolysis and apoptotic mitochondrial changes. In our re‐analysis of GSE69818, gene expression networks in relation to emphysema deepen insights into the molecular mechanism of COPD and also identify some promising therapeutic targets.     

伴HBV感染性肝癌调控枢纽基因筛选与初步鉴定

慢性乙型肝炎病毒(Hepatitis B virus,HBV)感染引起的原发性肝癌涉及多种基因、转录本和蛋白质的相互作用及调控。从单个基因的角度来看,某个基因的表达量的改变只能对肝癌发生发展的局部作出解释而无法从整体行为进行深入和全面的探索,无法满足高度复杂性的调控研究需要。筛选乙肝相关性肝癌的基因芯片数据获取差异表达基因后,应用加权基因共表达网络分析算法构建基因共表达网络,识别与肝癌发生相关的模块,利用可视化筛选枢纽基因,并针对枢纽基因进行基因本体富集分析和初步验证。富集分析和文献挖掘一致发现,某些枢纽基因确实与多种癌症的发生与发展存在显著的关联。权重基因共表达网络分析方法被证明是一个高效的系统生物学方法,应用该方法发现了新的HBV相关性肝癌枢纽基因。经实验验证,发现枢纽基因SHARPIN促进细胞迁移。该研究对肝癌发生的调控机制以及发现HBV慢性感染导致肝癌的新型诊断标志物和(或)药物作用靶点提供了新的视野。     

LncRNAs related key pathways and genes in ischemic stroke by weighted gene co-expression network analysis (WGCNA)

BackgroundIschemic stroke (IS) was a significant public health concern and long-chain noncoding RNAs (lncRNAs) were gaining particular importance in stroke biology, however, the potential mechanism of lncRNAs in IS was not fully understood.MethodsIn this study, three diagnosed patients with IS and three controls were selected to establish the lncRNA library. Weighted gene co-expression network analysis (WGCNA) was applied to screen key lncRNA modules associated with IS. The key lncRNAs were identified by module membership (MM) and gene significance (GS). The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was used to identify the key pathways and protein-protein interaction (PPI) network method was used to identify the key genes.ResultsA total of 3627 lncRNAs were investigated, followed by an analysis of 17 modules, and only one module was highly associated with the IS. The top 10 lncRNAs were identified based on GS and MM. KEGG pathways analysis revealed the top two pathways of the Human T cell Lymphotropic Virus-1 (HTLV-1) infection and the mTOR signaling pathway might influence the progress of IS. Further, genes meeting the top two degree (AKT1 and MAPK14) were selected as the hub genes in the PPI network.ConclusionTo summarize, this study identified the key pathways and genes, which might serve as biomarkers and targets for precise diagnosis and treatment of IS in the future.     

IRF2-ferroptosis related gene is associated with prognosis and EMT in gliomas

Ferroptosis is a new type of programmed cell death that has excellent anti-tumor potential in different tumors. However, the research on ferroptosis in glioma is still incomplete. In this study, we aimed to revealed the relationship between ferroptosis-related genes (FRGs) and glioma. We collected gene expression profiles of glioma patients from the TCGA and CGGA databases. All glioma samples were classified into five subtypes using the R software ConsensusClusterPlus. Subsequently, we performed single sample gene set enrichment analysis (ssGSEA) to explore the correlation between different subtypes and immune status and ferroptosis. Then co-expression modules were constructed via weighted gene co-expression network analysis (WGCNA). A Gene Ontology (GO) analysis was conducted to analyze the potential biological functions of the genes in the modules. Finally, we identified 10 hub genes using the PPI network. The in vitro experiments were used to validate our predictions. We found that the expression level of IRF2 is positively correlated with the grade of glioma. The overexpression of IRF2 could protect glioma cells from ferroptosis and enhance the invasive and migratory abilities. Silence of IRF2 had the opposite effect. In conclusion, we demonstrated a novel ferroptosis-related signature for predicting prognosis, and IRF2 could be a potential biomarker for diagnosis and treatment in glioma.     

Identification of gene coexpression modules,hub genes,and pathways related to spinal cord injury using integrated bioinformatics methods

Spinal cord injury (SCI) is characterized by dramatic neurons loss and axonal regeneration suppression. The underlying mechanism associated with SCI-induced immune suppression is still unclear. Weighted gene coexpression network analysis (WGCNA) is now widely applied for the identification of the coexpressed modules, hub genes, and pathways associated with clinic traits of diseases. We performed this study to identify hub genes associated with SCI development. Gene Expression Omnibus (GEO) data sets GSE45006 and GSE20907 were downloaded and the significant correlativity and connectivity between them were detected using WGCNA. Three significant consensus modules, including 567 eigengenes, were identified from the master GSE45006 data following the preconditions of approximate scale-free topology for WGCNA. Further bioinformatics analysis showed these eigengenes were involved in inflammatory and immune responses in SCI. Three hub genes Rac2, Itgb2, and Tyrobp and one pathway “natural killer cell-mediated cytotoxicity” were identified following short time-series expression miner, protein-protein interaction network, and functional enrichment analysis. Gradually upregulated expression patterns of Rac2, Itgb2, and Tyrobp genes at 0, 3, 7, and 14 days after SCI were confirmed based on GSE45006 and GSE20907 data set. Finally, we found that Rac2, Itgb2, and Tyrobp genes might take crucial roles in SCI development through the “natural killer cell–mediated cytotoxicity” pathway.     

基于网络药理学探讨“黄芪-白术-熟地黄”组方防治肾病综合征的作用机制

本研究旨在通过网络药理学方法和分子对接技术探讨黄芪-白术-熟地黄组方(HBS)治疗肾病综合征的作用机制.通过多个数据库获取肾病综合征基因并进行功能模块分解,找出肾病综合征基因参与的主要生物学过程.通过文献以及数据库查找HBS活性成分和基因靶点,筛选出HBS治疗肾病综合征的有效靶点.通过有效靶点的KEGG和GO富集分析,...     

基于网络药理学和分子对接探讨银杏叶治疗高血压病潜在作用机制

通过网络药理学和分子对接技术探讨银杏叶治疗高血压的潜在作用机制.首先,通过TCMSP、Swiss Target Prediction、Uniprot等数据库获取银杏叶的化学成分与对应靶点;运用OMIM、DrugBank及Gencards疾病数据库搜索高血压相关靶点.然后,取银杏叶对应靶点与高血压相关靶点的交集即可得到银...     

埃及伊蚊不同组织的基因共表达模式分析

【目的】利用权重基因共表达网络分析(weighted gene co-expression network analysis,WGCNA)探索埃及伊蚊Aedes aegypti不同组织基因共表达模式。【方法】从NCBI SRA数据库中选择埃及伊蚊不同组织的转录组数据中具代表性的9种组织(雌雄成蚊的触角和脑,雌蚊的喙、下颚须和卵巢,雄成蚊的前足、中足、后足和腹部末端)的双端测序数据;经过缺失值移除以及方差计算后,筛选出方差最大的5 000个基因,利用R软件中WGCNA包建立埃及伊蚊成蚊不同组织的基因共表达网络并划分模块;然后利用clusterProfiler包对组织特异性模块内的基因进行GO(Gene Ontology)和KEGG(Kyoto Encyclopediaof Genes and Genomes)富集分析,并用Cytoscape软件中的CytoHubba插件筛选共表达模块内的hub基因。【结果】从埃及伊蚊成蚊不同组织中共鉴定出11个基因共表达模块,在雌蚊触角、喙、卵巢、下颚须以及雄蚊脑、腹部末端组织中各鉴定出1个特异性表达模块,雄蚊前足、中足和后足组织中无特异性表达模块。6个组织特异性表达模块内基因功能注释到组织生物学功能;其中,雌蚊触角特异性green模块内基因具有气味结合和嗅觉受体活性等功能;雌蚊喙特异性purple模块内基因具有丝氨酸型肽链内切酶活性和丝氨酸水解酶活性等功能;雄蚊脑特异性blue模块内基因在生物学过程调节、信号转导和神经系统过程等生物学过程中发挥主要作用。利用CytoHubba进一步鉴定出所选组织特异性共表达模块中具有高连通性的hub基因,包括AAEL010426, AAEL002896, AAEL002600, AAEL000961, AAEL007784和AAEL006429。【结论】本研究依据埃及伊蚊不同组织转录组数据,利用WGCNA方法发现了许多重要的基因共表达模块。本研究的结果为蚊虫基因共表达模式分析提供新思路和方法基础,对探究蚊虫不同组织特有的基因资源信息以及功能基因生物信息学研究有参考价值。     

Identification of gene expression profiles and key genes in subchondral bone of osteoarthritis using weighted gene coexpression network analysis

Osteoarthritis (OA) significantly influences the quality life of people around the world. It is urgent to find an effective way to understand the genetic etiology of OA. We used weighted gene coexpression network analysis (WGCNA) to explore the key genes involved in the subchondral bone pathological process of OA. Fifty gene expression profiles of GSE51588 were downloaded from the Gene Expression Omnibus database. The OA‐associated genes and gene ontologies were acquired from JuniorDoc. Weighted gene coexpression network analysis was used to find disease‐related networks based on 21756 gene expression correlation coefficients, hub‐genes with the highest connectivity in each module were selected, and the correlation between module eigengene and clinical traits was calculated. The genes in the traits‐related gene coexpression modules were subject to functional annotation and pathway enrichment analysis using ClusterProfiler. A total of 73 gene modules were identified, of which, 12 modules were found with high connectivity with clinical traits. Five modules were found with enriched OA‐associated genes. Moreover, 310 OA‐associated genes were found, and 34 of them were among hub‐genes in each module. Consequently, enrichment results indicated some key metabolic pathways, such as extracellular matrix (ECM)‐receptor interaction (hsa04512), focal adhesion (hsa04510), the phosphatidylinositol 3'‐kinase (PI3K)‐Akt signaling pathway (PI3K‐AKT) (hsa04151), transforming growth factor beta pathway, and Wnt pathway. We intended to identify some core genes, collagen (COL)6A3, COL6A1, ITGA11, BAMBI, and HCK, which could influence downstream signaling pathways once they were activated. In this study, we identified important genes within key coexpression modules, which associate with a pathological process of subchondral bone in OA. Functional analysis results could provide important information to understand the mechanism of OA.     

Construction of miRNA‐lncRNA‐mRNA co‐expression network affecting EMT‐mediated cisplatin resistance in ovarian cancer

Platinum resistance is one of the major concerns in ovarian cancer treatment. Recent evidence shows the critical role of epithelial–mesenchymal transition (EMT) in this resistance. Epithelial‐like ovarian cancer cells show decreased sensitivity to cisplatin after cisplatin treatment. Our study prospected the association between epithelial phenotype and response to cisplatin in ovarian cancer. Microarray dataset {"type":"entrez-geo","attrs":{"text":"GSE47856","term_id":"47856"}}GSE47856 was acquired from the GEO database. After identifying differentially expressed genes (DEGs) between epithelial‐like and mesenchymal‐like cells, the module identification analysis was performed using weighted gene co‐expression network analysis (WGCNA). The gene ontology (GO) and pathway analyses of the most considerable modules were performed. The protein–protein interaction network was also constructed. The hub genes were specified using Cytoscape plugins MCODE and cytoHubba, followed by the survival analysis and data validation. Finally, the co‐expression of miRNA‐lncRNA‐TF with the hub genes was reconstructed. The co‐expression network analysis suggests 20 modules relating to the Epithelial phenotype. The antiquewhite4, brown and darkmagenta modules are the most significant non‐preserved modules in the Epithelial phenotype and contain the most differentially expressed genes. GO, and KEGG pathway enrichment analyses on these modules divulge that these genes were primarily enriched in the focal adhesion, DNA replication pathways and stress response processes. ROC curve and overall survival rate analysis show that the co‐expression pattern of the brown module''s hub genes could be a potential prognostic biomarker for ovarian cancer cisplatin resistance.