K-Module Algorithm: An Additional Step to Improve the Clustering Results of WGCNA Co-Expression Networks

Network medicine approaches for identification of novel prognostic systems biomarkers and drug candidates for papillary thyroid carcinoma

Papillary thyroid carcinoma (PTC) is one of the most common endocrine carcinomas worldwide and the aetiology of this cancer is still not well understood. Therefore, it remains important to understand the disease mechanism and find prognostic biomarkers and/or drug candidates for PTC. Compared with approaches based on single-gene assessment, network medicine analysis offers great promise to address this need. Accordingly, in the present study, we performed differential co-expressed network analysis using five transcriptome datasets in patients with PTC and healthy controls. Following meta-analysis of the transcriptome datasets, we uncovered common differentially expressed genes (DEGs) for PTC and, using these genes as proxies, found a highly clustered differentially expressed co-expressed module: a 'PTC-module'. Using independent data, we demonstrated the high prognostic capacity of the PTC-module and designated this module as a prognostic systems biomarker. In addition, using the nodes of the PTC-module, we performed drug repurposing and text mining analyzes to identify novel drug candidates for the disease. We performed molecular docking simulations, and identified: 4-demethoxydaunorubicin hydrochloride, AS605240, BRD-A60245366, ER 27319 maleate, sinensetin, and TWS119 as novel drug candidates whose efficacy was also confirmed by in silico analyzes. Consequently, we have highlighted here the need for differential co-expression analysis to gain a systems-level understanding of a complex disease, and we provide candidate prognostic systems biomarker and novel drugs for PTC.

Construction of a prognostic prediction model for renal clear cell carcinoma combining clinical traits

Kidney renal clear cell carcinoma (KIRC) is one of the common malignant tumors of the urinary system. Patients with different risk levels are other in terms of disease progression patterns and disease regression. The poorer prognosis for high-risk patients compared to low-risk patients. Therefore, it is essential to accurately high-risk screen patients and gives accurate and timely treatment. Differential gene analysis, weighted correlation network analysis, Protein-protein interaction network, and univariate Cox analysis were performed sequentially on the train set. Next, the KIRC prognostic model was constructed using the least absolute shrinkage and selection operator (LASSO), and the Cancer Genome Atlas (TCGA) test set and the Gene Expression Omnibus dataset verified the model's validity. Finally, the constructed models were analyzed; including gene set enrichment analysis (GSEA) and immune analysis. The differences in pathways and immune functions between the high-risk and low-risk groups were observed to provide a reference for clinical treatment and diagnosis. A four-step key gene screen resulted in 17 key factors associated with disease prognosis, including 14 genes and 3 clinical features. The LASSO regression algorithm selected the seven most critical key factors to construct the model: age, grade, stage, GDF3, CASR, CLDN10, and COL9A2. In the training set, the accuracy of the model in predicting 1-, 2- and 3-year survival rates was 0.883, 0.819, and 0.830, respectively. The accuracy of the TCGA dataset was 0.831, 0.801, and 0.791, and the accuracy of the GSE29609 dataset was 0.812, 0.809, and 0.851 in the test set. Model scoring divided the sample into a high-risk group and a low-risk group. There were significant differences in disease progression and risk scores between the two groups. GSEA analysis revealed that the enriched pathways in the high-risk group mainly included proteasome and primary immunodeficiency. Immunological analysis showed that CD8 (+) T cells, M1 macrophages, PDCD1, and CTLA4 were upregulated in the high-risk group. In contrast, antigen-presenting cell stimulation and T-cell co-suppression were more active in the high-risk group. This study added clinical characteristics to constructing the KIRC prognostic model to improve prediction accuracy. It provides help to assess the risk of patients more accurately. The differences in pathways and immunity between high and low-risk groups were also analyzed to provide ideas for treating KIRC patients.

A functional gene module identification algorithm in gene expression data based on genetic algorithm and gene ontology

Since genes do not function individually, the gene module is considered an important tool for interpreting gene expression profiles. In order to consider both functional similarity and expression similarity in module identification, GMIGAGO, a functional Gene Module Identification algorithm based on Genetic Algorithm and Gene Ontology, was proposed in this work. GMIGAGO is an overlapping gene module identification algorithm, which mainly includes two stages: In the first stage (initial identification of gene modules), Improved Partitioning Around Medoids Based on Genetic Algorithm (PAM-GA) is used for the initial clustering on gene expression profiling, and traditional gene co-expression modules can be obtained. Only similarity of expression levels is considered at this stage. In the second stage (optimization of functional similarity within gene modules), Genetic Algorithm for Functional Similarity Optimization (FSO-GA) is used to optimize gene modules based on gene ontology, and functional similarity within gene modules can be improved. Without loss of generality, we compared GMIGAGO with state-of-the-art gene module identification methods on six gene expression datasets, and GMIGAGO identified the gene modules with the highest functional similarity (much higher than state-of-the-art algorithms). GMIGAGO was applied in BRCA, THCA, HNSC, COVID-19, Stem, and Radiation datasets, and it identified some interesting modules which performed important biological functions. The hub genes in these modules could be used as potential targets for diseases or radiation protection. In summary, GMIGAGO has excellent performance in mining molecular mechanisms, and it can also identify potential biomarkers for individual precision therapy.

Integrated bioinformatic analyses investigate macrophage-M1-related biomarkers and tuberculosis therapeutic drugs

Tuberculosis (TB) is a common infectious disease linked to host genetics and the innate immune response. It is vital to investigate new molecular mechanisms and efficient biomarkers for Tuberculosis because the pathophysiology of the disease is still unclear, and there aren't any precise diagnostic tools. This study downloaded three blood datasets from the GEO database, two of which (GSE19435 and 83456) were used to build a weighted gene co-expression network for searching hub genes associated with macrophage M1 by the CIBERSORT and WGCNA algorithms. Furthermore, 994 differentially expressed genes (DEGs) were extracted from healthy and TB samples, four of which were associated with macrophage M1, naming RTP4, CXCL10, CD38, and IFI44. They were confirmed as upregulation in TB samples by external dataset validation (GSE34608) and quantitative real-time PCR analysis (qRT-PCR). CMap was used to predict potential therapeutic compounds for tuberculosis using 300 differentially expressed genes (150 downregulated and 150 upregulated genes), and six small molecules (RWJ-21757, phenamil, benzanthrone, TG-101348, metyrapone, and WT-161) with a higher confidence value were extracted. We used in-depth bioinformatics analysis to investigate significant macrophage M1-related genes and promising anti-Tuberculosis therapeutic compounds. However, more clinical trials were necessary to determine their effect on Tuberculosis.

KISL: knowledge-injected semi-supervised learning for biological co-expression network modules

The exploration of important biomarkers associated with cancer development is crucial for diagnosing cancer, designing therapeutic interventions, and predicting prognoses. The analysis of gene co-expression provides a systemic perspective on gene networks and can be a valuable tool for mining biomarkers. The main objective of co-expression network analysis is to discover highly synergistic sets of genes, and the most widely used method is weighted gene co-expression network analysis (WGCNA). With the Pearson correlation coefficient, WGCNA measures gene correlation, and uses hierarchical clustering to identify gene modules. The Pearson correlation coefficient reflects only the linear dependence between variables, and the main drawback of hierarchical clustering is that once two objects are clustered together, the process cannot be reversed. Hence, readjusting inappropriate cluster divisions is not possible. Existing co-expression network analysis methods rely on unsupervised methods that do not utilize prior biological knowledge for module delineation. Here we present a method for identification of outstanding modules in a co-expression network using a knowledge-injected semi-supervised learning approach (KISL), which utilizes apriori biological knowledge and a semi-supervised clustering method to address the issue existing in the current GCN-based clustering methods. To measure the linear and non-linear dependence between genes, we introduce a distance correlation due to the complexity of the gene-gene relationship. Eight RNA-seq datasets of cancer samples are used to validate its effectiveness. In all eight datasets, the KISL algorithm outperformed WGCNA when comparing the silhouette coefficient, Calinski-Harabasz index and Davies-Bouldin index evaluation metrics. According to the results, KISL clusters had better cluster evaluation values and better gene module aggregation. Enrichment analysis of the recognition modules demonstrated their effectiveness in discovering modular structures in biological co-expression networks. In addition, as a general method, KISL can be applied to various co-expression network analyses based on similarity metrics. Source codes for the KISL and the related scripts are available online at https://github.com/Mowonhoo/KISL.git.

Identification of biomarkers related to neutrophils and two molecular subtypes of systemic lupus erythematosus

Background

Systemic lupus erythematosus (SLE), an autoimmune disease with complex pathogenesis, poses a considerable threat to women’s health. Increasing evidence indicates that neutrophils play an important role in the development and progression of lupus.

Methods

Weighted correlation network analysis and single-sample gene set enrichment analysis (GSEA) were used to analyse SLE expression data from a comprehensive gene expression database and identify modules associated with neutrophils. Thereafter, the biomarkers most closely related to neutrophils were identified. We reclassified SLE into two molecular subtypes based on the aforementioned biomarkers and evaluated cell infiltration, molecular mechanisms, and signature pathways in each subtype.

Results

The results showed significant differences in immunological characteristics between the two molecular subtypes of SLE. Hub genes were significantly upregulated in the NEUT-H subtype, and they may be associated with lupus activity. The GSEA revealed associations between our biomarkers and key metabolic pathways.

Conclusions

Our study provides not only a classification for patients with SLE but also new cell and gene targets for immunotherapy, as well as a new experimental paradigm to explore immunotherapy for other autoimmune diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01306-9.

Molecular Regulatory Mechanisms Drive Emergent Pathogenetic Properties of Neisseria gonorrhoeae

Neisseria gonorrhoeae is the causative agent of the sexually transmitted infection (STI) gonorrhea, with an estimated 87 million annual cases worldwide. N. gonorrhoeae predominantly colonizes the male and female genital tract (FGT). In the FGT, N. gonorrhoeae confronts fluctuating levels of nutrients and oxidative and non-oxidative antimicrobial defenses of the immune system, as well as the resident microbiome. One mechanism utilized by N. gonorrhoeae to adapt to this dynamic FGT niche is to modulate gene expression primarily through DNA-binding transcriptional regulators. Here, we describe the major N. gonorrhoeae transcriptional regulators, genes under their control, and how these regulatory processes lead to pathogenic properties of N. gonorrhoeae during natural infection. We also discuss the current knowledge of the structure, function, and diversity of the FGT microbiome and its influence on gonococcal survival and transcriptional responses orchestrated by its DNA-binding regulators. We conclude with recent multi-omics data and modeling tools and their application to FGT microbiome dynamics. Understanding the strategies utilized by N. gonorrhoeae to regulate gene expression and their impact on the emergent characteristics of this pathogen during infection has the potential to identify new effective strategies to both treat and prevent gonorrhea.

Distance correlation application to gene co-expression network analysis

BACKGROUND

To construct gene co-expression networks, it is necessary to evaluate the correlation between different gene expression profiles. However, commonly used correlation metrics, including both linear (such as Pearson's correlation) and monotonic (such as Spearman's correlation) dependence metrics, are not enough to observe the nature of real biological systems. Hence, introducing a more informative correlation metric when constructing gene co-expression networks is still an interesting topic.

RESULTS

In this paper, we test distance correlation, a correlation metric integrating both linear and non-linear dependence, with other three typical metrics (Pearson's correlation, Spearman's correlation, and maximal information coefficient) on four different arrays (macrophage and liver) and RNA-seq (cervical cancer and pancreatic cancer) datasets. Among all the metrics, distance correlation is distribution free and can provide better performance on complex relationships and anti-outlier. Furthermore, distance correlation is applied to Weighted Gene Co-expression Network Analysis (WGCNA) for constructing a gene co-expression network analysis method which we named Distance Correlation-based Weighted Gene Co-expression Network Analysis (DC-WGCNA). Compared with traditional WGCNA, DC-WGCNA can enhance the result of enrichment analysis and improve the module stability.

CONCLUSIONS

Distance correlation is better at revealing complex biological relationships between gene profiles compared with other correlation metrics, which contribute to more meaningful modules when analyzing gene co-expression networks. However, due to the high time complexity of distance correlation, the implementation requires more computer memory.

A Novel Calibration Step in Gene Co-Expression Network Construction

High-throughput technologies such as DNA microarrays and RNA-sequencing are used to measure the expression levels of large numbers of genes simultaneously. To support the extraction of biological knowledge, individual gene expression levels are transformed to Gene Co-expression Networks (GCNs). In a GCN, nodes correspond to genes, and the weight of the connection between two nodes is a measure of similarity in the expression behavior of the two genes. In general, GCN construction and analysis includes three steps; 1) calculating a similarity value for each pair of genes 2) using these similarity values to construct a fully connected weighted network 3) finding clusters of genes in the network, commonly called modules. The specific implementation of these three steps can significantly impact the final output and the downstream biological analysis. GCN construction is a well-studied topic. Existing algorithms rely on relatively simple statistical and mathematical tools to implement these steps. Currently, software package WGCNA appears to be the most widely accepted standard. We hypothesize that the raw features provided by sequencing data can be leveraged to extract modules of higher quality. A novel preprocessing step of the gene expression data set is introduced that in effect calibrates the expression levels of individual genes, before computing pairwise similarities. Further, the similarity is computed as an inner-product of positive vectors. In experiments, this provides a significant improvement over WGCNA, as measured by aggregate p-values of the gene ontology term enrichment of the computed modules.