Metastatic Infiltration of Nervous Tissue and Periosteal Nerve Sprouting in Multiple Myeloma-Induced Bone Pain in Mice and Human.

Multiple myeloma (MM) is a neoplasia of B plasma cells that often induces bone pain. However, the mechanisms underlying myeloma-induced bone pain (MIBP) are mostly unknown. Using a syngeneic MM mouse model, we show that periosteal nerve sprouting of calcitonin gene-related peptide (CGRP+) and growth associated protein 43 (GAP43+) fibers occurs concurrent to the onset of nociception and its blockade provides transient pain relief. MM patient samples also showed increased periosteal innervation. Mechanistically, we investigated MM induced gene expression changes in the dorsal root ganglia (DRG) innervating the MM-bearing bone of male mice and found alterations in pathways associated with cell cycle, immune response and neuronal signaling. The MM transcriptional signature was consistent with metastatic MM infiltration to the DRG, a never-before described feature of the disease that we further demonstrated histologically. In the DRG, MM cells caused loss of vascularization and neuronal injury, which may contribute to late-stage MIBP. Interestingly, the transcriptional signature of a MM patient was consistent with MM cell infiltration to the DRG. Overall, our results suggest that MM induces a plethora of peripheral nervous system alterations that may contribute to the failure of current analgesics and suggest neuroprotective drugs as appropriate strategies to treat early onset MIBP.SIGNIFICANCE STATEMENT Multiple myeloma (MM) is a painful bone marrow cancer that significantly impairs the quality of life of the patients. Analgesic therapies for myeloma-induced bone pain (MIBP) are limited and often ineffective, and the mechanisms of MIBP remain unknown. In this manuscript, we describe cancer-induced periosteal nerve sprouting in a mouse model of MIBP, where we also encounter metastasis to the dorsal root ganglia (DRG), a never-before described feature of the disease. Concomitant to myeloma infiltration, the lumbar DRGs presented blood vessel damage and transcriptional alterations, which may mediate MIBP. Explorative studies on human tissue support our preclinical findings. Understanding the mechanisms of MIBP is crucial to develop targeted analgesic with better efficacy and fewer side effects for this patient population.

Identifying the genes impacted by cell proliferation in proteomics and transcriptomics studies.

Hypothesis-free high-throughput profiling allows relative quantification of thousands of proteins or transcripts across samples and thereby identification of differentially expressed genes. It is used in many biological contexts to characterize differences between cell lines and tissues, identify drug mode of action or drivers of drug resistance, among others. Changes in gene expression can also be due to confounding factors that were not accounted for in the experimental plan, such as change in cell proliferation. We combined the analysis of 1,076 and 1,040 cell lines in five proteomics and three transcriptomics data sets to identify 157 genes that correlate with cell proliferation rates. These include actors in DNA replication and mitosis, and genes periodically expressed during the cell cycle. This signature of cell proliferation is a valuable resource when analyzing high-throughput data showing changes in proliferation across conditions. We show how to use this resource to help in interpretation of in vitro drug screens and tumor samples. It informs on differences of cell proliferation rates between conditions where such information is not directly available. The signature genes also highlight which hits in a screen may be due to proliferation changes; this can either contribute to biological interpretation or help focus on experiment-specific regulation events otherwise buried in the statistical analysis.

A time-resolved multi-omics atlas of Acanthamoeba castellanii encystment.

Encystment is a common stress response of most protists, including free-living amoebae. Cyst formation protects the amoebae from eradication and can increase virulence of the bacteria they harbor. Here, we mapped the global molecular changes that occur in the facultatively pathogenic amoeba Acanthamoeba castellanii during the early steps of the poorly understood process of encystment. By performing transcriptomic, proteomic, and phosphoproteomic experiments during encystment, we identified more than 150,000 previously undescribed transcripts and thousands of protein sequences absent from the reference genome. These results provide molecular details to the regulation of expected biological processes, such as cell proliferation shutdown, and reveal new insights such as a rapid phospho-regulation of sites involved in cytoskeleton remodeling and translation regulation. This work constitutes the first time-resolved molecular atlas of an encysting organism and a useful resource for further investigation of amoebae encystment to allow for a better control of pathogenic amoebae.

Differentially Expressed miRNAs in Ulcerative Colitis and Crohn's Disease.

Differential microRNA (miRNA or miR) regulation is linked to the development and progress of many diseases, including inflammatory bowel disease (IBD). It is well-established that miRNAs are involved in the differentiation, maturation, and functional control of immune cells. miRNAs modulate inflammatory cascades and affect the extracellular matrix, tight junctions, cellular hemostasis, and microbiota. This review summarizes current knowledge of differentially expressed miRNAs in mucosal tissues and peripheral blood of patients with ulcerative colitis and Crohn's disease. We combined comprehensive literature curation with computational meta-analysis of publicly available high-throughput datasets to obtain a consensus set of miRNAs consistently differentially expressed in mucosal tissues. We further describe the role of the most relevant differentially expressed miRNAs in IBD, extract their potential targets involved in IBD, and highlight their diagnostic and therapeutic potential for future investigations.

Homology-directed repair protects the replicating genome from metabolic assaults.

Homology-directed repair (HDR) safeguards DNA integrity under various forms of stress, but how HDR protects replicating genomes under extensive metabolic alterations remains unclear. Here, we report that besides stalling replication forks, inhibition of ribonucleotide reductase (RNR) triggers metabolic imbalance manifested by the accumulation of increased reactive oxygen species (ROS) in cell nuclei. This leads to a redox-sensitive activation of the ATM kinase followed by phosphorylation of the MRE11 nuclease, which in HDR-deficient settings degrades stalled replication forks. Intriguingly, nascent DNA degradation by the ROS-ATM-MRE11 cascade is also triggered by hypoxia, which elevates signaling-competent ROS and attenuates functional HDR without arresting replication forks. Under these conditions, MRE11 degrades daughter-strand DNA gaps, which accumulate behind active replisomes and attract error-prone DNA polymerases to escalate mutation rates. Thus, HDR safeguards replicating genomes against metabolic assaults by restraining mutagenic repair at aberrantly processed nascent DNA. These findings have implications for cancer evolution and tumor therapy.

Identification of hyper-rewired genomic stress non-oncogene addiction genes across 15 cancer types.

Non-oncogene addiction (NOA) genes are essential for supporting the stress-burdened phenotype of tumours and thus vital for their survival. Although NOA genes are acknowledged to be potential drug targets, there has been no large-scale attempt to identify and characterise them as a group across cancer types. Here we provide the first method for the identification of conditional NOA genes and their rewired neighbours using a systems approach. Using copy number data and expression profiles from The Cancer Genome Atlas (TCGA) we performed comparative analyses between high and low genomic stress tumours for 15 cancer types. We identified 101 condition-specific differential coexpression modules, mapped to a high-confidence human interactome, comprising 133 candidate NOA rewiring hub genes. We observe that most modules lose coexpression in the high-stress state and that activated stress modules and hubs take part in homoeostasis maintenance processes such as chromosome segregation, oxireductase activity, mitotic checkpoint (PLK1 signalling), DNA replication initiation and synaptic signalling. We furthermore show that candidate NOA rewiring hubs are unique for each cancer type, but that their respective rewired neighbour genes largely are shared across cancer types.

ProtFus: A Comprehensive Method Characterizing Protein-Protein Interactions of Fusion Proteins.

Tailored therapy aims to cure cancer patients effectively and safely, based on the complex interactions between patients' genomic features, disease pathology and drug metabolism. Thus, the continual increase in scientific literature drives the need for efficient methods of data mining to improve the extraction of useful information from texts based on patients' genomic features. An important application of text mining to tailored therapy in cancer encompasses the use of mutations and cancer fusion genes as moieties that change patients' cellular networks to develop cancer, and also affect drug metabolism. Fusion proteins, which are derived from the slippage of two parental genes, are produced in cancer by chromosomal aberrations and trans-splicing. Given that the two parental proteins for predicted fusion proteins are known, we used our previously developed method for identifying chimeric protein-protein interactions (ChiPPIs) associated with the fusion proteins. Here, we present a validation approach that receives fusion proteins of interest, predicts their cellular network alterations by ChiPPI and validates them by our new method, ProtFus, using an online literature search. This process resulted in a set of 358 fusion proteins and their corresponding protein interactions, as a training set for a Naïve Bayes classifier, to identify predicted fusion proteins that have reliable evidence in the literature and that were confirmed experimentally. Next, for a test group of 1817 fusion proteins, we were able to identify from the literature 2908 PPIs in total, across 18 cancer types. The described method, ProtFus, can be used for screening the literature to identify unique cases of fusion proteins and their PPIs, as means of studying alterations of protein networks in cancers. Availability: http://protfus.md.biu.ac.il/.

Specifying RNA-Binding Regions in Proteins by Peptide Cross-Linking and Affinity Purification.

RNA-binding proteins (RBPs) allow cells to carry out pre-RNA processing and post-transcriptional regulation of gene expression, and aberrations in RBP functions have been linked to many diseases, including neurological disorders and cancer. Human cells encode thousands of RNA-binding proteins with unique RNA-binding properties. These properties are regulated through modularity of a large variety of RNA-binding domains, rendering RNA-protein interactions difficult to study. Recently, the introduction of proteomics methods has provided novel insights into RNA-binding proteins at a systems level. However, determining the exact protein sequence regions that interact with RNA remains challenging and laborious, especially considering that many RBPs lack canonical RNA-binding domains. Here we describe a streamlined proteomic workflow called peptide cross-linking and affinity purification (pCLAP) that allows rapid characterization of RNA-binding regions in proteins. pCLAP is based upon the combined use of UV cross-linking and enzymatic digestion of RNA-bound proteins followed by single-shot mass spectrometric analysis. To benchmark our method, we identified the binding regions for polyadenylated RNA-binding proteins in HEK293 cells, allowing us to map the mRNA interaction regions of more than 1000 RBPs with very high reproducibility from replicate single-shot analyses. Our results show specific enrichment of many known RNA-binding regions on many known RNA-binding proteins, confirming the specificity of our approach.

Mass spectrometry of human leukocyte antigen class I peptidomes reveals strong effects of protein abundance and turnover on antigen presentation.

HLA class I molecules reflect the health state of cells to cytotoxic T cells by presenting a repertoire of endogenously derived peptides. However, the extent to which the proteome shapes the peptidome is still largely unknown. Here we present a high-throughput mass-spectrometry-based workflow that allows stringent and accurate identification of thousands of such peptides and direct determination of binding motifs. Applying the workflow to seven cancer cell lines and primary cells, yielded more than 22,000 unique HLA peptides across different allelic binding specificities. By computing a score representing the HLA-I sampling density, we show a strong link between protein abundance and HLA-presentation (p < 0.0001). When analyzing overpresented proteins - those with at least fivefold higher density score than expected for their abundance - we noticed that they are degraded almost 3 h faster than similar but nonpresented proteins (top 20% abundance class; median half-life 20.8h versus 23.6h, p < 0.0001). This validates protein degradation as an important factor for HLA presentation. Ribosomal, mitochondrial respiratory chain, and nucleosomal proteins are particularly well presented. Taking a set of proteins associated with cancer, we compared the predicted immunogenicity of previously validated T-cell epitopes with other peptides from these proteins in our data set. The validated epitopes indeed tend to have higher immunogenic scores than the other detected HLA peptides. Remarkably, we identified five mutated peptides from a human colon cancer cell line, which have very recently been predicted to be HLA-I binders. Altogether, we demonstrate the usefulness of combining MS-analysis with immunogenesis prediction for identifying, ranking, and selecting peptides for therapeutic use.

DISEASES: text mining and data integration of disease-gene associations.

Text mining is a flexible technology that can be applied to numerous different tasks in biology and medicine. We present a system for extracting disease-gene associations from biomedical abstracts. The system consists of a highly efficient dictionary-based tagger for named entity recognition of human genes and diseases, which we combine with a scoring scheme that takes into account co-occurrences both within and between sentences. We show that this approach is able to extract half of all manually curated associations with a false positive rate of only 0.16%. Nonetheless, text mining should not stand alone, but be combined with other types of evidence. For this reason, we have developed the DISEASES resource, which integrates the results from text mining with manually curated disease-gene associations, cancer mutation data, and genome-wide association studies from existing databases. The DISEASES resource is accessible through a web interface at http://diseases.jensenlab.org/, where the text-mining software and all associations are also freely available for download.

Temporal disease trajectories condensed from population-wide registry data covering 6.2 million patients.

A key prerequisite for precision medicine is the estimation of disease progression from the current patient state. Disease correlations and temporal disease progression (trajectories) have mainly been analysed with focus on a small number of diseases or using large-scale approaches without time consideration, exceeding a few years. So far, no large-scale studies have focused on defining a comprehensive set of disease trajectories. Here we present a discovery-driven analysis of temporal disease progression patterns using data from an electronic health registry covering the whole population of Denmark. We use the entire spectrum of diseases and convert 14.9 years of registry data on 6.2 million patients into 1,171 significant trajectories. We group these into patterns centred on a small number of key diagnoses such as chronic obstructive pulmonary disease (COPD) and gout, which are central to disease progression and hence important to diagnose early to mitigate the risk of adverse outcomes. We suggest such trajectory analyses may be useful for predicting and preventing future diseases of individual patients.

Protein-driven inference of miRNA-disease associations.

MOTIVATION: MicroRNAs (miRNAs) are a highly abundant class of non-coding RNA genes involved in cellular regulation and thus also diseases. Despite miRNAs being important disease factors, miRNA-disease associations remain low in number and of variable reliability. Furthermore, existing databases and prediction methods do not explicitly facilitate forming hypotheses about the possible molecular causes of the association, thereby making the path to experimental follow-up longer. RESULTS: Here we present miRPD in which miRNA-Protein-Disease associations are explicitly inferred. Besides linking miRNAs to diseases, it directly suggests the underlying proteins involved, which can be used to form hypotheses that can be experimentally tested. The inference of miRNAs and diseases is made by coupling known and predicted miRNA-protein associations with protein-disease associations text mined from the literature. We present scoring schemes that allow us to rank miRNA-disease associations inferred from both curated and predicted miRNA targets by reliability and thereby to create high- and medium-confidence sets of associations. Analyzing these, we find statistically significant enrichment for proteins involved in pathways related to cancer and type I diabetes mellitus, suggesting either a literature bias or a genuine biological trend. We show by example how the associations can be used to extract proteins for disease hypothesis. AVAILABILITY AND IMPLEMENTATION: All datasets, software and a searchable Web site are available at http://mirpd.jensenlab.org.

Drug-induced regulation of target expression.

Drug perturbations of human cells lead to complex responses upon target binding. One of the known mechanisms is a (positive or negative) feedback loop that adjusts the expression level of the respective target protein. To quantify this mechanism systems-wide in an unbiased way, drug-induced differential expression of drug target mRNA was examined in three cell lines using the Connectivity Map. To overcome various biases in this valuable resource, we have developed a computational normalization and scoring procedure that is applicable to gene expression recording upon heterogeneous drug treatments. In 1290 drug-target relations, corresponding to 466 drugs acting on 167 drug targets studied, 8% of the targets are subject to regulation at the mRNA level. We confirmed systematically that in particular G-protein coupled receptors, when serving as known targets, are regulated upon drug treatment. We further newly identified drug-induced differential regulation of Lanosterol 14-alpha demethylase, Endoplasmin, DNA topoisomerase 2-alpha and Calmodulin 1. The feedback regulation in these and other targets is likely to be relevant for the success or failure of the molecular intervention.

Linear motif atlas for phosphorylation-dependent signaling.

Systematic and quantitative analysis of protein phosphorylation is revealing dynamic regulatory networks underlying cellular responses to environmental cues. However, matching these sites to the kinases that phosphorylate them and the phosphorylation-dependent binding domains that may subsequently bind to them remains a challenge. NetPhorest is an atlas of consensus sequence motifs that covers 179 kinases and 104 phosphorylation-dependent binding domains [Src homology 2 (SH2), phosphotyrosine binding (PTB), BRCA1 C-terminal (BRCT), WW, and 14-3-3]. The atlas reveals new aspects of signaling systems, including the observation that tyrosine kinases mutated in cancer have lower specificity than their non-oncogenic relatives. The resource is maintained by an automated pipeline, which uses phylogenetic trees to structure the currently available in vivo and in vitro data to derive probabilistic sequence models of linear motifs. The atlas is available as a community resource (http://netphorest.info).

Systematic discovery of in vivo phosphorylation networks.

Protein kinases control cellular decision processes by phosphorylating specific substrates. Thousands of in vivo phosphorylation sites have been identified, mostly by proteome-wide mapping. However, systematically matching these sites to specific kinases is presently infeasible, due to limited specificity of consensus motifs, and the influence of contextual factors, such as protein scaffolds, localization, and expression, on cellular substrate specificity. We have developed an approach (NetworKIN) that augments motif-based predictions with the network context of kinases and phosphoproteins. The latter provides 60%-80% of the computational capability to assign in vivo substrate specificity. NetworKIN pinpoints kinases responsible for specific phosphorylations and yields a 2.5-fold improvement in the accuracy with which phosphorylation networks can be constructed. Applying this approach to DNA damage signaling, we show that 53BP1 and Rad50 are phosphorylated by CDK1 and ATM, respectively. We describe a scalable strategy to evaluate predictions, which suggests that BCLAF1 is a GSK-3 substrate.