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

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.     

Identifying significant genes and functionally enriched pathways in familial hypercholesterolemia using integrated gene co-expression network analysis

Familial hypercholesterolemia (FH) is a monogenic lipid disorder which promotes atherosclerosis and cardiovascular diseases. Owing to the lack of sufficient published information, this study aims to identify the potential genetic biomarkers for FH by studying the global gene expression profile of blood cells. The microarray expression data of FH patients and controls was analyzed by different computational biology methods like differential expression analysis, protein network mapping, hub gene identification, functional enrichment of biological pathways, and immune cell restriction analysis. Our results showed the dysregulated expression of 115 genes connected to lipid homeostasis, immune responses, cell adhesion molecules, canonical Wnt signaling, mucin type O-glycan biosynthesis pathways in FH patients. The findings from expanded protein interaction network construction with known FH genes and subsequent Gene Ontology (GO) annotations have also supported the above findings, in addition to identifying the involvement of dysregulated thyroid hormone and ErbB signaling pathways in FH patients. The genes like CSNK1A1, JAK3, PLCG2, RALA, and ZEB2 were found to be enriched under all GO annotation categories. The subsequent phenotype ontology results have revealed JAK3I, PLCG2, and ZEB2 as key hub genes contributing to the inflammation underlying cardiovascular and immune response related phenotypes. Immune cell restriction findings show that above three genes are highly expressed by T-follicular helper CD4⁺ T cells, naïve B cells, and monocytes, respectively. These findings not only provide a theoretical basis to understand the role of immune dysregulations underlying the atherosclerosis among FH patients but may also pave the way to develop genomic medicine for cardiovascular diseases.     

Identification of biomarkers of chromophobe renal cell carcinoma by weighted gene co-expression network analysis

Background

Chromophobe renal cell carcinoma (ChRCC) is the second common subtype of non-clear cell renal cell carcinoma (nccRCC), which accounting for 4–5% of renal cell carcinoma (RCC). However, there is no effective bio-marker to predict clinical outcomes of this malignant disease. Bioinformatic methods may provide a feasible potential to solve this problem.

Methods

In this study, differentially expressed genes (DEGs) of ChRCC samples on The Cancer Genome Atlas database were filtered out to construct co-expression modules by weighted gene co-expression network analysis and the key module were identified by calculating module-trait correlations. Functional analysis was performed on the key module and candidate hub genes were screened out by co-expression and MCODE analysis. Afterwards, real hub genes were filter out in an independent dataset GSE15641 and validated by survival analysis.

Results

Overall 2215 DEGs were screened out to construct eight co-expression modules. Brown module was identified as the key module for the highest correlations with pathologic stage, neoplasm status and survival status. 29 candidate hub genes were identified. GO and KEGG analysis demonstrated most candidate genes were enriched in mitotic cell cycle. Three real hub genes (SKA1, ERCC6L, GTSE-1) were selected out after mapping candidate genes to GSE15641 and two of them (SKA1, ERCC6L) were significantly related to overall survivals of ChRCC patients.

Conclusions

In summary, our findings identified molecular markers correlated with progression and prognosis of ChRCC, which might provide new implications for improving risk evaluation, therapeutic intervention, and prognosis prediction in ChRCC patients.

     

Mining the key genes for ventilator-induced lung injury using co-expression network analysis

Mechanical ventilation is extensively adopted in general anesthesia and respiratory failure management, but it can also induce ventilator-induced lung injury (VILI). Therefore, it is of great urgency to explore the mechanisms involved in the VILI pathogenesis, which might contribute to its future prevention and treatment. Four microarray datasets from the GEO database were selected in our investigation, and were subjected to the Weighted Gene Co-Expression Network Analysis (WGCNA) to identify the VILI-correlated gene modules. The limma package in R software was used to identify the differentially expressed genes (DEGs) between the VILI and control groups. WGCNA was constructed by merging the {"type":"entrez-geo","attrs":{"text":"GSE9314","term_id":"9314"}}GSE9314, {"type":"entrez-geo","attrs":{"text":"GSE9368","term_id":"9368"}}GSE9368, {"type":"entrez-geo","attrs":{"text":"GSE11434","term_id":"11434"}}GSE11434 and {"type":"entrez-geo","attrs":{"text":"GSE11662","term_id":"11662"}}GSE11662 datasets. A total of 49 co-expression network modules were determined as associated with VILI. The intersected genes between hub genes screened from DEGs for VILI and those identified using WGCNA were as follows: Tlr2, Hmox1, Serpine1, Mmp9, Il6, Il1b, Ptgs2, Fos and Atf3, which were determined to be key genes for VILI. Those key genes were validated by {"type":"entrez-geo","attrs":{"text":"GSE86229","term_id":"86229"}}GSE86229 and quantitative PCR (qPCR) experiment to have significantly statistical difference in their expression between the VILI and control groups. In a nutshell, nine key genes with expression differences in VILI were screened by WGCNA by integrating multiple datasets.     

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.     

Bioinformatics role of the WGCNA analysis and co-expression network identifies of prognostic marker in lung cancer

Lung cancer is the most talked about cancer in the world. It is also one of the cancers that currently has a high mortality rate. The aim of our research is to find more effective therapeutic targets and prognostic markers for human lung cancer. First, we download gene expression data from the GEO database. We performed weighted co-expression network analysis on the selected genes, we then constructed scale-free networks and topological overlap matrices, and performed correlation modular analysis with the cancer group. We screened the 200 genes with the highest correlation in the cyan module for functional enrichment analysis and protein interaction network construction, found that most of them focused on cell division, tumor necrosis factor-mediated signaling pathways, cellular redox homeostasis, reactive oxygen species biosynthesis, and other processes, and were related to the cell cycle, apoptosis, HIF-1 signaling pathway, p53 signaling pathway, NF-κB signaling pathway, and several cancer disease pathways are involved. Finally, we used the GEPIA website data to perform survival analysis on some of the genes with GS > 0.6 in the cyan module. CBX3, AHCY, MRPL12, TPGB, TUBG1, KIF11, LRRC59, MRPL17, TMEM106B, ZWINT, TRIP13, and HMMR was identified as an important prognostic factor for lung cancer patients. In summary, we identified 12 mRNAs associated with lung cancer prognosis. Our study contributes to a deeper understanding of the molecular mechanisms of lung cancer and provides new insights into drug use and prognosis.     

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

【目的】分析广叶绣球菌（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个核心基因,其可能参与了基质降解及相关基因的表达调控。【结论】不同木质纤维素类型显著影响了广叶绣球菌木质纤维素降解基因的差异表达轮廓,这种差异反映了广叶绣球菌对不同木质纤维素特异的降解策略。     

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 reveals potential candidate genes affecting drip loss in pork

Drip loss is an essential evaluation indicator for pork quality. It is closely related to other meat quality indicators, including water-holding capacity, water loss rate and pH value at 45 min (pH₁) and 24 h post-mortem (pH₂), and is influenced by environmental and genetic factors and their interactions. We previously conducted differentially expressed gene analysis to identify candidate genes affecting drip loss using eight individuals with extremely high- and low-drip loss selected from 28 purebred Duroc pigs. Using 28 identical samples, in the present study, we performed weighted gene co-expression network analysis with drip loss and drip loss-related traits, including water-holding capacity, water loss rate, pH₁ and pH₂. A total of 25 modules were identified, and five of them correlated with at least two drip loss or drip loss-related traits. After functional enrichment analysis of genes in the five modules, three modules were found to be critical, as their genes were significantly involved in amino acid metabolism, immune response and apoptosis, which have potential relationships with drip loss. Furthermore, we identified five candidate genes affecting drip loss in one critical module, AASS, BCKDHB, ALDH6A1, MUT and MCCC1, as they overlapped with differentially expressed genes detected in our previous study, exhibited protein–protein interactions and had potential biological functions in affecting drip loss according to the literature. The outcomes of the present study enhance our understanding of the molecular mechanisms underlying drip loss and will aid in improving the pork quality.