Computational drug prediction in hepatoblastoma by integrating pan-cancer transcriptomics with pharmacological response.

BACKGROUND AND AIMS: Hepatoblastoma (HB) is the predominant form of pediatric liver cancer, though it remains exceptionally rare. While treatment outcomes for children with HB have improved, patients with advanced tumors face limited therapeutic choices. Additionally, survivors often suffer from long-term adverse effects due to treatment, including ototoxicity, cardiotoxicity, delayed growth, and secondary tumors. Consequently, there is a pressing need to identify new and effective therapeutic strategies for patients with HB. Computational methods to predict drug sensitivity from a tumor's transcriptome have been successfully applied for some common adult malignancies, but specific efforts in pediatric cancers are lacking because of the paucity of data. APPROACH AND RESULTS: In this study, we used DrugSense to assess drug efficacy in patients with HB, particularly those with the aggressive C2 subtype associated with poor clinical outcomes. Our method relied on publicly available collections of pan-cancer transcriptional profiles and drug responses across 36 tumor types and 495 compounds. The drugs predicted to be most effective were experimentally validated using patient-derived xenograft models of HB grown in vitro and in vivo. We thus identified 2 cyclin-dependent kinase 9 inhibitors, alvocidib and dinaciclib as potent HB growth inhibitors for the high-risk C2 molecular subtype. We also found that in a cohort of 46 patients with HB, high cyclin-dependent kinase 9 tumor expression was significantly associated with poor prognosis. CONCLUSIONS: Our work proves the usefulness of computational methods trained on pan-cancer data sets to reposition drugs in rare pediatric cancers such as HB, and to help clinicians in choosing the best treatment options for their patients.

A yeast synthetic network for in vivo assessment of reverse-engineering and modeling approaches.

Systems biology approaches are extensively used to model and reverse engineer gene regulatory networks from experimental data. Conversely, synthetic biology allows "de novo" construction of a regulatory network to seed new functions in the cell. At present, the usefulness and predictive ability of modeling and reverse engineering cannot be assessed and compared rigorously. We built in the yeast Saccharomyces cerevisiae a synthetic network, IRMA, for in vivo "benchmarking" of reverse-engineering and modeling approaches. The network is composed of five genes regulating each other through a variety of regulatory interactions; it is negligibly affected by endogenous genes, and it is responsive to small molecules. We measured time series and steady-state expression data after multiple perturbations. These data were used to assess state-of-the-art modeling and reverse-engineering techniques. A semiquantitative model was able to capture and predict the behavior of the network. Reverse engineering based on differential equations and Bayesian networks correctly inferred regulatory interactions from the experimental data.

Definition of the transcriptional units of inherited retinal disease genes by meta-analysis of human retinal transcriptome data.

BACKGROUND: Inherited retinal diseases (IRD) are genetically heterogeneous disorders that cause the dysfunction or loss of photoreceptor cells and ultimately lead to blindness. To date, next-generation sequencing procedures fail to detect pathogenic sequence variants in coding regions of known IRD disease genes in about 30-40% of patients. One of the possible explanations for this missing heritability is the presence of yet unidentified transcripts of known IRD genes. Here, we aimed to define the transcript composition of IRD genes in the human retina by a meta-analysis of publicly available RNA-seq datasets using an ad-hoc designed pipeline. RESULTS: We analysed 218 IRD genes and identified 5,054 transcripts, 3,367 of which were not previously reported. We assessed their putative expression levels and focused our attention on 435 transcripts predicted to account for at least 5% of the expression of the corresponding gene. We looked at the possible impact of the newly identified transcripts at the protein level and experimentally validated a subset of them. CONCLUSIONS: This study provides an unprecedented, detailed overview of the complexity of the human retinal transcriptome that can be instrumental in contributing to the resolution of some cases of missing heritability in IRD patients.

Pharmacoproteomics pinpoints HSP70 interaction for correction of the most frequent Wilson disease-causing mutant of ATP7B.

Pathogenic mutations in the copper transporter ATP7B have been hypothesized to affect its protein interaction landscape contributing to loss of function and, thereby, to hepatic copper toxicosis in Wilson disease. Although targeting mutant interactomes was proposed as a therapeutic strategy, druggable interactors for rescue of ATP7B mutants remain elusive. Using proteomics, we found that the frequent H1069Q substitution promotes ATP7B interaction with HSP70, thus accelerating endoplasmic reticulum (ER) degradation of the mutant protein and consequent copper accumulation in hepatic cells. This prompted us to use an HSP70 inhibitor as bait in a bioinformatics search for structurally similar Food and Drug Administration-approved drugs. Among the hits, domperidone emerged as an effective corrector that recovered trafficking and function of ATP7B-H1069Q by impairing its exposure to the HSP70 proteostatic network. Our findings suggest that HSP70-mediated degradation can be safely targeted with domperidone to rescue ER-retained ATP7B mutants and, hence, to counter the onset of Wilson disease.

A Siamese neural network model for the prioritization of metabolic disorders by integrating real and simulated data.

MOTIVATION: Untargeted metabolomic approaches hold a great promise as a diagnostic tool for inborn errors of metabolisms (IEMs) in the near future. However, the complexity of the involved data makes its application difficult and time consuming. Computational approaches, such as metabolic network simulations and machine learning, could significantly help to exploit metabolomic data to aid the diagnostic process. While the former suffers from limited predictive accuracy, the latter is normally able to generalize only to IEMs for which sufficient data are available. Here, we propose a hybrid approach that exploits the best of both worlds by building a mapping between simulated and real metabolic data through a novel method based on Siamese neural networks (SNN). RESULTS: The proposed SNN model is able to perform disease prioritization for the metabolic profiles of IEM patients even for diseases that it was not trained to identify. To the best of our knowledge, this has not been attempted before. The developed model is able to significantly outperform a baseline model that relies on metabolic simulations only. The prioritization performances demonstrate the feasibility of the method, suggesting that the integration of metabolic models and data could significantly aid the IEM diagnosis process in the near future. AVAILABILITY AND IMPLEMENTATION: Metabolic datasets used in this study are publicly available from the cited sources. The original data produced in this study, including the trained models and the simulated metabolic profiles, are also publicly available (Messa et al., 2020).

Drug Repurposing in Rare Diseases: An Integrative Study of Drug Screening and Transcriptomic Analysis in Nephropathic Cystinosis.

Diagnosis and cure for rare diseases represent a great challenge for the scientific community who often comes up against the complexity and heterogeneity of clinical picture associated to a high cost and time-consuming drug development processes. Here we show a drug repurposing strategy applied to nephropathic cystinosis, a rare inherited disorder belonging to the lysosomal storage diseases. This approach consists in combining mechanism-based and cell-based screenings, coupled with an affordable computational analysis, which could result very useful to predict therapeutic responses at both molecular and system levels. Then, we identified potential drugs and metabolic pathways relevant for the pathophysiology of nephropathic cystinosis by comparing gene-expression signature of drugs that share common mechanisms of action or that involve similar pathways with the disease gene-expression signature achieved with RNA-seq.

gep2pep: a Bioconductor package for the creation and analysis of pathway-based expression profiles.

SUMMARY: Pathway-based expression profiles allow for high-level interpretation of transcriptomic data and systematic comparison of dysregulated cellular programs. We have previously demonstrated the efficacy of pathway-based approaches with two different applications: the Drug Set Enrichment Analysis and the Gene2drug analysis. Here we present a software tool that allows to easily convert gene-based profiles to pathway-based profiles and analyze them within the popular R framework. We also provide pre-computed profiles derived from the original Connectivity Map and its next generation release, i.e. the LINCS database. AVAILABILITY AND IMPLEMENTATION: the tool is implemented as the R/Bioconductor package gep2pep and can be freely downloaded from https://bioconductor.org/packages/gep2pep. SUPPLEMENTARY INFORMATION: Supplementary data are available at http://dsea.tigem.it/lincs.

TRPV4 and purinergic receptor signalling pathways are separately linked in airway epithelia to CFTR and TMEM16A chloride channels.

KEY POINTS: Eact is a putative pharmacological activator of TMEM16A. Eact is strongly effective in recombinant Fischer rat thyroid (FRT) cells but not in airway epithelial cells with endogenous TMEM16A expression. Transcriptomic analysis, gene silencing and functional studies in FRT cells reveal that Eact is actually an activator of the Ca2+ -permeable TRPV4 channel. In airway epithelial cells TRPV4 and TMEM16A are expressed in separate cell types. Intracellular Ca2+ elevation by TRPV4 stimulation leads to CFTR channel activation. ABSTRACT: TMEM16A is a Ca2+ -activated Cl- channel expressed in airway epithelial cells, particularly under conditions of mucus hypersecretion. To investigate the role of TMEM16A, we used Eact, a putative TMEM16A pharmacological activator. However, in contrast to purinergic stimulation, we found little effect of Eact on bronchial epithelial cells under conditions of high TMEM16A expression. We hypothesized that Eact is an indirect activator of TMEM16A. By a combination of approaches, including short-circuit current recordings, bulk and single cell RNA sequencing, intracellular Ca2+ imaging and RNA interference, we found that Eact is actually an activator of the Ca2+ -permeable TRPV4 channel and that the modest effect of this compound in bronchial epithelial cells is due to a separate expression of TMEM16A and TRPV4 in different cell types. Importantly, we found that TRPV4 stimulation induced activation of the CFTR Cl- channel. Our study reveals the existence of separate Ca2+ signalling pathways linked to different Cl- secretory processes.

Comprehensive kinome NGS targeted expression profiling by KING-REX.

BACKGROUND: Protein kinases are enzymes controlling different cellular functions. Genetic alterations often result in kinase dysregulation, making kinases a very attractive class of druggable targets in several human diseases. Existing approved drugs still target a very limited portion of the human 'kinome', demanding a broader functional knowledge of individual and co-expressed kinase patterns in physiologic and pathologic settings. The development of novel rapid and cost-effective methods for kinome screening is therefore highly desirable, potentially leading to the identification of novel kinase drug targets. RESULTS: In this work, we describe the development of KING-REX (KINase Gene RNA EXpression), a comprehensive kinome RNA targeted custom assay-based panel designed for Next Generation Sequencing analysis, coupled with a dedicated data analysis pipeline. We have conceived KING-REX for the gene expression analysis of 512 human kinases; for 319 kinases, paired assays and custom analysis pipeline features allow the evaluation of 3'- and 5'-end transcript imbalances as readout for the prediction of gene rearrangements. Validation tests on cell line models harboring known gene fusions demonstrated a comparable accuracy of KING-REX gene expression assessment as in whole transcriptome analyses, together with a robust detection of transcript portion imbalances in rearranged kinases, even in complex RNA mixtures or in degraded RNA. CONCLUSIONS: These results support the use of KING-REX as a rapid and cost effective kinome investigation tool in the field of kinase target identification for applications in cancer biology and other human diseases.

gene2drug: a computational tool for pathway-based rational drug repositioning.

Motivation: Drug repositioning has been proposed as an effective shortcut to drug discovery. The availability of large collections of transcriptional responses to drugs enables computational approaches to drug repositioning directly based on measured molecular effects. Results: We introduce a novel computational methodology for rational drug repositioning, which exploits the transcriptional responses following treatment with small molecule. Specifically, given a therapeutic target gene, a prioritization of potential effective drugs is obtained by assessing their impact on the transcription of genes in the pathway(s) including the target. We performed in silico validation and comparison with a state-of-art technique based on similar principles. We next performed experimental validation in two different real-case drug repositioning scenarios: (i) upregulation of the glutamate-pyruvate transaminase (GPT), which has been shown to induce reduction of oxalate levels in a mouse model of primary hyperoxaluria, and (ii) activation of the transcription factor TFEB, a master regulator of lysosomal biogenesis and autophagy, whose modulation may be beneficial in neurodegenerative disorders. Availability and implementation: A web tool for Gene2drug is freely available at http://gene2drug.tigem.it. An R package is under development and can be obtained from https://github.com/franapoli/gep2pep. Contact: dibernardo@tigem.it. Supplementary information: Supplementary data are available at Bioinformatics online.

High-resolution analysis of the human retina miRNome reveals isomiR variations and novel microRNAs.

MicroRNAs play a fundamental role in retinal development and function. To characterise the miRNome of the human retina, we carried out deep sequencing analysis on sixteen individuals. We established the catalogue of retina-expressed miRNAs, determined their relative abundance and found that a small number of miRNAs accounts for almost 90% of the retina miRNome. We discovered more than 3000 miRNA variants (isomiRs), encompassing a wide range of sequence variations, which include seed modifications that are predicted to have an impact on miRNA action. We demonstrated that a seed-modifying isomiR of the retina-enriched miR-124-3p was endowed with different targeting properties with respect to the corresponding canonical form. Moreover, we identified 51 putative novel, retina-specific miRNAs and experimentally validated the expression for nine of them. Finally, a parallel analysis of the human Retinal Pigment Epithelium (RPE)/choroid, two tissues that are known to be crucial for retina homeostasis, yielded notably distinct miRNA enrichment patterns compared to the retina. The generated data are accessible through an ad hoc database. This study is the first to reveal the complexity of the human retina miRNome at nucleotide resolution and constitutes a unique resource to assess the contribution of miRNAs to the pathophysiology of the human retina.

An atlas of gene expression and gene co-regulation in the human retina.

The human retina is a specialized tissue involved in light stimulus transduction. Despite its unique biology, an accurate reference transcriptome is still missing. Here, we performed gene expression analysis (RNA-seq) of 50 retinal samples from non-visually impaired post-mortem donors. We identified novel transcripts with high confidence (Observed Transcriptome (ObsT)) and quantified the expression level of known transcripts (Reference Transcriptome (RefT)). The ObsT included 77 623 transcripts (23 960 genes) covering 137 Mb (35 Mb new transcribed genome). Most of the transcripts (92%) were multi-exonic: 81% with known isoforms, 16% with new isoforms and 3% belonging to new genes. The RefT included 13 792 genes across 94 521 known transcripts. Mitochondrial genes were among the most highly expressed, accounting for about 10% of the reads. Of all the protein-coding genes in Gencode, 65% are expressed in the retina. We exploited inter-individual variability in gene expression to infer a gene co-expression network and to identify genes specifically expressed in photoreceptor cells. We experimentally validated the photoreceptors localization of three genes in human retina that had not been previously reported. RNA-seq data and the gene co-expression network are available online (http://retina.tigem.it).

Drug-set enrichment analysis: a novel tool to investigate drug mode of action.

MOTIVATION: Automated screening approaches are able to rapidly identify a set of small molecules inducing a desired phenotype from large small-molecule libraries. However, the resulting set of candidate molecules is usually very diverse pharmacologically, thus little insight on the shared mechanism of action (MoA) underlying their efficacy can be gained. RESULTS: We introduce a computational method (Drug-Set Enrichment Analysis-DSEA) based on drug-induced gene expression profiles, which is able to identify the molecular pathways that are targeted by most of the drugs in the set. By diluting drug-specific effects unrelated to the phenotype of interest, DSEA is able to highlight phenotype-specific pathways, thus helping to formulate hypotheses on the MoA shared by the drugs in the set. We validated the method by analysing five different drug-sets related to well-known pharmacological classes. We then applied DSEA to identify the MoA shared by drugs known to be partially effective in rescuing mutant cystic fibrosis transmembrane conductance regulator (CFTR) gene function in Cystic Fibrosis. AVAILABILITY AND IMPLEMENTATION: The method is implemented as an online web tool publicly available at http://dsea.tigem.it. CONTACT: dibernardo@tigem.it SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Deficiency of multidrug resistance 2 contributes to cell transformation through oxidative stress.

Multidrug resistance 2 (Mdr2), also called adenosine triphosphate-binding cassette B4 (ABCB4), is the transporter of phosphatidylcholine (PC) at the canalicular membrane of mouse hepatocytes, which plays an essential role for bile formation. Mutations in human homologue MDR3 are associated with several liver diseases. Knockout of Mdr2 results in hepatic inflammation, liver fibrosis and hepatocellular carcinoma (HCC). Whereas the pathogenesis in Mdr2 (-/-) mice has been largely attributed to the toxicity of bile acids due to the absence of PC in the bile, the question of whether Mdr2 deficiency per se perturbs biological functions in the cell has been poorly addressed. As Mdr2 is expressed in many cell types, we used mouse embryonic fibroblasts (MEF) derived from Mdr2 (-/-) embryos to show that deficiency of Mdr2 increases reactive oxygen species accumulation, lipid peroxidation and DNA damage. We found that Mdr2 (-/-) MEFs undergo spontaneous transformation and that Mdr2 (-/-) mice are more susceptible to chemical carcinogen-induced intestinal tumorigenesis. Microarray analysis in Mdr2-/- MEFs and cap analysis of gene expression in Mdr2 (-/-) HCCs revealed extensively deregulated genes involved in oxidation reduction, fatty acid metabolism and lipid biosynthesis. Our findings imply a close link between Mdr2 (-/-) -associated tumorigenesis and perturbation of these biological processes and suggest potential extrahepatic functions of Mdr2/MDR3.

A reverse-engineering approach to dissect post-translational modulators of transcription factor's activity from transcriptional data.

BACKGROUND: Transcription factors (TFs) act downstream of the major signalling pathways functioning as master regulators of cell fate. Their activity is tightly regulated at the transcriptional, post-transcriptional and post-translational level. Proteins modifying TF activity are not easily identified by experimental high-throughput methods. RESULTS: We developed a computational strategy, called Differential Multi-Information (DMI), to infer post-translational modulators of a transcription factor from a compendium of gene expression profiles (GEPs). DMI is built on the hypothesis that the modulator of a TF (i.e. kinase/phosphatases), when expressed in the cell, will cause the TF target genes to be co-expressed. On the contrary, when the modulator is not expressed, the TF will be inactive resulting in a loss of co-regulation across its target genes. DMI detects the occurrence of changes in target gene co-regulation for each candidate modulator, using a measure called Multi-Information. We validated the DMI approach on a compendium of 5,372 GEPs showing its predictive ability in correctly identifying kinases regulating the activity of 14 different transcription factors. CONCLUSIONS: DMI can be used in combination with experimental approaches as high-throughput screening to efficiently improve both pathway and target discovery. An on-line web-tool enabling the user to use DMI to identify post-transcriptional modulators of a transcription factor of interest che be found at http://dmi.tigem.it.

The phytoestrogen genistein modulates lysosomal metabolism and transcription factor EB (TFEB) activation.

Genistein (5,7-dihydroxy-3-(4-hydroxyphenyl)-4H-1-benzopyran-4-one) has been previously proposed as a potential drug for use in substrate reduction therapy for mucopolysaccharidoses, a group of inherited metabolic diseases caused by mutations leading to inefficient degradation of glycosaminoglycans (GAGs) in lysosomes. It was demonstrated that this isoflavone can cross the blood-brain barrier, making it an especially desirable potential drug for the treatment of neurological symptoms present in most lysosomal storage diseases. So far, no comprehensive genomic analyses have been performed to elucidate the molecular mechanisms underlying the effect elicited by genistein. Therefore, the aim of this work was to identify the genistein-modulated gene network regulating GAG biosynthesis and degradation, taking into consideration the entire lysosomal metabolism. Our analyses identified over 60 genes with known roles in lysosomal biogenesis and/or function whose expression was enhanced by genistein. Moreover, 19 genes whose products are involved in both GAG synthesis and degradation pathways were found to be remarkably differentially regulated by genistein treatment. We found a regulatory network linking genistein-mediated control of transcription factor EB (TFEB) gene expression, TFEB nuclear translocation, and activation of TFEB-dependent lysosome biogenesis to lysosomal metabolism. Our data indicate that the molecular mechanism of genistein action involves not only impairment of GAG synthesis but more importantly lysosomal enhancement via TFEB. These findings contribute to explaining the beneficial effects of genistein in lysosomal storage diseases as well as envisage new therapeutic approaches to treat these devastating diseases.

Blood transcriptomics of drug-naïve sporadic Parkinson's disease patients.

BACKGROUND: Parkinson's disease (PD) is a chronic progressive neurodegenerative disorder that is clinically defined in terms of motor symptoms. These are preceded by prodromal non-motor manifestations that prove the systemic nature of the disease. Identifying genes and pathways altered in living patients provide new information on the diagnosis and pathogenesis of sporadic PD. METHODS: Changes in gene expression in the blood of 40 sporadic PD patients and 20 healthy controls ("Discovery set") were analyzed by taking advantage of the Affymetrix platform. Patients were at the onset of motor symptoms and before initiating any pharmacological treatment. Data analysis was performed by applying Ranking-Principal Component Analysis, PUMA and Significance Analysis of Microarrays. Functional annotations were assigned using GO, DAVID, GSEA to unveil significant enriched biological processes in the differentially expressed genes. The expressions of selected genes were validated using RT-qPCR and samples from an independent cohort of 12 patients and controls ("Validation set"). RESULTS: Gene expression profiling of blood samples discriminates PD patients from healthy controls and identifies differentially expressed genes in blood. The majority of these are also present in dopaminergic neurons of the Substantia Nigra, the key site of neurodegeneration. Together with neuronal apoptosis, lymphocyte activation and mitochondrial dysfunction, already found in previous analysis of PD blood and post-mortem brains, we unveiled transcriptome changes enriched in biological terms related to epigenetic modifications including chromatin remodeling and methylation. Candidate transcripts as CBX5, TCF3, MAN1C1 and DOCK10 were validated by RT-qPCR. CONCLUSIONS: Our data support the use of blood transcriptomics to study neurodegenerative diseases. It identifies changes in crucial components of chromatin remodeling and methylation machineries as early events in sporadic PD suggesting epigenetics as target for therapeutic intervention.

Mantra 2.0: an online collaborative resource for drug mode of action and repurposing by network analysis.

SUMMARY: Elucidation of molecular targets of a compound [mode of action (MoA)] and its off-targets is a crucial step in drug development. We developed an online collaborative resource (MANTRA 2.0) that supports this process by exploiting similarities between drug-induced transcriptional profiles. Drugs are organized in a network of nodes (drugs) and edges (similarities) highlighting 'communities' of drugs sharing a similar MoA. A user can upload gene expression profiles before and after drug treatment in one or multiple cell types. An automated processing pipeline transforms the gene expression profiles into a unique drug 'node' embedded in the drug-network. Visual inspection of the neighbouring drugs and communities helps in revealing its MoA and to suggest new applications of known drugs (drug repurposing). MANTRA 2.0 allows storing and sharing user-generated network nodes, thus making MANTRA 2.0 a collaborative ever-growing resource. AVAILABILITY AND IMPLEMENTATION: The web tool is freely available for academic use at http://mantra.tigem.it.

In-vivo real-time control of protein expression from endogenous and synthetic gene networks.

We describe an innovative experimental and computational approach to control the expression of a protein in a population of yeast cells. We designed a simple control algorithm to automatically regulate the administration of inducer molecules to the cells by comparing the actual protein expression level in the cell population with the desired expression level. We then built an automated platform based on a microfluidic device, a time-lapse microscopy apparatus, and a set of motorized syringes, all controlled by a computer. We tested the platform to force yeast cells to express a desired fixed, or time-varying, amount of a reporter protein over thousands of minutes. The computer automatically switched the type of sugar administered to the cells, its concentration and its duration, according to the control algorithm. Our approach can be used to control expression of any protein, fused to a fluorescent reporter, provided that an external molecule known to (indirectly) affect its promoter activity is available.

Reverse engineering a mouse embryonic stem cell-specific transcriptional network reveals a new modulator of neuronal differentiation.

Gene expression profiles can be used to infer previously unknown transcriptional regulatory interaction among thousands of genes, via systems biology 'reverse engineering' approaches. We 'reverse engineered' an embryonic stem (ES)-specific transcriptional network from 171 gene expression profiles, measured in ES cells, to identify master regulators of gene expression ('hubs'). We discovered that E130012A19Rik (E13), highly expressed in mouse ES cells as compared with differentiated cells, was a central 'hub' of the network. We demonstrated that E13 is a protein-coding gene implicated in regulating the commitment towards the different neuronal subtypes and glia cells. The overexpression and knock-down of E13 in ES cell lines, undergoing differentiation into neurons and glia cells, caused a strong up-regulation of the glutamatergic neurons marker Vglut2 and a strong down-regulation of the GABAergic neurons marker GAD65 and of the radial glia marker Blbp. We confirmed E13 expression in the cerebral cortex of adult mice and during development. By immuno-based affinity purification, we characterized protein partners of E13, involved in the Polycomb complex. Our results suggest a role of E13 in regulating the division between glutamatergic projection neurons and GABAergic interneurons and glia cells possibly by epigenetic-mediated transcriptional regulation.