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

Integrative Multiomics in the Lung Reveals a Protective Role of Asporin in Pulmonary Arterial Hypertension

BACKGROUND

Integrative multiomics can elucidate pulmonary arterial hypertension (PAH) pathobiology, but procuring human PAH lung samples is rare.

METHODS

We leveraged transcriptomic profiling and deep phenotyping of the largest multicenter PAH lung biobank to date (96 disease and 52 control) by integration with clinicopathologic data, genome-wide association studies, Bayesian regulatory networks, single-cell transcriptomics, and pharmacotranscriptomics.

RESULTS

We identified 2 potentially protective gene network modules associated with vascular cells, and we validated ASPN, coding for asporin, as a key hub gene that is upregulated as a compensatory response to counteract PAH. We found that asporin is upregulated in lungs and plasma of multiple independent PAH cohorts and correlates with reduced PAH severity. We show that asporin inhibits proliferation and transforming growth factor-β/phosphorylated SMAD2/3 signaling in pulmonary artery smooth muscle cells from PAH lungs. We demonstrate in Sugen-hypoxia rats that ASPN knockdown exacerbated PAH and recombinant asporin attenuated PAH.

CONCLUSIONS

Our integrative systems biology approach to dissect the PAH lung transcriptome uncovered asporin as a novel protective target with therapeutic potential in PAH.

Rhizobium etli CFN42 and Sinorhizobium meliloti 1021 bioinformatic transcriptional regulatory networks from culture and symbiosis

Rhizobium etli CFN42 proteome-transcriptome mixed data of exponential growth and nitrogen-fixing bacteroids, as well as Sinorhizobium meliloti 1021 transcriptome data of growth and nitrogen-fixing bacteroids, were integrated into transcriptional regulatory networks (TRNs). The one-step construction network consisted of a matrix-clustering analysis of matrices of the gene profile and all matrices of the transcription factors (TFs) of their genome. The networks were constructed with the prediction of regulatory network application of the RhizoBindingSites database (http://rhizobindingsites.ccg.unam.mx/). The deduced free-living Rhizobium etli network contained 1,146 genes, including 380 TFs and 12 sigma factors. In addition, the bacteroid R. etli CFN42 network contained 884 genes, where 364 were TFs, and 12 were sigma factors, whereas the deduced free-living Sinorhizobium meliloti 1021 network contained 643 genes, where 259 were TFs and seven were sigma factors, and the bacteroid Sinorhizobium meliloti 1021 network contained 357 genes, where 210 were TFs and six were sigma factors. The similarity of these deduced condition-dependent networks and the biological E. coli and B. subtilis independent condition networks segregates from the random Erdös-Rényi networks. Deduced networks showed a low average clustering coefficient. They were not scale-free, showing a gradually diminishing hierarchy of TFs in contrast to the hierarchy role of the sigma factor rpoD in the E. coli K12 network. For rhizobia networks, partitioning the genome in the chromosome, chromids, and plasmids, where essential genes are distributed, and the symbiotic ability that is mostly coded in plasmids, may alter the structure of these deduced condition-dependent networks. It provides potential TF gen-target relationship data for constructing regulons, which are the basic units of a TRN.

Phenomics in Livestock Research: Bottlenecks and Promises of Digital Phenotyping and Other Quantification Techniques on a Global Scale

Bottlenecks in moving genomics to real-life applications also include phenomics. This is true not only for genomics medicine and public health genomics but also in ecology and livestock phenomics. This expert narrative review explores the intricate relationship between genetic makeup and observable phenotypic traits across various biological levels in the context of livestock research. We unpack and emphasize the significance of precise phenotypic data in selective breeding outcomes and examine the multifaceted applications of phenomics, ranging from improvement to assessing welfare, reproductive traits, and environmental adaptation in livestock. As phenotypic traits exhibit strong correlations, their measurement alongside specific biological outcomes provides insights into performance, overall health, and clinical endpoints like morbidity and disease. In addition, automated assessment of livestock holds potential for monitoring the dynamic phenotypic traits across various species, facilitating a deeper comprehension of how they adapt to their environment and attendant stressors. A key challenge in genetic improvement in livestock is predicting individuals with optimal fitness without direct measurement. Temporal predictions from unmanned aerial systems can surpass genomic predictions, offering in-depth data on livestock. In the near future, digital phenotyping and digital biomarkers may further unravel the genetic intricacies of stress tolerance, adaptation and welfare aspects of animals enabling the selection of climate-resilient and productive livestock. This expert review thus delves into challenges associated with phenotyping and discusses technological advancements shaping the future of biological research concerning livestock.

Transcriptomic Analysis of Hub Genes Reveals Associated Inflammatory Pathways in Estrogen-Dependent Gynecological Diseases

Gynecological diseases are triggered by aberrant molecular pathways that alter gene expression, hormonal balance, and cellular signaling pathways, which may lead to long-term physiological consequences. This study was able to identify highly preserved modules and key hub genes that are mainly associated with gynecological diseases, represented by endometriosis (EM), ovarian cancer (OC), cervical cancer (CC), and endometrial cancer (EC), through the weighted gene co-expression network analysis (WGCNA) of microarray datasets sourced from the Gene Expression Omnibus (GEO) database. Five highly preserved modules were observed across the EM (GSE51981), OC (GSE63885), CC (GSE63514), and EC (GSE17025) datasets. The functional annotation and pathway enrichment analysis revealed that the highly preserved modules were heavily involved in several inflammatory pathways that are associated with transcription dysregulation, such as NF-kB signaling, JAK-STAT signaling, MAPK-ERK signaling, and mTOR signaling pathways. Furthermore, the results also include pathways that are relevant in gynecological disease prognosis through viral infections. Mutations in the ESR1 gene that encodes for ERα, which were shown to also affect signaling pathways involved in inflammation, further indicate its importance in gynecological disease prognosis. Potential drugs were screened through the Drug Repurposing Encyclopedia (DRE) based on the up-and downregulated hub genes, wherein a bacterial ribosomal subunit inhibitor and a benzodiazepine receptor agonist were the top candidates. Other drug candidates include a dihydrofolate reductase inhibitor, glucocorticoid receptor agonists, cholinergic receptor agonists, selective serotonin reuptake inhibitors, sterol demethylase inhibitors, a bacterial antifolate, and serotonin receptor antagonist drugs which have known anti-inflammatory effects, demonstrating that the gene network highlights specific inflammatory pathways as a therapeutic avenue in designing drug candidates for gynecological diseases.

Construction and evaluation of endometriosis diagnostic prediction model and immune infiltration based on efferocytosis-related genes

Background: Endometriosis (EM) is a long-lasting inflammatory disease that is difficult to treat and prevent. Existing research indicates the significance of immune infiltration in the progression of EM. Efferocytosis has an important immunomodulatory function. However, research on the identification and clinical significance of efferocytosis-related genes (EFRGs) in EM is sparse. Methods: The EFRDEGs (differentially expressed efferocytosis-related genes) linked to datasets associated with endometriosis were thoroughly examined utilizing the Gene Expression Omnibus (GEO) and GeneCards databases. The construction of the protein-protein interaction (PPI) and transcription factor (TF) regulatory network of EFRDEGs ensued. Subsequently, machine learning techniques including Univariate logistic regression, LASSO, and SVM classification were applied to filter and pinpoint diagnostic biomarkers. To establish and assess the diagnostic model, ROC analysis, multivariate regression analysis, nomogram, and calibration curve were employed. The CIBERSORT algorithm and single-cell RNA sequencing (scRNA-seq) were employed to explore immune cell infiltration, while the Comparative Toxicogenomics Database (CTD) was utilized for the identification of potential therapeutic drugs for endometriosis. Finally, immunohistochemistry (IHC) and reverse transcription quantitative polymerase chain reaction (RT-qPCR) were utilized to quantify the expression levels of biomarkers in clinical samples of endometriosis. Results: Our findings revealed 13 EFRDEGs associated with EM, and the LASSO and SVM regression model identified six hub genes (ARG2, GAS6, C3, PROS1, CLU, and FGL2). Among these, ARG2, GAS6, and C3 were confirmed as diagnostic biomarkers through multivariate logistic regression analysis. The ROC curve analysis of GSE37837 (AUC = 0.627) and GSE6374 (AUC = 0.635), along with calibration and DCA curve assessments, demonstrated that the nomogram built on these three biomarkers exhibited a commendable predictive capacity for the disease. Notably, the ratio of nine immune cell types exhibited significant differences between eutopic and ectopic endometrial samples, with scRNA-seq highlighting M0 Macrophages, Fibroblasts, and CD8 Tex cells as the cell populations undergoing the most substantial changes in the three biomarkers. Additionally, our study predicted seven potential medications for EM. Finally, the expression levels of the three biomarkers in clinical samples were validated through RT-qPCR and IHC, consistently aligning with the results obtained from the public database. Conclusion: we identified three biomarkers and constructed a diagnostic model for EM in this study, these findings provide valuable insights for subsequent mechanistic research and clinical applications in the field of endometriosis.

CXADR promote epithelial-mesenchymal transition in endometriosis by modulating AKT/GSK-3β signaling

Endometriosis is a benign high prevalent disease exhibiting malignant features. However, the underlying pathogenesis and key molecules of endometriosis remain unclear. By integrating and analysis of existing expression profile datasets, we identified coxsackie and adenovirus receptor (CXADR), as a novel key gene in endometriosis. Based on the results of immunohistochemistry (IHC), we confirmed significant down-regulation of CXADR in ectopic endometrial tissues obtained from women with endometriosis compared with healthy controls. Further in vitro investigation indicated that CXADR regulated the stability and function of the phosphatases and AKT inhibitors PHLPP2 (pleckstrin homology domain and leucine-rich repeat protein phosphatase 2) and PTEN (phosphatase and tensin homolog). Loss of CXADR led to phosphorylation of AKT and glycogen synthase kinase-3β (GSK-3β), which resulted in stabilization of an epithelial-mesenchymal transition (EMT) factor, SNAIL1 (snail family transcriptional repressor 1). Therefore, EMT processs was induced, and the proliferation, migration and invasion of Ishikawa cells were enhanced. Over-expression of CXADR showed opposite effects. These findings suggest a previously undefined role of AKT/GSK-3β signaling axis in regulating EMT and reveal the involvement of a CXADR-induced EMT, in pathogenic progression of endometriosis.

Identifying Immune Cell Infiltration and Hub Genes Related to M2 Macrophages in Endometriosis by Bioinformatics Analysis

Endometriosis (EM) is a chronic, estrogen-dependent inflammatory disease. Presently, the pathophysiology of EM is still unclear, and numerous studies have established that the immune system plays a major role in the pathophysiology of EM. Six microarray datasets were downloaded from the GEO public database. A total of 151 endometrial samples (72 ectopic endometria and 79 controls) were included in this study. CIBERSORT and ssGSEA were applied to calculate the immune infiltration of EM and control samples. Moreover, we validated four different correlation analyses to explore immune microenvironment of EM and finally identified M2 macrophage-related hub genes and further conducted the specific immunologic signaling pathway analysis by GSEA. The logistic regression model was investigated by ROC and further validated by two external datasets. From the results of the two immune infiltration assays, we concluded that M2 macrophages, regulatory T cells (Tregs), M1 macrophages, activated B cells, T follicular helper cells, activated dendritic cells, and resting NK cells have a significant difference between control and EM tissues. Through multidimensional correlation analysis, we found that macrophages play an important central role in cell-to-cell interactions, especially M2 macrophages. Four immune-related hub genes, namely FN1, CCL2, ESR1, and OCLN, are closely related to M2 macrophages and play a crucial role in the occurrence and immune microenvironment of endometriosis. The combined AUC of ROC prediction model in test and validation sets were 0.9815 and 0.8206, respectively. We conclude that M2 macrophages play a central role in the immune-infiltrating microenvironment of EM.

Co-Expression Networks in Sunflower: Harnessing the Power of Multi-Study Transcriptomic Public Data to Identify and Categorize Candidate Genes for Fungal Resistance

Fungal plant diseases are a major threat to food security worldwide. Current efforts to identify and list loci involved in different biological processes are more complicated than originally thought, even when complete genome assemblies are available. Despite numerous experimental and computational efforts to characterize gene functions in plants, about ~40% of protein-coding genes in the model plant Arabidopsis thaliana L. are still not categorized in the Gene Ontology (GO) Biological Process (BP) annotation. In non-model organisms, such as sunflower (Helianthus annuus L.), the number of BP term annotations is far fewer, ~22%. In the current study, we performed gene co-expression network analysis using eight terabytes of public transcriptome datasets and expression-based functional prediction to categorize and identify loci involved in the response to fungal pathogens. We were able to construct a reference gene network of healthy green tissue (GreenGCN) and a gene network of healthy and stressed root tissues (RootGCN). Both networks achieved robust, high-quality scores on the metrics of guilt-by-association and selective constraints versus gene connectivity. We were able to identify eight modules enriched in defense functions, of which two out of the three modules in the RootGCN were also conserved in the GreenGCN, suggesting similar defense-related expression patterns. We identified 16 WRKY genes involved in defense related functions and 65 previously uncharacterized loci now linked to defense response. In addition, we identified and classified 122 loci previously identified within QTLs or near candidate loci reported in GWAS studies of disease resistance in sunflower linked to defense response. All in all, we have implemented a valuable strategy to better describe genes within specific biological processes.

Gene Co-Expression Network Analysis Reveals Key Regulatory Genes in Metisa plana Hormone Pathways

Metisa plana Walker (Lepidoptera: Psychidae) is a major oil palm pest species distributed across Southeast Asia. M. plana outbreaks are regarded as serious ongoing threats to the oil palm industry due to their ability to significantly reduce fruit yield and subsequent productivity. Currently, conventional pesticide overuses may harm non-target organisms and severely pollute the environment. This study aims to identify key regulatory genes involved in hormone pathways during the third instar larvae stage of M. plana gene co-expression network analysis. A weighted gene co-expression network analysis (WGCNA) was conducted on the M. plana transcriptomes to construct a gene co-expression network. The transcriptome datasets were obtained from different development stages of M. plana, i.e., egg, third instar larvae, pupa, and adult. The network was clustered using the DPClusO algorithm and validated using Fisher's exact test and receiver operating characteristic (ROC) analysis. The clustering analysis was performed on the network and 20 potential regulatory genes (such as MTA1-like, Nub, Grn, and Usp) were identified from ten top-most significant clusters. Pathway enrichment analysis was performed to identify hormone signalling pathways and these pathways were identified, i.e., hormone-mediated signalling, steroid hormone-mediated signalling, and intracellular steroid hormone receptor signalling as well as six regulatory genes Hnf4, Hr4, MED14, Usp, Tai, and Trr. These key regulatory genes have a potential as important targets in future upstream applications and validation studies in the development of biorational pesticides against M. plana and the RNA interference (RNAi) gene silencing method.

Potential miRNA-gene interactions determining progression of various ATLL cancer subtypes after infection by HTLV-1 oncovirus

BACKGROUND

Adult T-cell Leukemia/Lymphoma (ATLL) is a rapidly progressing type of T-cell non-Hodgkin lymphoma that is developed after the infection by human T-cell leukemia virus type 1 (HTLV-1). It could be categorized into four major subtypes, acute, lymphoma, chronic, and smoldering. These different subtypes have some shared clinical manifestations, and there are no trustworthy biomarkers for diagnosis of them.

METHODS

We applied weighted-gene co-expression network analysis to find the potential gene and miRNA biomarkers for various ATLL subtypes. Afterward, we found reliable miRNA-gene interactions by identifying the experimentally validated-target genes of miRNAs.

RESULTS

The outcomes disclosed the interactions of miR-29b-2-5p and miR-342-3p with LSAMP in ATLL_acute, miR-575 with UBN2, miR-342-3p with ZNF280B, and miR-342-5p with FOXRED2 in ATLL_chronic, miR-940 and miR-423-3p with C6orf141, miR-940 and miR-1225-3p with CDCP1, and miR-324-3p with COL14A1 in ATLL_smoldering. These miRNA-gene interactions determine the molecular factors involved in the pathogenesis of each ATLL subtype and the unique ones could be considered biomarkers.

CONCLUSION

The above-mentioned miRNAs-genes interactions are suggested as diagnostic biomarkers for different ATLL subtypes.

Bioinformatic analysis and machine learning to identify the diagnostic biomarkers and immune infiltration in adenomyosis

Background: Adenomyosis is a hormone-dependent benign gynecological disease characterized by the invasion of the endometrium into the myometrium. Women with adenomyosis can suffer from abnormal uterine bleeding, severe pelvic pain, and subfertility or infertility, which can interfere with their quality of life. However, effective diagnostic biomarkers for adenomyosis are currently lacking. The aim of this study is to explore the mechanism of adenomyosis by identifying biomarkers and potential therapeutic targets for adenomyosis and analyzing their correlation with immune infiltration in adenomyosis. Methods: Two datasets, GSE78851 and GSE68870, were downloaded and merged for differential expression analysis and functional enrichment analysis using R software. Weighted gene co-expression network analysis (WGCNA), the least absolute shrinkage and selection operator (LASSO), and support vector machine-recursive feature elimination (SVE-RFE) were combined to explore candidate genes. Quantitative reverse transcriptase PCR (qRT-PCR) was conducted to verify the biomarkers and receiver operating characteristic curve analysis was used to assess the diagnostic value of each biomarker. Single-sample Gene Set Enrichment Analysis (ssGSEA) and CIBERSORT were used to explore immune cell infiltration in adenomyosis and the correlation between diagnostic biomarkers and immune cells. Results: A total of 318 genes were differentially expressed. Through the analysis of differentially expressed genes and WGCNA, we obtained 189 adenomyosis-related genes. After utilizing the LASSO and SVM-RFE algorithms, four hub genes, namely, six-transmembrane epithelial antigen of the prostate-1 (STEAP1), translocase of outer mitochondrial membrane 20 (TOMM20), glycosyltransferase eight domain-containing 2 (GLT8D2), and NME/NM23 family member 5 (NME5) expressed in nucleoside-diphosphate kinase, were identified and verified by qRT-PCR. Immune infiltration analysis indicated that T helper 17 cells, CD56dim natural killer cells, monocytes, and memory B-cell may be associated with the occurrence of adenomyosis. There were significant correlations between the diagnostic biomarkers and immune cells. Conclusion: STEAP1, TOMM20, GLT8D2, and NME5 were identified as potential biomarkers and therapeutic targets for adenomyosis. Immune infiltration may contribute to the onset and progression of adenomyosis.

Construction and analysis of gene co-expression network in the pathogenic fungus Ustilago maydis

Introduction

Biological systems respond to environmental disturbances and a wide range of compounds through complex gene interaction networks. The enormous growth of experimental information obtained using large-scale genomic techniques such as microarrays and RNA sequencing led to the construction of a wide variety of gene co-expression networks in recent years. These networks allow the discovery of clusters of co-expressed genes that potentially work in the same process linking them to biological processes often of interest to industrial, medicinal, and academic research.

Methods

In this study, we built the gene co-expression network of Ustilago maydis from the gene expression data of 168 samples belonging to 19 series, which correspond to the GPL3681 platform deposited in the NCBI using WGCNA software. This network was analyzed to identify clusters of co-expressed genes, gene hubs and Gene Ontology terms. Additionally, we identified relevant modules through a hypergeometric approach based on a predicted set of transcription factors and virulence genes.

Results and Discussion

We identified 13 modules in the gene co-expression network of U. maydis. The TFs enriched in the modules of interest belong to the superfamilies of Nucleic acid-binding proteins, Winged helix DNA-binding, and Zn2/Cys6 DNA-binding. On the other hand, the modules enriched with virulence genes were classified into diseases related to corn smut, Invasive candidiasis, among others. Finally, a large number of hypothetical, a large number of hypothetical genes were identified as highly co-expressed with virulence genes, making them possible experimental targets.

In-depth systems biological evaluation of bovine alveolar macrophages suggests novel insights into molecular mechanisms underlying Mycobacterium bovis infection

Objective

Bovine tuberculosis (bTB) is a chronic respiratory infectious disease of domestic livestock caused by intracellular Mycobacterium bovis infection, which causes ~$3 billion in annual losses to global agriculture. Providing novel tools for bTB managements requires a comprehensive understanding of the molecular regulatory mechanisms underlying the M. bovis infection. Nevertheless, a combination of different bioinformatics and systems biology methods was used in this study in order to clearly understand the molecular regulatory mechanisms of bTB, especially the immunomodulatory mechanisms of M. bovis infection.

Methods

RNA-seq data were retrieved and processed from 78 (39 non-infected control vs. 39 M. bovis-infected samples) bovine alveolar macrophages (bAMs). Next, weighted gene co-expression network analysis (WGCNA) was performed to identify the co-expression modules in non-infected control bAMs as reference set. The WGCNA module preservation approach was then used to identify non-preserved modules between non-infected controls and M. bovis-infected samples (test set). Additionally, functional enrichment analysis was used to investigate the biological behavior of the non-preserved modules and to identify bTB-specific non-preserved modules. Co-expressed hub genes were identified based on module membership (MM) criteria of WGCNA in the non-preserved modules and then integrated with protein-protein interaction (PPI) networks to identify co-expressed hub genes/transcription factors (TFs) with the highest maximal clique centrality (MCC) score (hub-central genes).

Results

As result, WGCNA analysis led to the identification of 21 modules in the non-infected control bAMs (reference set), among which the topological properties of 14 modules were altered in the M. bovis-infected bAMs (test set). Interestingly, 7 of the 14 non-preserved modules were directly related to the molecular mechanisms underlying the host immune response, immunosuppressive mechanisms of M. bovis, and bTB development. Moreover, among the co-expressed hub genes and TFs of the bTB-specific non-preserved modules, 260 genes/TFs had double centrality in both co-expression and PPI networks and played a crucial role in bAMs-M. bovis interactions. Some of these hub-central genes/TFs, including PSMC4, SRC, BCL2L1, VPS11, MDM2, IRF1, CDKN1A, NLRP3, TLR2, MMP9, ZAP70, LCK, TNF, CCL4, MMP1, CTLA4, ITK, IL6, IL1A, IL1B, CCL20, CD3E, NFKB1, EDN1, STAT1, TIMP1, PTGS2, TNFAIP3, BIRC3, MAPK8, VEGFA, VPS18, ICAM1, TBK1, CTSS, IL10, ACAA1, VPS33B, and HIF1A, had potential targets for inducing immunomodulatory mechanisms by M. bovis to evade the host defense response.

Conclusion

The present study provides an in-depth insight into the molecular regulatory mechanisms behind M. bovis infection through biological investigation of the candidate non-preserved modules directly related to bTB development. Furthermore, several hub-central genes/TFs were identified that were significant in determining the fate of M. bovis infection and could be promising targets for developing novel anti-bTB therapies and diagnosis strategies.

Gene expression analysis in endometriosis: Immunopathology insights, transcription factors and therapeutic targets

Background

Endometriosis is recognized as an estrogen-dependent inflammation disorder, estimated to affect 8%-15% of women of childbearing age. Currently, the etiology and pathogenesis of endometriosis are not completely clear. Underlying mechanism for endometriosis is still under debate and needs further exploration. The involvement of transcription factors and immune mediations may be involved in the pathophysiological process of endometriosis, but the specific mechanism remains to be explored. This study aims to investigate the underlying molecular mechanisms in endometriosis.

Methods

The gene expression profile of endometriosis was obtained from the gene expression omnibus (GEO) database. Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were applied to the endometriosis GSE7305 datasets. Cibersort and MCP-counter were used to explore the immune response gene sets, immune response pathway, and immune environment. Differentially expressed genes (DEGs) were identified and screened. Common biological pathways were being investigated using the kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis. Transcription factors were from The Human Transcription Factors. The least absolute shrinkage and selection operator (Lasso) model identified four differential expressions of transcription factors (AEBP1, HOXB6, KLF2, and RORB). Their diagnostic value was calculated by receiver operating characteristic (ROC) curve analysis and validated in the validation cohort (GSE11691, GSE23339). By constructing the interaction network of crucial transcription factors, weighted gene coexpression network analysis (WGCNA) was used to search for key module genes. Metascape was used for enrichment analysis of essential module genes and obtained HOXB6, KLF2. The HOXB6 and KLF2 were further verified as the only two intersection genes according to Support Vector Machine Recursive Feature Elimination (SVM-RFE) and random forest models. We constructed ceRNA (lncRNA-miRNA-mRNA) networks with four potential transcription factors. Finally, we performed molecular docking for goserelin and dienogest with four transcription factors (AEBP1, HOXB6, KLF2, and RORB) to screen potential drug targets.

Results

Immune and metabolic pathways were enriched in GSVA and GSEA. In single sample gene set enrichment analysis (ssGSEA), most immune infiltrating cells, immune response gene sets, and immune response pathways are differentially expressed between endometriosis and non-endometriosis. Twenty-seven transcription factors were screened from differentially expressed genes. Most of the twenty-seven transcription factors were correlated with immune infiltrating cells, immune response gene sets and immune response pathways. Furthermore, Adipocyte enhancer binding protein 1 (AEBP1), Homeobox B6 (HOXB6), Kruppel Like Factor 2 (KLF2) and RAR Related Orphan Receptor B (RORB) were selected out from twenty-seven transcription factors. ROC analysis showed that the four genes had a high diagnostic value for endometriosis. In addition, KLF2 and HOXB6 were found to play particularly important roles in multiple modules (String, WGCNA, SVM-RFE, random forest) on the gene interaction network. Using the ceRNA network, we found that NEAT1 may regulate the expressions of AEBP1, HOXB6 and RORB, while X Inactive Specific Transcript (XIST) may control the expressions of HOXB6, RORB and KLF2. Finally, we found that goserelin and dienogest may be potential drugs to regulate AEBP1, HOXB6, KLF2 and RORB through molecular docking.

Conclusions

AEBP1, HOXB6, KLF2, and RORB may be potential biomarkers for endometriosis. Two of them, KLF2 and HOXB6, are critical molecules in the gene interaction network of endometriosis. Discovered by molecular docking, AEBP1, HOXB6, KLF2, and RORB are targets for goserelin and dienogest.

Comparative transcriptome profiling and weighted gene co-expression network analysis to identify core genes in maize (Zea mays L.) silks infected by multiple fungi

Maize (Zea mays L.) is the third most popular Poaceae crop after wheat and rice and used in feed and pharmaceutical sectors. The maize silk contains bioactive components explored by traditional Chinese herbal medicine for various pharmacological activities. However, Fusarium graminearum, Fusarium verticillioides, Trichoderma atroviride, and Ustilago maydis can infect the maize, produce mycotoxins, hamper the quantity and quality of silk production, and further harm the primary consumer's health. However, the defense mechanism is not fully understood in multiple fungal infections in the silk of Z. mays. In this study, we applied bioinformatics approaches to use the publicly available transcriptome data of Z. mays silk affected by multiple fungal flora to identify core genes involved in combatting disease response. Differentially expressed genes (DEGs) were identified among intra- and inter-transcriptome data sets of control versus infected Z. mays silks. Upon further comparison between up- and downregulated genes within the control of datasets, 4,519 upregulated and 5,125 downregulated genes were found. The DEGs have been compared with genes in the modules of weighted gene co-expression network analysis to relevant specific traits towards identifying core genes. The expression pattern of transcription factors, carbohydrate-active enzymes (CAZyme), and resistance genes was analyzed. The present investigation is supportive of our findings that the gene ontology, immunity stimulus, and resistance genes are upregulated, but physical and metabolic processes such as cell wall organizations and pectin synthesis were downregulated respectively. Our results are indicative that terpene synthase TPS6 and TPS11 are involved in the defense mechanism against fungal infections in maize silk.

eQTLs as causal instruments for the reconstruction of hormone linked gene networks

Hormones act within in highly dynamic systems and much of the phenotypic response to variation in hormone levels is mediated by changes in gene expression. The increase in the number and power of large genetic association studies has led to the identification of hormone linked genetic variants. However, the biological mechanisms underpinning the majority of these loci are poorly understood. The advent of affordable, high throughput next generation sequencing and readily available transcriptomic databases has shown that many of these genetic variants also associate with variation in gene expression levels as expression Quantitative Trait Loci (eQTLs). In addition to further dissecting complex genetic variation, eQTLs have been applied as tools for causal inference. Many hormone networks are driven by transcription factors, and many of these genes can be linked to eQTLs. In this mini-review, we demonstrate how causal inference and gene networks can be used to describe the impact of hormone linked genetic variation upon the transcriptome within an endocrinology context.

A framework for non-preserved consensus gene module detection in Johne's disease

Johne's disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) is a major concern in dairy industry. Since, the pathogenesis of the disease is not clearly known, it is necessary to develop an approach to discover molecular mechanisms behind this disease with high confidence. Biological studies often suffer from issues with reproducibility. Lack of a method to find stable modules in co-expression networks from different datasets related to Johne's disease motivated us to present a computational pipeline to identify non-preserved consensus modules. Two RNA-Seq datasets related to MAP infection were analyzed, and consensus modules were detected and were subjected to the preservation analysis. The non-preserved consensus modules in both datasets were determined as they are modules whose connectivity and density are affected by the disease. Long non-coding RNAs (lncRNAs) and TF genes in the non-preserved consensus modules were identified to construct integrated networks of lncRNA-mRNA-TF. These networks were confirmed by protein-protein interactions (PPIs) networks. Also, the overlapped hub genes between two datasets were considered hub genes of the consensus modules. Out of 66 consensus modules, 21 modules were non-preserved consensus modules, which were common in both datasets and 619 hub genes were members of these modules. Moreover, 34 lncRNA and 152 TF genes were identified in 12 and 19 non-preserved consensus modules, respectively. The predicted PPIs in 17 non-preserved consensus modules were significant, and 283 hub genes were commonly identified in both co-expression and PPIs networks. Functional enrichment analysis revealed that eight out of 21 modules were significantly enriched for biological processes associated with Johne's disease including “inflammatory response,” “interleukin-1-mediated signaling pathway”, “type I interferon signaling pathway,” “cytokine-mediated signaling pathway,” “regulation of interferon-beta production,” and “response to interferon-gamma.” Moreover, some genes (hub mRNA, TF, and lncRNA) were introduced as potential candidates for Johne's disease pathogenesis such as TLR2, NFKB1, IRF1, ATF3, TREM1, CDH26, HMGB1, STAT1, ISG15, CASP3. This study expanded our knowledge of molecular mechanisms involved in Johne's disease, and the presented pipeline enabled us to achieve more valid results.

Identifying key genes in milk fat metabolism by weighted gene co-expression network analysis

Milk fat is the most important and energy-rich substance in milk, and its content and composition are important reference elements in the evaluation of milk quality. However, the current identification of valuable candidate genes affecting milk fat is limited. IlluminaPE150 was used to sequence bovine mammary epithelial cells (BMECs) with high and low milk fat rates (MFP), the weighted gene co-expression network (WGCNA) was used to analyze mRNA expression profile data in this study. As a result, a total of 10,310 genes were used to construct WGCNA, and the genes were classified into 18 modules. Among them, violet (r = 0.74), yellow (r = 0.75) and darkolivegreen (r =  − 0.79) modules were significantly associated with MFP, and 39, 181, 75 hub genes were identified, respectively. Combining enrichment analysis and differential genes (DEs), we screened five key candidate DEs related to lipid metabolism, namely PI4K2A, SLC16A1, ATP8A2, VEGFD and ID1, respectively. Relative to the small intestine, liver, kidney, heart, ovary and uterus, the gene expression of PI4K2A is the highest in mammary gland, and is significantly enriched in GO terms and pathways related to milk fat metabolism, such as monocarboxylic acid transport, phospholipid transport, phosphatidylinositol signaling system, inositol phosphate metabolism and MAPK signaling pathway. This study uses WGCNA to form an overall view of MFP, providing a theoretical basis for identifying potential pathways and hub genes that may be involved in milk fat synthesis.

Time-Course Transcriptome Analysis of Lungs From Mice Infected With Hypervirulent Klebsiella pneumoniae via Aerosolized Intratracheal Inoculation

Hypervirulent Klebsiella pneumoniae (hvKp) can cause life-threatening community-acquired infections among healthy young individuals and is thus of concern for global dissemination. In this study, a mouse model of acute primary hvKp pneumonia was established via aerosolized intratracheal (i.t.) inoculation, laying the foundation for conducting extensive studies related to hvKp. Subsequently, a time-course transcriptional profile was created of the lungs from the mouse model at 0, 12, 24, 48 and 60 hours post-infection (hpi) using RNA Sequencing (RNA-Seq). RNA-Seq data were analyzed with the use of Mfuzz time clustering, weighted gene co-expression network analysis (WGCNA) and Immune Cell Abundance Identifier for mouse (ImmuCellAI-mouse). A gradual change in the transcriptional profile of the lungs was observed that reflected expected disease progression. At 12 hpi, genes related to acute phase inflammatory response increased in expression and lipid metabolism appeared to have a pro-inflammatory effect. At 24 hpi, exacerbation of inflammation was observed and active IFN-γ suggested that signaling promoted activation and recruitment of macrophages occurred. Genes related to maintaining the structural integrity of lung tissues showed a sustained decrease in expression after infection and the decrease was especially marked at 48 hpi. TNF, IL-17, MAPK and NF-kB signaling pathways may play key roles in the immunopathogenesis mechanism at all stages of infection. Natural killer (NK) cells consistently decreased in abundance after infection, which has rarely been reported in hvKp infection and could provide a new target for treatment. Genes Saa1 and Slpi were significantly upregulated during infection. Both Saa1, which is associated with lipopolysaccharide (LPS) that elicits host inflammatory response, and Slpi, which encodes an antimicrobial protein, have not previously been reported in hvKp infections and could be important targets for subsequent studies. To t our knowledge, this paper represents the first study to investigate the pulmonary transcriptional response to hvKp infection. The results provide new insights into the molecular mechanisms underlying the pathogenesis of hvKp pulmonary infection that can contribute to the development of therapies to reduce hvKp pneumonia.

Blood and brain transcriptome analysis reveals APOE genotype-mediated and immune-related pathways involved in Alzheimer disease

BACKGROUND

While Alzheimer disease (AD) is generally considered as a brain disorder, blood biomarkers may be useful for the diagnosis and prediction of AD brain pathology. The APOE ε4 allele has shown cerebrovascular effects including acceleration of blood-brain barrier (BBB) breakdown.

METHODS

We evaluated the differential expression of previously established AD genes in brains from 344 pathologically confirmed AD cases and 232 controls and in blood from 112 pathologically confirmed AD cases and 67 controls from the Religious Orders Study and Memory and Aging Project. Differential gene expression between AD cases and controls was analyzed in the blood and brain jointly using a multivariate approach in the total sample and within APOE genotype groups. Gene set enrichment analysis was performed within APOE genotype groups using the results from the combined blood and brain analyses to identify biologically important pathways. Gene co-expression networks in brain and blood samples were investigated using weighted correlation network analysis. Top-ranked genes from networks and pathways were further evaluated with vascular injury traits.

RESULTS

We observed differentially expressed genes with P < 0.05 in both brain and blood for established AD genes INPP5D (upregulated) and HLA-DQA1 (downregulated). PIGHP1 and FRAS1 were differentially expressed at the transcriptome-wide level (P < 3.3 × 10-6) within ε2/ε3 and ε3/ε4 groups, respectively. Gene set enrichment analysis revealed 21 significant pathways (false discovery rate P < 0.05) in at least one APOE genotype group. Ten pathways were significantly enriched in the ε3/ε4 group, and six of these were unique to these subjects. Four pathways (allograft rejection, interferon gamma response, peroxisome, and TNFA signaling via NFKB) were enriched for AD upregulated genes in the ε3/ε4 group and AD downregulated genes in subjects lacking ε4. We identified a co-expressed gene network in the brain that reproduced in blood and showed higher average expression in ε4 carriers. Twenty-three genes from pathway and network analyses were significantly associated with at least one vascular injury trait.

CONCLUSION

These results suggest that the APOE genotype contributes to unique expression network profiles in both blood and brain. Several genes in these networks are associated with measures of vascular injury and potentially contribute to ε4's effect on the BBB.

MRI radiogenomics for intelligent diagnosis of breast tumors and accurate prediction of neoadjuvant chemotherapy responses-a review

BACKGROUND AND OBJECTIVE

This paper aims to overview multidimensional mining algorithms in relation to Magnetic Resonance Imaging (MRI) radiogenomics for computer aided detection and diagnosis of breast tumours. The work also aims to address a new problem in radiogenomics mining: how to combine structural radiomics information with non-structural genomics information for improving the accuracy and efficacy of Neoadjuvant Chemotherapy (NAC).

METHODS

This requires the automated extraction of parameters from non-structural breast radiomics data, and finding feature vectors with diagnostic value, which then are combined with genomics data. In order to address the problem of weakly labelled tumour images, a Generative Adiversarial Networks (GAN) based deep learning strategy is proposed for the classification of tumour types; this has significant potential for providing accurate real-time identification of tumorous regions from MRI scans. In order to efficiently integrate in a deep learning framework different features from radiogenomics datasets at multiple spatio-temporal resolutions, pyramid structured and multi-scale densely connected U-Nets are proposed. A bidirectional gated recurrent unit (BiGRU) combined with an attention based deep learning approach is also proposed.

RESULTS

The aim is to accurately predict NAC responses by combining imaging and genomic datasets. The approaches discussed incorporate some of the latest developments in of current signal processing and artificial intelligence and have significant potential in advancing and provide a development platform for future cutting-edge biomedical radiogenomics analysis.

CONCLUSIONS

The association of genotypic and phenotypic features is at the core of the emergent field of Precision Medicine. It makes use of advances in biomedical big data analysis, which enables the correlation between disease-associated phenotypic characteristics, genetics polymorphism and gene activation to be revealed.

Gene Co-Expression Network Analysis Identifies Vitamin D-Associated Gene Modules in Adult Normal Rectal Epithelium Following Supplementation

Colorectal cancer (CRC) is a common, multifactorial disease. While observational studies have identified an association between lower vitamin D and higher CRC risk, supplementation trials have been inconclusive and the mechanisms by which vitamin D may modulate CRC risk are not well understood. We sought to perform a weighted gene co-expression network analysis (WGCNA) to identify modules present after vitamin D supplementation (when plasma vitamin D level was sufficient) which were absent before supplementation, and then to identify influential genes in those modules. The transcriptome from normal rectal mucosa biopsies of 49 individuals free from CRC were assessed before and after 12 weeks of 3200IU/day vitamin D (Fultium-D3) supplementation using paired-end total RNAseq. While the effects on expression patterns following vitamin D supplementation were subtle, WGCNA identified highly correlated genes forming gene modules. Four of the 17 modules identified in the post-vitamin D network were not preserved in the pre-vitamin D network, shedding new light on the biochemical impact of supplementation. These modules were enriched for GO terms related to the immune system, hormone metabolism, cell growth and RNA metabolism. Across the four treatment-associated modules, 51 hub genes were identified, with enrichment of 40 different transcription factor motifs in promoter regions of those genes, including VDR:RXR. Six of the hub genes were nominally differentially expressed in studies of vitamin D effects on adult normal mucosa organoids: LCN2, HLA-C, AIF1L, PTPRU, PDE4B and IFI6. By taking a gene-correlation network approach, we have described vitamin D induced changes to gene modules in normal human rectal epithelium in vivo, the target tissue from which CRC develops.

Predictive Biomarkers for Postmyocardial Infarction Heart Failure Using Machine Learning: A Secondary Analysis of a Cohort Study

Background

There are few biomarkers with an excellent predictive value for postacute myocardial infarction (MI) patients who developed heart failure (HF). This study aimed to screen candidate biomarkers to predict post-MI HF.

Methods

This is a secondary analysis of a single-center cohort study including nine post-MI HF patients and eight post-MI patients who remained HF-free over a 6-month follow-up. Transcriptional profiling was analyzed using the whole blood samples collected at admission, discharge, and 1-month follow-up. We screened differentially expressed genes and identified key modules using weighted gene coexpression network analysis. We confirmed the candidate biomarkers using the developed external datasets on post-MI HF. The receiver operating characteristic curves were created to evaluate the predictive value of these candidate biomarkers.

Results

A total of 6,778, 1,136, and 1,974 genes (dataset 1) were differently expressed at admission, discharge, and 1-month follow-up, respectively. The white and royal blue modules were most significantly correlated with post-MI HF (dataset 2). After overlapping dataset 1, dataset 2, and external datasets (dataset 3), we identified five candidate biomarkers, including FCGR2A, GSDMB, MIR330, MED1, and SQSTM1. When GSDMB and SQSTM1 were combined, the area under the curve achieved 1.00, 0.85, and 0.89 in admission, discharge, and 1-month follow-up, respectively.

Conclusions

This study demonstrates that FCGR2A, GSDMB, MIR330, MED1, and SQSTM1 are the candidate predictive biomarker genes for post-MI HF, and the combination of GSDMB and SQSTM1 has a high predictive value.

Weighted Gene Co-Expression Analysis Network-Based Analysis on the Candidate Pathways and Hub Genes in Eggplant Bacterial Wilt-Resistance: A Plant Research Study

Solanum melongena L. (eggplant) bacterial wilt is a severe soil borne disease. Here, this study aimed to explore the regulation mechanism of eggplant bacterial wilt-resistance by transcriptomics with weighted gene co-expression analysis network (WGCNA). The different expression genes (DEGs) of roots and stems were divided into 21 modules. The module of interest (root: indianred4, stem: coral3) with the highest correlation with the target traits was selected to elucidate resistance genes and pathways. The selected module of roots and stems co-enriched the pathways of MAPK signalling pathway, plant pathogen interaction, and glutathione metabolism. Each top 30 hub genes of the roots and stems co-enriched a large number of receptor kinase genes. A total of 14 interesting resistance-related genes were selected and verified with quantitative polymerase chain reaction (qPCR). The qPCR results were consistent with those of WGCNA. The hub gene of EGP00814 (namely SmRPP13L4) was further functionally verified; SmRPP13L4 positively regulated the resistance of eggplant to bacterial wilt by qPCR and virus-induced gene silencing (VIGS). Our study provides a reference for the interaction between eggplants and bacterial wilt and the breeding of broad-spectrum and specific eggplant varieties that are bacterial wilt-resistant.

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.

An integrated analysis of mRNAs, lncRNAs, and miRNAs based on weighted gene co-expression network analysis involved in bovine endometritis

In dairy cattle, endometritis is a severe infectious disease that occurs following parturition. It is clear that genetic factors are involved in the etiology of endometritis, however, the molecular pathogenesis of endometritis is not entirely understood. In this study, a system biology approach was used to better understand the molecular mechanisms underlying the development of endometritis. Forty transcriptomic datasets comprising of 20 RNA-Seq (GSE66825) and 20 miRNA-Seq (GSE66826) were obtained from the GEO database. Next, the co-expressed modules were constructed based on RNA-Seq (Rb-modules) and miRNA-Seq (mb-modules) data, separately, using a weighted gene co-expression network analysis (WGCNA) approach. Preservation analysis was used to find the non-preserved Rb-modules in endometritis samples. Afterward, the non-preserved Rb-modules were assigned to the mb-modules to construct the integrated regulatory networks. Just highly connected genes (hubs) in the networks were considered and functional enrichment analysis was used to identify the biological pathways associated with the development of the disease. Furthermore, additional bioinformatic analysis including protein-protein interactions network and miRNA target prediction were applied to enhance the reliability of the results. Thirty-five Rb-modules and 10 mb-modules were identified and 19 and 10 modules were non-preserved, respectively, which were enriched in biological pathways related to endometritis like inflammation and ciliogenesis. Two non-preserved Rb-modules were significantly assigned to three mb-modules and three and two important sub-networks in the Rb-modules were identified, respectively, including important mRNAs, lncRNAs and miRNAs genes like IRAK1, CASP3, CCDC40, CCDC39, ZMYND10, FOXJ1, TLR4, IL10, STAT3, FN1, AKT1, CD68, ENSBTAG00000049936, ENSBTAG00000050527, ENSBTAG00000051242, ENSBTAG00000049287, bta-miR-449, bta-miR-484, bta-miR-149, bta-miR-30b and bta-miR-423. The potential roles of these genes have been previously demonstrated in endometritis or related pathways, which reinforced putative functions of the suggested integrated regulatory networks in the endometritis pathogenesis. These findings may help further elucidate the underlying mechanisms of bovine endometritis.

Hypermethylation and Downregulation of UTP6 Are Associated With Stemness Properties, Chemoradiotherapy Resistance, and Prognosis in Rectal Cancer: A Co-expression Network Analysis

Background

To identify the hub genes associated with chemoradiotherapy resistance in rectal cancer and explore the potential mechanism.

Methods

Weighted gene co-expression network analysis (WGCNA) was performed to identify the gene modules correlated with the chemoradiotherapy resistance of rectal cancer.

Results

The mRNA expression of 31 rectal cancer patients receiving preoperative chemoradiotherapy was described in our previous study. Through WGCNA, we demonstrated that the chemoradiotherapy resistance modules were enriched for translation, DNA replication, and the androgen receptor signaling pathway. Additionally, we identified and validated UTP6 as a new effective predictor for chemoradiotherapy sensitivity and a prognostic factor for the survival of colorectal cancer patients using our data and the GSE35452 dataset. Low UTP6 expression was correlated with significantly worse disease-free survival (DFS), overall survival (OS), and event- and relapse-free survival both in our data and the R2 Platform. Moreover, we verified the UTP6 expression in 125 locally advanced rectal cancer (LARC) patients samples by immunohistochemical analysis. The results demonstrated that low UTP6 expression was associated with worse DFS and OS by Kaplan-Meier and COX regression model analyses. Gene set enrichment and co-expression analyses showed that the mechanism of the UTP6-mediated chemoradiotherapy resistance may involve the regulation of FOXK2 expression by transcription factor pathways.

Conclusion

Low expression of the UTP6 was found to be associated with chemoradiotherapy resistance and the prognosis of colorectal cancer possibly via regulating FOXK2 expression by transcription factor pathways.

Identification and validation of potential novel biomarkers to predict distant metastasis in differentiated thyroid cancer

Background

Distant metastasis (DM) is not common in differentiated thyroid cancer (DTC). However, it is associated with a significantly poor prognosis. Early detection of high-risk DTC patients is difficult, and the molecular mechanism is still unclear. Therefore, the present study aims to establish a novel predictive model based on clinicopathological parameters and DM-related gene signatures to provide guidelines for clinicians in decision making.

Methods

Weighted gene co-expression network analysis (WGCNA) was performed to discover co-expressed gene modules and hub genes associated with DM. Univariate and multivariate analyses were carried out to identify independent clinicopathological risk factors based on The Cancer Genome Atlas (TCGA) database. An integrated nomogram prediction model was established. Finally, real hub genes were validated using the GSE60542 database and various thyroid cell lines.

Results

The midnightblue module was most significantly positively correlated with DM (R=0.56, P=9e-06) by as per WGCNA. DLX5 (AUC: 0.769), COX6B2 (AUC: 0.764), and LYPD1 (AUC: 0.760) were determined to be the real hub genes that play a crucial role in predicting DM. Meanwhile, univariate and multivariate analyses demonstrated that T-stage (OR, 15.03; 95% CI, 1.75-319.40; and P=0.024), histologic subtype (OR, 0.17; 95% CI, 0.03-0.92; and P=0.042) were the independent predictors of DM. Subsequently, a nomogram model was constructed based on gene signatures and independent clinical risk factors exhibited good performance. Additionally, the mRNA expressions of real hub genes in the GSE60542 dataset were consistent with TCGA.

Conclusions

The present study has provided a reliable model to predict DM in patients with DTC. This model is likely to serve as an individual risk assessment tool in therapeutic decision-making.

Identification of Key Genes and Pathways associated with Endometriosis by Weighted Gene Co-expression Network Analysis

Background: Endometriosis is a common gynecological disorder with high rates of infertility and pelvic pain. However, its pathogenesis and diagnostic biomarkers remain unclear. This study aimed to elucidate potential hub genes and key pathways associated with endometriosis in ectopic endometrium (EC) and eutopic endometrium (EU). Material and Method: EC and EU-associated microarray datasets were obtained from the gene expression omnibus (GEO) database. Gene set enrichment analysis was performed to obtain further biological insight into the EU and EC-associated genes. Weighted gene co-expression network analysis (WGCNA) was performed to find clinically significant modules of highly-correlated genes. The hub genes that belong to both the weighted gene co-expression network and protein-protein interaction (PPI) network were identified using a Venn diagram. Results: We obtained EC and EU-associated microarray datasets GSE7305 and GSE120103. Genes in the EC were mainly enriched in the immune response and immune cell trafficking, and genes in the EU were mainly enriched in stress response and steroid hormone biosynthesis. PPI networks and weighted gene co-expression networks were constructed. An EC-associated blue module and an EU-associated magenta module were identified, and their function annotations revealed that hormone receptor signaling or inflammatory microenvironments may promote EU passing through the oviducts and migrating to the ovarian surfaces, and adhesion and immune correlated genes may induce the successful ectopic implantation of the endometrium (EC). Twelve hub genes in the EC and sixteen hub genes in the EU were recognized and further validated in independent datasets. Conclusion: Our study identified, for the first time, the hub genes and enrichment pathways in the EC and EU using WGCNA, which may provide a comprehensive understanding of the pathogenesis of endometriosis and have important clinical implications for the treatment and diagnosis of endometriosis.

Expression of 5-methylcytosine regulators is highly associated with the clinical phenotypes of prostate cancer and DNMTs expression predicts biochemical recurrence

In patients with prostate cancer (PCa), there is a high rate of overdiagnosis and frequent overtreatment. Therefore, there is an urgent need for more accurate prediction of biochemical recurrence (BCR). DNA methylation regulation patterns play crucial roles in tumorigenicity, progression, and treatment efficacy in PCa. However, the global relationship between epigenetic alterations, changes in mRNA levels, and pathologic phenotypes of PCa remain largely undefined. Here, we conducted a systematic analysis to identify global coexpression and comethylation modules in PCa. We identified coregulated methylation and expression modules and the relationships between epigenetic modifications, tumor progression, and the corresponding immune microenvironment in PCa. Our results show that DNA methyltransferases (DNMTs) are strongly associated with pathologic phenotypes and immune infiltration patterns in PCa. We built a two-factor predictive model using the expression features of DNMT3B and DNMT1. The model was used to predict the BCR status of patients with PCa and achieved area under the receiver operating characteristic curve values of 0.70 and 0.88 in the training and independent testing datasets, respectively.

Identification and Exploration of Novel Macrophage M2-Related Biomarkers and Potential Therapeutic Agents in Endometriosis

Endometriosis (EM) is a chronic neuroinflammatory disorder that is associated with pain and infertility that affects ∼10% of reproductive-age women. The pathophysiology and etiology of EM remain poorly understood, and diagnostic delays are common. Exploration of the underlying molecular mechanism, as well as novel diagnostic biomarkers and therapeutic targets, is urgently needed. Inflammation is known to play a key role in the development of lesions, which are a defining feature of the disorder. In our research, the CIBERSORT and WGCNA algorithms were used to establish a weighted gene co-expression network and to identify macrophage-related hub genes using data downloaded from the GEO database (GSE11691, 7305). The analysis identified 1,157 differentially expressed genes (DEGs) in EM lesions, of which five were identified as being related to M2 macrophages and were validated as differentially expressed by qRT-PCR and immunohistochemistry (IHC). Of these putative novel biomarker genes, bridging integrator 2 (BIN2), chemokine receptor 5 (CCR5), and macrophage mannose receptor 1 (MRC1) were upregulated, while spleen tyrosine kinase (SYK) and metalloproteinase 12 (ADAM12) were downregulated in ectopic endometria vs. normal endometria. Meanwhile, 23 potentially therapeutic small molecules for EM were obtained from the cMAP database, among which topiramate, isoflupredone, adiphenine, dexverapamil, MS-275, and celastrol were the top six molecules with the highest absolute enrichment values. This is our first attempt to use the CIBERSORT and WGCNA algorithms for the identification of novel Mϕ2 macrophage-related biomarkers of EM. Our findings provide novel insights into the impact of immune cells on the etiology of EM; nevertheless, further investigation of these key genes and therapeutic drugs is needed to validate their effects on EM.

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.

Endometrial DNA damage response is modulated in endometriosis

STUDY QUESTION

Is the DNA damage response (DDR) dysregulated in the eutopic endometrium of women with endometriosis?

SUMMARY ANSWER

Endometrial expression of genes involved in DDR is modulated in women with endometriosis, compared to those without the disease.

WHAT IS KNOWN ALREADY

Ectopic endometriotic lesions are reported to harbour somatic mutations, thereby hinting at dysregulation of DDR and DNA repair pathways. However, it remains inconclusive whether the eutopic endometrium also manifests dysregulated DDR in endometriosis.

STUDY DESIGN, SIZE, DURATION

For this case-control study conducted between 2015 and 2019, eutopic endometrial (E) samples (EE- from women with endometriosis, CE- from women without endometriosis) were collected in either mid-proliferative (EE-MP, n = 23; CE-MP, n = 17) or mid-secretory (EE-MS, n = 17; CE-MS, n = 9) phases of the menstrual cycle. This study compares: (i) DNA damage marker localization, (ii) expression of DDR genes and (iii) expression of DNA repair genes in eutopic endometrial samples from women with and without endometriosis.

PARTICIPANTS/MATERIALS, SETTING, METHODS

The study included (i) 40 women (aged 31.9 ± 0.81 years) with endometriosis and (ii) 26 control women (aged 31.4 ± 1.02 years) without endometriosis. Eutopic endometrial samples from the two groups were divided into different parts for histological analysis, immunohistochemistry, RNA extraction, protein extraction and comet assays. Eighty-four genes of relevance in the DNA damage signalling pathway were evaluated for their expression in eutopic endometrial samples, using RT2 Profiler PCR arrays. Validations of the expression of two GADD (Growth Arrest DNA Damage Inducible) proteins - GADD45A and GADD45G were carried out by immunoblotting. DNA damage was assessed by immunohistochemical localization of γ-H2AFX (a phosphorylated variant of histone H2AX) and 8-OHdG (8-hydroxy-2'-deoxyguanosine). RNA sequencing data from mid-proliferative (EE-MP, n = 4; CE-MP, n = 3) and mid-secretory phase (EE-MS and CE-MS, n = 4 each) endometrial samples were scanned to compare the expression status of all the genes implicated in human DNA repair. PCNA (Proliferating Cell Nuclear Antigen) expression was determined to assess endometrial proliferation. Residual DNA damage in primary endometrial cells was checked by comet assays. Public datasets were also scanned for the expression of DDR and DNA repair genes as our RNASeq data were limited by small sample size. All the comparisons were made between phase-matched endometrial samples from women with and without endometriosis.

MAIN RESULTS AND THE ROLE OF CHANCE

Endometrial expression of DDR genes and intensity of immunolocalized γ-H2AFX were significantly (P < 0.05) higher in EE, compared to CE samples. DDR proteins, especially those belonging to the GADD family, were found to be differentially abundant in EE, as compared to CE. These patterns were evident in both mid-proliferative and mid-secretory phases. Intriguingly, higher DDR was associated with increased cell proliferation in EE-MP, compared to CE-MP. Furthermore, among the differentially expressed transcripts (DETs) encoded by DNA repair genes, the majority showed up-regulation in EE-MP, compared to CE-MP. Interestingly, CE-MP and EE-MP had a comparable percentage (P > 0.05) of cells with residual DNA damage. However, unlike the mid-proliferative phase data, many DETs encoded by DNA repair genes were down-regulated in EE-MS, compared to CE-MS. An analysis of the phase-matched control and endometriosis samples included in the GSE51981 dataset available in the Gene Expression Omnibus database also revealed significant (P < 0.05) alterations in the expression of DDR and DNA repair genes in EE, compared to CE.

LARGE-SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

The study was conducted on a limited number of endometrial samples. Also, the study does not reveal the causes underlying dysregulated DDR in the eutopic endometrium of women with endometriosis.

WIDER IMPLICATIONS OF THE FINDINGS

Alterations in the expression of DDR and DNA repair genes indirectly suggest that eutopic endometrium, as compared to its healthy counterpart, encounters DNA damage-inducing stimuli, either of higher strength or for longer duration in endometriosis. It will be worthwhile to identify the nature of such stimuli and also explore the role of higher genomic insults and dysregulated DDR/DNA repair in the origin and/or progression of endometriosis.

STUDY FUNDING/COMPETING INTEREST(S)

The study was supported by the Department of Biotechnology and Indian Council of Medical Research, Government of India. No conflict of interest is declared.

Integrative transcriptomic analysis of tissue-specific metabolic crosstalk after myocardial infarction

Myocardial infarction (MI) promotes a range of systemic effects, many of which are unknown. Here, we investigated the alterations associated with MI progression in heart and other metabolically active tissues (liver, skeletal muscle, and adipose) in a mouse model of MI (induced by ligating the left ascending coronary artery) and sham-operated mice. We performed a genome-wide transcriptomic analysis on tissue samples obtained 6- and 24 hr post MI or sham operation. By generating tissue-specific biological networks, we observed: (1) dysregulation in multiple biological processes (including immune system, mitochondrial dysfunction, fatty-acid beta-oxidation, and RNA and protein processing) across multiple tissues post MI and (2) tissue-specific dysregulation in biological processes in liver and heart post MI. Finally, we validated our findings in two independent MI cohorts. Overall, our integrative analysis highlighted both common and specific biological responses to MI across a range of metabolically active tissues.

The human body is like a state-of-the-art car, where each part must work together with all the others. When a car breaks down, most of the time the problem is not isolated to only one part, as it is an interconnected system. Diseases in the human body can also have systemic effects, so it is important to study their implications throughout the body. Most studies of heart attacks focus on the direct impact on the heart and the cardiovascular system. Learning more about how heart attacks affect rest of the body may help scientists identify heart attacks early or create improved treatments.

Arif and Klevstig et al. show that heart attacks affect the metabolism throughout the body. In the experiments, mice underwent a procedure that mimics either a heart attack or a fake procedure. Then, Arif and Klevstig et al. compared the activity of genes in the heart, muscle, liver and fat tissue of the two groups of mice 6- and 24-hours after the operations. This revealed disruptions in the immune system, metabolism and the production of proteins. The experiments also showed that changes in the activity of four important genes are key to these changes. This suggests that this pattern of changes could be used as a way to identify heart attacks.

The experiments show that heart attacks have important effects throughout the body, especially on metabolism. These discoveries may help scientists learn more about the underlying biological processes and develop new treatments that prevent the harmful systemic effects of heart attacks and boost recovery.

Candidate pharmacological treatments for substance use disorder and suicide identified by gene co-expression network-based drug repositioning

Patients with substance use disorders (SUD) are at high risk to die by suicide. So far, the neurobiology of the suicide-SUD association has not been elucidated. This study aimed to identify potential pharmacological targets among hub genes from brain gene co-expression networks of individuals with SUD in a suicidal and non-suicidal context. Post-mortem samples from the prefrontal cortex of 79 individuals were analyzed. Individuals were classified into the following groups: suicides with SUD (n = 28), suicides without SUD (n = 23), nonsuicides with SUD (n = 9), nonsuicides without SUD (n = 19). Gene expression profiles were evaluated with the Illumina HumanHT-12 v4 array. Co-expression networks were constructed in WGCNA using the differentially expressed genes found in the comparisons: (a) suicides with and without SUD and (b) nonsuicides with and without SUD. Hub genes were selected for drug-gene interaction testing in the DGIdb database. Among drugs interacting with hub genes in suicides we found MAOA inhibitors and dextromethorphan. In the nonsuicide individuals, we found interactions with eglumegad and antipsychotics (olanzapine, clozapine, loxapine). Modafinil was found to interact with genes in both suicides and nonsuicides. These drugs represent possible candidate treatments for patients with SUD with and without suicidal behavior and their study in each context is encouraged.

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 GSE9314, GSE9368, GSE11434 and 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 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.

Target identification and drug discovery by data-driven hypothesis and experimental validation in ovarian endometriosis

OBJECTIVE

To identify targets and discover drugs for ovarian endometriosis (OE) DESIGN: A basic study based on a data-driven hypothesis and experimental validation SETTING: Center for Reproductive Medicine PATIENT(S)/ANIMAL(S): Fourteen patients with OE and 7 healthy donors were recruited, and 15 female C57/BL6 mice were involved.

INTERVENTION(S)

Samples of OE lesions and normal endometrium were obtained. The ITPR1-knockdowned ectopic human endometrial stromal cells (HESCs) were subjected to ribonucleic acid (RNA) sequencing, cell-counting kit-8 (CCK-8) assay, 5-ethynyl-2'-deoxyuridine (EdU) staining, and flow cytometry. Camptothecin was administered to HESCs and in an OE mouse model.

MAIN OUTCOME MEASURE(S)

ITPR1 expression in OE lesions and normal endometrium, cell proliferation and apoptosis of HESCs with ITPR1 knockdown or camptothecin treatment, and autograft volume in the OE mouse model RESULT(S): Two significant OE-relevant gene modules were identified and involved the PI3K/Akt and aging-relevant pathways. Fifteen hub genes were identified and confirmed, among which the most significant gene, ITPR1, was robustly elevated in OE lesions. RNA sequencing revealed that ITPR1 was highly relevant to cell proliferation and apoptosis, which was further confirmed by CCK-8 assay, EdU staining, and flow cytometry analysis. ITPR1 knockdown inhibited cell proliferation and induced HESC apoptosis. The candidate drugs targeting these modules were screened, among which camptothecin and irinotecan were identified as promising drugs. Both compounds suppressed HESC proliferation and induced apoptosis; ITPR1 expression was suppressed by camptothecin. The therapeutic effect of camptothecin was also validated in the OE mouse model.

CONCLUSION(S)

This study identified the therapeutic targets and promising drugs for OE and shed light on the use of camptothecin in OE treatment.

CoExp: A Web Tool for the Exploitation of Co-expression Networks

Gene co-expression networks are a powerful type of analysis to construct gene groupings based on transcriptomic profiling. Co-expression networks make it possible to discover modules of genes whose mRNA levels are highly correlated across samples. Subsequent annotation of modules often reveals biological functions and/or evidence of cellular specificity for cell types implicated in the tissue being studied. There are multiple ways to perform such analyses with weighted gene co-expression network analysis (WGCNA) amongst one of the most widely used R packages. While managing a few network models can be done manually, it is often more advantageous to study a wider set of models derived from multiple independently generated transcriptomic data sets (e.g., multiple networks built from many transcriptomic sources). However, there is no software tool available that allows this to be easily achieved. Furthermore, the visual nature of co-expression networks in combination with the coding skills required to explore networks, makes the construction of a web-based platform for their management highly desirable. Here, we present the CoExp Web application, a user-friendly online tool that allows the exploitation of the full collection of 109 co-expression networks provided by the CoExpNets suite of R packages. We describe the usage of CoExp, including its contents and the functionality available through the family of CoExpNets packages. All the tools presented, including the web front- and back-ends are available for the research community so any research group can build its own suite of networks and make them accessible through their own CoExp Web application. Therefore, this paper is of interest to both researchers wishing to annotate their genes of interest across different brain network models and specialists interested in the creation of GCNs looking for a tool to appropriately manage, use, publish, and share their networks in a consistent and productive manner.

Identification of LncRNAs Associated With FOLFOX Chemoresistance in mCRC and Construction of a Predictive Model

Oxaliplatin, fluorouracil plus leucovorin (FOLFOX) regimen is the first-line chemotherapy of patients with metastatic colorectal cancer (mCRC). However, studies are limited regarding long non-coding RNAs (lncRNAs) associated with FOLFOX chemotherapy response and prognosis. This study aimed to identify lncRNAs associated with FOLFOX chemotherapy response and prognosis in mCRC patients and to construct a predictive model. We analyzed lncRNA expression in 11 mCRC patients treated with FOLFOX chemotherapy before surgery (four sensitive, seven resistant) by Gene Array Chip. The top eight lncRNAs (AC007193.8, CTD-2008N3.1, FLJ36777, RP11-509J21.4, RP3-508I15.20, LOC100130950, RP5-1042K10.13, and LINC00476) for chemotherapy response were identified according to weighted correlation network analysis (WGCNA). A competitive endogenous RNA (ceRNA) network was then constructed. The crucial functions of the eight lncRNAs enriched in chemotherapy resistance were mitogen-activated protein kinase (MAPK) and proteoglycans signaling pathway. Receiver operating characteristic (ROC) analysis demonstrated that the eight lncRNAs were potent predictors for chemotherapy resistance of mCRC patients. To further identify a signature model lncRNA chemotherapy response and prognosis, the validation set consisted of 196 CRC patients from our center was used to validate lncRNAs expression and prognosis by quantitative PCR (qPCR). The expression of the eight lncRNAs expression between CRC cancerous and adjacent non-cancerous tissues was also verified in the validation data set to determine the prognostic value. A generalized linear model was established to predict the probability of chemotherapy resistance and survival. Our findings showed that the eight-lncRNA signature may be a novel biomarker for the prediction of FOLFOX chemotherapy response and prognosis of mCRC patients.

Weighted gene co-expression network analysis of the salt-responsive transcriptomes reveals novel hub genes in green halophytic microalgae Dunaliella salina

Despite responses to salinity stress in Dunaliella salina, a unicellular halotolerant green alga, being subject to extensive study, but the underlying molecular mechanism remains unknown. Here, Empirical Bayes method was applied to identify the common differentially expressed genes (DEGs) between hypersaline and normal conditions. Then, using weighted gene co-expression network analysis (WGCNA), which takes advantage of a graph theoretical approach, highly correlated genes were clustered as a module. Subsequently, connectivity patterns of the identified modules in two conditions were surveyed to define preserved and non-preserved modules by combining the Zsummary and medianRank measures. Finally, common and specific hub genes in non-preserved modules were determined using Eigengene-based module connectivity or module membership (k_ME) measures and validation was performed by using leave-one-out cross-validation (LOOCV). In this study, the power of beta = 12 (scale-free R2 = 0.8) was selected as the soft-thresholding to ensure a scale-free network, which led to the identification of 15 co-expression modules. Results also indicate that green, blue, brown, and yellow modules are non-preserved in salinity stress conditions. Examples of enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in non-preserved modules are Sulfur metabolism, Oxidative phosphorylation, Porphyrin and chlorophyll metabolism, Vitamin B6 metabolism. Moreover, the systems biology approach was applied here, proposed some salinity specific hub genes, such as radical-induced cell death1 protein (RCD1), mitogen-activated protein kinase kinase kinase 13 (MAP3K13), long-chain acyl-CoA synthetase (ACSL), acetyl-CoA carboxylase, biotin carboxylase subunit (AccC), and fructose-bisphosphate aldolase (ALDO), for the development of metabolites accumulating strains in D. salina.

Systematic Identification of Key Functional Modules and Genes in Gastric Cancer

Gastric cancer (GC) is associated with high incidence and mortality rates worldwide. Differentially expressed gene (DEG) analysis and weighted gene coexpression network analysis (WGCNA) are important bioinformatic methods for screening core genes. In our study, DEG analysis and WGCNA were combined to screen the hub genes, and pathway enrichment analyses were performed on the DEGs. SBNO2 was identified as the hub gene based on the intersection between the DEGs and the purple module in WGCNA. The expression and prognostic value of SBNO2 were verified in UALCAN, GEPIA2, Human Cancer Metastasis Database, Kaplan–Meier plotter, and TIMER. We identified 1974 DEGs, and 28 modules were uncovered via WGCNA. The purple module was identified as the hub module in WGCNA. SBNO2 was identified as the hub gene, which was upregulated in tumour tissues. Moreover, patients with GC and higher SBNO2 expression had worse prognoses. In addition, SBNO2 was suggested to play an important role in immune cell infiltration. In summary, based on DEGs and key modules related to GC, we identified SBNO2 as a hub gene, thereby offering novel insights into the development and treatment of GC.

Integrative Analysis of MAPK14 as a Potential Biomarker for Cardioembolic Stroke

The aim of this study was to obtain the candidate genes and biomarkers that are significantly related to cardioembolic stroke (CS) by applying bioinformatics analysis. In accordance with the results of the weighted gene coexpression network analysis (WGCNA) in the GSE58294 dataset, 11 CS-related coexpression network modules were identified in this study. Correlation analysis showed that the black and pink modules are significantly associated with CS. A total of 18 core genes in the black module and one core gene in the pink module were determined. We then identified differentially expressed genes (DEGs) of CS at 3 h, 5 h, and 24 h postonset. After performing intersection, it was found that 311 genes were coexpressed at these three time points. These genes were majorly enriched in positive regulation of transferase activity and regulation of peptidase activity. The abovementioned coexpressed DEGs were subjected to protein-protein interaction analysis and subnetwork module analysis. Subsequently, we used cytoHubba to obtain 11 key genes from DEGs. The intersection of the core genes screened from WGCNA and the key genes selected from DEGs yielded the MAPK14 gene. The expression level of MAPK14 on the receiver operating characteristic (ROC) curves of CS at 3 h, 5 h, and 24 h showed that the area under the ROC curve (AUC) was 0.923, 0.934, and 0.941, respectively. In a nutshell, MAPK14 screened out by using WGCNA showed differential expression in CS. We conclude that MAPK14 can be used as a potential biological marker of CS and exhibits potential to predict the physiopathological condition of CS patients.

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.

Analysis of exosomal lncRNA, miRNA and mRNA expression profiles and ceRNA network construction in endometriosis

Aim: To investigate exosomal RNAs (long noncoding RNAs (lncRNAs), microRNAs (miRNAs) and messenger RNAs (mRNAs)) profiling and their related networks in endometriosis (EMs). Materials & methods: RNA sequence was performed in exosomes from ovarian endometriomas (EC), eutopic endometria (EU) and normal endometria (Control) stromal cells. The bioinformatics algorithms evaluated competing endogenous RNA (ceRNA) networks. The top-ranked ceRNA networks were confirmed by RT-PCR. Results: Overlapped differentially expressed 938 lncRNAs, 39 miRNAs and 1449 mRNAs were identified. 13 co-expression modules and 61 ceRNA networks were constructed. Conclusion: This study for the first time shows exosomal RNA biomarkers and lncRNA-related networks in EMs, which reveals a novel molecular mechanism of EMs and provides new resources for EM diagnosis and treatment.

FBXW4 Acts as a Protector of FOLFOX-Based Chemotherapy in Metastatic Colorectal Cancer Identified by Co-Expression Network Analysis

Background

FOLFOX chemotherapy is one of the most commonly used treatments for colorectal cancer (CRC) patients. However, the efficacy and tolerance of FOLFOX therapy varies between patients. The purpose of this study was to explore hub genes associated with primary chemotherapy-resistance and to explore the possible mechanisms involved from non-European patients.

Method

A weighted gene co-expression network was constructed to identify gene modules associated with chemotherapy resistance in mCRC from China.

Results

A Gene Array Chip was used to detect mRNA expression in 11 mCRC patients receiving preoperative FOLFOX chemotherapy. The immune response was associated with chemotherapy-resistance in microarray data. Through the use of WGCNA, we demonstrated that the crucial functions enriched in chemotherapy-resistance modules were cell proliferation, MAPK signaling pathways, and PI3K signaling pathways. Additionally, we identified and validated FBXW4 as a new effective predictor for chemotherapy sensitivity and a prognostic factor for survival of CRC patients by using our own data and GSE69657. Furthermore, a meta-analysis of 15 Gene Expression Omnibus–sourced datasets showed that FBXW4 messenger RNA levels were significantly lower in CRC tissues than in normal colon tissues. An analysis of the data from the R2: Genomics Analysis and Visualization Platform showed that low FBXW4 expression was correlated with a significantly worse event- and relapse-free survival. Gene set enrichment analysis showed that the mechanism of FBXW4-mediated chemotherapy resistance may involve the DNA replication signal pathway and the cell cycle.

Conclusion

FBXW4 is associated with chemotherapy resistance and prognosis of CRC probably by regulating DNA replication signaling pathways and the cell cycle.

Systems biology approach identifies functional modules and regulatory hubs related to secondary metabolites accumulation after transition from autotrophic to heterotrophic growth condition in microalgae

Heterotrophic growth mode is among the most promising strategies put forth to overcome the low biomass and secondary metabolites productivity challenge. To shedding light on the underlying molecular mechanisms, transcriptome meta-analysis was integrated with weighted gene co-expression network analysis (WGCNA), connectivity analysis, functional enrichment, and hubs identification. Meta-analysis and Functional enrichment analysis demonstrated that most of the biological processes are up-regulated at heterotrophic growth condition, which leads to change of genetic architectures and phenotypic outcomes. WGNCA analysis of meta-genes also resulted four significant functional modules across logarithmic (LG), transition (TR), and production peak (PR) phases. The expression pattern and connectivity characteristics of the brown module as a non-preserved module vary across LG, TR, and PR phases. Functional analysis identified Carotenoid biosynthesis, Fatty acid metabolism and Methane metabolism as enriched pathways in the non-preserved module. Our integrated approach was applied here, identified some hubs, such as a serine hydroxymethyltransferase (SHMT1), which is the best candidate for development of metabolites accumulating strains in microalgae. Current study provided a new insight into underlying metabolite accumulation mechanisms and opens new avenue for the future applied studies in the microalgae field.

Identification of Modules With Similar Gene Regulation and Metabolic Functions Based on Co-expression Data

Biological systems respond to environmental perturbations and to a large diversity of compounds through gene interactions, and these genetic factors comprise complex networks. In particular, a wide variety of gene co-expression networks have been constructed in recent years thanks to the dramatic increase of experimental information obtained with techniques, such as microarrays and RNA sequencing. These networks allow the identification of groups of co-expressed genes that can function in the same process and, in turn, these networks may be related to biological functions of industrial, medical and academic interest. In this study, gene co-expression networks for 17 bacterial organisms from the COLOMBOS database were analyzed via weighted gene co-expression network analysis and clustered into modules of genes with similar expression patterns for each species. These networks were analyzed to determine relevant modules through a hypergeometric approach based on a set of transcription factors and enzymes for each genome. The richest modules were characterized using PFAM families and KEGG metabolic maps. Additionally, we conducted a Gene Ontology analysis for enrichment of biological functions. Finally, we identified modules that shared similarity through all the studied organisms by using comparative genomics.

Cross-disorder analysis of endometriosis and its comorbid diseases reveals shared genes and molecular pathways and proposes putative biomarkers of endometriosis

RESEARCH QUESTION

Women with endometriosis are considered to be at higher risk of several chronic diseases, such as autoimmune disorders, gynaecological cancers, asthma/atopic diseases and cardiovascular and inflammatory bowel diseases. Could the study of endometriosis-associated comorbidities help to identify potential biomarkers and target pathways of endometriosis?

DESIGN

A systematic review was performed to identify all possible endometriosis-associated comorbid conditions. Next, this list of disorders was coded into MeSH terms, and the gene expression profiles were downloaded from the Phenopedia database and subsequently analysed following a systems biology approach.

RESULTS

The results identified a group of 127 candidate genes that were recurrently expressed in endometriosis and its closest comorbidities and that were defined as 'endometriosis sibling disorders' (ESD). The enrichment analysis showed that these candidate genes are principally involved in immune and drug responses, hormone metabolism and cell proliferation, which are well-known hallmarks of endometriosis. The expression of ESD genes was then validated on independent sample cohorts (n = 207 samples), in which the involvement of 16 genes (AGTR1, BDNF, C3, CCL2, CD40, CYP17A1, ESR1, IGF1, IGF2, IL10, MMP1, MMP7, MMP9, PGR, SERPINE1 and TIMP2) in endometriosis was confirmed. Several of these genes harbour polymorphisms that associate to either endometriosis or its comorbid conditions.

CONCLUSIONS

The study results highlight the molecular processes underlying the aetiopathogenesis of endometriosis and its comorbid conditions, and identify putative endometriosis biomarkers.