Drug-gene Network Research Articles

The role of m5C RNA methylation regulators in the diagnosis and immune microenvironment of osteoarthritis

The 5-methylcytosine (m5C) is a common post-transcriptional RNA methylation modification and is involved in the pathological process of many diseases. However, little is known about the role of m5C in osteoarthritis (OA). OA gene data and the corresponding information were downloaded from the Gene Expression Omnibus database. Based on 36 m5C regulators, we constructed the landscape and diagnostic model for OA. Later, two m5C modification patterns were identified, and functional analyses were performed to evaluate whether these patterns were related to endoplasmic reticulum (ER) stress and mitochondrial autophagy. We further comprehensively analyzed the immune cell infiltration characteristics in different modification patterns in OA. We also established the post-transcriptional regulatory networks and drug-gene networks. Our findings suggested that m5C regulators were differentially expressed between OA and normal samples and could serve as novel biomarkers for the diagnosis of OA. Besides, m5C regulators may be involved in regulating ER stress, mitochondrial autophagy, and immune infiltration in OA. The m5C modification can influence the sensitivity to drugs and the potential post-transcriptional regulatory mechanisms might provide promising targets.

Molecular characterization of PANoptosis-related genes in chronic kidney disease.

Chronic kidney disease (CKD) is characterized by fibrosis and inflammation in renal tissues. Several types of cell death have been implicated in CKD onset and progression. Unlike traditional forms of cell death, PANoptosis is characterized by the crosstalk among programmed cell death pathways. However, the interaction between PANoptosis and CKD remains unclear. Here, we used bioinformatics methods to identify differentially expressed genes and differentially expressed PANoptosis-related genes (DE-PRGs) using data from the GSE37171 dataset. Following this, we further performed gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, and gene set enrichment analysis using the data. We adopted a combined approach to select hub genes, using the STRING database and CytoHubba plug-in, and we used the GSE66494 as a validation dataset. In addition, we constructed ceRNA, transcription factor (TF)-gene, and drug-gene networks using Cytoscape. Lastly, we conducted immunohistochemical analysis and western blotting to validate the hub genes. We identified 57 PANoptosis-associated genes as DE-PRGs. We screened nine hub genes from the 57 DE-PRGs. We identified two hub genes (FOS and PTGS2) using the GSE66494 database, Nephroseq, immunohistochemistry, and western blotting. A common miRNA (Hsa-miR-101-3p) and three TFs (CREB1, E2F1, and RELA) may play a crucial role in the onset and progression of PANoptosis-related CKD. In our analysis of the drug-gene network, we identified eight drugs targeting FOS and 52 drugs targeting PTGS2.

Molecular mechanism underlying the protective effects of ischemic preconditioning in total knee arthroplasty

ProposeTo investigate the molecular mechanisms underlying the protective effects of ischemic preconditioning (IPC) in patients undergoing total knee arthroplasty. MethodsGSE21164 was extracted from an online database, followed by an investigation of differentially expressed genes (DEGs) between IPC treatment samples at 2 time points (T0T and T1T). Function and pathway enrichment analyses were performed on the DEGs. A protein-protein interaction network was constructed to identify hub genes according to 5 different algorithms, followed by enrichment analysis. In addition, long noncoding RNAs (lncRNAs) were identified between the T0T and T1T samples. Furthermore, a competing endogenous RNA network was predicted based on the identified lncRNA-messenger RNA (mRNA), lncRNA-microRNA (miRNA), and mRNA-miRNA relationships revealed in this study. Finally, a drug-gene network was investigated. Statistical analyses were performed using GraphPad Prism 8.0. Differences between groups were determined using an unpaired t-test. p < 0.05 was considered significant. ResultsA total of 343 DEGs at T0 and 10 DEGs at T1 were identified and compared with their respective control groups, followed by 100 DEGs between T0T and T1T. Based on these 100 DEGs, protein-protein interaction network analysis revealed 9 hub genes, mainly with mitochondria-related functions and the carbon metabolism pathway. Six differentially expressed lncRNAs were investigated between T0T and T1T. A competing endogenous RNA network was constructed using 259 lncRNA–miRNA–mRNA interactions, including alpha-2-macroglobulin antisense RNA 1-miR-7161-5p-iron-sulfur cluster scaffold. Finally, 13 chemical drugs associated with the hub genes were explored. ConclusionIron-sulfur cluster scaffold may promote IPC-induced ischemic tolerance mediated by alpha-2-macroglobulin antisense RNA 1-miR-7161-5p axis. Moreover, IPC may induce a protective response after total knee arthroplasty via mitochondria-related functions and the carbon metabolism pathway, which should be further validated in the near future.

Identifying Key Genes and Approved Medications Associated with Major Depressive Disorder Using Network Analysis and Systems Biology

Objective: Major depressive disorder (MDD) stands as one of the serious psychiatric conditions that detrimentally affect patients' quality of life and leads to a significant part of disability worldwide. Due to the limited understanding of the basic molecular mechanisms of depression and antidepressant medications, a clear understanding of the onset and development of MDD is unavailable. This study aims to figure out the pivotal genes and pathways implicated in the MDD development and identify medications that can potentially improve MDD treatment based on their relation with the key genes. Method: Symbols of human coding genes were retrieved from the HUGO Gene Nomenclature Committee database. These symbols were then queried for MDD-related associations using a Python script in PubMed. Subsequently, genes with two or more related articles to MDD were selected. A union of our search data and MDD-related genes in the DisGeNET database was found. The gene interaction network was generated and analyzed utilizing the STRING and Cytoscape, respectively. Finally, a drug-gene network was constructed and medications that can affect multiple genes were selected. Results: The union of our search data and DisGeNET data contained 1734 genes. Based on network analysis, TNF, IL1B, IL6, STAT1, and STAT3 were identified as the key genes in the MDD pathogenesis. Eleven drugs that affect more than one gene were detected through a drug-gene network. These medications include Acitretin, Adalimumab, Alteplase, Cisplatin, Digoxin, Etanercept, Infliximab, Insulin, Omeprazole, Pentoxifylline, and Rabeprazole. Conclusion: In summary, our findings identified five genes as key genes in MDD development, as well as medications related to key genes. This study provides a new vision of the pathogenesis and treatment of MDD. However, further experimental and clinical studies are necessary.

Identification of ibuprofen targeting CXCR family members to alleviate metabolic disturbance in lipodystrophy based on bioinformatics and in vivo experimental verification.

Lipodystrophy is a rare disease that is poorly diagnosed due to its low prevalence and frequent phenotypic heterogeneity. The main therapeutic measures for patients with clinical lipodystrophy are aimed at improving general metabolic complications such as diabetes mellitus, insulin resistance, and hypertriglyceridemia. Therefore, there is an urgent need to find new biomarkers to aid in the diagnosis and targeted treatment of patients with congenital generalized lipodystrophy (CGL). Dataset GSE159337 was obtained via the Gene Expression Omnibus database. First, differentially expressed genes (DEGs) between CGL and control samples were yielded via differential expression analysis and were analyzed for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment to explore the functional pathways. Next, protein-protein interaction analysis and the MCC algorithm were implemented to yield candidate genes, which were then subjected to receiver operating characteristic (ROC) analysis to identify biomarkers with an area under the curve value exceeding 0.8. Moreover, random forest (RF), logistic regression, and support vector machine (SVM) analyses were carried out to assess the diagnostic ability of biomarkers for CGL. Finally, the small-molecule drugs targeting biomarkers were predicted, and ibuprofen was further validated in lipodystrophy mice. A total of 71 DEGs in GSE159337 were sifted out and were involved in immune receptor activity, immune response-regulating signaling pathway, and secretory granule membrane. Moreover, CXCR2, TNFSF10, NLRC4, CCR2, CEACAM3, TLR10, TNFAIP3, and JUN were considered as biomarkers by performing ROC analysis on 10 candidate genes. Meanwhile, RF, logistic regression, and SVM analyses further described that those biomarkers had an excellent diagnosis capability for CGL. Eventually, the drug-gene network included ibuprofen-CXCR1, ibuprofen-CXCR1, cenicriviroc-CCR2, fenofibrate-JUN, and other relationship pairs. Ibuprofen treatment was also validated to downregulate CXCR1 and CXCR2 in peripheral blood mononuclear cells (PBMCs) and improve glucose tolerance, hypertriglyceridemia, hepatic steatosis, and liver inflammation in lipodystrophy mice. Eight biomarkers, namely, CXCR2, TNFSF10, NLRC4, CCR2, CEACAM3, TLR10, TNFAIP3, and JUN, were identified through bioinformatic analyses, and ibuprofen targeting CXCR1 and CXCR2 in PBMCs was shown to improve metabolic disturbance in lipodystrophy, contributing to studies related to the diagnosis and treatment of lipodystrophy.

CDK4/6 inhibitors in drug-induced liver injury: a pharmacovigilance study of the FAERS database and analysis of the drug-gene interaction network.

The aim of this study was to investigate the potential risk of drug-induced liver injury (DILI) caused by the CDK4/6 inhibitors (CDK4/6is abemaciclib, ribociclib, and palbociclib by comprehensively analyzing the FDA Adverse Event Reporting System (FAERS) database. Moreover, potential toxicological mechanisms of CDK4/6is-related liver injury were explored via drug-gene network analysis. In this retrospective observational study, we collected reports of DILI associated with CDK4/6i use from the FAERS dated January 2014 to March 2023. We conducted disproportionality analyses using the reporting odds ratio (ROR) with a 95% confidence interval (CI). Pathway enrichment analysis and drug-gene network analyses were subsequently performed to determine the potential mechanisms underlying CDK4/6i-induced liver injury. We found positive signals for DILI with ribociclib (ROR = 2.60) and abemaciclib (ROR = 2.37). DILIs associated with liver-related investigations, signs, and symptoms were confirmed in all three reports of CDK4/6is. Moreover, ascites was identified as an unlisted hepatic adverse effect of palbociclib. We isolated 189 interactive target genes linking CDK4/6 inhibitors to hepatic injury. Several key genes, such as STAT3, HSP90AA1, and EP300, were revealed via protein-protein analysis, emphasizing their central roles within the network. KEGG pathway enrichment of these genes highlighted multiple pathways. Our study revealed variations in hepatobiliary toxicity among the different CDK4/6 inhibitors, with ribociclib showing the highest risk of liver injury, followed by abemaciclib, while palbociclib appeared relatively safe. Our findings emphasize the need for cautious use of CDK4/6 inhibitors, and regular liver function monitoring is recommended for long-term CDK4/6 inhibitor use.

Identification of m6A-associated autophagy genes in non-alcoholic fatty liver.

Studies had shown that autophagy was closely related to nonalcoholic fat liver disease (NAFLD), while N6-methyladenosine (m6A) was involved in the regulation of autophagy. However, the mechanism of m6A related autophagy in NAFLD was unclear. The NAFLD related datasets were gained via the Gene Expression Omnibus (GEO) database, and we also extracted 232 autophagy-related genes (ARGs) and 37 m6A. First, differentially expressed ARGs (DE-ARGs) and differentially expressed m6A (DE-m6A) were screened out by differential expression analysis. DE-ARGs associated with m6A were sifted out by Pearson correlation analysis, and the m6A-ARGs relationship pairs were acquired. Then, autophagic genes in m6A-ARGs pairs were analyzed for machine learning algorithms to obtain feature genes. Further, we validated the relationship between feature genes and NAFLD through quantitative real-time polymerase chain reaction (qRT-PCR), Western blot (WB). Finally, the immuno-infiltration analysis was implement, and we also constructed the TF-mRNA and drug-gene networks. There were 19 DE-ARGs and four DE-m6A between NAFLD and normal samples. The three m6A genes and five AGRs formed the m6A-ARGs relationship pairs. Afterwards, genes obtained from machine learning algorithms were intersected to yield three feature genes (TBK1, RAB1A, and GOPC), which showed significant positive correlation with astrocytes, macrophages, smooth muscle, and showed significant negative correlation with epithelial cells, and endothelial cells. Besides, qRT-PCR and WB indicate that TBK1, RAB1A and GOPC significantly upregulated in NAFLD. Ultimately, we found that the TF-mRNA network included FOXP1-GOPC, ATF1-RAB1A and other relationship pairs, and eight therapeutic agents such as R-406 and adavosertib were predicted based on the TBK1. The study investigated the potential molecular mechanisms of m6A related autophagy feature genes (TBK1, RAB1A, and GOPC) in NAFLD through bioinformatic analyses and animal model validation. However, it is critical to note that these findings, although consequential, demonstrate correlations rather than cause-and-effect relationships. As such, more research is required to fully elucidate the underlying mechanisms and validate the clinical relevance of these feature genes.

Time-series analysis of hematopoietic stem cells.

This study aimed to investigate the molecular mechanisms underlying the aging of hematopoietic stem cells (HSCs). Gene expression profile GSE32719 was downloaded from the Gene Expression Omnibus database, including 14 young, 5 middle, and 8 old HSCs. Differential expression analysis, short time-series expression miner analysis, and weighted co-expression network analysis were conducted to screen for hub genes whose expression changed over time during HSC aging. Subsequently, functional enrichment and multiple regulatory network analyses of the hub genes were performed. A total of 124 intersecting time-dependent differentially expressed and module genes were obtained, which were considered hub genes whose expression changed over time during HSC aging. Hub genes were significantly enriched in pathways such as the Hippo and AMP-activated protein kinase (AMPK) signaling pathways. Moreover, AP-1 Transcription Factor Subunit (FOS) and sirtuin 1 (SIRT1) had higher degrees in the protein-protein interaction network, were regulated by more transcription factors (TFs), such as Sp1 transcription factor (SP1) and BRCA1 DNA repair-associated (BRCA1), in the TF-mRNA-miRNA network, were associated with more diseases in the disease-gene network, and could be targeted by more drugs in the drug-gene network. Furthermore, SIRT1 was targeted by miR-9-5p in the TF-mRNA-miRNA network. Hub genes such as FOS and SIRT1 and key pathways such as the Hippo and AMPK signaling pathways may play crucial roles in HSC aging. Moreover, FOS and SIRT1 were regulated by SP1 and BRCA1, respectively, during HSC aging. Furthermore, miR-9-5p may modulate HSC aging by targeting SIRT1. Thus, FOS and SIRT1 may be potential therapeutic targets for age-related hematopoietic dysfunction.

To Develop Biomarkers for Diabetic Nephropathy Based on Genes Related to Fibrosis and Propionate Metabolism and Their Functional Validation.

Propionate metabolism is important in the development of diabetes, and fibrosis plays an important role in diabetic nephropathy (DN). However, there are no studies on biomarkers related to fibrosis and propionate metabolism in DN. Hence, the current research is aimed at evaluating biomarkers associated with fibrosis and propionate metabolism and to explore their effect on DN progression. The GSE96804 (DN : control = 41 : 20) and GSE104948 (DN : control = 7 : 18) DN-related datasets and 924 propionate metabolism-related genes (PMRGs) and 656 fibrosis-related genes (FRGs) were acquired from the public database. First, DN differentially expressed genes (DN-DEGs) between the DN and control samples were sifted out via differential expression analysis. The PMRG scores of the DN samples were calculated based on PMRGs. The samples were divided into the high and low PMRG score groups according to the median scores. The PM-DEGs between the two groups were screened out. Second, the intersection of DN-DEGs, PM-DEGs, and FRGs was taken to yield intersected genes. Random forest (RF) and recursive feature elimination (RFE) analyses of the intersected genes were performed to sift out biomarkers. Then, single gene set enrichment analysis was conducted. Finally, immunoinfiltrative analysis was performed, and the transcription factor (TF)-microRNA (miRNA)-mRNA regulatory network and the drug-gene interaction network were constructed. There were 2633 DN-DEGs between the DN and control samples and 515 PM-DEGs between the high and low PMRG score groups. In total, 10 intersected genes were gained after taking the intersection of DN-DEGs, PM-DEGs, and FRGs. Seven biomarkers, namely, SLC37A4, ACOX2, GPD1, angiotensin-converting enzyme 2 (ACE2), SLC9A3, AGT, and PLG, were acquired via RF and RFE analyses, and they were found to be involved in various mechanisms such as glomerulus development, fatty acid metabolism, and peroxisome. The seven biomarkers were positively correlated with neutrophils. Moreover, 8 TFs, 60 miRNAs, and 7 mRNAs formed the TF-miRNA-mRNA regulatory network, including USF1-hsa-mir-1296-5p-AGT and HIF1A-hsa-mir-449a-5p-ACE2. The drug-gene network contained UROKINASE-PLG, ATENOLOL-AGT, and other interaction relationship pairs. Via bioinformatic analyses, the risk of fibrosis and propionate metabolism-related biomarkers in DN were explored, thereby providing novel ideas for research related to DN diagnosis and treatment.

Screening of Biomarkers Associated with Osteoarthritis Aging Genes and Immune Correlation Studies.

Osteoarthritis (OA) is a joint disease with a long and slow course, which is one of the major causes of disability in middle and old-aged people. This study was dedicated to excavating the cellular senescence-associated biomarkers of OA. The Gene Expression Omnibus (GEO) database was searched and five datasets pertaining to OA were obtained. After removing the batch effect, the GSE55235, GSE55457, GSE82107, and GSE12021 datasets were integrated together for screening of the candidate genes by differential analysis and weighted gene co-expression network analysis (WGCNA). Next, those genes were further filtered by machine learning algorithms to obtain cellular senescence-associated biomarkers of OA. Subsequently, enrichment analyses based on those biomarkers were conducted, and we profiled the infiltration levels of 22 types immune cells with the ERSORT algorithm. A lncRNA-miRNA-mRNA regulatory and drug-gene network were constructed. Finally, we validated the senescence-associated biomarkers at both in vivo and in vitro levels. Five genes (BCL6, MCL1, SLC16A7, PIM1, and EPHA3) were authenticated as cellular senescence-associated biomarkers in OA. ROC curves demonstrated the reliable capacity of the five genes as a whole to discriminate OA samples from normal samples. The nomogram diagnostic model based on 5 genes proved to be a reliable predictor of OA. Single-gene GSEA results pointed to the involvement of the five biomarkers in immune-related pathways and oxidative phosphorylation in the development of OA. Immune infiltration analysis manifested that the five genes were significantly correlated with differential immune cells. Subsequently, a lncRNA-miRNA-mRNA network and gene-drug network containing were generated based on five cellular senescence-associated biomarkers in OA. A foundation for understanding the pathophysiology of OA and new insights into OA diagnosis and treatment were provided by the identification of five genes, namely BCL6, MCL1, SLC16A7, PIM1, and EPHA3, as biomarkers associated with cellular senescence in OA.

Authentication and validation of key genes in the treatment of atopic dermatitis with Runfuzhiyang powder: combined RNA-seq, bioinformatics analysis, and experimental research.

Atopic dermatitis (AD) is inflammatory disease. So far, therapeutic mechanism of Runfuzhiyang powder on AD remains to be studied. This study aimed to mine key biomarkers to explore potential molecular mechanism for AD incidence and Runfuzhiyang powder treatment. The control group, AD group, treat group (AD mice treated with Runfuzhiyang powder were utilized for studying. Differentially expressed AD-related genes were acquired by intersecting of key module genes related to control group, AD group and treatment group which were screened by WGCNA and AD-related differentially expressed genes (DEGs). KEGG and GO analyses were further carried out. Next, LASSO regression analysis was utilized to screen feature genes. The ROC curves were applied to validate the diagnostic ability of feature genes to obtain AD-related biomarkers. Then protein-protein interaction (PPI) network, immune infiltration analysis and single-gene gene set enrichment analysis (GSEA) were presented. Finally, TF-mRNA-lncRNA and drug-gene networks of biomarkers were constructed. 4 AD-related biomarkers (Ddit4, Sbf2, Senp8 and Zfp777) were identified in AD groups compared with control group and treat group by LASSO regression analysis. The ROC curves revealed that four biomarkers had good distinguishing ability between AD group and control group, as well as AD group and treatment group. Next, GSEA revealed that pathways of E2F targets, KRAS signaling up and inflammatory response were associated with 4 biomarkers. Then, we found that Ddit4, Sbf2 and Zfp777 were significantly positively correlated with M0 Macrophage, and were significantly negatively relevant to Resting NK. Senp8 was the opposite. Finally, a TF-mRNA-lncRNA network including 200 nodes and 592 edges was generated, and 20 drugs targeting SENP8 were predicted. 4 AD-related and Runfuzhiyang powder treatment-related biomarkers (Ddit4, Sbf2, Senp8 and Zfp777) were identified, which could provide a new idea for targeted treatment and diagnosis of AD.

PathDIP 5: improving coverage and making enrichment analysis more biologically meaningful.

Pathway Data Integration Portal (PathDIP) is an integrated pathway database that was developed to increase functional gene annotation coverage and reduce bias in pathway enrichment analysis. PathDIP 5 provides multiple improvements to enable more interpretable analysis: users can perform enrichment analysis using all sources, separate sources or by combining specific pathway subsets; they can select the types of sources to use or the types of pathways for the analysis, reducing the number of resulting generic pathways or pathways not related to users' research question; users can use API. All pathways have been mapped to seven representative types. The results of pathway enrichment can be summarized through knowledge-based pathway consolidation. All curated pathways were mapped to 53 pathway ontology-based categories. In addition to genes, pathDIP 5 now includes metabolites. We updated existing databases, included two new sources, PathBank and MetabolicAtlas, and removed outdated databases. We enable users to analyse their results using Drugst.One, where a drug-gene network is created using only the user's genes in a specific pathway. Interpreting the results of any analysis is now improved by multiple charts on all the results pages. PathDIP 5 is freely available at https://ophid.utoronto.ca/pathDIP.

Cellular Heterogeneity of Activated Primary Human Macrophages and Associated Drug-Gene Networks: From Biology to Precision Therapeutics.

Interferon-γ (IFNγ) signaling plays a complex role in atherogenesis. IFNγ stimulation of macrophages permits in vitro exploration of proinflammatory mechanisms and the development of novel immune therapies. We hypothesized that the study of macrophage subpopulations could lead to anti-inflammatory interventions. Primary human macrophages activated by IFNγ (M(IFNγ)) underwent analyses by single-cell RNA sequencing, time-course cell-cluster proteomics, metabolite consumption, immunoassays, and functional tests (phagocytic, efferocytotic, and chemotactic). RNA-sequencing data were analyzed in LINCS (Library of Integrated Network-Based Cellular Signatures) to identify compounds targeting M(IFNγ) subpopulations. The effect of compound BI-2536 was tested in human macrophages in vitro and in a murine model of atherosclerosis. Single-cell RNA sequencing identified 2 major clusters in M(IFNγ): inflammatory (M(IFNγ)i) and phagocytic (M(IFNγ)p). M(IFNγ)i had elevated expression of inflammatory chemokines and higher amino acid consumption compared with M(IFNγ)p. M(IFNγ)p were more phagocytotic and chemotactic with higher Krebs cycle activity and less glycolysis than M(IFNγ)i. Human carotid atherosclerotic plaques contained 2 such macrophage clusters. Bioinformatic LINCS analysis using our RNA-sequencing data identified BI-2536 as a potential compound to decrease the M(IFNγ)i subpopulation. BI-2536 in vitro decreased inflammatory chemokine expression and secretion in M(IFNγ) by shrinking the M(IFNγ)i subpopulation while expanding the M(IFNγ)p subpopulation. BI-2536 in vivo shifted the phenotype of macrophages, modulated inflammation, and decreased atherosclerosis and calcification. We characterized 2 clusters of macrophages in atherosclerosis and combined our cellular data with a cell-signature drug library to identify a novel compound that targets a subset of macrophages in atherosclerosis. Our approach is a precision medicine strategy to identify new drugs that target atherosclerosis and other inflammatory diseases.

Identification of novel potential drugs and miRNAs biomarkers in lung cancer based on gene co-expression network analysis.

Non-small cell lung cancer (NSCLC) is an important cause of cancer-associated deaths worldwide. Therefore, the exact molecular mechanisms of NSCLC are unidentified. The present investigation aims to identify the miRNAs with predictive value in NSCLC. The two datasets were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed miRNAs (DEmiRNA) and mRNAs (DEmRNA) were selected from the normalized data. Next, miRNA-mRNA interactions were determined. Then, co-expression network analysis was completed using the WGCNA package in R software. The co-expression network between DEmiRNAs and DEmRNAs was calculated to prioritize the miRNAs. Next, the enrichment analysis was performed for DEmiRNA and DEmRNA. Finally, the drug-gene interaction network was constructed by importing the gene list to dgidb database. A total of 3,033 differentially expressed genes and 58 DE miRNA were recognized from two datasets. The co-expression network analysis was utilized to build a gene co-expression network. Next, four modules were selected based on the Zsummary score. In the next step, a bipartite miRNA-gene network was constructed and hub miRNAs (let-7a-2-3p, let-7d-5p, let-7b-5p, let-7a-5p, and let-7b-3p) were selected. Finally, a drug-gene network was constructed while SUNITINIB, MEDROXYPROGESTERONE ACETATE, DOFETILIDE, HALOPERIDOL, and CALCITRIOL drugs were recognized as a beneficial drug in NSCLC. The hub miRNAs and repurposed drugs may act a vital role in NSCLC progression and treatment, respectively; however, these results must validate in further clinical and experimental assessments.

COVID-19: A novel holistic systems biology approach to predict its molecular mechanisms (in vitro) and repurpose drugs.

COVID-19 strangely kills some youth with no history of physical weakness, and in addition to the lungs, it may even directly harm other organs. Its complex mechanism has led to the loss of any significantly effective drug, and some patients with severe forms still die daily. Common methods for identifying disease mechanisms and drug design are often time-consuming or reductionist. Here, we use a novel holistic systems biology approach to predict its molecular mechanisms (in vitro), significant molecular relations with SARS, and repurpose drugs. We have utilized its relative phylogenic similarity to SARS. Using the available omics data for SARS and the fewer data for COVID-19 to decode the mechanisms and their significant relations, We applied the Cytoscape analyzer, MCODE, STRING, and DAVID tools to predict the topographically crucial molecules, clusters, protein interaction mappings, and functional analysis. We also applied a novel approach to identify the significant relations between the two infections using the Fischer exact test for MCODE clusters. We then constructed and analyzed a drug-gene network using PharmGKB and DrugBank (retrieved using the dgidb). Some of the shared identified crucial molecules, BPs and pathways included Kaposi sarcoma-associated herpesvirus infection, Influenza A, and NOD-like receptor signaling pathways. Besides, our identified crucial molecules specific to host response against SARS-CoV-2 included FGA, BMP4, PRPF40A, and IFI16. We also introduced seven new repurposed candidate drugs based on the drug-gene network analysis for the identified crucial molecules. Therefore, we suggest that our newly recommended repurposed drugs be further investigated in Vitro and in Vivo against COVID-19.

Identification and validation of differentially expressed genes for targeted therapy in NSCLC using integrated bioinformatics analysis.

Despite the high prevalence of lung cancer, with a five-year survival rate of only 23%, the underlying molecular mechanisms of non-small cell lung cancer (NSCLC) remain unknown. There is a great need to identify reliable candidate biomarker genes for early diagnosis and targeted therapeutic strategies to prevent cancer progression. In this study, four datasets obtained from the Gene Expression Omnibus were evaluated for NSCLC- associated differentially expressed genes (DEGs) using bioinformatics analysis. About 10 common significant DEGs were shortlisted based on their p-value and FDR (DOCK4, ID2, SASH1, NPR1, GJA4, TBX2, CD24, HBEGF, GATA3, and DDR1). The expression of significant genes was validated using experimental data obtained from TCGA and the Human Protein Atlas database. The human proteomic data for post- translational modifications was used to interpret the mutations in these genes. Validation of DEGs revealed a significant difference in the expression of hub genes in normal and tumor tissues. Mutation analysis revealed 22.69%, 48.95%, and 47.21% sequence predicted disordered regions of DOCK4, GJA4, and HBEGF, respectively. The gene-gene and drug-gene network analysis revealed important interactions between genes and chemicals suggesting they could act as probable drug targets. The system-level network showed important interactions between these genes, and the drug interaction network showed that these genes are affected by several types of chemicals that could serve as potential drug targets. The study demonstrates the importance of systemic genetics in identifying potential drug- targeted therapies for NSCLC. The integrative system- level approach should contribute to a better understanding of disease etiology and may accelerate drug discovery for many cancer types.

Identification of key signaling pathways and hub genes related to immune infiltration in Kawasaki disease with resistance to intravenous immunoglobulin based on weighted gene co-expression network analysis.

Background: Kawasaki disease (KD) is an acute vasculitis, that is, the leading cause of acquired heart disease in children, with approximately 10%-20% of patients with KD suffering intravenous immunoglobulin (IVIG) resistance. Although the underlying mechanism of this phenomenon remains unclear, recent studies have revealed that immune cell infiltration may associate with its occurrence. Methods: In this study, we downloaded the expression profiles from the GSE48498 and GSE16797 datasets in the Gene Expression Omnibus database, analyzed differentially expressed genes (DEGs), and intersected the DEGs with the immune-related genes downloaded from the ImmPort database to obtain differentially expressed immune-related genes (DEIGs). Then CIBERSORT algorithm was used to calculate the immune cell compositions, followed by the WGCNA analysis to identify the module genes associated with immune cell infiltration. Next, we took the intersection of the selected module genes and DEIGs, then performed GO and KEGG enrichment analysis. Moreover, ROC curve validation, Spearman analysis with immune cells, TF, and miRNA regulation network, and potential drug prediction were implemented for the finally obtained hub genes. Results: The CIBERSORT algorithm showed that neutrophil expression was significantly higher in IVIG-resistant patients compared to IVIG-responsive patients. Next, we got differentially expressed neutrophil-related genes by intersecting DEIGs with neutrophil-related module genes obtained by WGCNA, for further analysis. Enrichment analysis revealed that these genes were associated with immune pathways, such as cytokine-cytokine receptor interaction and neutrophil extracellular trap formation. Then we combined the PPI network in the STRING database with the MCODE plugin in Cytoscape and identified 6 hub genes (TLR8, AQP9, CXCR1, FPR2, HCK, and IL1R2), which had good diagnostic performance in IVIG resistance according to ROC analysis. Furthermore, Spearman's correlation analysis confirmed that these genes were closely related to neutrophils. Finally, TFs, miRNAs, and potential drugs targeting the hub genes were predicted, and TF-, miRNA-, and drug-gene networks were constructed. Conclusion: This study found that the 6 hub genes (TLR8, AQP9, CXCR1, FPR2, HCK, and IL1R2) were significantly associated with neutrophil cell infiltration, which played an important role in IVIG resistance. In a word, this work rendered potential diagnostic biomarkers and prospective therapeutic targets for IVIG-resistant patients.

Novel potential drugs for the treatment of primary open-angle glaucoma using protein-protein interaction network analysis.

Glaucoma is the second leading cause of irreversible blindness, and primary open-angle glaucoma (POAG) is the most common type. Due to inadequate diagnosis, treatment is often not administered until symptoms occur. Hence, approaches enabling earlier prediction or diagnosis of POAG are necessary. We aimed to identify novel drugs for glaucoma through bioinformatics and network analysis. Data from 36 samples, obtained from the trabecular meshwork of healthy individuals and patients with POAG, were acquired from a dataset. Next, differentially expressed genes (DEGs) were identified to construct a protein-protein interaction (PPI) network. In both stages, the genes were enriched by studying the critical biological processes and pathways related to POAG. Finally, a drug-gene network was constructed, and novel drugs for POAG treatment were proposed. Genes with p < 0.01 and |log fold change| > 0.3 (1,350 genes) were considered DEGs and utilized to construct a PPI network. Enrichment analysis yielded several key pathways that were upregulated or downregulated. For example, extracellular matrix organization, the immune system, neutrophil degranulation, and cytokine signaling were upregulated among immune pathways, while signal transduction, the immune system, extracellular matrix organization, and receptor tyrosine kinase signaling were downregulated. Finally, novel drugs including metformin hydrochloride, ixazomib citrate, and cisplatin warrant further analysis of their potential roles in POAG treatment. The candidate drugs identified in this computational analysis require in vitro and in vivo validation to confirm their effectiveness in POAG treatment. This may pave the way for understanding life-threatening disorders such as cancer.

Molecular mechanism and diagnostic marker investigation of endoplasmic reticulum stress on periodontitis

PurposeThe aim of this study was to reveal the biological function of endoplasmic reticulum stress (ERS)-related genes (ERSGs) in periodontitis, and provide potential ERS diagnostic markers for clinical therapy of periodontitis.MethodsThe differentially expressed ERSGs (DE-ERSGs) were reveled based on periodontitis-related microarray dataset in Gene Expression Omnibus (GEO) database and 295 ERS in previous study, followed by a protein–protein interaction network construction. Then, the subtypes of periodontitis were explored, followed by validation with immune cell infiltration and gene set enrichment. Two machine learning algorithms were used to reveal potential ERS diagnostic markers of periodontitis. The diagnostic effect, target drug and immune correlation of these markers were further evaluated. Finally, a microRNA(miRNA)-gene interaction network was constructed.ResultsA total of 34 DE-ERSGs were revealed between periodontitis samples and control, followed by two subtypes investigated. There was a significant difference of ERS score, immune infiltration and Hallmark enrichment between two subtypes. Then, totally 7 ERS diagnostic markers including FCGR2B, XBP1, EDEM2, ATP2A3, ERLEC1, HYOU1 and YOD1 were explored, and the v the time-dependent ROC analysis showed a reliable result. In addition, a drug-gene network was constructed with 4 up-regulated ERS diagnostic markers and 24 drugs. Finally, based on 32 interactions, 5 diagnostic markers and 20 miRNAs, a miRNA-target network was constructed.ConclusionsUp-regulated miR-671-5p might take part in the progression of periodontitis via stimulating the expression of ATP2A3. ERSGs including XBP1 and FCGR2B might be novel diagnostic marker for periodontitis.

SARS-CoV-2 potential drugs, drug targets, and biomarkers: a viral-host interaction network-based analysis

COVID-19 is a global pandemic impacting the daily living of millions. As variants of the virus evolve, a complete comprehension of the disease and drug targets becomes a decisive duty. The Omicron variant, for example, has a notably high transmission rate verified in 155 countries. We performed integrative transcriptomic and network analyses to identify drug targets and diagnostic biomarkers and repurpose FDA-approved drugs for SARS-CoV-2. Upon the enrichment of 464 differentially expressed genes, pathways regulating the host cell cycle were significant. Regulatory and interaction networks featured hsa-mir-93-5p and hsa-mir-17-5p as blood biomarkers while hsa-mir-15b-5p as an antiviral agent. MYB, RRM2, ERG, CENPF, CIT, and TOP2A are potential drug targets for treatment. HMOX1 is suggested as a prognostic biomarker. Enhancing HMOX1 expression by neem plant extract might be a therapeutic alternative. We constructed a drug-gene network for FDA-approved drugs to be repurposed against the infection. The key drugs retrieved were members of anthracyclines, mitotic inhibitors, anti-tumor antibiotics, and CDK1 inhibitors. Additionally, hydroxyquinone and digitoxin are potent TOP2A inhibitors. Hydroxyurea, cytarabine, gemcitabine, sotalol, and amiodarone can also be redirected against COVID-19. The analysis enforced the repositioning of fluorouracil and doxorubicin, especially that they have multiple drug targets, hence less probability of resistance.