Systematic analysis on multiple Gene Expression Omnibus data sets reveals fierce immune response in hepatitis B virus‐related acute liver failure

Acute liver failure (ALF) caused by hepatitis B virus (HBV) is common type of liver failure in the world, with high morbidity and mortality rates. However, the prevalence, genetic background and factors determining the development of HBV‐related ALF are rarely studied. In this study, we examined three Gene Expression Omnibus (GEO) data sets by bioinformatics analysis to identify differentially expressed genes (DEGs), key biological processes and pathways. Immune infiltration analysis showed high immune cells infiltration in HBV‐related ALF tissue. We then confirmed natural killer cells and macrophages infiltration in clinical samples by immunohistochemistry assay, implying these cells play a significant role in HBV‐ALF. We found 1277 genes were co‐up‐regulated and that 1082 genes were co‐down‐regulated in the 3 data sets. Inflammation‐related pathways were enriched in the co‐up‐regulated genes and synthetic metabolic pathways were enriched in the co‐down‐regulated genes. WGCNA also revealed a key module enriching in immune inflammation response and identified 10 hub genes, differentially expressed in an independent data set. In conclusion, we identified fierce immune inflammatory response to elucidate the immune‐driven mechanism of HBV‐ALF and 10 hub genes based on gene expression profiles.     

Genome‐wide detection of CNVs associated with beak deformity in chickens using high‐density 600K SNP arrays

Beak deformity (crossed beaks) is found in several indigenous chicken breeds including Beijing‐You studied here. Birds with deformed beaks have reduced feed intake and poor production performance. Recently, copy number variation (CNV) has been examined in many species and is recognized as a source of genetic variation, especially for disease phenotypes. In this study, to unravel the genetic mechanisms underlying beak deformity, we performed genome‐wide CNV detection using Affymetrix chicken high‐density 600K data on 48 deformed‐beak and 48 normal birds using penncnv . As a result, two and eight CNV regions (CNVRs) covering 0.32 and 2.45 Mb respectively on autosomes were identified in deformed‐beak and normal birds respectively. Further RT‐qPCR studies validated nine of the 10 CNVRs. The ratios of six CNVRs were significantly different between deformed‐beak and normal birds (P < 0.01). Within these six regions, three and 21 known genes were identified in deformed‐beak and normal birds respectively. Bioinformatics analysis showed that these genes were enriched in six GO terms and one KEGG pathway. Five candidate genes in the CNVRs were further validated using RT‐qPCR. The expression of LRIG2 (leucine rich repeats and immunoglobulin like domains 2) was lower in birds with deformed beaks (P < 0.01). Therefore, the LRIG2 gene could be considered a key factor in view of its known functions and its potential roles in beak deformity. Overall, our results will be helpful for future investigations of the genomic structural variations underlying beak deformity in chickens.     

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.     

Construction of co-expression modules related to survival by WGCNA and identification of potential prognostic biomarkers in glioblastoma

Glioblastoma (GBM) is a malignant brain tumour with poor prognosis. The potential pathogenesis and therapeutic target are still need to be explored. Herein, TCGA expression profile data and clinical information were downloaded, and the WGCNA was conducted. Hub genes which closely related to poor prognosis of GBM were obtained. Further, the relationship between the genes of interest and prognosis of GBM, and immune microenvironment were analysed. Patients from TCGA were divided into high- and low-risk group. WGCNA was applied to the high- and low-risk group and the black module with the lowest preservation was identified which could distinguish the prognosis level of these two groups. The top 10 hub genes which were closely related to poor prognosis of patients were obtained. GO analysis showed the biological process of these genes mainly enriched in: Cell cycle, Progesterone-mediated oocyte maturation and Oocyte meiosis. CDCA5 and CDCA8 were screened out as the genes of interest. We found that their expression levels were closely related to overall survival. The difference analysis resulted from the TCGA database proved both CDCA5 and CDCA8 were highly expressed in GBM. After transfection of U87-MG cells with small interfering RNA, it revealed that knockdown of the CDCA5 and CDCA8 could influence the biological behaviours of proliferation, clonogenicity and apoptosis of GBM cells. Then, single-gene analysis was performed. CDCA5 and CDCA8 both had good correlations with genes that regulate cell cycle in the p53 signalling pathway. Moreover, it revealed that high amplification of CDCA5 was correlated with CD8⁺ T cells while CDCA8 with CD4⁺ T cells in GBM. These results might provide new molecular targets and intervention strategy for GBM.     

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.     

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

Gene expression and methylation profiles identified CXCL3 and CXCL8 as key genes for diagnosis and prognosis of colon adenocarcinoma

Colon adenocarcinoma (COAD) is one of the most common malignant tumors with high morbidity and mortality rates worldwide. Due to the poor clinical outcomes, it is indispensable to investigate novel biomarkers for the diagnosis and prognosis of COAD. The aim of this study is to explore key genes as potential biomarkers for the diagnosis and prognosis of COAD for clinical utility. Gene expression profiles (GSE44076 and GSE44861) and gene methylation profile (GSE29490) were analyzed to identify the aberrantly methylated-differentially expressed genes by R language and Perl software. Function enrichments were performed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. Moreover, hub genes were identified through protein–protein interaction (PPI) network. Besides, key genes were found by the module analysis and The Cancer Genome Atlas (TCGA) survival analysis. Finally, TCGA data and quantitative real-time polymerase chain reaction (RT-qPCR) was used to validate the key genes involved in COAD. Our study found two hypomethylation-high-expression genes (CXCL3 and CXCL8) in COAD tissues compared with the adjacent normal tissues. These results were also confirmed by RT-qPCR with 25 pairs of COAD and adjacent normal tissues. Meanwhile, low expression of the two genes was associated with poor survival in patients with COAD. CXCL3 and CXCL8 may serve as key genes in the diagnosis and prognosis for COAD.     

Expression and functional analysis of lncRNAs in the hippocampus of immature rats with status epilepticus

Long non‐coding RNAs (lncRNAs) have been implicated in the regulation of gene expression at various levels. However, to date, the expression profile of lncRNAs in status epilepticus (SE) was unclear. In our study, the expression profile of lncRNAs was investigated by high‐throughput sequencing based on a lithium/pilocarpine‐induced SE model in immature rats. Furthermore, weighted correlation network analysis (WGCNA), gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed to construct co‐expression networks and establish functions of the identified hub lncRNAs in SE. The functional role of a hub lncRNA (NONRATT010788.2) in SE was investigated in an in vitro model. Our results indicated that 7082 lncRNAs (3522 up‐regulated and 3560 down‐regulated), which are involved in cell proliferation, inflammatory responses, angiogenesis and autophagy, were dysregulated in the hippocampus of immature rats with SE. Additionally, WGCNA identified 667 up‐regulated hub lncRNAs in turquoise module that were involved in apoptosis, inflammatory responses and angiogenesis via regulation of HIF‐1, p53 and chemokine signalling pathways and via inflammatory mediator regulation of TRP channels. Knockdown of an identified hub lncRNA (NONRATT010788.2) inhibited neuronal apoptosis in vitro. Taken together, our study is the first to demonstrate the expression profile and potential function of lncRNAs in the hippocampus of immature rats with SE. The defined hub lncRNAs may participate in the pathogenesis of SE via lncRNA‐miRNA‐mRNA network.     

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.     

Weighted gene co‐expression network analysis can sort cancer‐associated fibroblast‐specific markers promoting bladder cancer progression

The role of cancer‐associated fibroblasts (CAFs) has been thoroughly investigated in tumour microenvironments but not in bladder urothelial carcinoma (BLCA). The cell fraction of CAFs gradually increased with BLCA progression. Weighted gene co‐expression network analysis (WGCNA) revealed a specific gene expression module of CAFs that are relevant to cancer progression and survival status. Fifteen key genes of the module were consistent with a fibroblast signature in single‐cell RNA sequencing, functionally related to the extracellular matrix, and significant in survival analysis and tumour staging. A comparison of the luminal‐infiltrated versus luminal‐papillary subtypes and fibroblast versus urothelial carcinoma cell lines and immunohistochemical data analysis demonstrated that the key genes were specifically expressed in CAFs. Moreover, these genes are highly correlated with previously reported CAF markers. In summary, CAFs play a major role in the progression of BLCA, and the 15 key genes act as BLCA‐specific CAF markers and can predict CAF changes. WGCNA can, therefore, be used to sort CAF‐specific gene set in cancer tissues.     

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.     

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.     

Variation in assembly stoichiometry in non‐metazoan homologs of the hub domain of Ca2+/calmodulin‐dependent protein kinase II

The multi‐subunit Ca²⁺/calmodulin‐dependent protein kinase II (CaMKII) holoenzyme plays a critical role in animal learning and memory. The kinase domain of CaMKII is connected by a flexible linker to a C‐terminal hub domain that assembles into a 12‐ or 14‐subunit scaffold that displays the kinase domains around it. Studies on CaMKII suggest that the stoichiometry and dynamic assembly/disassembly of hub oligomers may be important for CaMKII regulation. Although CaMKII is a metazoan protein, genes encoding predicted CaMKII‐like hub domains, without associated kinase domains, are found in the genomes of some green plants and bacteria. We show that the hub domains encoded by three related green algae, Chlamydomonas reinhardtii, Volvox carteri f. nagarensis, and Gonium pectoral, assemble into 16‐, 18‐, and 20‐subunit oligomers, as assayed by native protein mass spectrometry. These are the largest known CaMKII hub domain assemblies. A crystal structure of the hub domain from C. reinhardtii reveals an 18‐subunit organization. We identified four intra‐subunit hydrogen bonds in the core of the fold that are present in the Chlamydomonas hub domain, but not in metazoan hubs. When six point mutations designed to recapitulate these hydrogen bonds were introduced into the human CaMKII‐α hub domain, the mutant protein formed assemblies with 14 and 16 subunits, instead of the normal 12‐ and 14‐subunit assemblies. Our results show that the stoichiometric balance of CaMKII hub assemblies can be shifted readily by small changes in sequence.     

Identification of hub genes and pathways in adrenocortical carcinoma by integrated bioinformatic analysis

Adrenocortical carcinoma (ACC), a rare malignant neoplasm originating from adrenal cortical cells, has high malignancy and few treatments. Therefore, it is necessary to explore the molecular mechanism of tumorigenesis, screen and verify potential biomarkers, which will provide new clues for the treatment and diagnosis of ACC. In this paper, three gene expression profiles (GSE10927, GSE12368 and GSE90713) were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were obtained using the Limma package. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched by DAVID. Protein‐protein interaction (PPI) network was evaluated by STRING database, and PPI network was constructed by Cytoscape. Finally, GEPIA was used to validate hub genes’ expression. Compared with normal adrenal tissues, 74 up‐regulated DEGs and 126 down‐regulated DEGs were found in ACC samples; GO analysis showed that up‐regulated DEGs were enriched in organelle fission, nuclear division, spindle, et al, while down‐regulated DEGs were enriched in angiogenesis, proteinaceous extracellular matrix and growth factor activity; KEGG pathway analysis showed that up‐regulated DEGs were significantly enriched in cell cycle, cellular senescence and progesterone‐mediated oocyte maturation; Nine hub genes (CCNB1, CDK1, TOP2A, CCNA2, CDKN3, MAD2L1, RACGAP1, BUB1 and CCNB2) were identified by PPI network; ACC patients with high expression of 9 hub genes were all associated with worse overall survival (OS). These hub genes and pathways might be involved in the tumorigenesis, which will offer the opportunities to develop the new therapeutic targets of ACC.     

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.