SIGNOR 3.0, the SIGnaling network open resource 3.0: 2022 update

Comprehensive molecular interaction map of TGFβ induced epithelial to mesenchymal transition in breast cancer

Breast cancer is one of the prevailing cancers globally, with a high mortality rate. Metastatic breast cancer (MBC) is an advanced stage of cancer, characterised by a highly nonlinear, heterogeneous process involving numerous singling pathways and regulatory interactions. Epithelial-mesenchymal transition (EMT) emerges as a key mechanism exploited by cancer cells. Transforming Growth Factor-β (TGFβ)-dependent signalling is attributed to promote EMT in advanced stages of breast cancer. A comprehensive regulatory map of TGFβ induced EMT was developed through an extensive literature survey. The network assembled comprises of 312 distinct species (proteins, genes, RNAs, complexes), and 426 reactions (state transitions, nuclear translocations, complex associations, and dissociations). The map was developed by following Systems Biology Graphical Notation (SBGN) using Cell Designer and made publicly available using MINERVA ( http://35.174.227.105:8080/minerva/?id=Metastatic_Breast_Cancer_1 ). While the complete molecular mechanism of MBC is still not known, the map captures the elaborate signalling interplay of TGFβ induced EMT-promoting MBC. Subsequently, the disease map assembled was translated into a Boolean model utilising CaSQ and analysed using Cell Collective. Simulations of these have captured the known experimental outcomes of TGFβ induced EMT in MBC. Hub regulators of the assembled map were identified, and their transcriptome-based analysis confirmed their role in cancer metastasis. Elaborate analysis of this map may help in gaining additional insights into the development and progression of metastatic breast cancer.

The landscape of cancer-rewired GPCR signaling axes

We explored the dysregulation of G-protein-coupled receptor (GPCR) ligand systems in cancer transcriptomics datasets to uncover new therapeutics opportunities in oncology. We derived an interaction network of receptors with ligands and their biosynthetic enzymes. Multiple GPCRs are differentially regulated together with their upstream partners across cancer subtypes and are associated to specific transcriptional programs and to patient survival patterns. The expression of both receptor-ligand (or enzymes) partners improved patient stratification, suggesting a synergistic role for the activation of GPCR networks in modulating cancer phenotypes. Remarkably, we identified many such axes across several cancer molecular subtypes, including many involving receptor-biosynthetic enzymes for neurotransmitters. We found that GPCRs from these actionable axes, including, e.g., muscarinic, adenosine, 5-hydroxytryptamine, and chemokine receptors, are the targets of multiple drugs displaying anti-growth effects in large-scale, cancer cell drug screens, which we further validated. We have made the results generated in this study freely available through a webapp (gpcrcanceraxes.bioinfolab.sns.it).

Inference of drug off-target effects on cellular signaling using interactome-based deep learning

Many diseases emerge from dysregulated cellular signaling, and drugs are often designed to target specific signaling proteins. Off-target effects are, however, common and may ultimately result in failed clinical trials. Here we develop a computer model of the cell's transcriptional response to drugs for improved understanding of their mechanisms of action. The model is based on ensembles of artificial neural networks and simultaneously infers drug-target interactions and their downstream effects on intracellular signaling. With this, it predicts transcription factors' activities, while recovering known drug-target interactions and inferring many new ones, which we validate with an independent dataset. As a case study, we analyze the effects of the drug Lestaurtinib on downstream signaling. Alongside its intended target, FLT3, the model predicts an inhibition of CDK2 that enhances the downregulation of the cell cycle-critical transcription factor FOXM1. Our approach can therefore enhance our understanding of drug signaling for therapeutic design.

Unveiling the signaling network of FLT3-ITD AML improves drug sensitivity prediction

Currently, the identification of patient-specific therapies in cancer is mainly informed by personalized genomic analysis. In the setting of acute myeloid leukemia (AML), patient-drug treatment matching fails in a subset of patients harboring atypical internal tandem duplications (ITDs) in the tyrosine kinase domain of the FLT3 gene. To address this unmet medical need, here we develop a systems-based strategy that integrates multiparametric analysis of crucial signaling pathways, and patient-specific genomic and transcriptomic data with a prior knowledge signaling network using a Boolean-based formalism. By this approach, we derive personalized predictive models describing the signaling landscape of AML FLT3-ITD positive cell lines and patients. These models enable us to derive mechanistic insight into drug resistance mechanisms and suggest novel opportunities for combinatorial treatments. Interestingly, our analysis reveals that the JNK kinase pathway plays a crucial role in the tyrosine kinase inhibitor response of FLT3-ITD cells through cell cycle regulation. Finally, our work shows that patient-specific logic models have the potential to inform precision medicine approaches.

Exploration of functional relations among differentially co-expressed genes identifies regulators in glioblastoma

The conventional computational approaches to investigating a disease confront inherent constraints as they often need to improve in delving beyond protein functional associations and grasping their deeper contextual significance within the disease framework. Such context-specificity can be explored using clinical data by evaluating the change in interaction between the biological entities in different conditions by investigating the differential co-expression relationships. We believe that the integration and analysis of differential co-expression and the functional relationships, primarily focusing on the source nodes, will open novel insights about disease progression as the source proteins could trigger signaling cascades, mostly because they are transcription factors, cell surface receptors, or enzymes that respond instantly to a particular stimulus. A thorough contextual investigation of these nodes could lead to a helpful beginning point for identifying potential causal linkages and guiding subsequent scientific investigations to uncover mechanisms underlying observed associations. Our methodology includes functional protein-protein Interaction (PPI) data and co-expression information and filters functional linkages through a series of critical steps, culminating in the identification of a robust set of regulators. Our analysis identified eleven key regulators-AKT1, BRCA1, CAMK2G, CUL1, FGFR3, KIF3A, NUP210, PRKACB, RAB8A, RPS6KA2 and TGFB3-in glioblastoma. These regulators play a pivotal role in disease classification, cell growth control, and patient survivability and exhibit associations with immune infiltrations and disease hallmarks. This underscores the importance of assessing correlation towards causality in unraveling complex biological insights.

Transomics2cytoscape: an automated software for interpretable 2.5-dimensional visualization of trans-omic networks

Biochemical network visualization is one of the essential technologies for mechanistic interpretation of omics data. In particular, recent advances in multi-omics measurement and analysis require the development of visualization methods that encompass multiple omics data. Visualization in 2.5 dimension (2.5D visualization), which is an isometric view of stacked X-Y planes, is a convenient way to interpret multi-omics/trans-omics data in the context of the conventional layouts of biochemical networks drawn on each of the stacked omics layers. However, 2.5D visualization of trans-omics networks is a state-of-the-art method that primarily relies on time-consuming human efforts involving manual drawing. Here, we present an R Bioconductor package 'transomics2cytoscape' for automated visualization of 2.5D trans-omics networks. We confirmed that transomics2cytoscape could be used for rapid visualization of trans-omics networks presented in published papers within a few minutes. Transomics2cytoscape allows for frequent update/redrawing of trans-omics networks in line with the progress in multi-omics/trans-omics data analysis, thereby enabling network-based interpretation of multi-omics data at each research step. The transomics2cytoscape source code is available at https://github.com/ecell/transomics2cytoscape .

Camptothecin structure simplification elaborated new imidazo[2,1-b]quinazoline derivative as a human topoisomerase I inhibitor with efficacy against bone cancer cells and colon adenocarcinoma

Camptothecin is a pentacyclic natural alkaloid that inhibits the hTop1 enzyme involved in DNA transcription and cancer cell growth. Camptothecin structure pitfalls prompted us to design new congeners using a structure simplification strategy to reduce the ring extension number from pentacyclic to tetracyclic while maintaining potential stacking of the new compounds with the DNA base pairs at the Top1-mediated cleavage complex and aqueous solubility, as well as minimizing compound-liver toxicity. The principal axis of this study was the verification of hTop1 inhibiting activity as a possible mechanism of action and the elaboration of new simplified inhibitors with improved pharmacodynamic and pharmacokinetic profiling using three structure panels (A-C) of (isoquinolinoimidazoquinazoline), (imidazoquinazoline), and (imidazoisoquinoline), respectively. DNA relaxation assay identified five compounds as hTop1 inhibitors belonging to the imidazoisoquinolines 3a,b, the imidazoquinazolines 12, and the isoquinolinoimidazoquinazolines 7a,b. In an MTT cytotoxicity assay against different cancer cell lines, compound 12 was the most potent against HOS bone cancer cells (IC50 = 1.47 μM). At the same time, the other inhibitors had no detectable activity against any cancer cell type. Compound (12) demonstrated great penetrating power in the HOS cancer cells' 3D-multicellular tumor spheroid model. Bioinformatics research of the hTop1 gene revealed that the TP53 cell proliferative gene is in the network of hTop1. The finding is confirmed empirically using the gene expression assay that proved the increase in p53 expression. The impact of structure simplification on compound 12 profile, characterized by the absence of acute oral liver toxicity when compared to Doxorubicin as a standard inhibitor, the lethal dose measured on Swiss Albino female mice and reported at LD50 = 250 mg/kg, and therapeutic significance in reducing colon adenocarcinoma tumor volume by 75.36 % after five weeks of treatment with compound 12. The molecular docking solutions of the active CPT-based derivative 12 and the inactive congener 14 into the active site of hTop1 and the activity cliffing of such MMP directed us to recommend the addition of HBD and HBA variables to compound 12 imidazoquinazoline core scaffold to enhance the potency via hydrogen bond formation with the major groove amino acids (Asp533, Lys532) as well as maintaining the hydrogen bond with the minor groove amino acid Arg364.

CardiOmics signatures reveal therapeutically actionable targets and drugs for cardiovascular diseases

Cardiovascular diseases are the leading cause of death worldwide, with heart failure being a complex condition that affects millions of individuals. Single-nucleus RNA sequencing has recently emerged as a powerful tool for unraveling the molecular mechanisms behind cardiovascular diseases. This cutting-edge technology enables the identification of molecular signatures, intracellular networks, and spatial relationships among cardiac cells, including cardiomyocytes, mast cells, lymphocytes, macrophages, lymphatic endothelial cells, endocardial cells, endothelial cells, epicardial cells, adipocytes, fibroblasts, neuronal cells, pericytes, and vascular smooth muscle cells. Despite these advancements, the discovery of essential therapeutic targets and drugs for precision cardiology remains a challenge. To bridge this gap, we conducted comprehensive in silico analyses of single-nucleus RNA sequencing data, functional enrichment, protein interactome network, and identification of the shortest pathways to physiological phenotypes. This integrated multi-omics analysis generated CardiOmics signatures, which allowed us to pinpoint three therapeutically actionable targets (ADRA1A1, PPARG, and ROCK2) and 15 effective drugs, including adrenergic receptor agonists, adrenergic receptor antagonists, norepinephrine precursors, PPAR receptor agonists, and Rho-associated kinase inhibitors, involved in late-stage cardiovascular disease clinical trials.

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

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

omicSynth: An open multi-omic community resource for identifying druggable targets across neurodegenerative diseases

Treatments for neurodegenerative disorders remain rare, but recent FDA approvals, such as lecanemab and aducanumab for Alzheimer disease (MIM: 607822), highlight the importance of the underlying biological mechanisms in driving discovery and creating disease modifying therapies. The global population is aging, driving an urgent need for therapeutics that stop disease progression and eliminate symptoms. In this study, we create an open framework and resource for evidence-based identification of therapeutic targets for neurodegenerative disease. We use summary-data-based Mendelian randomization to identify genetic targets for drug discovery and repurposing. In parallel, we provide mechanistic insights into disease processes and potential network-level consequences of gene-based therapeutics. We identify 116 Alzheimer disease, 3 amyotrophic lateral sclerosis (MIM: 105400), 5 Lewy body dementia (MIM: 127750), 46 Parkinson disease (MIM: 605909), and 9 progressive supranuclear palsy (MIM: 601104) target genes passing multiple test corrections (pSMR_multi < 2.95 × 10-6 and pHEIDI > 0.01). We created a therapeutic scheme to classify our identified target genes into strata based on druggability and approved therapeutics, classifying 41 novel targets, 3 known targets, and 115 difficult targets (of these, 69.8% are expressed in the disease-relevant cell type from single-nucleus experiments). Our novel class of genes provides a springboard for new opportunities in drug discovery, development, and repurposing in the pre-competitive space. In addition, looking at drug-gene interaction networks, we identify previous trials that may require further follow-up such as riluzole in Alzheimer disease. We also provide a user-friendly web platform to help users explore potential therapeutic targets for neurodegenerative diseases, decreasing activation energy for the community.

The landscape of SETBP1 gene expression and transcription factor activity across human tissues

The SET binding protein 1 (SETBP1) gene encodes a transcription factor (TF) involved in various cellular processes. Variants in SETBP1 can result in three different diseases determined by the introduction (germline vs. somatic) and location of the variant. Germline variants cause the ultra-rare pediatric Schinzel Giedion Syndrome (SGS) and SETBP1 haploinsufficiency disorder (SETBP1-HD), characterized by severe multisystemic abnormalities with neurodegeneration or a less severe brain phenotype accompanied by hypotonia and strabismus, respectively. Somatic variants in SETBP1 are associated with hematological malignancies and cancer development in other tissues in adults. To better understand the tissue-specific mechanisms involving SETBP1, we analyzed publicly available RNA-sequencing (RNA-seq) data from the Genotype-Tissue Expression (GTEx) project. We found SETBP1 and its known target genes were widely expressed across 31 adult human tissues. K-means clustering identified three distinct expression patterns of SETBP1 targets across tissues. Functional enrichment analysis (FEA) of each cluster revealed gene sets related to transcriptional regulation, DNA binding, and mitochondrial function. TF activity analysis of SETBP1 and its target TFs revealed tissue-specific TF activity, underscoring the role of tissue context-driven regulation and suggesting its impact in SETBP1-associated disease. In addition to uncovering tissue-specific molecular signatures of SETBP1 expression and TF activity, we provide a Shiny web application to facilitate exploring TF activity across human tissues for 758 TFs. This study provides insight into the landscape of SETBP1 expression and TF activity across 31 non-diseased human tissues and reveals tissue-specific expression and activity of SETBP1 and its targets. In conjunction with the web application we constructed, our framework enables researchers to generate hypotheses related to the role tissue backgrounds play with respect to gene expression and TF activity in different disease contexts.

Curation of causal interactions mediated by genes associated with autism accelerates the understanding of gene-phenotype relationships underlying neurodevelopmental disorders

Autism spectrum disorder (ASD) comprises a large group of neurodevelopmental conditions featuring, over a wide range of severity and combinations, a core set of manifestations (restricted sociality, stereotyped behavior and language impairment) alongside various comorbidities. Common and rare variants in several hundreds of genes and regulatory regions have been implicated in the molecular pathogenesis of ASD along a range of causation evidence strength. Despite significant progress in elucidating the impact of few paradigmatic individual loci, such sheer complexity in the genetic architecture underlying ASD as a whole has hampered the identification of convergent actionable hubs hypothesized to relay between the vastness of risk alleles and the core phenotypes. In turn this has limited the development of strategies that can revert or ameliorate this condition, calling for a systems-level approach to probe the cross-talk of cooperating genes in terms of causal interaction networks in order to make convergences experimentally tractable and reveal their clinical actionability. As a first step in this direction, we have captured from the scientific literature information on the causal links between the genes whose variants have been associated with ASD and the whole human proteome. This information has been annotated in a computer readable format in the SIGNOR database and is made freely available in the resource website. To link this information to cell functions and phenotypes, we have developed graph algorithms that estimate the functional distance of any protein in the SIGNOR causal interactome to phenotypes and pathways. The main novelty of our approach resides in the possibility to explore the mechanistic links connecting the suggested gene-phenotype relations.

Knowledge Graphs and Their Applications in Drug Discovery

Knowledge graphs represent information in the form of entities and relationships between those entities. Such a representation has multiple potential applications in drug discovery, including democratizing access to biomedical data, contextualizing or visualizing that data, and generating novel insights through the application of machine learning approaches. Knowledge graphs put data into context and therefore offer the opportunity to generate explainable predictions, which is a key topic in contemporary artificial intelligence. In this chapter, we outline some of the factors that need to be considered when constructing biomedical knowledge graphs, examine recent advances in mining such systems to gain insights for drug discovery, and identify potential future areas for further development.

Specifying cellular context of transcription factor regulons for exploring context-specific gene regulation programs

Understanding the role of transcription and transcription factors in cellular identity and disease, such as cancer and autoimmunity, is essential. However, comprehensive data resources for cell line-specific transcription factor-to-target gene annotations are currently limited. To address this, we developed a straightforward method to define regulons that capture the cell-specific aspects of TF binding and transcript expression levels. By integrating cellular transcriptome and transcription factor binding data, we generated regulons for four common cell lines comprising both proximal and distal cell line-specific regulatory events. Through systematic benchmarking involving transcription factor knockout experiments, we demonstrated performance on par with state-of-the-art methods, with our method being easily applicable to other cell types of interest. We present case studies using three cancer single-cell datasets to showcase the utility of these cell-type-specific regulons in exploring transcriptional dysregulation. In summary, this study provides a valuable tool and a resource for systematically exploring cell line-specific transcriptional regulations, emphasizing the utility of network analysis in deciphering disease mechanisms.

Precise Therapeutic Targeting of Distinct NRXN1+/- Mutations

As genetic studies continue to identify risk loci that are significantly associated with risk for neuropsychiatric disease, a critical unanswered question is the extent to which diverse mutations--sometimes impacting the same gene-- will require tailored therapeutic strategies. Here we consider this in the context of rare neuropsychiatric disorder-associated copy number variants (2p16.3) resulting in heterozygous deletions in NRXN1, a pre-synaptic cell adhesion protein that serves as a critical synaptic organizer in the brain. Complex patterns of NRXN1 alternative splicing are fundamental to establishing diverse neurocircuitry, vary between the cell types of the brain, and are differentially impacted by unique (non-recurrent) deletions. We contrast the cell-type-specific impact of patient-specific mutations in NRXN1 using human induced pluripotent stem cells, finding that perturbations in NRXN1 splicing result in divergent cell-type-specific synaptic outcomes. Via distinct loss-of-function (LOF) and gain-of-function (GOF) mechanisms, NRXN1+/- deletions cause decreased synaptic activity in glutamatergic neurons, yet increased synaptic activity in GABAergic neurons. Stratification of patients by LOF and GOF mechanisms will facilitate individualized restoration of NRXN1 isoform repertoires; towards this, antisense oligonucleotides knockdown mutant isoform expression and alters synaptic transcriptional signatures, while treatment with β-estradiol rescues synaptic function in glutamatergic neurons. Given the increasing number of mutations predicted to engender both LOF and GOF mechanisms in brain disease, our findings add nuance to future considerations of precision medicine.

The effects of pathogenic and likely pathogenic variants for inherited hemostasis disorders in 140 214 UK Biobank participants

Rare genetic diseases affect millions, and identifying causal DNA variants is essential for patient care. Therefore, it is imperative to estimate the effect of each independent variant and improve their pathogenicity classification. Our study of 140 214 unrelated UK Biobank (UKB) participants found that each of them carries a median of 7 variants previously reported as pathogenic or likely pathogenic. We focused on 967 diagnostic-grade gene (DGG) variants for rare bleeding, thrombotic, and platelet disorders (BTPDs) observed in 12 367 UKB participants. By association analysis, for a subset of these variants, we estimated effect sizes for platelet count and volume, and odds ratios for bleeding and thrombosis. Variants causal of some autosomal recessive platelet disorders revealed phenotypic consequences in carriers. Loss-of-function variants in MPL, which cause chronic amegakaryocytic thrombocytopenia if biallelic, were unexpectedly associated with increased platelet counts in carriers. We also demonstrated that common variants identified by genome-wide association studies (GWAS) for platelet count or thrombosis risk may influence the penetrance of rare variants in BTPD DGGs on their associated hemostasis disorders. Network-propagation analysis applied to an interactome of 18 410 nodes and 571 917 edges showed that GWAS variants with large effect sizes are enriched in DGGs and their first-order interactors. Finally, we illustrate the modifying effect of polygenic scores for platelet count and thrombosis risk on disease severity in participants carrying rare variants in TUBB1 or PROC and PROS1, respectively. Our findings demonstrate the power of association analyses using large population datasets in improving pathogenicity classifications of rare variants.

Prechordal structures act cooperatively in early trabeculae development of gnathostome skull

The vertebrate skull is formed by mesoderm and neural crest (NC) cells. The mesoderm contributes to the skull chordal domain, with the notochord playing an essential role in this process. The NC contributes to the skull prechordal domain, prompting investigation into the embryonic structures involved in prechordal neurocranium cartilage formation. The trabeculae cartilage, a structure of the prechordal neurocranium, arises at the convergence of prechordal plate (PCP), ventral midline (VM) cells of the diencephalon, and dorsal oral ectoderm. This study examines the molecular participation of these embryonic structures in gnathostome trabeculae development. PCP-secreted SHH induces its expression in VM cells of the diencephalon, initiating a positive feedback loop involving SIX3 and GLI1. SHH secreted by the VM cells of the diencephalon acts on the dorsal oral ectoderm, stimulating condensation of NC cells to form trabeculae. SHH from the prechordal region affects the expression of SOX9 in NC cells. BMP7 and SHH secreted by PCP induce NKX2.1 expression in VM cells of the diencephalon, but this does not impact trabeculae formation. Molecular cooperation between PCP, VM cells of the diencephalon, and dorsal oral ectoderm is crucial for craniofacial development by NC cells in the prechordal domain.

Monocytes subsets altered distribution and dysregulated plasma hsa-miR-21-5p and hsa-miR-155-5p in HCV-linked liver cirrhosis progression to hepatocellular carcinoma

PURPOSE

The authors aim to investigate the altered monocytes subsets distribution in liver cirrhosis (LC) and subsequent hepatocellular carcinoma (HCC) in association with the expression level of plasma Homo sapiens (has)-miR-21-5p and hsa-miR-155-5p. A step toward non-protein coding (nc) RNA precision medicine based on the immune perturbation manifested as altered monocytes distribution, on top of LC and HCC.

METHODS

Seventy-nine patients diagnosed with chronic hepatitis C virus (CHCV) infection with LC were enrolled in the current study. Patients were sub-classified into LC group without HCC (n = 40), LC with HCC (n = 39), and 15 apparently healthy controls. Monocyte subsets frequencies were assessed by flow cytometry. Real-time quantitative PCR was used to measure plasma hsa-miR-21-5p and hsa-miR-155-5p expression.

RESULTS

Hsa-miR-21-5p correlated with intermediate monocytes (r = 0.30, p = 0.007), while hsa-miR-155-5p negatively correlated with non-classical monocytes (r = - 0.316, p = 0.005). ROC curve analysis revealed that combining intermediate monocytes frequency and hsa-miR-21 yielded sensitivity = 79.5%, specificity = 75%, and AUC = 0.84. In comparison, AFP yielded a lower sensitivity = 69% and 100% specificity with AUC = 0.85. Logistic regression analysis proved that up-regulation of intermediate monocytes frequency and hsa-miR-21-5p were independent risk factors for LC progression to HCC, after adjustment for co-founders.

CONCLUSION

Monocyte subsets differentiation in HCC was linked to hsa-miR-21-5p and hsa-miR-155-5p. Combined up-regulation of intermediate monocytes frequency and hsa-miR-21-5p expression could be considered a sensitive indicator of LC progression to HCC. Circulating intermediate monocytes and hsa-miR-21-5p were independent risk factors for HCC evolution, clinically and in silico proved.

A framework for considering prior information in network-based approaches to omics data analysis

For decades, molecular biologists have been uncovering the mechanics of biological systems. Efforts to bring their findings together have led to the development of multiple databases and information systems that capture and present pathway information in a computable network format. Concurrently, the advent of modern omics technologies has empowered researchers to systematically profile cellular processes across different modalities. Numerous algorithms, methodologies, and tools have been developed to use prior knowledge networks (PKNs) in the analysis of omics datasets. Interestingly, it has been repeatedly demonstrated that the source of prior knowledge can greatly impact the results of a given analysis. For these methods to be successful it is paramount that their selection of PKNs is amenable to the data type and the computational task they aim to accomplish. Here we present a five-level framework that broadly describes network models in terms of their scope, level of detail, and ability to inform causal predictions. To contextualize this framework, we review a handful of network-based omics analysis methods at each level, while also describing the computational tasks they aim to accomplish.

The Cyclin-dependent kinase 1: more than a cell cycle regulator

The Cyclin-dependent kinase 1, as a serine/threonine protein kinase, is more than a cell cycle regulator as it was originally identified. During the last decade, it has been shown to carry out versatile functions during the last decade. From cell cycle control to gene expression regulation and apoptosis, CDK1 is intimately involved in many cellular events that are vital for cell survival. Here, we provide a comprehensive catalogue of the CDK1 upstream regulators and substrates, describing how this kinase is implicated in the control of key 'cell cycle-unrelated' biological processes. Finally, we describe how deregulation of CDK1 expression and activation has been closely associated with cancer progression and drug resistance.

The landscape of SETBP1 gene expression and transcription factor activity across human tissues

Background

The SET binding protein 1 (SETBP1) gene encodes a transcription factor (TF) involved in various cellular processes. Distinct SETBP1 variants have been linked to three different diseases. Germline variants cause the ultra-rare pediatric Schinzel Giedion Syndrome (SGS) and SETBP1 haploinsufficiency disorder (SETBP1-HD), characterized by severe multisystemic abnormalities with neurodegeneration or a less severe brain phenotype accompanied by hypotonia and strabismus, respectively. Somatic variants in SETBP1 are associated with hematological malignancies and cancer development in other tissues in adults.

Results

To better understand the tissue-specific mechanisms involving SETBP1, we analyzed publicly available RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project. We found SETBP1, and its known target genes were widely expressed across 31 adult human tissues. K-means clustering identified three distinct expression patterns of SETBP1 targets across tissues. Functional enrichment analysis (FEA) of each cluster revealed gene sets related to transcription regulation, DNA binding, and mitochondrial function. TF activity analysis of SETBP1 and its target TFs revealed tissue-specific TF activity, underscoring the role of tissue context-driven regulation and suggesting its impact in SETBP1-associated disease. In addition to uncovering tissue-specific molecular signatures of SETBP1 expression and TF activity, we provide a Shiny web application to facilitate exploring TF activity across human tissues for 758 TFs.

Conclusions

This study provides insight into the landscape of SETBP1 expression and TF activity across 31 non-diseased human tissues and reveals tissue-specific expression and activity of SETBP1 and its targets. In conjunction with the web application we constructed, our framework enables researchers to generate hypotheses related to the role tissue backgrounds play with respect to gene expression and TF activity in different disease contexts.

The landscape of cancer rewired GPCR signaling axes

We explored the dysregulation of GPCR ligand signaling systems in cancer transcriptomics datasets to uncover new therapeutics opportunities in oncology. We derived an interaction network of receptors with ligands and their biosynthetic enzymes, which revealed that multiple GPCRs are differentially regulated together with their upstream partners across cancer subtypes. We showed that biosynthetic pathway enrichment from enzyme expression recapitulated pathway activity signatures from metabolomics datasets, providing valuable surrogate information for GPCRs responding to organic ligands. We found that several GPCRs signaling components were significantly associated with patient survival in a cancer type-specific fashion. The expression of both receptor-ligand (or enzymes) partners improved patient stratification, suggesting a synergistic role for the activation of GPCR networks in modulating cancer phenotypes. Remarkably, we identified many such axes across several cancer molecular subtypes, including many pairs involving receptor-biosynthetic enzymes for neurotransmitters. We found that GPCRs from these actionable axes, including e.g., muscarinic, adenosine, 5-hydroxytryptamine and chemokine receptors, are the targets of multiple drugs displaying anti-growth effects in large-scale, cancer cell drug screens. We have made the results generated in this study freely available through a webapp (gpcrcanceraxes.bioinfolab.sns.it).

A Network of 17 Microtubule-Related Genes Highlights Functional Deregulations in Breast Cancer

A wide panel of microtubule-associated proteins and kinases is involved in coordinated regulation of the microtubule cytoskeleton and may thus represent valuable molecular markers contributing to major cellular pathways deregulated in cancer. We previously identified a panel of 17 microtubule-related (MT-Rel) genes that are differentially expressed in breast tumors showing resistance to taxane-based chemotherapy. In the present study, we evaluated the expression, prognostic value and functional impact of these genes in breast cancer. We show that 14 MT-Rel genes (KIF4A, ASPM, KIF20A, KIF14, TPX2, KIF18B, KIFC1, AURKB, KIF2C, GTSE1, KIF15, KIF11, RACGAP1, STMN1) are up-regulated in breast tumors compared with adjacent normal tissue. Six of them (KIF4A, ASPM, KIF20A, KIF14, TPX2, KIF18B) are overexpressed by more than 10-fold in tumor samples and four of them (KIF11, AURKB, TPX2 and KIFC1) are essential for cell survival. Overexpression of all 14 genes, and underexpression of 3 other MT-Rel genes (MAST4, MAPT and MTUS1) are associated with poor breast cancer patient survival. A Systems Biology approach highlighted three major functional networks connecting the 17 MT-Rel genes and their partners, which are centered on spindle assembly, chromosome segregation and cytokinesis. Our studies identified mitotic Aurora kinases and their substrates as major targets for therapeutic approaches against breast cancer.

Inferring single-cell gene regulatory network by non-redundant mutual information

Gene regulatory network plays a crucial role in controlling the biological processes of living creatures. Deciphering the complex gene regulatory networks from experimental data remains a major challenge in system biology. Recent advances in single-cell RNA sequencing technology bring massive high-resolution data, enabling computational inference of cell-specific gene regulatory networks (GRNs). Many relevant algorithms have been developed to achieve this goal in the past years. However, GRN inference is still less ideal due to the extra noises involved in pseudo-time information and large amounts of dropouts in datasets. Here, we present a novel GRN inference method named Normi, which is based on non-redundant mutual information. Normi manipulates these problems by employing a sliding size-fixed window approach on the entire trajectory and conducts average smoothing strategy on the gene expression of the cells in each window to obtain representative cells. To further alleviate the impact of dropouts, we utilize the mixed KSG estimator to quantify the high-order time-delayed mutual information among genes, then filter out the redundant edges by adopting Max-Relevance and Min Redundancy algorithm. Moreover, we determined the optimal time delay for each gene pair by distance correlation. Normi outperforms other state-of-the-art GRN inference methods on both simulated data and single-cell RNA sequencing (scRNA-seq) datasets, demonstrating its superiority in robustness. The performance of Normi in real scRNA-seq data further reveals its ability to identify the key regulators and crucial biological processes.

omicSynth: an Open Multi-omic Community Resource for Identifying Druggable Targets across Neurodegenerative Diseases

Treatments for neurodegenerative disorders remain rare, although recent FDA approvals, such as Lecanemab and Aducanumab for Alzheimer's Disease, highlight the importance of the underlying biological mechanisms in driving discovery and creating disease modifying therapies. The global population is aging, driving an urgent need for therapeutics that stop disease progression and eliminate symptoms. In this study, we create an open framework and resource for evidence-based identification of therapeutic targets for neurodegenerative disease. We use Summary-data-based Mendelian Randomization to identify genetic targets for drug discovery and repurposing. In parallel, we provide mechanistic insights into disease processes and potential network-level consequences of gene-based therapeutics. We identify 116 Alzheimer's disease, 3 amyotrophic lateral sclerosis, 5 Lewy body dementia, 46 Parkinson's disease, and 9 Progressive supranuclear palsy target genes passing multiple test corrections (pSMR_multi < 2.95×10-6 and pHEIDI > 0.01). We created a therapeutic scheme to classify our identified target genes into strata based on druggability and approved therapeutics - classifying 41 novel targets, 3 known targets, and 115 difficult targets (of these 69.8% are expressed in the disease relevant cell type from single nucleus experiments). Our novel class of genes provides a springboard for new opportunities in drug discovery, development and repurposing in the pre-competitive space. In addition, looking at drug-gene interaction networks, we identify previous trials that may require further follow-up such as Riluzole in AD. We also provide a user-friendly web platform to help users explore potential therapeutic targets for neurodegenerative diseases, decreasing activation energy for the community [https://nih-card-ndd-smr-home-syboky.streamlit.app/].

Identifying miRNA biomarkers for breast cancer and ovarian cancer: a text mining perspective

BACKGROUND

microRNA (miRNAs) are small, non-coding RNAs that mediate post-transcriptional gene silencing. Numerous studies have demonstrated the critical role of miRNAs in the development of breast cancer and ovarian cancer. To reduce potential bias from individual studies, a more comprehensive approach of exploring miRNAs in cancer research is essential. This study aims to explore the role of miRNAs in the development of breast cancer and ovarian cancer.

METHODS

Abstracts of the publications were tokenized and the biomedical terms (miRNA, gene, disease, species) were identified and extracted for vectorization. Predictive analyses were conducted with four machine learning models: K-Nearest Neighbors (KNN), Support Vector Machines (SVM), Random Forest (RF), and Naïve Bayes. Both holdout validation and cross-validation were utilized. Feature importance will be identified for miRNA-cancer networks construction.

RESULTS

We found that miR-182 is highly specific to female cancers. miR-182 targets different genes in regulating breast cancer and ovarian cancer. Naïve Bayes provided a promising prediction model for breast cancer and ovarian cancer with miRNAs and genes combination, with an accuracy score greater than 60%. Feature importance identified miR-155 and miR-199 are critical for breast cancer and ovarian cancer prediction, with miR-155 being highly related to breast cancer, whereas miR-199 being more associated with ovarian cancer.

CONCLUSION

Our approach effectively identified potential miRNA biomarkers associated with breast cancer and ovarian cancer, providing a solid foundation for generating novel research hypotheses and guiding future experimental studies.

Resources and tools for rare disease variant interpretation

Collectively, rare genetic disorders affect a substantial portion of the world's population. In most cases, those affected face difficulties in receiving a clinical diagnosis and genetic characterization. The understanding of the molecular mechanisms of these diseases and the development of therapeutic treatments for patients are also challenging. However, the application of recent advancements in genome sequencing/analysis technologies and computer-aided tools for predicting phenotype-genotype associations can bring significant benefits to this field. In this review, we highlight the most relevant online resources and computational tools for genome interpretation that can enhance the diagnosis, clinical management, and development of treatments for rare disorders. Our focus is on resources for interpreting single nucleotide variants. Additionally, we present use cases for interpreting genetic variants in clinical settings and review the limitations of these results and prediction tools. Finally, we have compiled a curated set of core resources and tools for analyzing rare disease genomes. Such resources and tools can be utilized to develop standardized protocols that will enhance the accuracy and effectiveness of rare disease diagnosis.

New trends in synthetic drugs and natural products targeting 20S proteasomes in cancers

Cancer is a global health challenge that remains to be a field of extensive research aiming to find new anticancer therapeutics. The 20S proteasome complex is one of the targets of anticancerdrugs, as it is correlated with several cancer types. Herein, we aim to discuss the 20S proteasome subunits and investigatethe currently studied proteasome inhibitors targeting the catalytically active proteasome subunits. In this review, we summarize the proteindegradation mechanism of the 20S proteasome complex and compareit with the 26S proteasome complex. Afterwards, the localization of the 20S proteasome is summarized as well as its use as a diagnosticandprognostic marker. The FDA-approved proteasome inhibitors (PIs) under clinical trials are summarized and their current limited use in solid tumors is also reviewed in addition to the expression of theβ5 subunit in differentcell lines. The review discusses in-silico analysis of the active subunit of the 20S proteasome complex. For development of new proteasome inhibitor drugs, the natural products inhibiting the 20S proteasome are summarized, as well as novel methodologies and challenges for the natural product discovery and current information about the biosynthetic gene clusters encoding them. We herein briefly summarize some resistancemechanismsto the proteasomeinhibitors. Additionally, we focus on the three main classes of proteasome inhibitors: 1] boronic acid, 2] beta-lactone and 3] epoxide inhibitor classes, as well as other PI classes, and their IC₅₀ values and their structure-activity relationship (SAR). Lastly,we summarize several future prospects of developing new proteasome inhibitors towards the treatment of tumors, especially solid tumors.

Recent advances in kinase signaling network profiling by mass spectrometry

Mass spectrometry-based phosphoproteomics is currently the leading methodology for the study of global kinase signaling. The scientific community is continuously releasing technological improvements for sensitive and fast identification of phosphopeptides, and their accurate quantification. To interpret large-scale phosphoproteomics data, numerous bioinformatic resources are available that help understanding kinase network functional role in biological systems upon perturbation. Some of these resources are databases of phosphorylation sites, protein kinases and phosphatases; others are bioinformatic algorithms to infer kinase activity, predict phosphosite functional relevance and visualize kinase signaling networks. In this review, we present the latest experimental and bioinformatic tools to profile protein kinase signaling networks and provide examples of their application in biomedicine.