The identification of drug targets is highly challenging, particularly for diseases of the brain. To address this problem, we developed and experimentally validated a general computational framework for drug target discovery that combines gene regulatory information with causal reasoning (“Causal Reasoning Analytical Framework for Target discovery”—CRAFT). Using a systems genetics approach and starting from gene expression data from the target tissue, CRAFT provides a predictive framework for identifying cell membrane receptors with a direction-specified influence over disease-related gene expression profiles. As proof of concept, we applied CRAFT to epilepsy and predicted the tyrosine kinase receptor Csf1R as a potential therapeutic target. The predicted effect of Csf1R blockade in attenuating epilepsy seizures was validated in three pre-clinical models of epilepsy. These results highlight CRAFT as a systems-level framework for target discovery and suggest Csf1R blockade as a novel therapeutic strategy in epilepsy. CRAFT is applicable to disease settings other than epilepsy.

BackgroundThe relationship between monogenic and polygenic forms of epilepsy is poorly understood and the extent to which the genetic and acquired epilepsies share common pathways is unclear. Here, we use an integrated systems-level analysis of brain gene expression data to identify molecular networks disrupted in epilepsy.ResultsWe identified a co-expression network of 320 genes (M30), which is significantly enriched for non-synonymous de novo mutations ascertained from patients with monogenic epilepsy and for common variants associated with polygenic epilepsy. The genes in the M30 network are expressed widely in the human brain under tight developmental control and encode physically interacting proteins involved in synaptic processes. The most highly connected proteins within the M30 network were preferentially disrupted by deleterious de novo mutations for monogenic epilepsy, in line with the centrality-lethality hypothesis. Analysis of M30 expression revealed consistent downregulation in the epileptic brain in heterogeneous forms of epilepsy including human temporal lobe epilepsy, a mouse model of acquired temporal lobe epilepsy, and a mouse model of monogenic Dravet (SCN1A) disease. These results suggest functional disruption of M30 via gene mutation or altered expression as a convergent mechanism regulating susceptibility to epilepsy broadly. Using the large collection of drug-induced gene expression data from Connectivity Map, several drugs were predicted to preferentially restore the downregulation of M30 in epilepsy toward health, most notably valproic acid, whose effect on M30 expression was replicated in neurons.ConclusionsTaken together, our results suggest targeting the expression of M30 as a potential new therapeutic strategy in epilepsy.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-1097-7) contains supplementary material, which is available to authorized users.

Glioblastomas are aggressive astrocytomas characterized by endothelial cell proliferation and abnormal vasculature, which can cause brain edema and increase patient morbidity. We identified the heparin-binding cytokine pleiotrophin as a driver of vascular abnormalization in glioma. Pleiotrophin abundance was greater in high-grade human astrocytomas and correlated with poor survival. Anaplastic lymphoma kinase (ALK), which is a receptor that is activated by pleiotrophin, was present in mural cells associated with abnormal vessels. Orthotopically implanted gliomas formed from GL261 cells that were engineered to produce pleiotrophin showed increased microvessel density and enhanced tumor growth compared with gliomas formed from control GL261 cells. The survival of mice with pleiotrophin-producing gliomas was shorter than that of mice with gliomas that did not produce pleiotrophin. Vessels in pleiotrophin-producing gliomas were poorly perfused and abnormal, a phenotype that was associated with increased deposition of vascular endothelial growth factor (VEGF) in direct proximity to the vasculature. The growth of pleiotrophin-producing GL261 gliomas was inhibited by treatment with the ALK inhibitor crizotinib, the ALK inhibitor ceritinib, or the VEGF receptor inhibitor cediranib, whereas control GL261 tumors did not respond to either inhibitor. Our findings link pleiotrophin abundance in gliomas with survival in humans and mice, and show that pleiotrophin promotes glioma progression through increased VEGF deposition and vascular abnormalization.

NMT is a powerful method for dissecting tumour microglia and macrophage subpopulations and can guide further investigation of BMDC subsets in glioma and neuro-inflammatory diseases.

Pleiotrophin (PTN) augments tumor growth by increasing proliferation of tumor cells and promoting vascular abnormalization, but its role in early gliomagenesis has not been evaluated. Through analysis of publically available datasets, we demonstrate that increased PTN mRNA expression is associated with amplification of chromosome 7, identified as one of the earliest steps in glioblastoma development. To elucidate the role of PTN in tumor initiation we employed the RCAS/tv-a model that allows glioma induction by RCAS-virus mediated expression of oncogenes in neural progenitor cells. Intracranial injection of RCAS-PTN did not induce glioma formation when administrated alone, but significantly enhanced RCAS-platelet derived growth factor (PDGF)B-induced gliomagenesis. PTN co-treatment augmented PDGFB-induced Akt activation in neural progenitor cells in vitro, and enhanced neural sphere size associated with increased proliferation. Our data indicates that PTN expression is associated with chromosome 7 gain, and that PTN enhances PDGFB-induced gliomagenesis by stimulating proliferation of neural progenitor cells.

Objective To identify a convergent, multitarget proliferation characteristic for astrocytoma transformation that could be targeted for therapy discovery. Methods Using an integrated functional genomics approach, we prioritized networks associated with astrocytoma progression using the following criteria: differential co‐expression between grade II and grade III IDH1‐mutated and 1p/19q euploid astrocytomas, preferential enrichment for genetic risk to cancer, association with patient survival and sample‐level genomic features. Drugs targeting the identified multitarget network characteristic for astrocytoma transformation were computationally predicted using drug transcriptional perturbation data and validated using primary human astrocytoma cells. Results A single network, M2, consisting of 177 genes, was associated with glioma progression on the basis of the above criteria. Functionally, M2 encoded physically interacting proteins regulating cell cycle processes and analysis of genome‐wide gene‐regulatory interactions using mutual information and DNA–protein interactions revealed the known regulators of cell cycle processes FoxM1, B‐Myb, and E2F2 as key regulators of M2. These results suggest functional disruption of M2 via gene mutation or altered expression as a convergent pathway regulating astrocytoma transformation. By considering M2 as a multitarget drug target regulating astrocytoma transformation, we identified several drugs that are predicted to restore M2 expression in anaplastic astrocytoma toward its low‐grade profile and of these, we validated the known antiproliferative drug resveratrol as down‐regulating multiple nodes of M2 including at nanomolar concentrations achievable in human cerebrospinal fluid by oral dosing. Interpretation Our results identify M2 as a multitarget network characteristic for astrocytoma progression and encourage M2‐based drug screening to identify new compounds for preventing glioma transformation.

Human bronchial epithelial cells (HBECs) have first-line contact with harmful substances during smoking, and changes in their metabolism most likely represent a defining factor in coping with the stress and development of airway diseases. This study was designed to determine the dynamics of metabolome changes in HBECs treated with cigarette smoke condensate (CSC), and to test whether normal metabolism can be restored by synthetic antioxidants. Principal component analysis, based on untargeted mass spectra, indicated that treatment of CSC-exposed HBECs with O-methyl-L-tyrosinyl-γ-L-glutamyl-L-cysteinylglycine (UPF1) acted faster than did N-acetylcysteine to revert the effect of CSC. The maximum effect of 10 μg/ml CSC itself on HBEC cell line, BEAS-2B, metabolism was seen at 2 hours after treatment, with return to the baseline level by 7 hours. In primary HBECs, the initial maximum effect was seen at 1 hour after CSC exposure. Certain metabolites associated with redox pathways and energy production were affected by CSC. Subsequent restoration of their content by UPF1 supports the hypothetical protective capacity of UPF1 against the oxidative stress and increased energy demand, respectively. Furthermore, UPF1 up-regulated the contents of phospholipid species identified as phosphatidylcholines and phosphatidylethanolamines in the CSC-exposed HBECs, indicating possible suppression of inflammatory processes along with an increase in spermidine as an endogenous cytoprotector. In conclusion, with this dynamic metabolomics study, we characterize the durability of the CSC-induced metabolic changes in BEAS-2B line cells and primary HBECs, and demonstrate the ability of UPF1 to significantly accelerate the recovery of HBECs from CSC insult.

Translating genome-wide association loci to therapies requires knowledge of the causal genes, their directionality of effect and the cell-types in which they act. To infer these relationships in the human brain, we implemented Mendelian randomisation using single cell-type expression quantitative trait loci (eQTLs) as genetic anchors. Expression QTLs were mapped across 8 major cell-types in brain tissue exclusively ascertained from donors with no history of brain disease. We report evidence for a causal association between the change in expression of 118 genes and one or more of 16 brain phenotypes, revealing candidate targets for risk mitigation and opportunities for shared preventative therapeutic strategies. We highlight key causal genes for neurodegenerative and neuropsychiatric disease and for each, we report its cellular context and the therapeutic directionality required for risk mitigation. Our use of control samples establishes a new resource for the causal interpretation of GWAS risk alleles for human brain phenotypes.