Integration of the Drug-Gene Interaction Database (DGIdb 4.0) with open crowdsource efforts

Identification of mitophagy-related hub genes during the progression of spinal cord injury by integrated multinomial bioinformatics analysis

Spinal cord injury (SCI) is a disturbance of peripheral and central nerve conduction that causes disability in sensory and motor function. Currently, there is no effective treatment for SCI. Mitophagy plays a vital role in mitochondrial quality control during various physiological and pathological processes. The study aimed to elucidate the role of mitophagy and identify potential mitophagy-related hub genes in SCI pathophysiology. Two datasets (GSE15878 and GSE138637) were analyzed. Firstly, the differentially expressed genes (DEGs) were identified and mitophagy-related genes were obtained from GeneCards, then the intersection between SCI and mitophagy-related genes was determined. Next, we performed gene set enrichment analysis (GSEA), weighted gene co-expression network analysis (WGCNA), protein-protein interaction network (PPI network), least absolute shrinkage and selection operator (LASSO), and cluster analysis to identify and define the hub genes in SCI. Finally, the link between hub genes and infiltrating immune cells was investigated and the potential transcriptional regulation/small molecular compounds to target hub genes were predicted. In total, SKP1 and BAP1 were identified as hub genes of mitophagy-related DEGs during SCI development and regulatory T cells (Tregs)/resting NK cells/activated mast cells may play an essential role in the progression of SCI. LINC00324 and SNHG16 may regulate SKP1 and BAP1, respectively, through miRNAs. Eleven and eight transcriptional factors (TFs) regulate SKP1 and BAP1, respectively, and six small molecular compounds target BAP1. Then, the mRNA expression levels of BAP1 and SKP1 were detected in the injured sites of spinal cord of SD rats at 6 h and 72 h after injury using RT-qPCR, and found that the level were decreased. Therefore, the pathways of mitophagy are downregulated during the pathophysiology of SCI, and SKP1 and BAP1 could be accessible targets for diagnosing and treating SCI.

Unraveling genetic threads: Identifying novel therapeutic targets for allergic rhinitis through Mendelian randomization

Background

Allergic rhinitis (AR) is a pervasive global health issue, and currently, there is a scarcity of targeted drug therapies available. This study aims to identify potential druggable target genes for AR using Mendelian randomization (MR) analysis.

Methods

MR analysis was conducted to assess the causal effect of expression quantitative trait loci (eQTL) in the blood on AR. Data on AR were collected from 2 datasets: FinnGen(R9) (11,009 cases and 359,149 controls) and UK Biobank (25,486 cases and 87,097 controls). Colocalization analysis was utilized to assess the common causal genetic variations between the identified drug target genes and AR. We also employed available genome-wide association studies (GWAS) data to gauge the impact of druggable genes on AR biomarkers and other allergic diseases.

Results

This study employs MR to analyze the relationship between 3410 druggable genes and AR. After Bonferroni correction, 10 genes were found to be significantly associated with AR risk (P < 0.05/3410). Colocalization analysis revealed a significant causal relationship between the expression variation of CFL1 and EFEMP2 genes and AR, sharing direct causal variants (colocalization probability PP.H3 + PP.H4 > 0.8), highlighting their importance as potential therapeutic targets for AR. The CFL1 gene showed a causal link with levels of thymic stromal lymphopoietin (TSLP), eosinophil count, and interleukin-13 (IL-13) (P = 0.016, 7.45E-16, 0.00091, respectively). EFEMP2 was also causally related to eosinophil count, IL-13, and interleukin-17 (IL-17) (P = 0.00012, 0.00091, 0.032, respectively). PheWAS analysis revealed significant associations of CFL1 with asthma, whereas EFEMP2 showed associations with both asthma and eczema. Protein-Protein Interaction (PPI) network analysis further unveiled the direct interactions of EFEMP2 and CFL1 with proteins related to immune regulation and inflammatory responses, with 77.64% of the network consisting of direct bindings, indicating their key roles in modulating AR-related immune and inflammatory responses. Notably, there was an 8.01% significant correlation between immune-related pathways and genes involved in inflammatory responses.

Conclusion

These genes present notable associations with AR biomarkers and other autoimmune diseases, offering valuable targets for developing new AR therapies.

Identification of the differentially expressed activated memory CD4+ T-cells-related genes and ceRNAs in oral lichen planus

Background

Oral lichen planus (OLP) is a common chronic oral mucosal disease with 1.4 % malignant transformation rate, and its etiology especially immune pathogenesis remains unclear. This study was aimed at investigating the immune cells related molecular underlying the pathophysiology of OLP through bioinformatics analysis.

Methods

The dataset GSE52130 obtained from the Gene Expression Omnibus (GEO) database was conducted a comprehensive analysis in this study. The CIBERSORTx was used for investigating immune cells infiltration. The gene set enrichment analysis (GSEA) and gene ontology (GO) enrichment were performed for exploring the biological functions and gene annotation. The protein-protein interactions (PPI) were constructed by STRING database and visualized by Cytoscape software. The cytohubba plugin was utilized for screening hub genes. The receiver operating characteristic (ROC) was performed for evaluating diagnostic value of hub genes. The miRNAs, lncRNAs and drugs were respectively predicted by NetworkAnalyst, miRTarbase, ENCORI, and DGIdb database.

Results

This study identified 595 differentially expressed genes (DEGs). The GSEA indicated keratinization, innate immune system and biological oxidation were involved in OLP. GO analysis showed extracellular matrix and keratinocyte were mainly enriched. And we found the activated memory CD4+ T cells were lowly infiltrated in OLP. We identified 101 activated memory CD4+ T-cells-related DEGs. Three hub genes (APP, IL1B, TF) were selected. APP and IL1B were significantly up-regulated, whereas TF was down-regulated in OLP. The three hub genes show high diagnostic value in OLP. Additionally, they were involved in MAPK signal, NF-kappaB signal and iron metabolism in OLP. What's more, NEAT1/XIST - miR - 15a - 5p/miR - 155-5p - APP/IL1B signal axis was focused in competing endogenous RNA (ceRNA) network. In addition, 35 drugs were predicted for OLP.

Conclusion

Three activated memory CD4+ T-cells-related DEGs were identified by integrative analysis. It may provide novel insight into the pathogenesis of OLP and suggest potential therapeutic targets for OLP.

Identification of diagnostic biomarkers and potential therapeutic targets for biliary atresia via WGCNA and machine learning methods

Biliary atresia (BA) is a severe and progressive biliary obstructive disease in infants that requires early diagnosis and new therapeutic targets. This study employed bioinformatics methods to identify diagnostic biomarkers and potential therapeutic targets for BA. Our analysis of mRNA expression from Gene Expression Omnibus datasets revealed 3,273 differentially expressed genes between patients with BA and those without BA (nBA). Weighted gene coexpression network analysis determined that the turquoise gene coexpression module, consisting of 298 genes, is predominantly associated with BA. The machine learning method then filtered out the top 2 important genes, CXCL8 and TMSB10, from the turquoise module. The area under receiver operating characteristic curves for TMSB10 and CXCL8 were 0.961 and 0.927 in the training group and 0.819 and 0.791 in the testing group, which indicated a high diagnostic value. Besides, combining TMSB10 and CXCL8, a nomogram with better diagnostic performance was built for clinical translation. Several studies have highlighted the potential of CXCL8 as a therapeutic target for BA, while TMSB10 has been shown to regulate cell polarity, which was related to BA progression. Our analysis with qRT PCR and immunohistochemistry also confirmed the upregulation of TMSB10 at mRNA and protein levels in BA liver samples. These findings highlight the sensitivity of CXCL8 and TMSB10 as diagnostic biomarkers and their potential as therapeutic targets for BA.

scRank infers drug-responsive cell types from untreated scRNA-seq data using a target-perturbed gene regulatory network

Cells respond divergently to drugs due to the heterogeneity among cell populations. Thus, it is crucial to identify drug-responsive cell populations in order to accurately elucidate the mechanism of drug action, which is still a great challenge. Here, we address this problem with scRank, which employs a target-perturbed gene regulatory network to rank drug-responsive cell populations via in silico drug perturbations using untreated single-cell transcriptomic data. We benchmark scRank on simulated and real datasets, which shows the superior performance of scRank over existing methods. When applied to medulloblastoma and major depressive disorder datasets, scRank identifies drug-responsive cell types that are consistent with the literature. Moreover, scRank accurately uncovers the macrophage subpopulation responsive to tanshinone IIA and its potential targets in myocardial infarction, with experimental validation. In conclusion, scRank enables the inference of drug-responsive cell types using untreated single-cell data, thus providing insights into the cellular-level impacts of therapeutic interventions.

Plasma proteome analysis implicates novel proteins as potential therapeutic targets for chronic kidney disease: A proteome-wide association study

Chronic kidney disease (CKD) is prevalent globally with limited therapeutic drugs available. To systemically identify novel proteins involved in the pathogenesis of CKD and possible therapeutic targets, we integrated human plasma proteomes with the genome-wide association studies (GWASs) of CKD, estimated glomerular filtration rate (eGFR) and blood urea nitrogen (BUN) to perform proteome-wide association study (PWAS), Mendelian Randomization and Bayesian colocalization analyses. The single-cell RNA sequencing data of healthy human and mouse kidneys were analyzed to explore the cell-type specificity of identified genes. Functional enrichment analysis was conducted to investigate the involved signaling pathways. The PWAS identified 22 plasma proteins significantly associated with CKD. Of them, the significant associations of three proteins (INHBC, LMAN2, and SNUPN) were replicated in the GWASs of eGFR, and BUN. Mendelian Randomization analyses showed that INHBC and SNUPN were causally associated with CKD, eGFR, and BUN. The Bayesian colocalization analysis identified shared causal variants for INHBC in CKD, eGFR, and BUN (all PP4 > 0.75). The single-cell RNA sequencing revealed that the INHBC gene was sparsely scattered within the kidney cells. This proteomic study revealed that INHBC, LMAN2, and SNUPN may be involved in the pathogenesis of CKD, which represent novel therapeutic targets and warrant further exploration in future research.

Application of angiogenesis-related genes associated with immune infiltration in the molecular typing and diagnosis of acute myocardial infarction

BACKGROUND

Angiogenesis has been discovered to be a critical factor in developing tumors and ischemic diseases. However, the role of angiogenesis-related genes (ARGs) in acute myocardial infarction (AMI) remains unclear.

METHODS

The GSE66360 dataset was used as the training cohort, and the GSE48060 dataset was used as the external validation cohort. The random forest (RF) algorithm was used to identify the signature genes. Consensus clustering analysis was used to identify robust molecular clusters associated with angiogenesis. The ssGSEA was used to analyze the correlation between ARGs and immune cell infiltration. In addition, we constructed miRNA-gene, transcription factor network, and targeted drug network of signature genes. RT-qPCR was used to verify the expression levels of signature genes.

RESULTS

Seven signature ARGs were identified based on the RF algorithm. Receiver operating characteristic curves confirmed the classification accuracy of the risk predictive model based on signature ARGs (area under the curve [AUC] = 0.9596 in the training cohort and AUC = 0.7773 in the external validation cohort). Subsequently, the ARG clusters were identified by consensus clustering. Cluster B had a more generalized high expression of ARGs and was significantly associated with immune infiltration. The miRNA and transcription factor network provided new ideas for finding potential upstream targets and biomarkers. Finally, the results of RT-qPCR were consistent with the bioinformatics analysis, further validating our results.

CONCLUSIONS

Angiogenesis is closely related to AMI, and characterizing the angiogenic features of patients with AMI can help to risk-stratify patients and provide personalized treatment.

Mining key circadian biomarkers for major depressive disorder by integrating bioinformatics and machine learning

OBJECTIVE

This study aimed to identify key clock genes closely associated with major depressive disorder (MDD) using bioinformatics and machine learning approaches.

METHODS

Gene expression data of 128 MDD patients and 64 healthy controls from blood samples were obtained. Differentially expressed were identified and weighted gene co-expression network analysis (WGCNA) was first performed to screen MDD-related key genes. These genes were then intersected with 1475 known circadian rhythm genes to identify circadian rhythm genes associated with MDD. Finally, multiple machine learning algorithms were applied for further selection, to determine the most critical 4 circadian rhythm biomarkers.

RESULTS

Four key circadian rhythm genes (ABCC2, APP, HK2 and RORA) were identified that could effectively distinguish MDD samples from controls. These genes were significantly enriched in circadian pathways and showed strong correlations with immune cell infiltration. Drug target prediction suggested that small molecules like melatonin and escitalopram may target these circadian rhythm proteins.

CONCLUSION

This study revealed discovered 4 key circadian rhythm genes closely associated with MDD, which may serve as diagnostic biomarkers and therapeutic targets. The findings highlight the important roles of circadian disruptions in the pathogenesis of MDD, providing new insights for precision diagnosis and targeted treatment of MDD.

Identifying therapeutic target genes for migraine by systematic druggable genome-wide Mendelian randomization

BACKGROUND

Currently, the treatment and prevention of migraine remain highly challenging. Mendelian randomization (MR) has been widely used to explore novel therapeutic targets. Therefore, we performed a systematic druggable genome-wide MR to explore the potential therapeutic targets for migraine.

METHODS

We obtained data on druggable genes and screened for genes within brain expression quantitative trait locis (eQTLs) and blood eQTLs, which were then subjected to two-sample MR analysis and colocalization analysis with migraine genome-wide association studies data to identify genes highly associated with migraine. In addition, phenome-wide research, enrichment analysis, protein network construction, drug prediction, and molecular docking were performed to provide valuable guidance for the development of more effective and targeted therapeutic drugs.

RESULTS

We identified 21 druggable genes significantly associated with migraine (BRPF3, CBFB, CDK4, CHD4, DDIT4, EP300, EPHA5, FGFRL1, FXN, HMGCR, HVCN1, KCNK5, MRGPRE, NLGN2, NR1D1, PLXNB1, TGFB1, TGFB3, THRA, TLN1 and TP53), two of which were significant in both blood and brain (HMGCR and TGFB3). The results of phenome-wide research showed that HMGCR was highly correlated with low-density lipoprotein, and TGFB3 was primarily associated with insulin-like growth factor 1 levels.

CONCLUSIONS

This study utilized MR and colocalization analysis to identify 21 potential drug targets for migraine, two of which were significant in both blood and brain. These findings provide promising leads for more effective migraine treatments, potentially reducing drug development costs.

Integrated Bioinformatics Analysis Revealed Immune Checkpoint Genes Relevant to Type 2 Diabetes

Objective

Chronic low-grade inflammation of the pancreatic islets is the characteristic of type 2 diabetes (T2D), and some of the immune checkpoints may play important roles in the pancreatic islet inflammation. Thus, we aim to explore the immune checkpoint genes (ICGs) associated with T2D, thereby revealing the role of ICGs in the pathogenesis of T2D based on bioinformatic analyses.

Methods

Differentially expressed genes (DEGs) and immune checkpoint genes (ICGs) of islets between T2D and control group were screened from datasets of the Gene Expression Omnibus (GEO). A risk model was built based on the coefficients of ICGs calculated by ridge regression. Functional enrichment analysis and immune cell infiltration estimation were conducted. Correlations between ICGs and hub genes, T2D-related disease genes, insulin secretion genes, and beta cell function-related genes were analyzed. Finally, we conducted RT-PCR to verify the expression of these ICGs.

Results

In total, pancreatic islets from 19 cases of T2D and 84 healthy subjects were included. We identified 458 DEGs. Six significantly upregulated ICGs (CD44, CD47, HAVCR2, SIRPA, TNFSF9, and VTCN1) in T2D were screened out. These ICGs were significantly correlated with several hub genes and T2D-related genes; furthermore, they were correlated with insulin secretion and β cell function-related genes. The analysis of immune infiltration showed that the concentrations of eosinophils, T cells CD4 naive, and T cells regulatory (Tregs) were significantly higher, but CD4 memory resting T cells and monocytes were lower in islets of T2D patients. The infiltrated immune cells in T2D pancreatic islet were associated with these six ICGs. Finally, the expression levels of four ICGs were confirmed by RT-PCR, and three ICGs were validated in another independent dataset.

Conclusion

In conclusion, the identified ICGs may play an important role in T2D. Identification of these differential genes may provide new clues for the diagnosis and treatment of T2D.

Elucidating sleep disorders: a comprehensive bioinformatics analysis of functional gene sets and hub genes

Background

Sleep disorders (SD) are known to have a profound impact on human health and quality of life although their exact pathogenic mechanisms remain poorly understood.

Methods

The study first accessed SD datasets from the GEO and identified DEGs. These DEGs were then subjected to gene set enrichment analysis. Several advanced techniques, including the RF, SVM-RFE, PPI networks, and LASSO methodologies, were utilized to identify hub genes closely associated with SD. Additionally, the ssGSEA approach was employed to analyze immune cell infiltration and functional gene set scores in SD. DEGs were also scrutinized in relation to miRNA, and the DGIdb database was used to explore potential pharmacological treatments for SD. Furthermore, in an SD murine model, the expression levels of these hub genes were confirmed through RT-qPCR and Western Blot analyses.

Results

The findings of the study indicate that DEGs are significantly enriched in functions and pathways related to immune cell activity, stress response, and neural system regulation. The analysis of immunoinfiltration demonstrated a marked elevation in the levels of Activated CD4+ T cells and CD8+ T cells in the SD cohort, accompanied by a notable rise in Central memory CD4 T cells, Central memory CD8 T cells, and Natural killer T cells. Using machine learning algorithms, the study also identified hub genes closely associated with SD, including IPO9, RAP2A, DDX17, MBNL2, PIK3AP1, and ZNF385A. Based on these genes, an SD diagnostic model was constructed and its efficacy validated across multiple datasets. In the SD murine model, the mRNA and protein expressions of these 6 hub genes were found to be consistent with the results of the bioinformatics analysis.

Conclusion

In conclusion, this study identified 6 genes closely linked to SD, which may play pivotal roles in neural system development, the immune microenvironment, and inflammatory responses. Additionally, the key gene-based SD diagnostic model constructed in this study, validated on multiple datasets showed a high degree of reliability and accuracy, predicting its wide potential for clinical applications. However, limited by the range of data sources and sample size, this may affect the generalizability of the results.

Mendelian randomization reveals interactions of the blood proteome and immunome in mitral valve prolapse

BACKGROUND

Mitral valve prolapse (MVP) is a common heart disorder characterized by an excessive production of proteoglycans and extracellular matrix in mitral valve leaflets. Large-scale genome-wide association study (GWAS) underlined that MVP is heritable. The molecular underpinnings of the disease remain largely unknown.

METHODS

We interrogated cross-modality data totaling more than 500,000 subjects including GWAS, 4809 molecules of the blood proteome, and genome-wide expression of mitral valves to identify candidate drivers of MVP. Data were investigated through Mendelian randomization, network analysis, ligand-receptor inference and digital cell quantification.

RESULTS

In this study, Mendelian randomization identify that 33 blood proteins, enriched in networks for immunity, are associated with the risk of MVP. MVP- associated blood proteins are enriched in ligands for which their cognate receptors are differentially expressed in mitral valve leaflets during MVP and enriched in cardiac endothelial cells and macrophages. MVP-associated blood proteins are involved in the renewal-polarization of macrophages and regulation of adaptive immune response. Cytokine activity profiling and digital cell quantification show in MVP a shift toward cytokine signature promoting M2 macrophage polarization. Assessment of druggability identify CSF1R, CX3CR1, CCR6, IL33, MMP8, ENPEP and angiotensin receptors as actionable targets in MVP.

CONCLUSIONS

Hence, integrative analysis identifies networks of candidate molecules and cells involved in immune control and remodeling of the extracellular matrix, which drive the risk of MVP.

A Diagnostic Model for Parkinson's Disease Based on Anoikis-Related Genes

Parkinson's disease (PD) is the second most prevalent neurodegenerative disease, and its pathological mechanisms are thought to be closely linked to apoptosis. Anoikis, a specific type of apoptosis, has recently been suggested to play a role in the progression of Parkinson's disease; however, the underlying mechanisms are not well understood. To explore the potential mechanisms involved in PD, we selected genes from the GSE28894 dataset and compared their expression in PD patients and healthy controls to identify differentially expressed genes (DEGs), and selected anoikis-related genes (ANRGs) from the DEGs. Furthermore, the least absolute shrinkage and selection operator (LASSO) regression approach and multivariate logistic regression highlighted five key genes-GSK3B, PCNA, CDC42, DAPK2, and SRC-as biomarker candidates. Subsequently, we developed a nomogram model incorporating these 5 genes along with age and sex to predict and diagnose PD. To evaluate the model's coherence, clinical applicability, and distinguishability, we utilized receiver operating characteristic (ROC) curves, the C-index, and calibration curves and validated it in both the GSE20295 dataset and our center's external clinical data. In addition, we confirmed the differential expression of the 5 model genes in human blood samples through qRT-PCR and Western blotting. Our constructed anoikis-related PD diagnostic model exhibits satisfactory predictive accuracy and offers novel insights into both diagnosis and treatment strategies for Parkinson's disease while facilitating its implementation in clinical practice.

A cloud-based precision oncology framework for whole genome sequence analysis

Cancer is one of the wide-ranging diseases which have a high mortality rate impacting globally. This scenario can be switched by early detection and correct precision treatment, a major concern for cancer patients. Clinicians can figure out the best-suited treatments for cancer patients by analyzing the patient's genome, which will treat the patient well and minimize the chances of side effects as well. Therefore, we have developed a fast, robust, and efficient solution as our precision oncology framework based on the whole genome sequencing of the individual's DNA. This platform can perform the entire genomic analysis, starting from the quality assessment of the input file to the variant annotation and functional prediction, followed by a certain level of interpretation. This analysis helps in the molecular profiling of the tumors for the identification of the targetable alterations. It takes in FASTQ or BAM file as an input and provides us with two output reports: a primary report, which consists of the patients' details, a summary of the analysis, and a secondary report, which is an elaborated report comprised of numerous results obtained from the analysis such as base changes, codon changes, amino acid changes, TMB analysis, MSI analysis, the variant frequency with its effects and impacts, affected biomarkers, etc. This framework can be effectively utilized for cancer treatment guidance, identification and validation of novel biomarkers, oncology research & development, genomic analysis, and gene manipulation.

Comprehensive Analysis of ceRNA Network and Immune Cell Infiltration Pattern of Autophagy-Related Genes in IgA Nephropathy

INTRODUCTION

IgA nephropathy (IgAN) is a prevalent worldwide glomerular disease with a complex pathophysiology that has significant economic implications. Despite the lack of successful research, this study aims to discover the potential competing endogenous RNA (ceRNA) network of autophagy-associated genes in IgAN and examine their correlation with immune cell infiltration.

METHODS

Autophagy-related hub genes were discovered by assessing the GSE116626 dataset and constructing a protein-protein interaction network. Nephroseq v5 analysis engine was used to analyze correlations between hub genes and proteinuria, glomerular filtration rate (GFR), and serum creatinine levels. Then, a ceRNA network construction and the CIBERSORT tool for immune cell infiltration analysis were also performed. Additionally, the differentially expressed autophagy-related genes were used to predict potential targeted medications for IgAN.

RESULTS

Overall, 1,396 differentially expressed genes were identified in IgAN along with 25 autophagy-related differentially expressed messenger RNAs. Enrichment analysis revealed significant involvement of autophagy and apoptosis in biological processes. Next, we evaluated the top hub nodes based on their highest degrees. The ability of IgAN discrimination was confirmed in the GSE35487 and GSE37460 datasets by validating the five hub genes: SIRT1, FOS, CCL2, CDKN1A, and MYC. In the Nephroseq v5 analysis engine, the clinical correlation of the five hub genes was confirmed. Furthermore, the ceRNA network identified 18 circular RNAs and 2 microRNAs associated with hub autophagy-related genes in IgAN. Our investigation identified hsa-miR-32-3p and hsa-let-7i-5p as having elevated expression levels and substantial diagnostic value. Finally, four distinctively infiltrated immune cells were found to be associated with the hub autophagy-related genes, and 67 drugs were identified as potential therapeutic options for IgAN.

CONCLUSION

This study sheds light on a novel ceRNA regulatory network mechanism associated with autophagy in IgAN development.

Identification of pleiotropic and specific therapeutic targets for cardio-cerebral diseases: A large-scale proteome-wide mendelian randomization and colocalization study

BACKGROUND

The cardiac-brain connection has been identified as the basis for multiple cardio-cerebral diseases. However, effective therapeutic targets for these diseases are still limited. Therefore, this study aimed to identify pleiotropic and specific therapeutic targets for cardio-cerebral diseases using Mendelian randomization (MR) and colocalization analyses.

METHODS

This study included two large protein quantitative trait loci studies with over 4,000 plasma proteins were included in the discovery and replication cohorts, respectively. We initially used MR to estimate the associations between protein and 20 cardio-cerebral diseases. Subsequently, Colocalization analysis was employed to enhance the credibility of the results. Protein target prioritization was based solely on including highly robust significant results from both the discovery and replication phases. Lastly, the Drug-Gene Interaction Database was utilized to investigate protein-gene-drug interactions further.

RESULTS

A total of 46 target proteins for cardio-cerebral diseases were identified as robust in the discovery and replication phases by MR, comprising 7 pleiotropic therapeutic proteins and 39 specific target proteins. Followed by colocalization analysis and prioritization of evidence grades for target protein, 6 of these protein-disease pairs have achieved the highly recommended level. For instance, the PILRA protein presents a pleiotropic effect on sick sinus syndrome and Alzheimer's disease, whereas GRN exerts specific effects on the latter. APOL3, LRP4, and F11, on the other hand, have specific effects on cardiomyopathy and ischemic stroke, respectively.

CONCLUSIONS

This study successfully identified important therapeutic targets for cardio-cerebral diseases, which benefits the development of preventive or therapeutic drugs.

Exploring cutting-edge strategies for drug repurposing in female cancers - An insight into the tools of the trade

Female cancers, which include breast and gynaecological cancers, represent a significant global health burden for women. Despite advancements in research pertinent to unearthing crucial pathological characteristics of these cancers, challenges persist in discovering potential therapeutic strategies. This is further exacerbated by economic burdens associated with de novo drug discovery and clinical intricacies such as development of drug resistance and metastasis. Drug repurposing, an innovative approach leveraging existing FDA-approved drugs for new indications, presents a promising avenue to expedite therapeutic development. Computational techniques, including virtual screening and analysis of drug-target-disease relationships, enable the identification of potential candidate drugs. Integration of diverse data types, such as omics and clinical information, enhances the precision and efficacy of drug repurposing strategies. Experimental approaches, including high-throughput screening assays, in vitro, and in vivo models, complement computational methods, facilitating the validation of repurposed drugs. This review highlights various target mining strategies based on analysis of differential gene expression, weighted gene co-expression, protein-protein interaction network, and host-pathogen interaction, among others. To unearth drug candidates, the technicalities of leveraging information from databases such as DrugBank, STITCH, LINCS, and ChEMBL, among others are discussed. Further in silico validation techniques encompassing molecular docking, pharmacophore modelling, molecular dynamic simulations, and ADMET analysis are elaborated. Overall, this review delves into the exploration of individual case studies to offer a wide perspective of the ever-evolving field of drug repurposing, emphasizing the multifaceted approaches and methodologies employed for the same to confront female cancers.

Genome-wide association study reveals a locus in ADARB2 for complete freedom from headache in Danish Blood Donors

Headache disorders are the most common disorders of the nervous system. The lifetime prevalence of headache disorders show that some individuals never experience headache. The etiology of complete freedom from headache is not known. To assess genetic variants associated with complete freedom from headache, we performed a genome-wide association study of individuals who have never experienced a headache. We included 63,992 individuals (2,998 individuals with complete freedom from headache and 60,994 controls) from the Danish Blood Donor Study Genomic Cohort. Participants were included in two rounds, from 2015 to 2018 and in 2020. We discovered a genome-wide significant association, with the lead variant rs7904615[G] in ADARB2 (EAF = 27%, OR = 1.20 [1.13-1.27], p = 3.92 × 10-9). The genomic locus was replicated in a non-overlapping cohort of 13,032 individuals (539 individuals with complete freedom from headache and 12,493 controls) from the Danish Blood Donor Study Genomic Cohort (p < 0.05, two-sided). Participants for the replication were included from 2015 to 2020. In conclusion, we show that complete freedom from headache has a genetic component, and we suggest that ADARB2 is involved in complete freedom from headache. The genomic locus was specific for complete freedom from headache and was not associated with any primary headache disorders.

Genetic associations in ankylosing spondylitis: circulating proteins as drug targets and biomarkers

Background

Ankylosing spondylitis (AS) is a complex condition with a significant genetic component. This study explored circulating proteins as potential genetic drug targets or biomarkers to prevent AS, addressing the need for innovative and safe treatments.

Methods

We analyzed extensive data from protein quantitative trait loci (pQTLs) with up to 1,949 instrumental variables (IVs) and selected the top single-nucleotide polymorphism (SNP) associated with AS risk. Utilizing a two-sample Mendelian randomization (MR) approach, we assessed the causal relationships between identified proteins and AS risk. Colocalization analysis, functional enrichment, and construction of protein-protein interaction networks further supported these findings. We utilized phenome-wide MR (phenMR) analysis for broader validation and repurposing of drugs targeting these proteins. The Drug-Gene Interaction database (DGIdb) was employed to corroborate drug associations with potential therapeutic targets. Additionally, molecular docking (MD) techniques were applied to evaluate the interaction between target protein and four potential AS drugs identified from the DGIdb.

Results

Our analysis identified 1,654 plasma proteins linked to AS, with 868 up-regulated and 786 down-regulated. 18 proteins (AGER, AIF1, ATF6B, C4A, CFB, CLIC1, COL11A2, ERAP1, HLA-DQA2, HSPA1L, IL23R, LILRB3, MAPK14, MICA, MICB, MPIG6B, TNXB, and VARS1) that show promise as therapeutic targets for AS or biomarkers, especially MAPK14, supported by evidence of colocalization. PhenMR analysis linked these proteins to AS and other diseases, while DGIdb analysis identified potential drugs related to MAPK14. MD analysis indicated strong binding affinities between MAPK14 and four potential AS drugs, suggesting effective target-drug interactions.

Conclusion

This study underscores the utility of MR analysis in AS research for identifying biomarkers and therapeutic drug targets. The involvement of Th17 cell differentiation-related proteins in AS pathogenesis is particularly notable. Clinical validation and further investigation are essential for future applications.

Identification and verification of a novel signature that combines cuproptosis-related genes with ferroptosis-related genes in osteoarthritis using bioinformatics analysis and experimental validation

BACKGROUND

Exploring the pathogenesis of osteoarthritis (OA) is important for its prevention, diagnosis, and treatment. Therefore, we aimed to construct novel signature genes (c-FRGs) combining cuproptosis-related genes (CRGs) with ferroptosis-related genes (FRGs) to explore the pathogenesis of OA and aid in its treatment.

MATERIALS AND METHODS

Differentially expressed c-FRGs (c-FDEGs) were obtained using R software. Enrichment analysis was performed and a protein-protein interaction (PPI) network was constructed based on these c-FDEGs. Then, seven hub genes were screened. Three machine learning methods and verification experiments were used to identify four signature biomarkers from c-FDEGs, after which gene set enrichment analysis, gene set variation analysis, single-sample gene set enrichment analysis, immune function analysis, drug prediction, and ceRNA network analysis were performed based on these signature biomarkers. Subsequently, a disease model of OA was constructed using these biomarkers and validated on the GSE82107 dataset. Finally, we analyzed the distribution of the expression of these c-FDEGs in various cell populations.

RESULTS

A total of 63 FRGs were found to be closely associated with 11 CRGs, and 40 c-FDEGs were identified. Bioenrichment analysis showed that they were mainly associated with inflammation, external cellular stimulation, and autophagy. CDKN1A, FZD7, GABARAPL2, and SLC39A14 were identified as OA signature biomarkers, and their corresponding miRNAs and lncRNAs were predicted. Finally, scRNA-seq data analysis showed that the differentially expressed c-FRGs had significantly different expression distributions across the cell populations.

CONCLUSION

Four genes, namely CDKN1A, FZD7, GABARAPL2, and SLC39A14, are excellent biomarkers and prospective therapeutic targets for OA.

IL-2RG as a possible immunotherapeutic target in CRC predicting poor prognosis and regulated by miR-7-5p and miR-26b-5p

Despite advances in treatment strategies, colorectal cancer (CRC) continues to cause significant morbidity and mortality, with mounting evidence a close link between immune system dysfunctions issued. Interleukin-2 receptor gamma (IL-2RG) plays a pivotal role as a common subunit receptor in the IL-2 family cytokines and activates the JAK-STAT pathway. This study delves into the role of Interleukin-2 receptor gamma (IL-2RG) within the tumor microenvironment and investigates potential microRNAs (miRNAs) that directly inhibit IL-2RG, aiming to discern their impact on CRC clinical outcomes. Bioinformatics analysis revealed a significant upregulation of IL-2RG mRNA in TCGA-COAD samples and showed strong correlations with the infiltration of various lymphocytes. Single-cell analysis corroborated these findings, highlighting IL-2RG expression in critical immune cell subsets. To explore miRNA involvement in IL-2RG dysregulation, mRNA was isolated from the tumor tissues and lymphocytes of 258 CRC patients and 30 healthy controls, and IL-2RG was cloned into the pcDNA3.1/CT-GFP-TOPO vector. Human embryonic kidney cell lines (HEK-293T) were transfected with this construct. Our research involved a comprehensive analysis of miRPathDB, miRWalk, and Targetscan databases to identify the miRNAs associated with the 3' UTR of human IL-2RG. The human microRNA (miRNA) molecules, hsa-miR-7-5p and hsa-miR-26b-5p, have been identified as potent suppressors of IL-2RG expression in CRC patients. Specifically, the downregulation of hsa-miR-7-5p and hsa-miR-26b-5p has been shown to result in the upregulation of IL-2RG mRNA expression in these patients. Prognostic evaluation of IL-2RG, hsa-miR-7-5p, and hsa-miR-26b-5p, using TCGA-COAD data and patient samples, established that higher IL-2RG expression and lower expression of both miRNAs were associated with poorer outcomes. Additionally, this study identified several long non-coding RNAs (LncRNAs), such as ZFAS1, SOX21-AS1, SNHG11, SNHG16, SNHG1, DLX6-AS1, GAS5, SNHG6, and MALAT1, which may act as competing endogenous RNA molecules for IL2RG by sequestering shared hsa-miR-7-5p and hsa-miR-26b-5p. In summary, this investigation underscores the potential utility of IL-2RG, hsa-miR-7-5p, and hsa-miR-26b-5p as serum and tissue biomarkers for predicting CRC patient prognosis while also offering promise as targets for immunotherapy in CRC management.

Circadian transcriptome of pancreatic adenocarcinoma unravels chronotherapeutic targets

Pancreatic ductal adenocarcinoma (PDA) is a lethal cancer characterized by a poor outcome and an increasing incidence. A significant majority (>80%) of newly diagnosed cases are deemed unresectable, leaving chemotherapy as the sole viable option, though with only moderate success. This necessitates the identification of improved therapeutic options for PDA. We hypothesized that there are temporal variations in cancer-relevant processes within PDA tumors, offering insights into the optimal timing of drug administration - a concept termed chronotherapy. In this study, we explored the presence of the circadian transcriptome in PDA using patient-derived organoids and validated these findings by comparing PDA data from The Cancer Genome Atlas with noncancerous healthy pancreas data from GTEx. Several PDA-associated pathways (cell cycle, stress response, Rho GTPase signaling) and cancer driver hub genes (EGFR and JUN) exhibited a cancer-specific rhythmic pattern intricately linked to the circadian clock. Through the integration of multiple functional measurements for rhythmic cancer driver genes, we identified top chronotherapy targets and validated key findings in molecularly divergent pancreatic cancer cell lines. Testing the chemotherapeutic efficacy of clinically relevant drugs further revealed temporal variations that correlated with drug-target cycling. Collectively, our study unravels the PDA circadian transcriptome and highlights a potential approach for optimizing chrono-chemotherapeutic efficacy.

Integrated analysis of public datasets for the discovery and validation of survival-associated genes in solid tumors

Identifying genes with prognostic significance that can act as biomarkers in solid tumors can help stratify patients and uncover novel therapy targets. Here, our goal was to expand our previous ranking analysis of survival-associated genes in various solid tumors to include colon cancer specimens with available transcriptomic and clinical data. A Gene Expression Omnibus search was performed to identify available datasets with clinical data and raw gene expression measurements. A combined database was set up and integrated into our Kaplan-Meier plotter, making it possible to identify genes with expression changes linked to altered survival. As a demonstration of the utility of the platform, the most powerful genes linked to overall survival in colon cancer were identified using uni- and multivariate Cox regression analysis. The combined colon cancer database includes 2,137 tumor samples from 17 independent cohorts. The most significant genes associated with relapse-free survival with a false discovery rate below 1% in colon cancer carcinoma were RBPMS (hazard rate [HR] = 2.52), TIMP1 (HR = 2.44), and COL4A2 (HR = 2.36). The three strongest genes associated with shorter survival in stage II colon cancer include CSF1R (HR = 2.86), FLNA (HR = 2.88), and TPBG (HR = 2.65). In summary, a new integrated database for colon cancer is presented. A colon cancer analysis subsystem was integrated into our Kaplan-Meier plotter that can be used to mine the entire database (https://www.kmplot.com). The portal has the potential to be employed for the identification and prioritization of promising biomarkers and therapeutic target candidates in multiple solid tumors including, among others, breast, lung, ovarian, gastric, pancreatic, and colon cancers.

DIRMC: a database of immunotherapy-related molecular characteristics

Cancer immunotherapy has brought about a revolutionary breakthrough in the field of cancer treatment. Immunotherapy has changed the treatment landscape for a variety of solid and hematologic malignancies. To assist researchers in efficiently uncovering valuable information related to cancer immunotherapy, we have presented a manually curated comprehensive database called DIRMC, which focuses on molecular features involved in cancer immunotherapy. All the content was collected manually from published literature, authoritative clinical trial data submitted by clinicians, some databases for drug target prediction such as DrugBank, and some experimentally confirmed high-throughput data sets for the characterization of immune-related molecular interactions in cancer, such as a curated database of T-cell receptor sequences with known antigen specificity (VDJdb), a pathology-associated TCR database (McPAS-TCR) et al. By constructing a fully connected functional network, ranging from cancer-related gene mutations to target genes to translated target proteins to protein regions or sites that may specifically affect protein function, we aim to comprehensively characterize molecular features related to cancer immunotherapy. We have developed the scoring criteria to assess the reliability of each MHC-peptide-T-cell receptor (TCR) interaction item to provide a reference for users. The database provides a user-friendly interface to browse and retrieve data by genes, target proteins, diseases and more. DIRMC also provides a download and submission page for researchers to access data of interest for further investigation or submit new interactions related to cancer immunotherapy targets. Furthermore, DIRMC provides a graphical interface to help users predict the binding affinity between their own peptide of interest and MHC or TCR. This database will provide researchers with a one-stop resource to understand cancer immunotherapy-related targets as well as data on MHC-peptide-TCR interactions. It aims to offer reliable molecular characteristics support for both the analysis of the current status of cancer immunotherapy and the development of new immunotherapy. DIRMC is available at http://www.dirmc.tech/. Database URL: http://www.dirmc.tech/.

EPHB2 as a key mediator of glioma progression: Insights from microenvironmental receptor ligand-related prognostic gene signature

Malignant gliomas, characterized by pronounced heterogeneity, a complex microenvironment, and a propensity for relapse and drug resistaniguree, pose significant challenges in oncology. This study aimed to investigate the prognostic value of Ligand and Receptor related genes (LRRGs) within the glioma microenvironment. An intersection of 71 ligand-related genes (LRGs) and 2628 receptor-related genes (RRGs) yielded a total of 69 LRRGs. Utilizing the least absolute shrinkage and selection operator (LASSO) regression analysis, a prognostic RiskScore model comprising 28 LRRGs was constructed. The model demonstrated robust prognostic value, further validated in the TCGA-GBMLGG dataset. Subsequent analyses included differential gene expression, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), gene set enrichment (GSEA), and gene set variation (GSVA) within RiskScore groups. Additionally, evaluations of PPI, mRNA-RBP, mRNA-TF, and mRNA-drug interaction networks were conducted. Four hub genes were identified through differential expression analysis of the 28 LRRGs across various GSE datasets. A multivariate Cox prognostic model was constructed for nomogram analysis, gene mutation analysis, and related expression distribution. This study underscores the role of LRRGs in intercellular communication within the glioma microenvironment and identifies four hub genes crucial for prognostic assessment in clinical glioma patients. These findings offer a potential evaluation framework for glioma patients, enhancing our understanding of the disease and informing future therapeutic strategies.

Exploiting meta-analysis of genome-wide interaction with serum 25-hydroxyvitamin D to identify novel genetic loci associated with pulmonary function

BACKGROUND

Higher 25-hydroxyvitamin D (25(OH)D) concentrations in serum has a positive association with pulmonary function. Investigating genome-wide interactions with 25(OH)D may reveal new biological insights into pulmonary function.

OBJECTIVES

We aimed to identify novel genetic variants associated with pulmonary function by accounting for 25(OH)D interactions.

METHODS

We included 211,264 participants from the observational United Kingdom Biobank study with pulmonary function tests (PFTs), genome-wide genotypes, and 25(OH)D concentrations from 4 ancestral backgrounds-European, African, East Asian, and South Asian. Among PFTs, we focused on forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) because both were previously associated with 25(OH)D. We performed genome-wide association study (GWAS) analyses that accounted for variant×25(OH)D interaction using the joint 2 degree-of-freedom (2df) method, stratified by participants' smoking history and ancestry, and meta-analyzed results. We evaluated interaction effects to determine how variant-PFT associations were modified by 25(OH)D concentrations and conducted pathway enrichment analysis to examine the biological relevance of our findings.

RESULTS

Our GWAS meta-analyses, accounting for interaction with 25(OH)D, revealed 30 genetic variants significantly associated with FEV1 or FVC (P2df <5.00×10-8) that were not previously reported for PFT-related traits. These novel variant signals were enriched in lung function-relevant pathways, including the p38 MAPK pathway. Among variants with genome-wide-significant 2df results, smoking-stratified meta-analyses identified 5 variants with 25(OH)D interactions that influenced FEV1 in both smoking groups (never smokers P1df interaction<2.65×10-4; ever smokers P1df interaction<1.71×10-5); rs3130553, rs2894186, rs79277477, and rs3130929 associations were only evident in never smokers, and the rs4678408 association was only found in ever smokers.

CONCLUSION

Genetic variant associations with lung function can be modified by 25(OH)D, and smoking history can further modify variant×25(OH)D interactions. These results expand the known genetic architecture of pulmonary function and add evidence that gene-environment interactions, including with 25(OH)D and smoking, influence lung function.

Genome-wide association analyses identify 95 risk loci and provide insights into the neurobiology of post-traumatic stress disorder

Post-traumatic stress disorder (PTSD) genetics are characterized by lower discoverability than most other psychiatric disorders. The contribution to biological understanding from previous genetic studies has thus been limited. We performed a multi-ancestry meta-analysis of genome-wide association studies across 1,222,882 individuals of European ancestry (137,136 cases) and 58,051 admixed individuals with African and Native American ancestry (13,624 cases). We identified 95 genome-wide significant loci (80 new). Convergent multi-omic approaches identified 43 potential causal genes, broadly classified as neurotransmitter and ion channel synaptic modulators (for example, GRIA1, GRM8 and CACNA1E), developmental, axon guidance and transcription factors (for example, FOXP2, EFNA5 and DCC), synaptic structure and function genes (for example, PCLO, NCAM1 and PDE4B) and endocrine or immune regulators (for example, ESR1, TRAF3 and TANK). Additional top genes influence stress, immune, fear and threat-related processes, previously hypothesized to underlie PTSD neurobiology. These findings strengthen our understanding of neurobiological systems relevant to PTSD pathophysiology, while also opening new areas for investigation.

DRSPRING: Graph convolutional network (GCN)-Based drug synergy prediction utilizing drug-induced gene expression profile

Great efforts have been made over the years to identify novel drug pairs with synergistic effects. Although numerous computational approaches have been proposed to analyze diverse types of biological big data, the pharmacogenomic profiles, presumably the most direct proxy of drug effects, have been rarely used due to the data sparsity problem. In this study, we developed a composite deep-learning-based model that predicts the drug synergy effect utilizing pharmacogenomic profiles as well as molecular properties. Graph convolutional network (GCN) was used to represent and integrate the chemical structure, genetic interactions, drug-target information, and gene expression profiles of cell lines. Insufficient amount of pharmacogenomic data, i.e., drug-induced expression profiles from the LINCS project, was resolved by augmenting the data with the predicted profiles. Our method learned and predicted the Loewe synergy score in the DrugComb database and achieved a better or comparable performance compared to other published methods in a benchmark test. We also investigated contribution of various input features, which highlighted the value of basal gene expression and pharmacogenomic profiles of each cell line. Importantly, DRSPRING (DRug Synergy PRediction by INtegrated GCN) can be applied to any drug pairs and any cell lines, greatly expanding its applicability compared to previous methods.

Effect of smoking cessation medications on intracranial aneurysm risk: A Mendelian randomization study

INTRODUCTION

We aim to assess the association between smoking behavior and intracranial aneurysms (IAs) and the effect of smoking cessation medications on IAs at the genetic level.

METHODS

Causal effects of four phenotypes: 1) age at initiation of regular smoking, 2) cigarettes smoked per day, 3) smoking cessation, and 4) smoking initiation on IAs, were analyzed using two-sample inverse-variance weighted Mendelian randomization analyses. The effects of genes interacting with the smoking cessation medications were analyzed using cis-expression quantitative trait loci genetic instruments on IAs using summary statistics-based Mendelian randomization analyses. Colocalization analyses were then used to test whether the genes shared causal variants with IAs. The role of confounding phenotypes as potential causative mechanisms of IAs at these gene loci was tested.

RESULTS

Cigarettes smoked per day (OR=2.89; 95% CI:1.85-4.51) and smoking initiation on IAs (OR=4.64; 95% CI: 2.64-8.15) were significantly associated with IA risk. However, age at initiation of regular smoking (OR=0.54; 95% CI: 0.10-2.8) and smoking cessation (OR=6.80; 95% CI: 0.01-4812) had no overall effect on IAs. Of 88 genes that interacted with smoking cessation medications, two had a causal effect on IA risk. Genetic variants affecting HYKK levels showed strong evidence of colocalization with IA risk. Higher HYKK levels in the blood were associated with a lower IA risk. Gene target analyses revealed that cigarettes/day could be a main mediator of HYKK's effect on IA risk.

CONCLUSIONS

This study provides evidence supporting that smoking initiation on IAs and cigarettes/day may increase IA risk. Increased HYKK gene expression may reduce IA risk. This can be explained by the increased number of cigarettes consumed daily. HYKK could also reduce IA risk due to the positive effect of continuous abstinence and varenicline therapy on smoking cessation.

Molecular subtype identification of cerebral ischemic stroke based on ferroptosis-related genes

Cerebral ischemic stroke (CIS) has the characteristics of a high incidence, disability, and mortality rate. Here, we aimed to explore the potential pathogenic mechanisms of ferroptosis-related genes (FRGs) in CIS. Three microarray datasets from the Gene Expression Omnibus (GEO) database were utilized to analyze differentially expressed genes (DEGs) between CIS and normal controls. FRGs were obtained from a literature report and the FerrDb database. Weighted gene co-expression network analysis (WGCNA) and protein-protein interaction (PPI) network were used to screen hub genes. The receiver operating characteristic (ROC) curve was adopted to evaluate the diagnostic value of key genes in CIS, followed by analysis of immune microenvironment, transcription factor (TF) regulatory network, drug prediction, and molecular docking. In total, 128 CIS samples were divided into 2 subgroups after clustering analysis. Compared with cluster A, 1560 DEGs were identified in cluster B. After the construction of the WGCNA and PPI network, 5 hub genes, including MAPK3, WAS, DNAJC5, PRKCD, and GRB2, were identified for CIS. Interestingly, MAPK3 was a FRG that differentially expressed between cluster A and cluster B. The expression levels of 5 hub genes were all specifically highly in cluster A subtype. It is noted that neutrophils were the most positively correlated with all 5 real hub genes. PRKCD was one of the target genes of FASUDIL. In conclusion, five real hub genes were identified as potential diagnostic markers, which can distinguish the two subtypes well.

Unveiling the link between lactate metabolism and rheumatoid arthritis through integration of bioinformatics and machine learning

Rheumatoid arthritis (RA) is a persistent autoimmune condition characterized by synovitis and joint damage. Recent findings suggest a potential link to abnormal lactate metabolism. This study aims to identify lactate metabolism-related genes (LMRGs) in RA and investigate their correlation with the molecular mechanisms of RA immunity. Data on the gene expression profiles of RA synovial tissue samples were acquired from the gene expression omnibus (GEO) database. The RA database was acquired by obtaining the common LMRDEGs, and selecting the gene collection through an SVM model. Conducting the functional enrichment analysis, followed by immuno-infiltration analysis and protein-protein interaction networks. The results revealed that as possible markers associated with lactate metabolism in RA, KCNN4 and SLC25A4 may be involved in regulating macrophage function in the immune response to RA, whereas GATA2 is involved in the immune mechanism of DC cells. In conclusion, this study utilized bioinformatics analysis and machine learning to identify biomarkers associated with lactate metabolism in RA and examined their relationship with immune cell infiltration. These findings offer novel perspectives on potential diagnostic and therapeutic targets for RA.

HPV DNA Integration at Actionable Cancer-Related Genes Loci in HPV-Associated Carcinomas

In HPV-associated carcinomas, some examples of cancer-related genes altered by viral insertion and corresponding to potential therapeutic targets have been described, but no quantitative assessment of these events, including poorly recurrent targets, has been reported to date. To document these occurrences, we built and analyzed a database comprised of 1455 cases, including HPV genotypes and tumor localizations. Host DNA sequences targeted by viral integration were classified as "non-recurrent" (one single reported case; 838 loci), "weakly recurrent" (two reported cases; 82 loci), and highly recurrent (≥3 cases; 43 loci). Whereas the overall rate of cancer-related target genes was 3.3% in the Gencode database, this rate increased to 6.5% in "non-recurrent", 11.4% in "weakly recurrent", and 40.1% in "highly recurrent" genes targeted by integration (p = 4.9 × 10-4). This rate was also significantly higher in tumors associated with high-risk HPV16/18/45 than other genotypes. Among the genes targeted by HPV insertion, 30.2% corresponded to direct or indirect druggable targets, a rate rising to 50% in "highly recurrent" targets. Using data from the literature and the DepMap 23Q4 release database, we found that genes targeted by viral insertion could be new candidates potentially involved in HPV-associated oncogenesis. A more systematic characterization of HPV/host fusion DNA sequences in HPV-associated cancers should provide a better knowledge of HPV-driven carcinogenesis and favor the development of personalize patient treatments.

Sparassis latifolia and exercise training as complementary medicine mitigated the 5-fluorouracil potent side effects in mice with colorectal cancer: bioinformatics approaches, novel monitoring pathological metrics, screening signatures, and innovative management tactic

BACKGROUND

Prompt identification and assessment of the disease are essential for reducing the death rate associated with colorectal cancer (COL). Identifying specific causal or sensitive components, such as coding RNA (cRNA) and non-coding RNAs (ncRNAs), may greatly aid in the early detection of colorectal cancer.

METHODS

For this purpose, we gave natural chemicals obtained from Sparassis latifolia (SLPs) either alone or in conjunction with chemotherapy (5-Fluorouracil to a mouse colorectal tumor model induced by AOM-DSS. The transcription profile of non-coding RNAs (ncRNAs) and their target hub genes was evaluated using qPCR Real-Time, and ELISA techniques.

RESULTS

MSX2, MMP7, ITIH4, and COL1A2 were identified as factors in inflammation and oxidative stress, leading to the development of COL. The hub genes listed, upstream regulatory factors such as lncRNA PVT1, NEAT1, KCNQ1OT1, SNHG16, and miR-132-3p have been discovered as biomarkers for prognosis and diagnosis of COL. The SLPs and exercise, effectively decreased the size and quantity of tumors.

CONCLUSIONS

This effect may be attributed to the modulation of gene expression levels, including MSX2, MMP7, ITIH4, COL1A2, PVT1, NEAT1, KCNQ1OT1, SNHG16, and miR-132-3p. Ultimately, SLPs and exercise have the capacity to be regarded as complementing and enhancing chemotherapy treatments, owing to their efficacious components.

Disulfidptosis-related signature elucidates the prognostic, immunologic, and therapeutic characteristics in ovarian cancer

Introduction

Ovarian cancer (OC) is the deadliest malignancy in gynecology, but the mechanism of its initiation and progression is poorly elucidated. Disulfidptosis is a novel discovered type of regulatory cell death. This study aimed to develop a novel disulfidptosis-related prognostic signature (DRPS) for OC and explore the effects and potential treatment by disulfidptosis-related risk stratification.

Methods

The disulfidptosis-related genes were first analyzed in bulk RNA-Seq and a prognostic nomogram was developed and validated by LASSO algorithm and multivariate cox regression. Then we systematically assessed the clinicopathological and mutational characteristics, pathway enrichment analysis, immune cell infiltration, single-cell-level expression, and drug sensitivity according to DRPS.

Results

The DRPS was established with 6 genes (MYL6, PDLIM1, ACTN4, FLNB, SLC7A11, and CD2AP) and the corresponding prognostic nomogram was constructed based on the DRPS, FIGO stage, grade, and residual disease. Stratified by the risk score derived from DRPS, patients in high-risk group tended to have worse prognosis, lower level of disulfidptosis, activated oncogenic pathways, inhibitory tumor immune microenvironment, and higher sensitivity to specific drugs including epirubicin, stauroporine, navitoclax, and tamoxifen. Single-cell transcriptomic analysis revealed the expression level of genes in the DRPS significantly varied in different cell types between tumor and normal tissues. The protein-level expression of genes in the DRPS was validated by the immunohistochemical staining analysis.

Conclusion

In this study, the DRPS and corresponding prognostic nomogram for OC were developed, which was important for OC prognostic assessment, tumor microenvironment modification, drug sensitivity prediction, and exploration of potential mechanisms in tumor development.

New treatment alternatives for primary and metastatic colorectal cancer by an integrated transcriptome and network analyses

Metastatic colorectal cancer (CRC) is still in need of effective treatments. This study applies a holistic approach to propose new targets for treatment of primary and liver metastatic CRC and investigates their therapeutic potential in-vitro. An integrative analysis of primary and metastatic CRC samples was implemented for alternative target and treatment proposals. Integrated microarray samples were grouped based on a co-expression network analysis. Significant gene modules correlated with primary CRC and metastatic phenotypes were identified. Network clustering and pathway enrichments were applied to gene modules to prioritize potential targets, which were shortlisted by independent validation. Finally, drug-target interaction search led to three agents for primary and liver metastatic CRC phenotypes. Hesperadin and BAY-1217389 suppress colony formation over a 14-day period, with Hesperadin showing additional efficacy in reducing cell viability within 48 h. As both candidates target the G2/M phase proteins NEK2 or TTK, we confirmed their anti-proliferative properties by Ki-67 staining. Hesperadinin particular arrested the cell cycle at the G2/M phase. IL-29A treatment reduced migration and invasion capacities of TGF-β induced metastatic cell lines. In addition, this anti-metastatic treatment attenuated TGF-β dependent mesenchymal transition. Network analysis suggests IL-29A induces the JAK/STAT pathway in a preventive manner.

Proteomics study the potential targets for Rifampicin-resistant spinal tuberculosis

Introduction: The escalating global surge in Rifampicin-resistant strains poses a formidable challenge to the worldwide campaign against tuberculosis (TB), particularly in developing countries. The frequent reports of suboptimal treatment outcomes, complications, and the absence of definitive treatment guidelines for Rifampicin-resistant spinal TB (DSTB) contribute significantly to the obstacles in its effective management. Consequently, there is an urgent need for innovative and efficacious drugs to address Rifampicin-resistant spinal tuberculosis, minimizing the duration of therapy sessions. This study aims to investigate potential targets for DSTB through comprehensive proteomic and pharmaco-transcriptomic analyses. Methods: Mass spectrometry-based proteomics analysis was employed to validate potential DSTB-related targets. PPI analysis confirmed by Immunohistochemistry (IHC) and Western blot analysis. Results: The proteomics analysis revealed 373 differentially expressed proteins (DEPs), with 137 upregulated and 236 downregulated proteins. Subsequent Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses delved into the DSTB-related pathways associated with these DEPs. In the context of network pharmacology analysis, five key targets-human leukocyte antigen A chain (HLAA), human leukocyte antigen C chain (HLA-C), HLA Class II Histocompatibility Antigen, DRB1 Beta Chain (HLA-DRB1), metalloproteinase 9 (MMP9), and Phospholipase C-like 1 (PLCL1)-were identified as pivotal players in pathways such as "Antigen processing and presentation" and "Phagosome," which are crucially enriched in DSTB. Moreover, pharmaco-transcriptomic analysis can confirm that 58 drug compounds can regulate the expression of the key targets. Discussion: This research confirms the presence of protein alterations during the Rifampicin-resistant process in DSTB patients, offering novel insights into the molecular mechanisms underpinning DSTB. The findings suggest a promising avenue for the development of targeted drugs to enhance the management of Rifampicin-resistant spinal tuberculosis.

Cullin-3 proteins be a novel biomarkers and therapeutic targets for hyperchloremia induced by oral poisoning

Oral poisoning can trigger diverse physiological reactions, determined by the toxic substance involved. One such consequence is hyperchloremia, characterized by an elevated level of chloride in the blood and leads to kidney damage and impairing chloride ion regulation. Here, we conducted a comprehensive genome-wide analysis to investigate genes or proteins linked to hyperchloremia. Our analysis included functional enrichment, protein-protein interactions, gene expression, exploration of molecular pathways, and the identification of potential shared genetic factors contributing to the development of hyperchloremia. Functional enrichment analysis revealed that oral poisoning owing hyperchloremia is associated with 4 proteins e.g. Kelch-like protein 3, Serine/threonine-protein kinase WNK4, Serine/threonine-protein kinase WNK1 and Cullin-3. The protein-protein interaction network revealed Cullin-3 as an exceptional protein, displaying a maximum connection of 18 nodes. Insufficient data from transcriptomic analysis indicates that there are lack of information having direct associations between these proteins and human-related functions to oral poisoning, hyperchloremia, or metabolic acidosis. The metabolic pathway of Cullin-3 protein revealed that the derivative is Sulfonamide which play role in, increasing urine output, and metabolic acidosis resulted in hypertension. Based on molecular docking results analysis it found that Cullin-3 proteins has the lowest binding energies score and being suitable proteins. Moreover, no major variations were observed in unbound Cullin-3 and all three peptide bound complexes shows that all systems remain compact during 50 ns simulations. The results of our study revealed Cullin-3 proteins be a strong foundation for the development of potential drug targets or biomarker for future studies.

Single-cell systems pharmacology identifies development-driven drug response and combination therapy in B cell acute lymphoblastic leukemia

Leukemia can arise at various stages of the hematopoietic differentiation hierarchy, but the impact of developmental arrest on drug sensitivity is unclear. Applying network-based analyses to single-cell transcriptomes of human B cells, we define genome-wide signaling circuitry for each B cell differentiation stage. Using this reference, we comprehensively map the developmental states of B cell acute lymphoblastic leukemia (B-ALL), revealing its strong correlation with sensitivity to asparaginase, a commonly used chemotherapeutic agent. Single-cell multi-omics analyses of primary B-ALL blasts reveal marked intra-leukemia heterogeneity in asparaginase response: resistance is linked to pre-pro-B-like cells, with sensitivity associated with the pro-B-like population. By targeting BCL2, a driver within the pre-pro-B-like cell signaling network, we find that venetoclax significantly potentiates asparaginase efficacy in vitro and in vivo. These findings demonstrate a single-cell systems pharmacology framework to predict effective combination therapies based on intra-leukemia heterogeneity in developmental state, with potentially broad applications beyond B-ALL.

A Multivariate Genome-Wide Association Study Reveals Neural Correlates and Common Biological Mechanisms of Psychopathology Spectra

There is considerable comorbidity across externalizing and internalizing behavior dimensions of psychopathology. We applied genomic structural equation modeling (gSEM) to genome-wide association study (GWAS) summary statistics to evaluate the factor structure of externalizing and internalizing psychopathology across 16 traits and disorders among European-ancestry individuals (n's = 16,400 to 1,074,629). We conducted GWAS on factors derived from well-fitting models. Downstream analyses served to identify biological mechanisms, explore drug repurposing targets, estimate genetic overlap between the externalizing and internalizing spectra, and evaluate causal effects of psychopathology liability on physical health. Both a correlated factors model, comprising two factors of externalizing and internalizing risk, and a higher-order single-factor model of genetic effects contributing to both spectra demonstrated acceptable t. GWAS identified 409 lead single nucleotide polymorphisms (SNPs) associated with externalizing and 85 lead SNPs associated with internalizing, while the second-order GWAS identified 256 lead SNPs contributing to broad psychopathology risk. In bivariate causal mixture models, nearly all externalizing and internalizing causal variants overlapped, despite a genetic correlation of only 0.37 (SE = 0.02) between them. Externalizing genes showed cell-type specific expression in GABAergic, cortical, and hippocampal neurons, and internalizing genes were associated with reduced subcallosal cortical volume, providing insight into the neurobiological underpinnings of psychopathology. Genetic liability for externalizing, internalizing, and broad psychopathology exerted causal effects on pain, general health, cardiovascular diseases, and chronic illnesses. These findings underscore the complex genetic architecture of psychopathology, identify potential biological pathways for the externalizing and internalizing spectra, and highlight the physical health burden of psychiatric comorbidity.

Exploring the cross-cancer effect of circulating proteins and discovering potential intervention targets for 13 site-specific cancers

BACKGROUND

The proteome is an important reservoir of potential therapeutic targets for cancer. This study aimed to examine the causal associations between plasma proteins and cancer risk and to identify proteins with cross-cancer effects.

METHODS

Genetic instruments for 3991 plasma proteins were extracted from a large-scale proteomic study. Summary-level data of 13 site-specific cancers were derived from publicly available datasets. Proteome-wide Mendelian randomization and colocalization analyses were used to investigate the causal effect of circulating proteins on cancers. Protein-protein interactions and druggability assessment were conducted to prioritize potential therapeutic targets. Finally, systematical Mendelian randomization analysis between healthy lifestyle factors and cancer-related proteins was conducted to identify which proteins could act as interventional targets by lifestyle changes.

RESULTS

Genetically determined circulating levels of 58 proteins were statistically significantly associated with 7 site-specific cancers. A total of 39 proteins were prioritized by colocalization, of them, 11 proteins (ADPGK, CD86, CLSTN3, CSF2RA, CXCL10, GZMM, IL6R, NCR3, SIGLEC5, SIGLEC14, and TAPBP) were observed to have cross-cancer effects. Notably, 5 of these identified proteins (CD86, CSF2RA, CXCL10, IL6R, and TAPBP) have been targeted for drug development in cancer therapy; 8 proteins (ADPGK, CD86, CXCL10, GZMM, IL6R, SIGLEC5, SIGLEC14, TAPBP) could be modulated by healthy lifestyles.

CONCLUSION

Our study identified 39 circulating protein biomarkers with convincing causal evidence for 7 site-specific cancers, with 11 proteins demonstrating cross-cancer effects, and prioritized the proteins as potential intervention targets by either drugs or lifestyle changes, which provided new insights into the etiology, prevention, and treatment of cancers.

Genetic determinants of plasma protein levels in the Estonian population

The proteome holds great potential as an intermediate layer between the genome and phenome. Previous protein quantitative trait locus studies have focused mainly on describing the effects of common genetic variations on the proteome. Here, we assessed the impact of the common and rare genetic variations as well as the copy number variants (CNVs) on 326 plasma proteins measured in up to 500 individuals. We identified 184 cis and 94 trans signals for 157 protein traits, which were further fine-mapped to credible sets for 101 cis and 87 trans signals for 151 proteins. Rare genetic variation contributed to the levels of 7 proteins, with 5 cis and 14 trans associations. CNVs were associated with the levels of 11 proteins (7 cis and 5 trans), examples including a 3q12.1 deletion acting as a hub for multiple trans associations; and a CNV overlapping NAIP, a sensor component of the NAIP-NLRC4 inflammasome which is affecting pro-inflammatory cytokine interleukin 18 levels. In summary, this work presents a comprehensive resource of genetic variation affecting the plasma protein levels and provides the interpretation of identified effects.

A novel hypergraph model for identifying and prioritizing personalized drivers in cancer

Cancer development is driven by an accumulation of a small number of driver genetic mutations that confer the selective growth advantage to the cell, while most passenger mutations do not contribute to tumor progression. The identification of these driver genes responsible for tumorigenesis is a crucial step in designing effective cancer treatments. Although many computational methods have been developed with this purpose, the majority of existing methods solely provided a single driver gene list for the entire cohort of patients, ignoring the high heterogeneity of driver events across patients. It remains challenging to identify the personalized driver genes. Here, we propose a novel method (PDRWH), which aims to prioritize the mutated genes of a single patient based on their impact on the abnormal expression of downstream genes across a group of patients who share the co-mutation genes and similar gene expression profiles. The wide experimental results on 16 cancer datasets from TCGA showed that PDRWH excels in identifying known general driver genes and tumor-specific drivers. In the comparative testing across five cancer types, PDRWH outperformed existing individual-level methods as well as cohort-level methods. Our results also demonstrated that PDRWH could identify both common and rare drivers. The personalized driver profiles could improve tumor stratification, providing new insights into understanding tumor heterogeneity and taking a further step toward personalized treatment. We also validated one of our predicted novel personalized driver genes on tumor cell proliferation by vitro cell-based assays, the promoting effect of the high expression of Low-density lipoprotein receptor-related protein 1 (LRP1) on tumor cell proliferation.

Cuproptosis-Related Genes MTF1 and LIPT1 as Novel Prognostic Biomarker in Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is a life-threatening hematologic malignant disease with high morbidity and mortality in both adults and children. Cuproptosis, a novel mode of cell death, plays an important role in tumor development, but the functional mechanisms of cuproptosis-related genes (CRGs) in AML are unclear. The differential expression of CRGs between tumors such as AML and normal tissues in UCSC XENA, TCGA and GTEx was verified using R (version: 3.6.3). Lasso regression, Cox regression and Nomogram were used to screen for prognostic biomarkers of AML and to construct corresponding prognostic models. Kaplan-Meier analysis, ROC analysis, clinical correlation analysis, immune infiltration analysis and enrichment analysis were used to further investigate the correlation and functional mechanisms of CRGs with AML. The ceRNA regulatory network was used to identify the mRNA-miRNA-lncRNA regulatory axis. Cuproptosis-related genes LIPT1, MTF1, GLS and CDKN2A were highly expressed in AML, while FDX1, LIAS, DLD, DLAT, PDHA1, SLC31A1 and ATP7B were lowly expressed in AML. Lasso regression, Cox regression, Nomogram and calibration curve finally identified MTF1 and LIPT1 as two novel prognostic biomarkers of AML and constructed the corresponding prognostic models. In addition, all 12 CRGs had predictive power for AML, with MTF1, LIAS, SLC31A1 and CDKN2A showing more reliable results. Further analysis showed that ATP7B was closely associated with mutation types such as FLT3, NPM1, RAS and IDH1 R140 in AML, while the expression of MTF1, LIAS and ATP7B in AML was closely associated with immune infiltration. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA) revealed that biological functions such as metal ion transmembrane transporter activity, haptoglobin binding and oxygen carrier activity, pathways such as interferon alpha response, coagulation, UV response DN, apoptosis, hypoxia and heme metabolism all play a role in the development of AML. The ceRNA regulatory network revealed that 6 lncRNAs such as MALAT1, interfere with MTF1 expression through 6 miRNAs such as hsa-miR-32-5p, which in turn affect the development and progression of AML. In addition, APTO-253 has the potential to become an AML-targeted drug. The cuproptosis-related genes MTF1 and LIPT1 can be used as prognostic biomarkers in AML. A total of six lncRNAs, including MALAT1, are involved in the expression and regulation of MTF1 in AML through six miRNAs such as hsa-miR-32-5p.

Links between melanoma germline risk loci, driver genes and comorbidities: insight from a tissue-specific multi-omic analysis

Genome-wide association studies (GWAS) have associated 76 loci with the risk of developing melanoma. However, understanding the molecular basis of such associations has remained a challenge because most of these loci are in non-coding regions of the genome. Here, we integrated data on epigenomic markers, three-dimensional (3D) genome organization, and expression quantitative trait loci (eQTL) from melanoma-relevant tissues and cell types to gain novel insights into the mechanisms underlying melanoma risk. This integrative approach revealed a total of 151 target genes, both near and far away from the risk loci in linear sequence, with known and novel roles in the etiology of melanoma. Using protein-protein interaction networks, we identified proteins that interact-directly or indirectly-with the products of the target genes. The interacting proteins were enriched for known melanoma driver genes. Further integration of these target genes into tissue-specific gene regulatory networks revealed patterns of gene regulation that connect melanoma to its comorbidities. Our study provides novel insights into the biological implications of genetic variants associated with melanoma risk.

Identification and analysis of the molecular targets of statins in colorectal cancer

Colorectal cancer (CRC) is the third most common cancer in the world. According to several types of research, statins may impact the development and treatment of CRC. This work aimed to use bioinformatics to discover the relationship between statin targets and differentially expressed genes (DEGs) in CRC patients and determine the possible molecular effect of statins on CRC suppression. We used CRC datasets from the GEO database to select CRC-related DEGs. DGIdb and STITCH databases were used to identify gene targets of subtypes of statin. Further, we identified the statin target of CRC DEGs hub genes by using a Venn diagram of CRC DEGs and statin targets. Funrich and enrichr databases were carried out for the KEGG pathway and gene ontology (GO) enrichment analysis, respectively. GSE74604 and GSE10950 were used to identify CRC DEGs. After analyzing datasets,1370 genes were identified as CRC DEGs, and 345 targets were found for statins. We found that 35 genes are CRC DEGs statin targets. We found that statin targets in CRC were enriched in the receptor and metallopeptidase activity for molecular function, cytoplasm and plasma membrane for cellular component, signal transduction, and cell communication for biological process genes were substantially enriched based on FunRich enrichment. Analysis of the KEGG pathways revealed that the overexpressed DEGs were enriched in the IL-17, PPAR, and Toll-like receptor signaling pathways. Finally, CCNB1, DNMT1, AURKB, RAC1, PPARGC1A, CDKN1A, CAV1, IL1B, and HSPD1 were identified as hub CRC DEGs statin targets. The genetic and molecular aspects of our findings reveal that statins might have a therapeutic effect on CRC.

Identified senescence endotypes in aged cartilage are reflected in the blood metabolome

Heterogeneous accumulation of senescent cells expressing the senescence-associated secretory phenotype (SASP) affects tissue homeostasis which leads to diseases, such as osteoarthritis (OA). In this study, we set out to characterize heterogeneity of cellular senescence within aged articular cartilage and explored the presence of corresponding metabolic profiles in blood that could function as representative biomarkers. Hereto, we set out to perform cluster analyses, using a gene-set of 131 senescence genes (N = 57) in a previously established RNA sequencing dataset of aged articular cartilage and a generated metabolic dataset in overlapping blood samples. Using unsupervised hierarchical clustering and pathway analysis, we identified two robust cellular senescent endotypes. Endotype-1 was enriched for cell proliferating pathways, expressing forkhead box protein O4 (FOXO4), RB transcriptional corepressor like 2 (RBL2), and cyclin-dependent kinase inhibitor 1B (CDKN1B); the FOXO mediated cell cycle was identified as possible target for endotype-1 patients. Endotype-2 showed enriched inflammation-associated pathways, expressed by interleukin 6 (IL6), matrix metallopeptidase (MMP)1/3, and vascular endothelial growth factor (VEGF)C and SASP pathways were identified as possible targets for endotype-2 patients. Notably, plasma-based metabolic profiles in overlapping blood samples (N = 21) showed two corresponding metabolic clusters in blood. These non-invasive metabolic profiles could function as biomarkers for patient-tailored targeting of senescence in OA.

Strength of spatial correlation between gray matter connectivity and patterns of proto-oncogene and neural network construction gene expression is associated with diffuse glioma survival

Introduction

Like other forms of neuropathology, gliomas appear to spread along neural pathways. Accordingly, our group and others have previously shown that brain network connectivity is highly predictive of glioma survival. In this study, we aimed to examine the molecular mechanisms of this relationship via imaging transcriptomics.

Methods

We retrospectively obtained presurgical, T1-weighted MRI datasets from 669 adult patients, newly diagnosed with diffuse glioma. We measured brain connectivity using gray matter networks and coregistered these data with a transcriptomic brain atlas to determine the spatial co-localization between brain connectivity and expression patterns for 14 proto-oncogenes and 3 neural network construction genes.

Results

We found that all 17 genes were significantly co-localized with brain connectivity (p < 0.03, corrected). The strength of co-localization was highly predictive of overall survival in a cross-validated Cox Proportional Hazards model (mean area under the curve, AUC = 0.68 +/- 0.01) and significantly (p < 0.001) more so for a random forest survival model (mean AUC = 0.97 +/- 0.06). Bayesian network analysis demonstrated direct and indirect causal relationships among gene-brain co-localizations and survival. Gene ontology analysis showed that metabolic processes were overexpressed when spatial co-localization between brain connectivity and gene transcription was highest (p < 0.001). Drug-gene interaction analysis identified 84 potential candidate therapies based on our findings.

Discussion

Our findings provide novel insights regarding how gene-brain connectivity interactions may affect glioma survival.

SGCLDGA: unveiling drug-gene associations through simple graph contrastive learning

Drug repurposing offers a viable strategy for discovering new drugs and therapeutic targets through the analysis of drug-gene interactions. However, traditional experimental methods are plagued by their costliness and inefficiency. Despite graph convolutional network (GCN)-based models' state-of-the-art performance in prediction, their reliance on supervised learning makes them vulnerable to data sparsity, a common challenge in drug discovery, further complicating model development. In this study, we propose SGCLDGA, a novel computational model leveraging graph neural networks and contrastive learning to predict unknown drug-gene associations. SGCLDGA employs GCNs to extract vector representations of drugs and genes from the original bipartite graph. Subsequently, singular value decomposition (SVD) is employed to enhance the graph and generate multiple views. The model performs contrastive learning across these views, optimizing vector representations through a contrastive loss function to better distinguish positive and negative samples. The final step involves utilizing inner product calculations to determine association scores between drugs and genes. Experimental results on the DGIdb4.0 dataset demonstrate SGCLDGA's superior performance compared with six state-of-the-art methods. Ablation studies and case analyses validate the significance of contrastive learning and SVD, highlighting SGCLDGA's potential in discovering new drug-gene associations. The code and dataset for SGCLDGA are freely available at https://github.com/one-melon/SGCLDGA.

Identifying and validating MMP family members (MMP2, MMP9, MMP12, and MMP16) as therapeutic targets and biomarkers in kidney renal clear cell carcinoma (KIRC)

Kidney Renal Clear Cell Carcinoma (KIRC) is a malignant tumor that carries a substantial risk of morbidity and mortality. The MMP family assumes a crucial role in tumor invasion and metastasis. This study aimed to uncover the mechanistic relevance of the MMP gene family as a therapeutic target and diagnostic biomarker in Kidney Renal Clear Cell Carcinoma (KIRC) through a comprehensive approach encompassing both computational and molecular analyses. STRING, Cytoscape, UALCAN, GEPIA, OncoDB, HPA, cBioPortal, GSEA, TIMER, ENCORI, DrugBank, targeted bisulfite sequencing (bisulfite-seq), conventional PCR, Sanger sequencing, and RT-qPCR based analyses were used in the present study to analyze MMP gene family members to accurately determine a few hub genes that can be utilized as both therapeutic targets and diagnostic biomarkers for KIRC. By performing STRING and Cytohubba analyses of the 24 MMP gene family members, MMP2 (matrix metallopeptidase 2), MMP9 (matrix metallopeptidase 9), MMP12 (matrix metallopeptidase 12), and MMP16 (matrix metallopeptidase 16) genes were denoted as hub genes having highest degree scores. After analyzing MMP2, MMP9, MMP12, and MMP16 via various TCGA databases and RT-qPCR technique across clinical samples and KIRC cell lines, interestingly, all these hub genes were found significantly overexpressed at mRNA and protein levels in KIRC samples relative to controls. The notable effect of the up-regulated MMP2, MMP9, MMP12, and MMP16 was also documented on the overall survival (OS) of the KIRC patients. Moreover, targeted bisulfite-sequencing (bisulfite-seq) analysis revealed that promoter hypomethylation pattern was associated with up-regulation of hub genes (MMP2, MMP9, MMP12, and MMP16). In addition to this, hub genes were involved in various diverse oncogenic pathways. The MMP gene family members (MMP2, MMP9, MMP12, and MMP16) may serve as therapeutic targets and prognostic biomarkers in KIRC.

scDrugPrio: a framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases

BACKGROUND

Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs.

METHODS

Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs.

RESULTS

scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment.

CONCLUSIONS

We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package ( https://github.com/SDTC-CPMed/scDrugPrio ).