Screening of mRNA markers in early bovine tuberculosis blood samples

Bovine tuberculosis (bTB) is a chronic zoonotic disease caused by Mycobacterium bovis. A large number of cattle are infected with bTB every year, resulting in huge economic losses. How to control bTB is an important issue in the current global livestock economy. In this study, the original transcriptome sequences related to this study were obtained from the dataset GSE192537 by searching the Gene Expression Omnibus (GEO) database. Our differential gene analysis showed that there were obvious biological activities related to immune activation and immune regulation in the early stage of bTB. Immune-related biological processes were more active in the early stage of bTB than in the late. There were obvious immune activation and immune cell recruitment in the early stage of bTB. Regulations in immune receptors are associated with pathophysiological processes of the early stage of bTB. A gene module consisting of 236 genes significantly related to the early stage of bTB was obtained by weighted gene co-expression network analysis, and 18 hub genes were further identified as potential biomarkers or therapeutic targets. Finally, by random forest algorithm and logistic regression modeling, FCRL1 was identified as a representative mRNA marker in early bTB blood. FCRL1 has the potential to be a diagnostic biomarker in early bTB.

Molecular mechanism of salinity and waterlogging tolerance in mangrove Kandelia obovata

Mangrove forests are colloquially referred to as "Earth's kidneys" and serve many important ecological and commercial functions. Salinity and waterlogging stress are the most important abiotic stressors restricting the growth and development of mangroves. Kandelia obovata (K. obovata) is the greatest latitudinally-distributed salt mangrove species in China.Here, morphology and transcriptomics were used to study the response of K. obovata to salt and waterlogging stress. In addition, weighted gene co-expression network analysis of the combined gene expression and phenotypic datasets was used to identify core salinity- and waterlogging-responsive modules. In this study, we observed that both high salinity and waterlogging significantly inhibited growth and development in K. obovata. Notably, growth was negatively correlated with salt concentration and positively correlated with waterlogging duration, and high salinity was significantly more inhibitive than waterlogging. A total of 7, 591 salt-responsive and 228 waterlogging-responsive differentially expressed genes were identified by RNA sequencing. Long-term salt stress was highly correlated with the measured physiological parameters while long-term waterlogging was poorly correlated with these traits. At the same time, 45 salinity-responsive and 16 waterlogging-responsive core genes were identified. All 61 core genes were mainly involved in metabolic and biosynthesis of secondary metabolites pathways. This study provides valuable insight into the molecular mechanisms of salinity and waterlogging tolerance in K. obovata, as well as a useful genetic resource for the improvement of mangrove stress tolerance using molecular breeding techniques.

Hub Gene Mining and Co-Expression Network Construction of Low-Temperature Response in Maize of Seedling by WGCNA

Weighted gene co-expression network analysis (WGCNA) is a research method in systematic biology. It is widely used to identify gene modules related to target traits in multi-sample transcriptome data. In order to further explore the molecular mechanism of maize response to low-temperature stress at the seedling stage, B144 (cold stress tolerant) and Q319 (cold stress sensitive) provided by the Maize Research Institute of Heilongjiang Academy of Agricultural Sciences were used as experimental materials, and both inbred lines were treated with 5 °C for 0 h, 12 h, and 24 h, with the untreated material as a control. Eighteen leaf samples were used for transcriptome sequencing, with three biological replicates. Based on the above transcriptome data, co-expression networks of weighted genes associated with low-temperature-tolerance traits were constructed by WGCNA. Twelve gene modules significantly related to low-temperature tolerance at the seedling stage were obtained, and a number of hub genes involved in low-temperature stress regulation pathways were discovered from the four modules with the highest correlation with target traits. These results provide clues for further study on the molecular genetic mechanisms of low-temperature tolerance in maize at the seedling stage.

Gene co-expression network identifies critical genes, pathways and regulatory motifs mediating the progression of rift valley fever in Bostaurus

Rift Valley Fever (RVF) is a mosquito-borne viral disease caused by the Rift Valley Fever Virus. The disease is a zoonosis that largely affects domestic animals, including sheep, goats, and cattle, resulting in severe morbidity and mortality marked by massive storm abortions. To halt human and livestock deaths due to RVF, the development of efficacious vaccines and therapeutics is a compelling and urgent priority. We sought to identify potential key modules (gene clusters), hub genes, and regulatory motifs involved in the pathogenesis of RVF in Bos taurus that are amenable to inhibition. We analyzed 39 Bos taurus RNA-Seq samples using the weighted gene co-expression network analysis (WGCNA) R package and uncovered significantly enriched modules containing genes with potential pivotal roles in RVF progression. Moreover, regulatory motif analysis conducted using the Multiple Expectation Maximization for Motif Elicitation (MEME) suite identified motifs that probably modulate vital biological processes. Gene ontology terms associated with identified motifs were inferred using the GoMo human database. The gene co-expression network constructed in WGCNA using 5000 genes contained seven (7) modules, out of which four were significantly enriched for terms associated with response to viruses, response to interferon-alpha, innate immune response, and viral defense. Additionally, several biological pathways implicated in developmental processes, anatomical structure development, and multicellular organism development were identified. Regulatory motifs analysis identified short, repeated motifs whose function(s) may be amenable to disruption by novel therapeutics. Predicted functions of identified motifs include tissue development, embryonic organ development, and organ morphogenesis. We have identified several hub genes in enriched co-expressed gene modules and regulatory motifs potentially involved in the pathogenesis of RVF in B. taurus that are likely viable targets for disruption by novel therapeutics.

An epithelial-mesenchymal transition-related prognostic model for colorectal cancer based on weighted gene co-expression network analysis

OBJECTIVE

To identify susceptibility modules and genes for colorectal cancer (CRC) using weighted gene co-expression network analysis (WGCNA).

METHODS

Four microarray datasets were downloaded from the Gene Expression Omnibus database. We divided the tumor samples into three subgroups based on consensus clustering of gene expression, and analyzed the correlations between the subgroups and clinical features. The genetic features of the subgroups were investigated by gene set enrichment analysis (GSEA). A gene expression network was constructed using WGCNA, and a protein-protein interaction (PPI) network was used to identify the key genes. Gene modules were annotated by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses.

RESULTS

We divided the cancer cases into three subgroups based on consensus clustering (subgroups I, II, III). The green module identified by WGCNA was correlated with clinical characteristics. Ten key genes were identified according to their degree of connectivity in the protein-protein interaction network: FYN, SEMA3A, AP2M1, L1CAM, NRP1, TLN1, VWF, ITGB3, ILK, and ACTN1.

CONCLUSION

We identified 10 hub genes as candidate biomarkers for CRC. These key genes may provide a theoretical basis for targeted therapy against CRC.

Integrated physiological and weighted gene co-expression network analysis reveals the hub genes engaged in nitrate-regulated alleviation of ammonium toxicity at the seedling stage in wheat (Triticum aestivum L.)

Wheat has a specific preference for NO₃ ^- and shows toxicity symptoms under high NH₄ ⁺ concentrations. Increasing the nitrate supply may alleviate ammonium stress. Nevertheless, the mechanisms underlying the nitrate regulation of wheat root growth to alleviate ammonium toxicity remain unclear. In this study, we integrated physiological and weighted gene co-expression network analysis (WGCNA) to identify the hub genes involved in nitrate alleviation of ammonium toxicity at the wheat seedling stage. Five NH₄ ⁺/NO₃ ^- ratio treatments, including 100/0 (N_a), 75/25 (N_r1), 50/50 (N_r2), 25/75 (N_r3), and 0/100 (N_n) were tested in this study. The results showed that sole ammonium treatment (N_a) increased the lateral root number but reduced root biomass. Increasing the nitrate supply significantly increased the root biomass. Increasing nitrate levels decreased abscisic acid (ABA) content and increased auxin (IAA) content. Furthermore, we identified two modules (blue and turquoise) using transcriptome data that were significantly related to root physiological growth indicators. TraesCS6A02G178000 and TraesCS2B02G056300 were identified as hub genes in the two modules which coded for plastidic ATP/ADP-transporter and WRKY62 transcription factors, respectively. Additionally, network analysis showed that in the blue module, TraesCS6A02G178000 interacts with downregulated genes that coded for indolin-2-one monooxygenase, SRG1, DETOXIFICATION, and wall-associated receptor kinase. In the turquoise module, TraesCS2B02G056300 was highly related to the genes that encoded ERD4, ERF109, CIGR2, and WD40 proteins, and transcription factors including WRKY24, WRKY22, MYB30, and JAMYB, which were all upregulated by increasing nitrate supply. These studies suggest that increasing the nitrate supply could improve root growth and alleviate ammonium toxicity through physiological and molecular regulation networks, including ROS, hormonal crosstalk, and transcription factors.

Screening of key genes associated with m6A methylation in diabetic nephropathy patients by CIBERSORT and weighted gene coexpression network analysis

Diabetic nephropathy (DN) is a common complication of diabetes. Due to its complex pathogenesis, there is no effective treatment. M6A is a newly discovered epigenetic mechanism that may be involved in the development of diabetic nephropathy. In this study, we analyzed differentially expressed genes (DEG) in the GEO database (GSE96804) and paid attention to genes with m6A methylation. 623 DEGs in glomerular tissue were identified by comparing diabetic nephropathy with normal. Correlation analysis with 21 genes involved in m6A modification showed that 492 genes were associated with m6A methylation. According to the CIBERSORT algorithm, the infiltration of M1 macrophages in DN patients was significantly higher than that in normal samples. Weighted gene coexpression network analysis (WGCNA) was used to screen for the modules most correlated with the clinical features of M1 macrophages. The genes in the selected modules and 492 m6A-related DEGs were intersected by a Venn diagram, and 43 key genes were obtained. GO and KEGG analyses showed that these genes were mainly related to the positive regulation of protein aggregation and the transforming growth factor β receptor signaling pathway. According to a literature review, among the top 10 genes, HSPA1A, HSPA1B, CHI3L1, TYRO3 and PTH1R are markers in diabetic nephropathy, and their abnormal expression is associated with renal hypertrophy, proteinuria and glomerulosclerosis. These findings may provide evidence for the diagnosis and treatment of diabetic nephropathy.

Identification of Transcriptional Pattern Related to Immune Cell Infiltration With Gene Co-Expression Network in Papillary Thyroid Cancer

BACKGROUND

A growing body of evidence suggests that immune cell infiltration in cancer is closely related to clinical outcomes. However, there is still a lack of research on papillary thyroid cancer (PTC).

METHODS

Based on single-sample gene set enrichment analysis (SSGSEA) algorithm and weighted gene co-expression network analysis (WGCNA) tool, the infiltration level of immune cell and key modules and genes associated with the level of immune cell infiltration were identified using PTC gene expression data from The Cancer Genome Atlas (TCGA) database. In addition, the co-expression network and protein-protein interactions network analysis were used to identify the hub genes. Moreover, the immunological and clinical characteristics of these hub genes were verified in TCGA and GSE35570 datasets and quantitative real-time polymerase chain reaction (PCR). Finally, receiver operating characteristic (ROC) curve analysis was used to evaluate the diagnostic value of hub genes.

RESULTS

Activated B cell, activated dendritic cell, CD56bright natural killer cell, CD56dim natural killer cell, Eosinophil, Gamma delta T cell, Immature dendritic cell, Macrophage, Mast cell, Monocyte, Natural killer cell, Neutrophil and Type 17 T helper cell were significantly changed between PTC and adjacent normal groups. WGCNA results showed that the black model had the highest correlation with the infiltration level of activated dendritic cells. We found 14 hub genes whose expression correlated to the infiltration level of activated dendritic cells in both TCGA and GSE35570 datasets. After validation in the TCGA dataset, the expression level of only 5 genes (C1QA, HCK, HLA-DRA, ITGB2 and TYROBP) in 14 hub genes were differentially expressed between PTC and adjacent normal groups. Meanwhile, the expression levels of these 5 hub genes were successfully validated in GSE35570 dataset. Quantitative real-time PCR results showed the expression of these 4 hub genes (except C1QA) was consistent with the results in TCGA and GSE35570 dataset. Finally, these 4 hub genes had diagnostic value to distinguish PTC and adjacent normal controls.

CONCLUSIONS

HCK, HLA-DRA, ITGB2 and TYROBP may be key diagnostic biomarkers and immunotherapy targets in PTC.

Optimizing weighted gene co-expression network analysis with a multi-threaded calculation of the topological overlap matrix

Biomolecular networks are often assumed to be scale-free hierarchical networks. The weighted gene co-expression network analysis (WGCNA) treats gene co-expression networks as undirected scale-free hierarchical weighted networks. The WGCNA R software package uses an Adjacency Matrix to store a network, next calculates the topological overlap matrix (TOM), and then identifies the modules (sub-networks), where each module is assumed to be associated with a certain biological function. The most time-consuming step of WGCNA is to calculate TOM from the Adjacency Matrix in a single thread. In this paper, the single-threaded algorithm of the TOM has been changed into a multi-threaded algorithm (the parameters are the default values of WGCNA). In the multi-threaded algorithm, Rcpp was used to make R call a C++ function, and then C++ used OpenMP to start multiple threads to calculate TOM from the Adjacency Matrix. On shared-memory MultiProcessor systems, the calculation time decreases as the number of CPU cores increases. The algorithm of this paper can promote the application of WGCNA on large data sets, and help other research fields to identify sub-networks in undirected scale-free hierarchical weighted networks. The source codes and usage are available at https://github.com/do-somethings-haha/multi-threaded_calculate_unsigned_TOM_from_unsigned_or_signed_Adjacency_Matrix_of_WGCNA.

Integrated Network Analysis to Identify Key Modules and Potential Hub Genes Involved in Bovine Respiratory Disease: A Systems Biology Approach

Background: Bovine respiratory disease (BRD) is the most common disease in the beef and dairy cattle industry. BRD is a multifactorial disease resulting from the interaction between environmental stressors and infectious agents. However, the molecular mechanisms underlying BRD are not fully understood yet. Therefore, this study aimed to use a systems biology approach to systematically evaluate this disorder to better understand the molecular mechanisms responsible for BRD. Methods: Previously published RNA-seq data from whole blood of 18 healthy and 25 BRD samples were downloaded from the Gene Expression Omnibus (GEO) and then analyzed. Next, two distinct methods of weighted gene coexpression network analysis (WGCNA), i.e., module-trait relationships (MTRs) and module preservation (MP) analysis were used to identify significant highly correlated modules with clinical traits of BRD and non-preserved modules between healthy and BRD samples, respectively. After identifying respective modules by the two mentioned methods of WGCNA, functional enrichment analysis was performed to extract the modules that are biologically related to BRD. Gene coexpression networks based on the hub genes from the candidate modules were then integrated with protein-protein interaction (PPI) networks to identify hub-hub genes and potential transcription factors (TFs). Results: Four significant highly correlated modules with clinical traits of BRD as well as 29 non-preserved modules were identified by MTRs and MP methods, respectively. Among them, two significant highly correlated modules (identified by MTRs) and six nonpreserved modules (identified by MP) were biologically associated with immune response, pulmonary inflammation, and pathogenesis of BRD. After aggregation of gene coexpression networks based on the hub genes with PPI networks, a total of 307 hub-hub genes were identified in the eight candidate modules. Interestingly, most of these hub-hub genes were reported to play an important role in the immune response and BRD pathogenesis. Among the eight candidate modules, the turquoise (identified by MTRs) and purple (identified by MP) modules were highly biologically enriched in BRD. Moreover, STAT1, STAT2, STAT3, IRF7, and IRF9 TFs were suggested to play an important role in the immune system during BRD by regulating the coexpressed genes of these modules. Additionally, a gene set containing several hub-hub genes was identified in the eight candidate modules, such as TLR2, TLR4, IL10, SOCS3, GZMB, ANXA1, ANXA5, PTEN, SGK1, IFI6, ISG15, MX1, MX2, OAS2, IFIH1, DDX58, DHX58, RSAD2, IFI44, IFI44L, EIF2AK2, ISG20, IFIT5, IFITM3, OAS1Y, HERC5, and PRF1, which are potentially critical during infection with agents of bovine respiratory disease complex (BRDC). Conclusion: This study not only helps us to better understand the molecular mechanisms responsible for BRD but also suggested eight candidate modules along with several promising hub-hub genes as diagnosis biomarkers and therapeutic targets for BRD.

Integrated Analysis of lncRNAs, mRNAs, and TFs to Identify Regulatory Networks Underlying MAP Infection in Cattle

Johne’s disease is a chronic infection of ruminants that burdens dairy herds with a significant economic loss. The pathogenesis of the disease has not been revealed clearly due to its complex nature. In order to achieve deeper biological insights into molecular mechanisms involved in MAP infection resulting in Johne’s disease, a system biology approach was used. As far as is known, this is the first study that considers lncRNAs, TFs, and mRNAs, simultaneously, to construct an integrated gene regulatory network involved in MAP infection. Weighted gene coexpression network analysis (WGCNA) and functional enrichment analysis were conducted to explore coexpression modules from which nonpreserved modules had altered connectivity patterns. After identification of hub and hub-hub genes as well as TFs and lncRNAs in the nonpreserved modules, integrated networks of lncRNA-mRNA-TF were constructed, and cis and trans targets of lncRNAs were identified. Both cis and trans targets of lncRNAs were found in eight nonpreserved modules. Twenty-one of 47 nonpreserved modules showed significant biological processes related to the immune system and MAP infection. Some of the MAP infection’s related pathways in the most important nonpreserved modules comprise “positive regulation of cytokine-mediated signaling pathway,” “negative regulation of leukocyte migration,” “T-cell differentiation,” “neutrophil activation,” and “defense response.” Furthermore, several genes were identified in these modules, including SLC11A1, MAPK8IP1, HMGCR, IFNGR1, CMPK2, CORO1A, IRF1, LDLR, BOLA-DMB, and BOLA-DMA, which are potentially associated with MAP pathogenesis. This study not only enhanced our knowledge of molecular mechanisms behind MAP infection but also highlighted several promising hub and hub-hub genes involved in macrophage-pathogen interaction.

Weighted Gene Co-Expression Network Analysis Reveals Key Genes and Potential Drugs in Abdominal Aortic Aneurysm

Abdominal aortic aneurysm (AAA) is a prevalent aortic disease that causes high mortality due to asymptomatic gradual expansion and sudden rupture. The underlying molecular mechanisms and effective pharmaceutical therapy for preventing AAA progression have not been fully identified. In this study, we identified the key modules and hub genes involved in AAA growth from the GSE17901 dataset in the Gene Expression Omnibus (GEO) database through the weighted gene co-expression network analysis (WGCNA). Key genes were further selected and validated in the mouse dataset (GSE12591) and human datasets (GSE7084, GSE47472, and GSE57691). Finally, we predicted drug candidates targeting key genes using the Drug-Gene Interaction database. Overall, we identified key modules enriched in the mitotic cell cycle, GTPase activity, and several metabolic processes. Seven key genes (CCR5, ADCY5, ADCY3, ACACB, LPIN1, ACSL1, UCP3) related to AAA progression were identified. A total of 35 drugs/compounds targeting the key genes were predicted, which may have the potential to prevent AAA progression.

Identification of Key Modules and Hub Genes of Annulus Fibrosus in Intervertebral Disc Degeneration

Background: Intervertebral disc degeneration impairs the quality of patients lives. Even though there has been development of many therapeutic strategies, most of them remain unsatisfactory due to the limited understanding of the mechanisms that underlie the intervertebral disc degeneration. Questions/purposes: This study is meant to identify the key modules and hub genes related to the annulus fibrosus in intervertebral disc degeneration (IDD) through: (1) constructing a weighted gene co-expression network; (2) identifying key modules and hub genes; (3) verifying the relationships of key modules and hub genes with IDD; and (4) confirming the expression pattern of hub genes in clinical samples. Methods: The Gene Expression Omnibus provided 24 sets of annulus fibrosus microarray data. Differentially expressed genes between the annulus fibrosus of degenerative and non-degenerative intervertebral disc samples have gone through the Gene Ontology (GO) and pathway analysis. The construction of a gene network and classification of genes into different modules were conducted through performing Weighted Gene Co-expression Network Analysis. The identification of modules and hub genes that were most related to intervertebral disc degeneration was proceeded. In order to verify the relationships of the module and hub genes with intervertebral disc degeneration, Ingenuity Pathway Analysis was operated. Clinical samples were adopted to help verify the hub gene expression profile. Results: One thousand one hundred ninety differentially expressed genes were identified. Terms and pathways associated with intervertebral disc degeneration were presented by GO and pathway analysis. The construction of a Weighted Gene Coexpression Network was completed and clustering differentially expressed genes into four modules was also achieved. The module with the lowest P-value and the highest absolute correlation coefficient was selected and its relationship with intervertebral disc degeneration was confirmed by Ingenuity Pathway Analysis. The identification of hub genes and the confirmation of their expression profile were also realized. Conclusions: This study generated a comprehensive overview of the gene networks underlying annulus fibrosus in intervertebral disc degeneration. Clinical Relevance: Modules and hub genes identified in this study are highly associated with intervertebral disc degeneration, and may serve as potential therapeutic targets for intervertebral disc degeneration.

Time Series Gene Expression Profiles Analysis Identified Several Potential Biomarkers for Sepsis

Sepsis is a life-threatening disorder and leads to organ dysfunction and death. Therefore, searching for more alternative biomarkers is of great significance for sepsis assessment and surveillance. In our study, the gene expression profiles of 163 samples from healthy controls and septic patients were analyzed and 8 gene co-expression modules were identified by constructing weighted gene co-expression network. The blue and yellow modules showed close correlations with the phenotypic trait "days postsepsis." Besides, differentially expressed genes (DEGs) over time in septic patients were screened using Short Time-series Expression Miner (STEM) program. The intersection of genes in the blue and yellow modules and DEGs, which were significantly enriched in "HTLV-1 infection" pathway, was analyzed with protein-protein interaction network. The logistic regression model based on these eight mRNAs was constructed to determine the type of the sample reliably. Eight vital genes CECR1, ANXA2, ELANE, CTSG, AZU1, PRTN3, LYZ, and DEFA4 presented high scores and may be associated with sepsis, which provided candidate biomarkers for sepsis.

Identification of Gene Modules and Hub Genes Involved in Mastitis Development Using a Systems Biology Approach

Objective

Mastitis is defined as the inflammation of the mammary gland, which impact directly on the production performance and welfare of dairy cattle. Since, mastitis is a multifactorial complex disease and the molecular pathways underlying this disorder have not been clearly understood yet, a system biology approach was used in this study to a better understanding of the molecular mechanisms behind mastitis.

Methods

Publicly available RNA-Seq data containing samples from milk of five infected and five healthy Holstein cows at five time points were retrieved. Gene Co-expression network analysis (WGCNA) approach and functional enrichment analysis were then applied with the aim to find the non-preserved module of genes that their connectivity were altered under infected condition. Hub genes were identified in the non-preserved modules and were subjected to protein-protein interactions (PPI) network construction.

Results

Among the 25 modules identified, eight modules were non-preserved and were also biologically associated with inflammation, immune response and mastitis development. Interestingly most of the hub genes in the eight modules were also densely connected in the PPI network. Of the hub genes, 250 genes were hubs in both co-expression and PPI networks and most of them were reported to play important roles in immune response or inflammatory pathways. The blue module was highly enriched in inflammatory responses and STAT1 was suggested to play an important role in mastitis development by regulating the immune related genes in this module. Moreover, a set of highly connected genes were identified such as BIRC3, PSMA6, FYN, F11R, NFKBIZ, NFKBIA, GRO1, PHB, CD3E, IL16, GSN, SOCS2, HCK, VAV1 and TLR6, which have been established to be critical for mastitis pathogenesis.

Conclusion

This study improved the understanding of the mechanisms underlying bovine mastitis and suggested eight non-preserved modules along with several most important genes with promising potential in etiology of mastitis.

Identification of DDX58 and CXCL10 as Potential Biomarkers in Acute Respiratory Distress Syndrome

Acute respiratory distress syndrome (ARDS) is a devastating condition of acute inflammatory lung injury and causes high morbidity and mortality. Therefore, investigations on the effective biomarkers will be significant for the understanding of ARDS. In our research, the gene expression profiles of 27 samples from ARDS patients (n = 18) and healthy controls (n = 9) were analyzed and eight gene co-expression modules were identified by constructing weighted gene co-expression network. The correlation analysis of modules with phenotypes showed that genes in the yellow and black modules, which were significantly enriched in the ARDS-related pathways, such as TNF signaling pathway, Toll-like receptor signaling pathway, and NF-kappa B signaling pathway, were associated with the phenotype "time postinfection." Genes DDX58 and CXCL10, which were highly expressed after infection and significantly enriched in ARDS-related pathways, presented high score in protein-protein interaction analysis, indicating that they may be associated with ARDS and providing novel biomarkers for its diagnosis, treatment, and surveillance.

Weighted Gene Co-Expression Network Analysis Identifies Critical Genes in the Development of Heart Failure After Acute Myocardial Infarction

Background: The development of heart failure (HF) remains a common complication following an acute myocardial infarction (AMI), and is associated with substantial adverse outcomes. However, the specific predictive biomarkers and candidate therapeutic targets for post-infarction HF have not been fully established. We sought to perform a weighted gene co-expression network analysis (WGCNA) to identify key modules, hub genes, and possible regulatory targets involved in the development of HF following AMI.

Methods: Genes exhibiting the most (top 50%) variation in expression levels across samples in a GSE59867 dataset were imported to the WGCNA. Gene Ontology and pathway enrichment analyses were performed on genes identified in the key module by Metascape. Gene regulatory networks were constructed using the microarray probe reannotation and bioinformatics database. Hub genes were screened out from the key module and validated using other datasets.

Results: A total of 10,265 most varied genes and six modules were identified between AMI patients who developed HF within 6 months of follow-up and those who did not. Specifically, the blue module was found to be the most significantly related to the development of post-infarction HF. Functional enrichment analysis revealed that the blue module was primarily associated with the inflammatory response, immune system, and apoptosis. Seven transcriptional factors, including SPI1, ZBTB7A, IRF8, PPARG, P65, KLF4, and Fos, were identified as potential regulators of the expression of genes identified in the blue module. Further, non-coding RNAs, including miR-142-3p and LINC00537, were identified as having close interactions with genes from the blue module. A total of six hub genes (BCL3, HCK, PPIF, S100A9, SERPINA1, and TBC1D9B) were identified and validated for their predictive value in identifying future HFs.

Conclusions: By using the WGCNA, we provide new insights into the underlying molecular mechanism and molecular markers correlated with HF development following an AMI, which may serve to improve risk stratification, therapeutic decisions, and prognosis prediction in AMI patients.

Identification of key genes involved in the development and progression of early-onset colorectal cancer by co-expression network analysis

A number of studies have revealed that there is an increasing incidence of early-onset colorectal cancer (CRC) in young adults (before the age of 50 years) and a progressive decline in CRC among older patients, after the age of 50 years (late-onset CRC). However, the etiology of early-onset CRC is not fully understood. The aim of the present study was to identify key genes associated with the development of early-onset CRC through weighted gene co-expression network analysis (WGCNA). The GSE39582 dataset was downloaded from the Gene Expression Omnibus database, and the data profiles of tissues from patients diagnosed before the age of 50 years were selected. The top 10,000 genes with the highest variability were used to construct the WGCNA. Hub genes were identified from the modules associated with clinical traits using gene significance >0.2 and module membership >0.8 as the cut-off criteria. Gene Ontology and pathway analyses were subsequently performed on the hub genes and a protein-protein interaction network (PPI) was constructed. The diagnostic value of module hub genes with a degree score >5 in the PPI network was verified in samples from patients with CRC diagnosed before the age of 50 years obtained from The Cancer Genome Atlas. Eight co-expressed gene modules were identified in the WGCNA and two modules (blue and turquoise) were associated with the tumor-node-metastasis stage. A total of 140 module hub genes were identified and found to be enriched in 'mitochondrial large ribosomal subunit', 'structural constituent of ribosome', 'poly (A) RNA binding', 'collagen binding', 'protein ubiquitination' and 'ribosome pathway'. Twenty-six module hub genes were found to have a degree score >5 in the PPI network, seven of which [secreted protein acidic and cysteine rich (SPARC), decorin (DCN), fibrillin 1 (FBN1), WW domain containing transcription regulator 1 (WWTR1), transgelin (TAGLN), DEAD-box helicase 28 (DDX28) and cold shock domain containing C2 (CSDC2)], had good prognostic values for patients with early-onset CRC, but not late-onset CRC. Therefore, SPARC, DCN, FBN1, WWTR1, TAGLN, DDX28 and CSDC2 may contribute to the development of early-onset CRC and may serve as potential diagnostic biomarkers.

Identification of hub genes and construction of transcriptional regulatory network for the progression of colon adenocarcinoma hub genes and TF regulatory network of colon adenocarcinoma

The aim of this study was to identify key genes related to the progression of colon adenocarcinoma (COAD), and to investigate the regulatory network of hub genes and transcription factors (TFs). Dataset GSE20916 including 44 normal colon, 55 adenoma, and 36 adenocarcinoma tissue samples was used to construct co‐expression networks via weighted gene co‐expression network. Gene Ontology annotation and the Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis for the objective module were performed using the online Database for Annotation, Visualization and Integrated Discovery. Hub genes were identified by taking the intersection of differentially expressed genes between dataset GSE20916 and GSE39582 and validated using The Cancer Genome Atlas (TCGA) database. The correlations between microRNA (miRNA) and hub genes were analyzed using the online website StarBase. Cytoscape was used to establish a regulatory network of TF‐miRNA‐target gene. We found that the orange module was a key module related to the tumor progression in COAD. In datasets GSE20916 and GSE39582, a total of eight genes (BGN, SULF1, COL1A1, FAP, THBS2, CTHRC1, COL5A2, and COL1A2) were selected, which were closely related with patients’ survivals in TCGA database and dataset GSE20916. COAD patients with higher expressions of each hub gene had a worse prognosis than those with lower expressions. A regulatory network of TF‐miRNA‐target gene with 144 TFs, 26 miRNAs, and 7 hub genes was established, including model KLF11‐miR149‐BGN, TCEAL6‐miR29B2‐COL1A1, and TCEAL6‐miR29B2‐COL1A2. In conclusion, during the progression of COAD, eight core genes (BGN, SULF1, COL1A1, FAP, THBS2, CTHRC1, COL5A2, and COL1A2) play vital roles. Regulatory networks of TF‐miRNA‐target gene can help to understand the disease progression and optimize treatment strategy.

Identification of an eight-lncRNA prognostic model for breast cancer using WGCNA network analysis and a Cox‑proportional hazards model based on L1-penalized estimation

An ever‑increasing number of long noncoding (lnc)RNAs has been identified in breast cancer. The present study aimed to establish an lncRNA signature for predicting survival in breast cancer. RNA expression profiling was performed using microarray gene expression data from the National Center for Biotechnology Information Gene Expression Omnibus, followed by the identification of breast cancer‑related preserved modules using weighted gene co‑expression network (WGCNA) network analysis. From the lncRNAs identified in these preserved modules, prognostic lncRNAs were selected using univariate Cox regression analysis in combination with the L1‑penalized (LASSO) Cox‑proportional Hazards (Cox‑PH) model. A risk score based on these prognostic lncRNAs was calculated and used for risk stratification. Differentially expressed RNAs (DERs) in breast cancer were identified using MetaDE. Gene Set Enrichment Analysis pathway enrichment analysis was conducted for these prognostic lncRNAs and the DERs related to the lncRNAs in the preserved modules. A total of five preserved modules comprising 73 lncRNAs were mined. An eight‑lncRNA signature (IGHA1, IGHGP, IGKV2‑28, IGLL3P, IGLV3‑10, AZGP1P1, LINC00472 and SLC16A6P1) was identified using the LASSO Cox‑PH model. Risk score based on these eight lncRNAs could classify breast cancer patients into two groups with significantly different survival times. The eight‑lncRNA signature was validated using three independent cohorts. These prognostic lncRNAs were significantly associated with the cell adhesion molecules pathway, JAK‑signal transducer and activator of transcription 5A pathway, and erbb pathway and are potentially involved in regulating angiotensin II receptor type 1, neuropeptide Y receptor Y1, KISS1 receptor, and C‑C motif chemokine ligand 5. The developed eight‑lncRNA signature may have clinical implications for predicting prognosis in breast cancer. Overall, this study provided possible molecular targets for the development of novel therapies against breast cancer.

Weighted Gene Co-expression Network Analysis of Endometriosis and Identification of Functional Modules Associated With Its Main Hallmarks

Although many genes have been identified using high throughput technologies in endometriosis (ES), only a small number of individual genes have been analyzed functionally. This is due to the complexity of the disease that has different stages and is affected by various genetic and environmental factors. Many genes are upregulated or downregulated at each stage of the disease, thus making it difficult to identify key genes. In addition, little is known about the differences between the different stages of the disease. We assumed that the study of the identified genes in ES at a system-level can help to better understand the molecular mechanism of the disease at different stages of the development. We used publicly available microarray data containing archived endometrial samples from women with minimal/mild endometriosis (MMES), mild/severe endometriosis (MSES) and without endometriosis. Using weighted gene co-expression analysis (WGCNA), functional modules were derived from normal endometrium (NEM) as the reference sample. Subsequently, we tested whether the topology or connectivity pattern of the modules was preserved in MMES and/or MSES. Common and specific hub genes were identified in non-preserved modules. Accordingly, hub genes were detected in the non-preserved modules at each stage. We identified sixteen co-expression modules. Of the 16 modules, nine were non-preserved in both MMES and MSES whereas five were preserved in NEM, MMES, and MSES. Importantly, two non-preserved modules were found in either MMES or MSES, highlighting differences between the two stages of the disease. Analyzing the hub genes in the non-preserved modules showed that they mostly lost or gained their centrality in NEM after developing the disease into MMES and MSES. The same scenario was observed, when the severeness of the disease switched from MMES to MSES. Interestingly, the expression analysis of the new selected gene candidates including CC2D2A, AEBP1, HOXB6, IER3, and STX18 as well as IGF-1, CYP11A1 and MMP-2 could validate such shifts between different stages. The overrepresented gene ontology (GO) terms were enriched in specific modules, such as genetic disposition, estrogen dependence, progesterone resistance and inflammation, which are known as endometriosis hallmarks. Some modules uncovered novel co-expressed gene clusters that were not previously discovered.

Co-expression network analysis identified CDH11 in association with progression and prognosis in gastric cancer

Background and aims

Gastric cancer (GC) is one of the most common cancers worldwide, and its pathogenesis is related to a complex network of gene interactions. The aims of our study were to find hub genes associated with the progression and prognosis of GC and illustrate the underlying mechanisms.

Methods

Weighted gene co-expression network analysis (WGCNA) was conducted using the microarray dataset and clinical data of GC patients from Gene Expression Omnibus (GEO) database to identify significant gene modules and hub genes associated with TNM stage in GC. Functional enrichment analysis and protein-protein interaction network analysis were performed using the significant module genes. We regarded the common hub genes in the co-expression network and protein-protein interaction (PPI) network as "real" hub genes for further analysis. Hub gene was validated in another independent dataset and The Cancer Genome Atlas (TCGA) dataset.

Results

In the significant purple module (R ²=0.35), a total of 12 network hub genes were identified, among which six were also hub nodes in the PPI network of the module genes. Functional annotation revealed that the genes in the purple module focused on the biological processes of system development, biological adhesion, extracellular structure organization and metabolic process. In terms of validation, CDH11 had a higher correlation with the TNM stage than other hub genes and was strongly correlated with biological adhesion based on GO functional enrichment analysis. Data obtained from the Gene Expression Profiling Interactive Analysis (GEPIA) showed that CDH11 expression had a strong positive correlation with GC stages (P<0.0001). In the testing set and Oncomine dataset, CDH11 was highly expressed in GC tissues (P<0.0001). Survival analysis indicated that samples with a high CDH11 expression showed a poor prognosis. Cox regression analysis demonstrated an independent predictor of CDH11 expression in GC prognosis (HR=1.482, 95% CI: 1.015-2.164). Furthermore, gene set enrichment analysis (GSEA) demonstrated that multiple tumor-related pathways, especially focal adhesion, were enriched in CDH11 highly expressed samples.

Conclusion

CDH11 was identified and validated in association with progression and prognosis in GC, probably by regulating biological adhesion and focal adhesion-related pathways.

Transcriptome sequencing identified hub genes for hepatocellular carcinoma by weighted-gene co-expression analysis

Hepatocellular carcinoma (HCC) is one of the most common malignancies worldwide, and it remains a challenge to understand the genetic mechanisms underlying hepatocarcinogenesis. A global gene network of differential expression profiles in HCC has yet to be fully characterized. In the present study, we performed transcriptome sequencing (mRNA and lncRNA) in liver cancer and cirrhotic tissues of nine HCC patients. We identified differentially expressed genes (DEGs) and constructed a weighted gene co-expression network for the DEGs. In total, 755 DEGs (747 mRNA and eight lncRNA) were identified, and several co-expression modules were significantly associated with HCC clinical traits, including tumor location, tumor grade, and the α-fetoprotein (AFP) level. Of note, we identified 15 hub genes in the module associated with AFP level, and three (SPX, AFP and ADGRE1) of four hub genes were validated in an independent HCC cohort (n=78). Identification of hub genes for HCC clinical traits has implications for further understanding of the molecular genetic basis of HCC.

Identification and evaluation of lncRNA and mRNA integrative modules in human peripheral blood mononuclear cells

AIM

To identify sets of functionally related long noncoding RNAs (lncRNAs) and mRNAs and to evaluate the importance of lncRNAs in an lncRNA-mRNA network.

METHODS

We carried out weighted gene co-expression network analysis and enrichment analyses to identify functional modules of co-expressed lncRNAs and mRNAs in human peripheral blood mononuclear cells of 43 females.

RESULTS

We identified seven modules and found hub lncRNAs in each module. Four of the seven modules had significant gene ontology enrichments. Some of the hub lncRNAs (e.g., SSX8, UCA1, HOXA-AS2, STARD4-AS1 and PCBP1-AS1) have known functions related with diseases such as cancers.

CONCLUSION

We identified seven biologically important lncRNA and mRNA integrative modules in females and showed that lncRNAs might play important roles in lncRNA-mRNA co-expression modules.