利用加权基因共表达网络分析方法挖掘香菇发育不同阶段的关键基因

香菇是世界产量第二大食用菌,栽培历史悠久.在木屑袋料栽培模式下,香菇发育可以分为菌丝生长期(G)、菌丝褐化期(B)、原基形成期(P)以及出菇期(FB)4个阶段.褐化期和原基形成期是香菇从营养生长期到生殖生长两个关键发育阶段,对香菇子实体产量和质量至关重要.本研究以3种不同栽培材料为重复,对香菇发育的前3个阶段进行了转录...     

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.     

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.     

Prognostic genes of hepatocellular carcinoma based on gene coexpression network analysis

Hepatocellular carcinoma (HCC) is the most common subtype in liver cancer whose prognosis is affected by malignant progression associated with complex gene interactions. However, there is currently no available biomarkers associated with HCC progression in clinical application. In our study, RNA sequencing expression data of 50 normal samples and 374 tumor samples was analyzed and 9225 differentially expressed genes were screened. Weighted gene coexpression network analysis was then conducted and the blue module we were interested was identified by calculating the correlations between 17 gene modules and clinical features. In the blue module, the calculation of topological overlap was applied to select the top 30 genes and these 30 genes were divided into the green group (11 genes) and the yellow group (19 genes) through searching whether these genes were validated by in vitro or in vivo experiments. The genes in the green group which had never been validated by any experiments were recognized as hub genes. These hub genes were subsequently validated by a new data set GSE76427 and KM Plotter Online Tool, and the results indicated that 10 genes (FBXO43, ARHGEF39, MXD3, VIPR1, DNASE1L3, PHLDA1, CSRNP1, ADR2B, C1RL, and CDC37L1) could act as prognosis and progression biomarkers of HCC. In summary, 10 genes who have never been mentioned in HCC were identified to be associated with malignant progression and prognosis of patients. These findings may contribute to the improvement of the therapeutic decision, risk stratification, and prognosis prediction for HCC patients.     

Key genes and functional coexpression modules involved in the pathogenesis of systemic lupus erythematosus

We performed a systematic review of genome‐wide gene expression datasets to identify key genes and functional modules involved in the pathogenesis of systemic lupus erythematosus (SLE) at a systems level. Genome‐wide gene expression datasets involving SLE patients were searched in Gene Expression Omnibus and ArrayExpress databases. Robust rank aggregation (RRA) analysis was used to integrate those public datasets and identify key genes associated with SLE. The weighted gene coexpression network analysis (WGCNA) was adapted to identify functional modules involved in SLE pathogenesis, and the gene ontology enrichment analysis was utilized to explore their functions. The aberrant expressions of several randomly selected key genes were further validated in SLE patients through quantitative real‐time polymerase chain reaction. Fifteen genome‐wide gene expression datasets were finally included, which involved a total of 1,778 SLE patients and 408 healthy controls. A large number of significantly upregulated or downregulated genes were identified through RRA analysis, and some of those genes were novel SLE gene signatures and their molecular roles in etiology of SLE remained vague. WGCNA further successfully identified six main functional modules involved in the pathogenesis of SLE. The most important functional module involved in SLE included 182 genes and mainly enriched in biological processes, including defense response to virus, interferon signaling pathway, and cytokine‐mediated signaling pathway. This study identifies a number of key genes and functional coexpression modules involved in SLE, which provides deepening insights into the molecular mechanism of SLE at a systems level and also provides some promising therapeutic targets.     

Identification of key gene modules and pathways of human glioma through coexpression network

Glioma causes great harm to people worldwide. Systemic coexpression analysis of this disease could be beneficial for the identification and development of new prognostic and predictive markers in the clinical management of glioma. In this study, we extracted data sets from the Gene Expression Omnibus data set by using “glioma” as the keyword. Then, a coexpression module was constructed with the help of Weighted Gene Coexpression Network Analysis software. Besides, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on the genes in these modules. As a result, the critical modules and target genes were identified. Eight coexpression modules were constructed using the 4,000 genes with a high expression value of the total 141 glioma samples. The result of the analysis of the interaction among these modules showed that there was a high scale independence degree among them. The GO and KEGG enrichment analyses showed that there was a significant difference in the enriched terms and degree among these eight modules, and module 5 was identified as the most important module. Besides, the pathways it was enriched in, hsa04510: Focal adhesion and hsa04610: Complement and coagulation cascades, were determined as the most important pathways. In summary, module 5 and the pathways it was enriched in, hsa04510: Focal adhesion and has 04610: Complement and coagulation cascades, have the potential to serve as biomarkers for patients with glioma.     

利用基因共表达网络分析蛋白质的相互作用

基于Hom in的基因共表达网络的比较分析,发现人类基因共表达网络和蛋白质相互作用数据之间存在一定的相关性。采用基因本体论对这两个网络重叠区域进行基因分类后发现,这些编码的蛋白质主要集中在对刺激物的应答途径之中。通过对该途径中的蛋白质相互作用网络作图,获得了两个独立的功能模块。通过对模块中的基因分类和关键基因分析得出两者分别对应于内外源刺激物的应答功能。本研究对于利用不断丰富的核酸公共数据信息挖掘蛋白质相互作用的研究具有积极的促进作用。     

Gene coexpression analysis offers important modules and pathway of human lung adenocarcinomas

Lung adenocarcinomas injured greatly on the people worldwide. Although clinic experiments and gene profiling analyses had been well performed, to our knowledge, systemic coexpression analysis of human genes for this cancer is still limited to date. Here, using the published data GSE75037, we built the coexpression modules of genes by Weighted Gene Co-Expression Network Analysis (WGCNA), and investigated function and protein–protein interaction network of coexpression genes by Database for Annotation, visualization, and Integrated Discovery (DAVID) and String database, respectively. First, 11 coexpression modules were conducted for 5,000 genes in the 83 samples recently. Number of genes for each module ranged from 90 to 1,260, with the mean of 454. Second, interaction relationships of hub-genes between pairwise modules showed great differences, suggesting relatively high scale independence of the modules. Third, functional enrichment of the coexpression modules showed great differences. We found that genes in modules 8 significantly enriched in the biological process and/or pathways of cell adhesion, extracellular matrix (ECM)–receptor interaction, focal adhesion, and PI3K-Akt signaling pathway, and so forth. It was inferred as the key module underlying lung adenocarcinomas. Furthermore, PPI analysis revealed that the genes COL1A1, COL1A2, COL3A1, CTGF, and BGN owned the largest number of adjacency genes, unveiling that they may functioned importantly during the occurrence of lung adenocarcinomas. To summary, genes involved in cell adhesion, ECM–receptor interaction, focal adhesion, and PI3K-Akt signaling pathway play crucial roles in human lung adenocarcinomas.     

应用GSEA和WGCNA方法分析细菌性败血症宿主关键差异基因及意义

【目的】采用生物信息学方法分析公共数据库来源的细菌性败血症患者全血转录组学表达谱,探讨细菌败血症相关的宿主关键差异基因及意义。【方法】基于GEO数据库中GSE80496和GSE72829全血转录组基因数据集,采用GEO2R、基因集富集分析(GSEA)联用加权基因共表达网络分析(WGCNA)筛选细菌性败血症患者相比健康人群显著改变的差异基因,通过R软件对交集基因进行GO功能分析和KEGG富集分析。同时,通过String 11.0和Cytoscape分析枢纽基因,验证枢纽基因在数据集GSE72809(Health组52例,Definedsepsis组52例)全血标本中的表达情况,并探讨婴儿性别、月(胎)龄、出生体重、是否接触抗生素等因素与靶基因表达谱间的关系。【结果】分析GSE80496和GSE72829数据集分别筛选得到932个基因和319个基因,联合WGCNA枢纽模块交集得到与细菌性败血症发病相关的10个枢纽基因(MMP9、ITGAM、CSTD、GAPDH、PGLYRP1、FOLR3、OSCAR、TLR5、IL1RN和TIMP1);GSEA分析获得关键通路(氨基酸糖类-核糖代谢、PPAR信号通路、聚糖生物合成通路、自噬调控通路、补体、凝血因子级联反应、尼古丁和烟酰胺代谢、不饱和脂肪酸生物合成和阿尔兹海默症通路)及生物学过程(类固醇激素分泌、腺苷酸环化酶的激活、细胞外基质降解和金属离子运输)。【结论】本项研究通过GEO2R、GSEA联用WGCNA分析,筛选出与细菌性败血症发病相关的2个枢纽模块、10个枢纽基因以及一些关键信号通路和生物学过程,可为后续深入研究细菌性败血症致病机制奠定理论依据。     

基于转录组和加权基因共表达网络的广叶绣球菌木质纤维素降解特性分析

【目的】分析广叶绣球菌（Sparassis latifolia）在不同木质纤维素诱导条件下基因表达差异,为广叶绣球菌木质纤维素降解关键基因和分子机制研究提供参考。【方法】以松木、杉木、甘蔗渣和天然堆积发酵后的杉木和发酵后的甘蔗渣为碳源,在液体培养条件下培养诱导广叶绣球菌,对其转录组进行测序研究,并对不同木质纤维素诱导样本进行加权基因共表达网络分析（weighted gene co-expression network analysis,WGCNA）。【结果】杉木培养与松木培养比较组差异表达基因最少（20个）,蔗渣培养与松木培养比较组差异表达基因最多（486个）。基因本体（gene ontology,GO）富集分析结果表明,差异表达基因主要涉及氧化还原酶活性、单加氧酶活性和铁离子结合活性等,京都基因和基因组百科全书（Kyoto encyclopedia of genes and genomes,KEGG）通路富集分析结果表明,差异表达基因主要涉及戊糖和葡萄糖醛酸转换、甲烷代谢和乙醛酸盐和二羧酸盐代谢等通路。发酵甘蔗渣为碳源培养时,纤维素和半纤维素降解相关的糖苷水解酶基因表达量总体上较高,而未发酵的松木、杉木和甘蔗渣为碳源培养时木质素降解或修饰相关的碳水化合物辅助酶基因表达量总体上较高。利用WGCNA共鉴定出10个共表达模块,其中green模块与未发酵蔗渣诱导显著正相关,blue模块与发酵甘蔗渣诱导显著正相关,magenta和turquoise模块与发酵杉木诱导显著正相关。GO富集分析结果表明,turquoise模块内基因显著富集到尿素跨膜转运子活性、甲基转移酶活性和单加酶活性等,blue模块基因显著富集到水解酶活性和β-甘露糖苷酶活性。KEGG通路富集分析结果表明,blue模块内基因显著富集的通路有半乳糖代谢、果糖和甘露糖代谢、苯丙氨酸代谢、精氨酸和脯氨酸代谢等。通过构建互作网络图挖掘到12个核心基因,其可能参与了基质降解及相关基因的表达调控。【结论】不同木质纤维素类型显著影响了广叶绣球菌木质纤维素降解基因的差异表达轮廓,这种差异反映了广叶绣球菌对不同木质纤维素特异的降解策略。     

Coexpression network analysis linked H2AFJ to chemoradiation resistance in colorectal cancer

Neoadjuvant chemoradiotherapy (CRT) resistance is a complex phenomenon and it remains a major problem for patients with a priori resistant tumor. Therefore, there is a strong need to investigate molecular biomarkers which may guide for treatment decision-making. In our study, weighted gene coexpression network analysis was applied to identify CRT-resistance hub modules in 12 colorectal cancer (CRC) cell lines with different CRT sensitivities from GSE20298 data set. The green module and purple module had the highest correlations with CRT resistance. Gene ontology enrichment analysis indicated that the function of these two modules focused on interferon-mediated signaling pathway, immune response, chromatin modulation, Rho GTPases activities, and regulation of apoptotic process. Then, 15 hub genes in both the coexpression and protein-protein interaction networks were selected. Among these hub genes, higher H2A histone family member J (H2AFJ) expression was independently validated in patient cohorts from two testing data sets of GSE46862 and GSE68204 to be related to CRT resistance. The receiver operating characteristic curve showed that H2AFJ could efficiently distinguish CRT-resistance cases from CRT-sensitive cases in another two testing data sets. Furthermore, meta-analysis of 12 Gene Expression Omnibus–sourced data sets showed that H2AFJ messenger RNA levels were significantly higher in CRC tissues than in normal colon tissues. High H2AFJ expression was correlated with a significant worse event- and relapse-free survival by analyzing the data from the R2: Genomics Analysis and Visualization Platform. Gene set enrichment analysis determined that the mechanism of H2AFJ-mediated CRT resistance might involve the ERK5 (MAPK7), human immunodeficiency virus Nef (HIV Nef), and inflammatory pathways. This study is the first, to the best of our knowledge, to implicate and verify H2AFJ as an effective new marker for CRT response prediction.     

Identification of hub genes in colorectal cancer based on weighted gene co-expression network analysis and clinical data from The Cancer Genome Atlas

Colorectal cancer (CRC) is one of the most common tumors worldwide and is associated with high mortality. Here we performed bioinformatics analysis, which we validated using immunohistochemistry in order to search for hub genes that might serve as biomarkers or therapeutic targets in CRC. Based on data from The Cancer Genome Atlas (TCGA), we identified 4832 genes differentially expressed between CRC and normal samples (1562 up-regulated and 3270 down-regulated in CRC). Gene ontology (GO) analysis showed that up-regulated genes were enriched mainly in organelle fission, cell cycle regulation, and DNA replication; down-regulated genes were enriched primarily in the regulation of ion transmembrane transport and ion homeostasis. Weighted gene co-expression network analysis (WGCNA) identified eight gene modules that were associated with clinical characteristics of CRC patients, including brown and blue modules that were associated with cancer onset. Analysis of the latter two hub modules revealed the following six hub genes: adhesion G protein-coupled receptor B3 (BAI3, also known as ADGRB3), cyclin F (CCNF), cytoskeleton-associated protein 2 like (CKAP2L), diaphanous-related formin 3 (DIAPH3), oxysterol binding protein-like 3 (OSBPL3), and RERG-like protein (RERGL). Expression levels of these hub genes were associated with prognosis, based on Kaplan–Meier survival analysis of data from the Gene Expression Profiling Interactive Analysis database. Immunohistochemistry of CRC tumor tissues confirmed that OSBPL3 is up-regulated in CRC. Our findings suggest that CCNF, DIAPH3, OSBPL3, and RERGL may be useful as therapeutic targets against CRC. BAI3 and CKAP2L may be novel biomarkers of the disease.