Systematic prediction of key genes for ovarian cancer by co-expression network analysis

Ovarian cancer (OC) is the most lethal gynaecological malignancy, characterized by high recurrence and mortality. However, the mechanisms of its pathogenesis remain largely unknown, hindering the investigation of the functional roles. This study sought to identify key hub genes that may serve as biomarkers correlated with prognosis. Here, we conduct an integrated analysis using the weighted gene co-expression network analysis (WGCNA) to explore the clinically significant gene sets and identify candidate hub genes associated with OC clinical phenotypes. The gene expression profiles were obtained from the MERAV database. Validations of candidate hub genes were performed with RNASeqV2 data and the corresponding clinical information available from The Cancer Genome Atlas (TCGA) database. In addition, we examined the candidate genes in ovarian cancer cells. Totally, 19 modules were identified and 26 hub genes were extracted from the most significant module (R² = .53) in clinical stages. Through the validation of TCGA data, we found that five hub genes (COL1A1, DCN, LUM, POSTN and THBS2) predicted poor prognosis. Receiver operating characteristic (ROC) curves demonstrated that these five genes exhibited diagnostic efficiency for early-stage and advanced-stage cancer. The protein expression of these five genes in tumour tissues was significantly higher than that in normal tissues. Besides, the expression of COL1A1 was associated with the TAX resistance of tumours and could be affected by the autophagy level in OC cell line. In conclusion, our findings identified five genes could serve as biomarkers related to the prognosis of OC and may be helpful for revealing pathogenic mechanism and developing further research.     

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.     

Analysis of gene co-expression networks and function modules at different developmental stages of chicken breast muscle

The development of poultry muscle fibers after hatching is closely related to meat quality and production efficiency. It is necessary to identify functional modules (groups of functionally related genes) related to muscle development at different developmental stages, and to investigate their relationships based on the weighted gene co-expression network analysis (WGCNA) methods. Accordingly, we investigated the co-expression associations between genes related to chicken breast muscle at four different developmental stages (between 2 and 14 weeks of age), and systematically analyzed the network topology in Jinmao Hua chicken. As a result, 2341 differentially expressed genes were identified and subjected to co-expression analysis. Four modules were identified to be related to a particular growth stage for the development of breast muscle. A series of genes with the highest connectivity were identified in the pink (2 weeks), yellow (6 weeks), green (10 weeks) and black modules (14 weeks), respectively, and visualized by Cytoscape. These hub genes (FGF, MAPKAPK5, NRG1, SCD, ACSL1, PPAR etc.) were mainly enriched in 15 pathways, such as MAPK signaling pathway, NRG/ErbB signaling pathway, and insulin signaling pathway. They shared biological functions related to development of breast muscle and adipogenesis. This is the first study of gene network with different stages of muscle development in Jinmao Hua chicken to observe co-expression patterns. It may contribute to the underlied molecular mechanisms of chicken breast muscle development.     

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.     

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

【目的】利用权重基因共表达网络分析(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方法发现了许多重要的基因共表达模块。本研究的结果为蚊虫基因共表达模式分析提供新思路和方法基础,对探究蚊虫不同组织特有的基因资源信息以及功能基因生物信息学研究有参考价值。     

刺盘孢菌的蛋白互作预测及侵染早期共表达模块分析

刺盘孢菌是一类重要的植物病原真菌,在全球范围内危害众多单双子叶植物。有许多研究对病菌的侵染模式进行了探索,但仍未阐明其确切的分子机制。本研究在预测病菌蛋白互作的基础上,结合表达谱数据对病菌在活体生存环境下的共表达模块进行挖掘和分析,以期为分子机制的研究提供新的线索。通过同源映射法和结构域法,预测得到刺盘孢菌的4 288个蛋白之间存在41 700个潜在互作,其中39 776个互作发生于异源蛋白之间,1 924个互作发生于同一蛋白内。将蛋白互作数据分别与4个表达谱数据进行整合,构建得到离体I、离体II、活体I和活体II 4个共表达互作组。对离体和活体互作组的共有基因的表达水平进行比较分析,结果表明,与离体互作组相比,活体互作组中与翻译、蛋白代谢等有关的基因表达水平下降,与离子转运、糖物质转运等有关的基因表达水平上升,暗示了物质转运在刺盘孢菌侵染早期的重要作用。进一步对活体互作组进行模块化分析,结果表明,活体I和活体II的特异模块分别与胁迫响应、肌动蛋白纤维长度调控有关,其中胁迫响应子网是以热激蛋白Hsp70为核心的互作簇,可能参与病菌对寄主的识别与对抗;肌动蛋白纤维长度调控子网则可能与病菌菌丝在寄主细胞间的延伸有关。     

Novel key genes in triple-negative breast cancer identified by weighted gene co-expression network analysis

Triple-negative breast cancer (TNBC) is a special subtype of breast cancer (BC) with poor prognosis. Although some molecular mechanisms of TNBC have been elucidated, the efficacy of current treatments is limited. Therefore, it is urgently demanded to screen for novel biomarkers and drug targets for TNBC. In this study, we obtained four independent data sets (GSE76250, GSE31448, GSE43358, and METABRIC) from the Gene Expression Omnibus (GEO) database and the cBioPortal website. In the GSE76250 data set, 890 differentially expressed genes were identified and weighted gene co-expression network analysis was performed based on them. Then, two preserved modules associated with the KI67 score were detected. Gene ontology and pathway enrichment analyses showed genes in the modules participated in some cancer-related biological processes or pathways. Non-SMC condensin I complex subunit G (NCAPG) and ATP-binding cassette subfamily A member 9 (ABCA9) were identified as hub genes of the modules, and the significance of hub genes was validated in the GSE43358 data set. Finally, their prognostic value was assessed by survival analysis. These findings suggested that NCAPG and ABCA9 may be the key genes of TNBC. Moreover, ABCA9 was first reported in TNBC. They deserved further studies.     

应用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个枢纽基因以及一些关键信号通路和生物学过程,可为后续深入研究细菌性败血症致病机制奠定理论依据。     

Identification of key lncRNAs in the carcinogenesis and progression of colon adenocarcinoma by co-expression network analysis

Colon adenocarcinoma (COAD) is one of the most common cancers, and its carcinogenesis and progression is influenced by multiple long non-coding RNAs (lncRNA), especially through the miRNA sponge effect. In this study, more than 4000 lncRNAs were re-annotated from the microarray datasets through probe sequence mapping to obtain reliable lncRNA expression profiles. As a systems biology method for describing the correlation patterns among genes across microarray samples, weighted gene co-expression network analysis was conducted to identify lncRNA modules associated with the five stepwise stages from normal colonic samples to COAD (n = 94). In the most relevant module (R² = −0.78, P = 4E-20), four hub lncRNAs were identified (CTD-2396E7.11, PCGF5, RP11-33O4.1, and RP11-164P12.5). Then, these four hub lncRNAs were validated using two other independent datasets including GSE20916 (n = 145) and GSE39582 (n = 552). The results indicated that all hub lncRNAs were significantly negatively correlated with the three-stage colonic carcinogenesis, as well as TNM stages in COAD (one-way analysis of variance P < 0.05). Kaplan-Meier survival curve showed that patients with higher expression of each hub lncRNA had a significantly higher overall survival rate and lower relapse risk (log-rank P < 0.05). In conclusion, through co-expression analysis, we identified and validated four key lncRNAs in association with the carcinogenesis and progression of COAD, and these lncRNAs might have important clinical implications for improving the risk stratification, therapeutic decision and prognosis prediction in COAD patients.     

circRNA-associated ceRNA network construction reveals the circRNAs involved in the progression and prognosis of breast cancer

Recently, increasing evidences show that circular RNAs (circRNAs) are important regulators of various diseases, especially cancer. However, the regulatory role and the potential mechanism of action of circRNAs in breast cancer remain largely unknown. In this study, weighted gene co-expression network analysis was conducted with the differentially expressed miRNAs and mRNAs in breast cancer from The Cancer Genome Atlas database to identify the key modules associated with the carcinogenesis of breast cancer. In the significant turquoise and brown modules, 22 miRNAs and 1877 mRNAs were identified, respectively. Then, We compared and predicted the target genes and performed survival analysis to identify the miRNAs and mRNAs related to the prognosis of breast cancer. A circRNA-related competitive endogenous RNA network was identified by database co-screening, and deleted in liver cancer 1 (DLC1) was identified as a key gene. Finally, to assess how genes in key modules and key genes contribute to the development of breast cancer, relevant pathway information was obtained through DAVID and Gene Set Enrichment Analysis. These data demonstrated that three circRNAs (hsa-circ-0083373, hsa-circ-0083374, and hsa-circ-0083375) that regulate DLC1 expression via hsa-mir-511 and are involved in the pathogenesis and development of breast cancer.     

前列腺癌多阶段相关枢纽蛋白的差异共表达模式

通过枢纽蛋白与其互作邻居之间的共表达关系在从前列腺正常组织到原位癌及原位癌到转移的两个发展阶段的变化,寻找在前列腺癌发生和发展过程中动态改变的枢纽蛋白,为探讨前列腺癌的发病与转移机理提供线索。基于前列腺正常组织、原位癌与癌转移后组织的基因表达谱数据,在前列腺癌的两个发展阶段分别检测到1578和2347个动态改变的枢纽蛋白。平均超过95%的在一个发展阶段动态改变的枢纽蛋白与在另一发展阶段动态改变的枢纽蛋白中的至少一个蛋白共享显著多的互作邻居。另外,两组动态改变的枢纽蛋白中有780个交叠蛋白,其中大约50%的蛋白与其互作邻居之间的共表达关系在两个发展阶段的强弱变化方向相反。结果提示,在前列腺癌的发生和发展过程中,动态改变的枢纽蛋白影响的功能通路高度一致,但是它们与其互作邻居之间有着不同的共表达模式。     

镉胁迫下地钱转录组的性别特异性响应机制

地钱(Marchantiapolymorpha)作为雌雄异株的苔纲植物,具有对重金属胁迫敏感的特性,且种群中雌雄比例差异较大,是研究重金属环境下植物性二态的理想对象之一。为探究雌雄地钱转录组对重金属镉的响应机制差异,利用高通量测序和加权基因共表达网络分析(WGCNA),鉴定了在镉胁迫条件下与不同性别地钱有关的模块和基因。结果表明,镉胁迫响应转录因子雌雄差异主要集中于bHLH、 ERF和MYB家族。关键差异响应基因包括MARPO＿0147s0038、MARPO＿1326s0001和MARPO＿0015s0058。差异响应通路雌性主要集中在信号转导通路,雄性则集中于类黄酮生物合成。研究结果有助于揭示雌雄异株地钱对镉胁迫的反应过程和响应机制,可为理解雌雄异株植物对重金属胁迫的性别特异性响应机制提供参考。     

Finding disease specific alterations in the co-expression of genes

MOTIVATION: Standard analysis routines for microarray data aim at differentially expressed genes. In this paper, we address the complementary problem of detecting sets of differentially co-expressed genes in two phenotypically distinct sets of expression profiles. RESULTS: We introduce a score for differential co-expression and suggest a computationally efficient algorithm for finding high scoring sets of genes. The use of our novel method is demonstrated in the context of simulations and on real expression data from a clinical study.     

Network-based approach to identify prognosis-related genes in tamoxifen-treated patients with estrogen receptor-positive breast cancer

Tamoxifen is an estrogen receptor (ER) antagonist that is most commonly used for the treatment of ER-positive breast cancer. However, tamoxifen resistance remains a major cause of cancer recurrence and progression. Here, we aimed to identify hub genes implicated in the progression and prognosis of ER-positive breast cancer following tamoxifen treatment. Microarray data ({"type":"entrez-geo","attrs":{"text":"GSE9893","term_id":"9893"}}GSE9893) for 155 tamoxifen-treated primary ER-positive breast cancer samples were obtained from the Gene Expression Omnibus database. In total, 1706 differentially expressed genes (DEGs), including 859 up-regulated and 847 down-regulated genes, were identified between relapse and relapse-free samples. Weighted correlation network analysis clustered genes into 13 modules, among which the tan and blue modules were the most significantly related to prognosis. From these two modules, we further identified and validated two prognosis-related hub genes (G-rich RNA sequence binding factor 1 (GRSF1) and microtubule-associated protein τ (MAPT)) via survival analysis based on several publicly available datasets. High expression of GRSF1 predicted poor prognosis, whereas MAPT indicated favorable outcomes in ER-positive breast cancer. Using breast cancer cell lines and tissue samples, we confirmed that GRSF1 was significantly up-regulated and MAPT was down-regulated in the tamoxifen-resistant group compared with the tamoxifen-sensitive group. The prognostic value of GRSF1 and MAPT was also verified in 48 tamoxifen-treated ER-positive breast cancer patients in our hospital. Gene set enrichment analysis (GSEA) suggested that GRSF1 was potentially involved in RNA degradation and cell cycle pathways, while MAPT was strongly linked to immune-related signaling pathways. Taken together, our findings established novel prognostic biomarkers to predict tamoxifen sensitivity, which may facilitate individualized management of breast cancer.