In silico Identification of 10 Hub Genes and an miRNA-mRNA Regulatory Network in Acute Kawasaki Disease

Identification of hub biomarkers and immune-related pathways participating in the progression of Kawasaki disease by integrated bioinformatics analysis

BACKGROUND

Kawasaki disease (KD) is a systemic vasculitis that commonly affects children and its etiology remains unknown. Growing evidence suggests that immune-mediated inflammation and immune cells in the peripheral blood play crucial roles in the pathophysiology of KD. The objective of this research was to find important biomarkers and immune-related mechanisms implicated in KD, along with their correlation with immune cells in the peripheral blood.

MATERIAL/METHODS

Gene microarray data from the Gene Expression Omnibus (GEO) was utilized in this study. Three datasets, namely GSE63881 (341 samples), GSE73463 (233 samples), and GSE73461 (279 samples), were obtained. To find intersecting genes, we employed differentially expressed genes (DEGs) analysis and weighted gene co-expression network analysis (WGCNA). Subsequently, functional annotation, construction of protein-protein interaction (PPI) networks, and Least Absolute Shrinkage and Selection Operator (LASSO) regression were performed to identify hub genes. The accuracy of these hub genes in identifying KD was evaluated using the receiver operating characteristic curve (ROC). Furthermore, Gene Set Variation Analysis (GSVA) was employed to explore the composition of circulating immune cells within the assessed datasets and their relationship with the hub gene markers.

RESULTS

WGCNA yielded eight co-expression modules, with one hub module (MEblue module) exhibiting the strongest association with acute KD. 425 distinct genes were identified. Integrating WGCNA and DEGs yielded a total of 277 intersecting genes. By conducting LASSO analysis, five hub genes (S100A12, MMP9, TLR2, NLRC4 and ARG1) were identified as potential biomarkers for KD. The diagnostic value of these five hub genes was demonstrated through ROC curve analysis, indicating their high accuracy in diagnosing KD. Analysis of the circulating immune cell composition within the assessed datasets revealed a significant association between KD and various immune cell types, including activated dendritic cells, neutrophils, immature dendritic cells, macrophages, and activated CD8 T cells. Importantly, all five hub genes exhibited strong correlations with immune cells.

CONCLUSION

Activated dendritic cells, neutrophils, and macrophages were closely associated with the pathogenesis of KD. Furthermore, the hub genes (S100A12, MMP9, TLR2, NLRC4, and ARG1) are likely to participate in the pathogenic mechanisms of KD through immune-related signaling pathways.

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.

Human microRNA-4433 (hsa-miR-4443) Targets 18 Genes to be a Risk Factor of Neurodegenerative Diseases

BACKGROUND

Neurodegenerative diseases, such as Alzheimer's disease patients (AD), Huntington's disease (HD) and Parkinson's disease (PD), are common causes of morbidity, mortality, and cognitive impairment in older adults.

OBJECTIVE

We aimed to understand the transcriptome characteristics of the cortex of neurodegenerative diseases and to provide an insight into the target genes of differently expressed microRNAs in the occurrence and development of neurodegenerative diseases.

METHODS

The Limma package of R software was used to analyze GSE33000, GSE157239, GSE64977 and GSE72962 datasets to identify the differentially expressed genes (DEGs) and microRNAs in the cortex of neurodegenerative diseases. Bioinformatics methods, such as GO enrichment analysis, KEGG enrichment analysis and gene interaction network analysis, were used to explore the biological functions of DEGs. Weighted gene co-expression network analysis (WGCNA) was used to cluster DEGs into modules. RNA22, miRDB, miRNet 2.0 and TargetScan7 databases were performed to predict the target genes of microRNAs.

RESULTS

Among 310 Alzheimer's disease (AD) patients, 157 Huntington's disease (HD) patients and 157 non-demented control (Con) individuals, 214 co-DEGs were identified. Those co-DEGs were filtered into 2 different interaction network complexes, representing immune-related genes and synapserelated genes. The WGCNA results identified five modules: yellow, blue, green, turquoise, and brown. Most of the co-DEGs were clustered into the turquoise module and blue module, which respectively regulated synapse-related function and immune-related function. In addition, human microRNA-4433 (hsa-miR-4443), which targets 18 co-DEGs, was the only 1 co-up-regulated microRNA identified in the cortex of neurodegenerative diseases.

CONCLUSION

214 DEGs and 5 modules regulate the immune-related and synapse-related function of the cortex in neurodegenerative diseases. Hsa-miR-4443 targets 18 co-DEGs and may be a potential molecular mechanism in neurodegenerative diseases' occurrence and development.

MicroRNAs in Kawasaki disease: An update on diagnosis, therapy and monitoring

Kawasaki disease (KD) is an acute autoimmune vascular disease featured with a long stage of febrile. It predominantly afflicts children under 5 years old and causes an increased risk of cardiovascular combinations. The onset and progression of KD are impacted by many aspects, including genetic susceptibility, infection, and immunity. In recent years, many studies revealed that miRNAs, a novel class of small non-coding RNAs, may play an indispensable role in the development of KD via differential expression and participation in the central pathogenesis of KD comprise of the modulation of immunity, inflammatory response and vascular dysregulation. Although specific diagnose criteria remains unclear up to date, accumulating clinical evidence indicated that miRNAs, as small molecules, could serve as potential diagnostic biomarkers and exhibit extraordinary specificity and sensitivity. Besides, miRNAs have gained attention in affecting therapies for Kawasaki disease and providing new insights into personalized treatment. Through consanguineous coordination with classical therapies, miRNAs could overcome the inevitable drug-resistance and poor prognosis problem in a novel point of view. In this review, we systematically reviewed the existing literature and summarized those findings to analyze the latest mechanism to explore the role of miRNAs in the treatment of KD from basic and clinical aspects retrospectively. Our discussion helps to better understand the pathogenesis of KD and may offer profound inspiration on KD diagnosis, treatment, and prognosis.

Integrated Microarray for Identifying the Hub mRNAs and Constructed MiRNA-mRNA Network in Coronary In-stent Restenosis

As a major complication after percutaneous coronary intervention (PCI) in patients who suffer from coronary artery disease, in-stent restenosis (ISR) poses a significant challenge for clinical management. A miRNA-mRNA regulatory network of ISR can be constructed to better reveal the occurrence of ISR. The relevant dataset from the Gene Expression Omnibus (GEO) database was downloaded, and 284 differentially expressed miRNAs (DE-miRNAs) and 849 differentially expressed mRNAs (DE-mRNAs) were identified. As predicted by online tools, 65 final functional genes (FmRNAs) were overlapping DE-mRNAs and DE-miRNAs target genes. In the biological process (BP) terms of Gene Ontology (GO) functional analysis, the FmRNAs were mainly enriched in cellular response to peptide, epithelial cell proliferation and response to peptide hormone. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the FmRNAs were mainly enriched in breast cancer, endocrine resistance and cushing syndrome. Jun Proto-Oncogene, AP-1 Transcription Factor Subunit (JUN), Insulin Like Growth Factor 1 Receptor (IGF1R), Member RAS Oncogene Family (RAB14), Specificity Protein 1 (SP1), Protein Tyrosine Phosphatase Non-Receptor Type1(PTPN1), DDB1 And CUL4 Associated Factor 10 (DCAF10), Retinoblastoma-Binding Protein 5 (RBBP5) and Eukaryotic Initiation Factor 4A-I (EIF4A1) were hub genes in the protein-protein interaction network (PPI network). The miRNA-mRNA network containing DE-miRNA and hub genes was built. Hsa-miR-139-5p-JUN, hsa-miR-324-5p-SP1 axis pairs were found in the miRNA-mRNA network, which could promote ISR development. The above results indicate that the miRNA-mRNA network constructed in ISR has a regulatory role in the development of ISR, and may provide new approaches for clinical treatment and experimental development.

MiR-144-3p induced by SP1 promotes IL-1β-induced pyroptosis in chondrocytes via PTEN/PINK1/Parkin axis

Rheumatoid arthritis (RA) often leads to functional disabilities and deformities. MiRNA plays a vital role in cell pyroptosis. Nevertheless, the function and underlying mechanism of miR-144-3p in pyroptosis during the progression of RA remains unclear. In this study, N1511 cells were stimulated with IL-1β to construct a RA model. 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide (MTT) assay was performed to assess the cell viability. Cell pyroptosis was detected by flow cytometry. The levels of inflammatory cytokines (TNF-α, IL-6, and IL-18) were assessed by enzyme-linked immunosorbent assay (ELISA). The relationship among specific protein 1 (SP1), microRNA-144-3p (miR-144-3p), and phosphatase and tensin homolog (PTEN) was explored by dual-luciferase reporter assay, RNA immunoprecipitation (RIP), and chromatin immunoprecipitation (ChIP), respectively. The level of miR-144-3p in N1511 cells was upregulated by IL-1β. MiR-144-3p knockdown inhibited IL-1β-induced pyroptosis in N1511 cells, and the expressions of NOD-like receptor family pyrin domain containing 3 (NLRP3), Cleaved caspase-1, Gasdermin D (GSDMD), and Cleaved caspase-3 in IL-1β-stimulated N1511 cells were increased. The levels of inflammatory cytokines in N1511 cells were increased by IL-1β, which were restored by miR-144-3p knockdown. MiR-144-3p knockdown abolished IL-1β-induced inactivation of putative kinase 1 (PINK1)/Parkin RBR E3 ubiquitin-protein (Parkin) signalling. Moreover, transcription factor SP1 could upregulate miR-144-3p expression and miR-144-3p negatively regulated PTEN expression. In summary, MiR-144-3p induced by SP1 could promote IL-1β-induced chondrocyte pyroptosis via inhibiting PTEN expression and suppressing the activation of PINK1/Parkin signalling, which provided a new strategy against RA.