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

k- Strong Inference Algorithm: A Hybrid Information Theory Based Gene Network Inference Algorithm

Gene networks allow researchers to understand the underlying mechanisms between diseases and genes while reducing the need for wet lab experiments. Numerous gene network inference (GNI) algorithms have been presented in the literature to infer accurate gene networks. We proposed a hybrid GNI algorithm, k-Strong Inference Algorithm (ksia), to infer more reliable and robust gene networks from omics datasets. To increase reliability, ksia integrates Pearson correlation coefficient (PCC) and Spearman rank correlation coefficient (SCC) scores to determine mutual information scores between molecules to increase diversity of relation predictions. To infer a more robust gene network, ksia applies three different elimination steps to remove redundant and spurious relations between genes. The performance of ksia was evaluated on microbe microarrays database in the overlap analysis with other GNI algorithms, namely ARACNE, C3NET, CLR, and MRNET. Ksia inferred less number of relations due to its strict elimination steps. However, ksia generally performed better on Escherichia coli (E.coli) and Saccharomyces cerevisiae (yeast) gene expression datasets due to F- measure and precision values. The integration of association estimator scores and three elimination stages slightly increases the performance of ksia based gene networks. Users can access ksia R package and user manual of package via https://github.com/ozgurcingiz/ksia .

Aortic Aneurysm: Finding the Right Target

This Special Issue of Biomedicines highlights many important scientific findings in aneurysm research [...].

Identification of key monocytes/macrophages related gene set of the early-stage abdominal aortic aneurysm by integrated bioinformatics analysis and experimental validation

Objective

Abdominal aortic aneurysm (AAA) is a lethal peripheral vascular disease. Inflammatory immune cell infiltration is a central part of the pathogenesis of AAA. It’s critical to investigate the molecular mechanisms underlying immune infiltration in early-stage AAA and look for a viable AAA marker.

Methods

In this study, we download several mRNA expression datasets and scRNA-seq datasets of the early-stage AAA models from the NCBI-GEO database. mMCP-counter and CIBERSORT were used to assess immune infiltration in early-stage experimental AAA. The scRNA-seq datasets were then utilized to analyze AAA-related gene modules of monocytes/macrophages infiltrated into the early-stage AAA by Weighted Correlation Network analysis (WGCNA). After that, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis for the module genes was performed by ClusterProfiler. The STRING database was used to create the protein-protein interaction (PPI) network. The Differentially Expressed Genes (DEGs) of the monocytes/macrophages were explored by Limma-Voom and the key gene set were identified. Then We further examined the expression of key genes in the human AAA dataset and built a logistic diagnostic model for distinguishing AAA patients and healthy people. Finally, real-time quantitative polymerase chain reaction (RT-qPCR) and Enzyme Linked Immunosorbent Assay (ELISA) were performed to validate the gene expression and serum protein level between the AAA and healthy donor samples in our cohort.

Results

Monocytes/macrophages were identified as the major immune cells infiltrating the early-stage experimental AAA. After pseudocell construction of monocytes/macrophages from scRNA-seq datasets and WGCNA analysis, four gene modules from two datasets were identified positively related to AAA, mainly enriched in Myeloid Leukocyte Migration, Collagen-Containing Extracellular matrix, and PI3K-Akt signaling pathway by functional enrichment analysis. Thbs1, Clec4e, and Il1b were identified as key genes among the hub genes in the modules, and the high expression of Clec4e, Il1b, and Thbs1 was confirmed in the other datasets. Then, in human AAA transcriptome datasets, the high expression of CLEC4E, IL1B was confirmed and a logistic regression model based on the two gene expressions was built, with an AUC of 0.9 in the train set and 0.79 in the validated set. Additionally, in our cohort, we confirmed the increased serum protein levels of IL-1β and CLEC4E in AAA patients as well as the increased expression of these two genes in AAA aorta samples.

Conclusion

This study identified monocytes/macrophages as the main immune cells infiltrated into the early-stage AAA and constructed a logistic regression model based on monocytes/macrophages related gene set. This study could aid in the early diagnostic of AAA.

Comprehensive RNA Expression Analysis Revealed Biological Functions of Key Gene Sets and Identified Disease-Associated Cell Types Involved in Rat Traumatic Brain Injury

Traumatic brain injury (TBI) is a worldwide public health concern without major therapeutic breakthroughs over the past decades. Developing effective treatment options and improving the prognosis of TBI depends on a better understanding of the mechanisms underlying TBI. This study performed a comprehensive analysis of 15 RNA expression datasets of rat TBIs from the GEO database. By integrating the results from the various analyses, this study investigated the biological processes, pathways, and cell types associated with TBI and explored the activity of these cells during various TBI phases. The results showed the response to cytokine, inflammatory response, bacteria-associated response, metabolic and biosynthetic processes, and pathways of neurodegeneration to be involved in the pathogenesis of TBI. The cellular abundance of microglia, perivascular macrophages (PM), and neurons were found to differ after TBI and at different times postinjury. In conclusion, immune- and inflammation-related pathways, as well as pathways of neurodegeneration, are closely related to TBI. Microglia, PM, and neurons are thought to play roles in TBI with different activities that vary by phase of TBI.

Identification of the Key Genes and Potential Therapeutic Compounds for Abdominal Aortic Aneurysm Based on a Weighted Correlation Network Analysis

Background: There is still an unmet need for therapeutic drugs for patients with an abdominal aortic aneurysm (AAA), especially for candidates unsuitable for surgical or interventional repair. Therefore, the purpose of this in silico study is to identify significant genes and regulatory mechanisms in AAA patients to predicate the potential therapeutic compounds for significant genes. Methods: The GSE57691 dataset was obtained from Gene Expression Omnibus (GEO) and used to identify the differentially expressed genes (DEGs) and weighted correlation network analysis (WGCNA). The biological function of DEGs was determined using gene ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). AAA-related genes were obtained from the Comparative Toxicogenomics Database (CTD) using the keywords: aortic aneurysm and abdominal. The hub genes in AAA were obtained by overlapping DEGs, WGCNA-based hub genes, and CTD-based genes. The diagnostic values of hub genes were determined using ROC curve analysis. Hereby, a TF-miRNA-hub gene network was constructed based on the miRnet database. Using these data, potential therapeutic compounds for the therapy of AAA were predicted based on the Drug Gene Interaction Database (DGIdb). Results: A total of 218 DEGs (17 upregulated and 201 downregulated) and their biological function were explored; 4093 AAA-related genes were derived by text mining. Three hub modules and 144 hub genes were identified by WGCNA. asparagine synthetase (ASNS), axin-related protein 2 (AXIN2), melanoma cell adhesion molecule (MCAM), and the testis-specific Y-encoded-like protein 1 (TSPYL1) were obtained as intersecting hub genes and the diagnostic values were confirmed with ROC curves. As potential compounds targeting the hub genes, asparaginase was identified as the target compound for ASNS. Prednisolone and abiraterone were identified as compounds targeting TSPYL1. For MCAM and TSPYL1, no potential therapeutic compound could be predicted. Conclusion: Using WGCNA analysis and text mining, pre-existing gene expression data were used to provide novel insight into potential AAA-related protein targets. For two of these targets, compounds could be predicted.