Identification of Hub Genes and MicroRNAs Associated With Idiopathic Pulmonary Arterial Hypertension by Integrated Bioinformatics Analyses

Molecular crosstalk between COVID-19 and Alzheimer's disease using microarray and RNA-seq datasets: A system biology approach

Objective

Coronavirus disease 2019 (COVID-19) is an infectious disease caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2). The clinical and epidemiological analysis reported the association between SARS-CoV-2 and neurological diseases. Among neurological diseases, Alzheimer's disease (AD) has developed as a crucial comorbidity of SARS-CoV-2. This study aimed to understand the common transcriptional signatures between SARS-CoV-2 and AD.

Materials and methods

System biology approaches were used to compare the datasets of AD and COVID-19 to identify the genetic association. For this, we have integrated three human whole transcriptomic datasets for COVID-19 and five microarray datasets for AD. We have identified differentially expressed genes for all the datasets and constructed a protein-protein interaction (PPI) network. Hub genes were identified from the PPI network, and hub genes-associated regulatory molecules (transcription factors and miRNAs) were identified for further validation.

Results

A total of 9,500 differentially expressed genes (DEGs) were identified for AD and 7,000 DEGs for COVID-19. Gene ontology analysis resulted in 37 molecular functions, 79 cellular components, and 129 biological processes were found to be commonly enriched in AD and COVID-19. We identified 26 hub genes which includes AKT1, ALB, BDNF, CD4, CDH1, DLG4, EGF, EGFR, FN1, GAPDH, INS, ITGB1, ACTB, SRC, TP53, CDC42, RUNX2, HSPA8, PSMD2, GFAP, VAMP2, MAPK8, CAV1, GNB1, RBX1, and ITGA2B. Specific miRNA targets associated with Alzheimer's disease and COVID-19 were identified through miRNA target prediction. In addition, we found hub genes-transcription factor and hub genes-drugs interaction. We also performed pathway analysis for the hub genes and found that several cell signaling pathways are enriched, such as PI3K-AKT, Neurotrophin, Rap1, Ras, and JAK-STAT.

Conclusion

Our results suggest that the identified hub genes could be diagnostic biomarkers and potential therapeutic drug targets for COVID-19 patients with AD comorbidity.

Identification of immune-related hub genes and analysis of infiltrated immune cells of idiopathic pulmonary artery hypertension

Objectives

Idiopathic pulmonary artery hypertension (IPAH) is a rare but life-threaten disease. However, the mechanism underlying IPAH is unclear. In this study, underlying mechanism, infiltration of immune cells, and immune-related hub genes of IPAH were analyzed via bioinformatics.

Methods

GSE15197, GSE48149, GSE113439, and GSE117261 were merged as lung dataset. Weighted gene correlation network analysis (WGCNA) was used to construct the co-expression gene networks of IPAH. Gene Ontology and pathway enrichment analysis were performed using DAVID, gene set enrichment analysis (GSEA), and gene set variation analysis (GSVA). Infiltration of immune cells in lung samples was analyzed using CIBERSORT. GSE22356 and GSE33463 were merged as peripheral blood mononuclear cells (PBMCs) dataset. Immune-related differentially expressed genes (IRDEGs) of lung and PBMCs dataset were analyzed. Based on the intersection between two sets of IRDEGs, hub genes were screened using machine learning algorithms and validated by RT-qPCR. Finally, competing endogenous RNA (ceRNA) networks of hub genes were constructed.

Results

The gray module was the most relevant module and genes in the module enriched in terms like inflammatory and immune responses. The results of GSEA and GSVA indicated that increasement in cytosolic calcium ion, and metabolism dysregulation play important roles in IPAH. The proportions of T cells CD4 memory resting and macrophage M1 were significantly greater in IPAH group, while the proportions of monocytes and neutrophils were significantly lower in IPAH group. IRDEGs of two datasets were analyzed and the intersection between two set of IRDEGs were identified as candidate hub genes. Predictive models for IPAH were constructed using data from PBMCs dataset with candidate hub genes as potential features via LASSO regression and XGBoost algorithm, respectively. CXCL10 and VIPR1 were identified as hub genes and ceRNA networks of CXCL10 was constructed.

Conclusion

Inflammatory response, increasement in cytosolic calcium ion, and metabolism dysregulation play important roles in IPAH. T cells CD4 memory resting and macrophage M1 were significantly infiltrated in lung samples from patients with IPAH. IRDEGs of lung dataset and PBMCs dataset were analyzed, and CXCL10 and VIPR1 were identified as hub genes.

Hub Genes and Immune Cell Infiltration in Hypoxia-Induced Pulmonary Hypertension: Bioinformatics Analysis and In Vivo Validation

BACKGROUND

Hypoxia-induced pulmonary hypertension (HPH) represents a severe pulmonary disorder with high morbidity and mortality, which necessitates identifying the critical molecular mechanisms underlying HPH pathogenesis.

METHODS

The mRNA expression microarray GSE15197 (containing 8 pulmonary tissues from HPH and 13 normal controls) was downloaded from Gene Expression Omnibus (GEO). Gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) were executed by RStudio software. The Protein-Protein Interaction (PPI) network was visualized and established using Cytoscape, and the cytoHubba app from Cytoscape was used to pick out the hub modules. The infiltration of immune cells in HPH was analyzed using the CIBERSORTx. To confirm the potential hub genes, real-time quantitative reverse transcription PCR (qRT-PCR) was conducted using lung tissues of rat HPH models and controls.

RESULTS

A total of 852 upregulated and 547 downregulated genes were identified. The top terms in biological processes were apoptosis, proliferation, and regulation of the MAPK cascade, including ERK1/2. Cytoplasm, cytosol, and membrane were enriched in cellular component groups. Molecular functions mainly focus on protein binding, protein serine/threonine kinase activity and identical protein binding. KEGG analysis identified pathways in cancer, regulation of actin cytoskeleton and rap1 signaling pathway. There was significantly different immune cell infiltration between HPH and normal control samples. High proportions of the memory subsets of B cells and CD4 cells, Macrophages M2 subtype, and resting Dendritic cells were found in HPH samples, while high proportions of naive CD4 cells and resting mast cells were found in normal control samples. The qRT-PCR results showed that among the ten identified hub modules, FBXL3, FBXL13 and XCL1 mRNA levels were upregulated, while NEDD4L, NPFFR2 and EDN3 were downregulated in HPH rats compared with control rats.

CONCLUSION

Our study revealed the key genes and the involvement of immune cell infiltration in HPH, thus providing new insight into the pathogenesis of HPH and potential treatment targets for patients with HPH.

The identification and verification of hub genes associated with pulmonary arterial hypertension using weighted gene co-expression network analysis

BACKGROUND

Pulmonary arterial hypertension (PAH) is characterized by a progressive increase in pulmonary vascular resistance and pulmonary arterial pressure, with complex etiology, difficult treatment and poor prognosis. The objective of this study was to investigate the potential biomarkers for PAH based on bioinformatics analysis.

METHODS

The GSE117261 datasets were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified by screening PAH patients and controls. Then the DEGs were analyzed using a Weighted Gene Co-expression Network Analysis (WGCNA) and the key modules were determined, and to further explore their potential biological functions via Gene Ontology analysis (GO), Kyoto Encyclopedia of Genes and Genomes Pathway analysis (KEGG), and Gene Set Enrichment Analysis (GSEA). Moreover, Protein-protein interaction (PPI) networks were constructed to identify hub gene candidates in the key modules. Finally, real-time quantitative polymerase chain reaction was supplied to detect the expressions of hub genes in human pulmonary arterial smooth cells treated with cobalt chloride (COCl₂) which was used to mimic hypoxia.

RESULTS

There were 2299 DEGs identified. WGCNA indicated that yellow module was the key one correlated with PAH. GO and KEGG analysis demonstrated that genes in the yellow module were mainly enriched in 'Pathways in cancer'. GSEA revealed that 'HALLMARK_MYC_TARGETS_V1' was remarkably enriched in PAH. Based on the PPI network, vascular endothelial growth factor A, proto-oncogene receptor tyrosine kinase (KIT), PNN interacting serine and arginine rich protein (PNISR) and heterogeneous nuclear ribonucleoprotein H1 (HNRNPH1) were identified as the hub genes. Additionally, the PCR indicated that the elevated expressions of PNISR and HNRNPH1 were in line with the bioinformatics analysis. ROC analysis determined that PNISR and HNRNPH1 may be potential biomarkers to provide better diagnosis of PAH.

CONCLUSION

PNISR and HNRNPH1 were potential biomarkers to diagnosis PAH. In summary, the identified DEGs, modules, pathways, and hub genes provide clues and shed light on the potential molecular mechanisms of PAH.

Connections for Matters of the Heart: Network Medicine in Cardiovascular Diseases

Cardiovascular diseases (CVD) are diverse disorders affecting the heart and vasculature in millions of people worldwide. Like other fields, CVD research has benefitted from the deluge of multiomics biomedical data. Current CVD research focuses on disease etiologies and mechanisms, identifying disease biomarkers, developing appropriate therapies and drugs, and stratifying patients into correct disease endotypes. Systems biology offers an alternative to traditional reductionist approaches and provides impetus for a comprehensive outlook toward diseases. As a focus area, network medicine specifically aids the translational aspect of in silico research. This review discusses the approach of network medicine and its application to CVD research.

Identifying Key Biomarkers and Immune Infiltration in Female Patients with Ischemic Stroke Based on Weighted Gene Co-Expression Network Analysis

Stroke is one of the leading causes of death and disability worldwide. Evidence shows that ischemic stroke (IS) accounts for nearly 80 percent of all strokes and that the etiology, risk factors, and prognosis of this disease differ by gender. Female patients may bear a greater burden than male patients. The immune system may play an important role in the pathophysiology of females with IS. Therefore, it is critical to investigate the key biomarkers and immune infiltration of female IS patients to develop effective treatment methods. Herein, we used weighted gene co-expression network analysis (WGCNA) to determine the key modules and core genes in female IS patients using the GSE22255, GSE37587, and GSE16561 datasets from the GEO database. Subsequently, we performed functional enrichment analysis and built a protein-protein interaction (PPI) network. Ten genes were selected as the true central genes for further investigation. After that, we explored the specific molecular and biological functions of these hub genes to gain a better understanding of the underlying pathogenesis of female IS patients. Moreover, the “Cell type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT)” was used to examine the distribution pattern of immune subtypes in female patients with IS and normal controls, revealing a new potential target for clinical treatment of the disease.

Stress Reactivity, Susceptibility to Hypertension, and Differential Expression of Genes in Hypertensive Compared to Normotensive Patients

Although half of hypertensive patients have hypertensive parents, known hypertension-related human loci identified by genome-wide analysis explain only 3% of hypertension heredity. Therefore, mainstream transcriptome profiling of hypertensive subjects addresses differentially expressed genes (DEGs) specific to gender, age, and comorbidities in accordance with predictive preventive personalized participatory medicine treating patients according to their symptoms, individual lifestyle, and genetic background. Within this mainstream paradigm, here, we determined whether, among the known hypertension-related DEGs that we could find, there is any genome-wide hypertension theranostic molecular marker applicable to everyone, everywhere, anytime. Therefore, we sequenced the hippocampal transcriptome of tame and aggressive rats, corresponding to low and high stress reactivity, an increase of which raises hypertensive risk; we identified stress-reactivity-related rat DEGs and compared them with their known homologous hypertension-related animal DEGs. This yielded significant correlations between stress reactivity-related and hypertension-related fold changes (log2 values) of these DEG homologs. We found principal components, PC1 and PC2, corresponding to a half-difference and half-sum of these log2 values. Using the DEGs of hypertensive versus normotensive patients (as the control), we verified the correlations and principal components. This analysis highlighted downregulation of β-protocadherins and hemoglobin as whole-genome hypertension theranostic molecular markers associated with a wide vascular inner diameter and low blood viscosity, respectively.