Melanoma 2.0. Skin cancer as a paradigm for emerging diagnostic technologies, computational modelling and artificial intelligence.

We live in an unprecedented time in oncology. We have accumulated samples and cases in cohorts larger and more complex than ever before. New technologies are available for quantifying solid or liquid samples at the molecular level. At the same time, we are now equipped with the computational power necessary to handle this enormous amount of quantitative data. Computational models are widely used helping us to substantiate and interpret data. Under the label of systems and precision medicine, we are putting all these developments together to improve and personalize the therapy of cancer. In this review, we use melanoma as a paradigm to present the successful application of these technologies but also to discuss possible future developments in patient care linked to them. Melanoma is a paradigmatic case for disruptive improvements in therapies, with a considerable number of metastatic melanoma patients benefiting from novel therapies. Nevertheless, a large proportion of patients does not respond to therapy or suffers from adverse events. Melanoma is an ideal case study to deploy advanced technologies not only due to the medical need but also to some intrinsic features of melanoma as a disease and the skin as an organ. From the perspective of data acquisition, the skin is the ideal organ due to its accessibility and suitability for many kinds of advanced imaging techniques. We put special emphasis on the necessity of computational strategies to integrate multiple sources of quantitative data describing the tumour at different scales and levels.

Monitoring the maturation of the sarcomere network: a super-resolution microscopy-based approach.

The in vitro generation of human cardiomyocytes derived from induced pluripotent stem cells (iPSC) is of great importance for cardiac disease modeling, drug-testing applications and for regenerative medicine. Despite the development of various cultivation strategies, a sufficiently high degree of maturation is still a decisive limiting factor for the successful application of these cardiac cells. The maturation process includes, among others, the proper formation of sarcomere structures, mediating the contraction of cardiomyocytes. To precisely monitor the maturation of the contractile machinery, we have established an imaging-based strategy that allows quantitative evaluation of important parameters, defining the quality of the sarcomere network. iPSC-derived cardiomyocytes were subjected to different culture conditions to improve sarcomere formation, including prolonged cultivation time and micro patterned surfaces. Fluorescent images of α-actinin were acquired using super-resolution microscopy. Subsequently, we determined cell morphology, sarcomere density, filament alignment, z-Disc thickness and sarcomere length of iPSC-derived cardiomyocytes. Cells from adult and neonatal heart tissue served as control. Our image analysis revealed a profound effect on sarcomere content and filament orientation when iPSC-derived cardiomyocytes were cultured on structured, line-shaped surfaces. Similarly, prolonged cultivation time had a beneficial effect on the structural maturation, leading to a more adult-like phenotype. Automatic evaluation of the sarcomere filaments by machine learning validated our data. Moreover, we successfully transferred this approach to skeletal muscle cells, showing an improved sarcomere formation cells over different differentiation periods. Overall, our image-based workflow can be used as a straight-forward tool to quantitatively estimate the structural maturation of contractile cells. As such, it can support the establishment of novel differentiation protocols to enhance sarcomere formation and maturity.

Objective assessment of stored blood quality by deep learning.

Stored red blood cells (RBCs) are needed for life-saving blood transfusions, but they undergo continuous degradation. RBC storage lesions are often assessed by microscopic examination or biochemical and biophysical assays, which are complex, time-consuming, and destructive to fragile cells. Here we demonstrate the use of label-free imaging flow cytometry and deep learning to characterize RBC lesions. Using brightfield images, a trained neural network achieved 76.7% agreement with experts in classifying seven clinically relevant RBC morphologies associated with storage lesions, comparable to 82.5% agreement between different experts. Given that human observation and classification may not optimally discern RBC quality, we went further and eliminated subjective human annotation in the training step by training a weakly supervised neural network using only storage duration times. The feature space extracted by this network revealed a chronological progression of morphological changes that better predicted blood quality, as measured by physiological hemolytic assay readouts, than the conventional expert-assessed morphology classification system. With further training and clinical testing across multiple sites, protocols, and instruments, deep learning and label-free imaging flow cytometry might be used to routinely and objectively assess RBC storage lesions. This would automate a complex protocol, minimize laboratory sample handling and preparation, and reduce the impact of procedural errors and discrepancies between facilities and blood donors. The chronology-based machine-learning approach may also improve upon humans' assessment of morphological changes in other biomedically important progressions, such as differentiation and metastasis.

Uncoupling protein 2 deficiency of non-cancerous tissues inhibits the progression of pancreatic cancer in mice.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is a disease of the elderly mostly because its development from preneoplastic lesions depends on the accumulation of gene mutations and epigenetic alterations over time. How aging of non-cancerous tissues of the host affects tumor progression, however, remains largely unknown. METHODS: We took advantage of a model of accelerated aging, uncoupling protein 2-deficient (Ucp2 knockout, Ucp2 KO) mice, to investigate the growth of orthotopically transplanted Ucp2 wild-type (WT) PDAC cells (cell lines Panc02 and 6606PDA) in vivo and to study strain-dependent differences of the PDAC microenvironment. RESULTS: Measurements of tumor weights and quantification of proliferating cells indicated a significant growth advantage of Panc02 and 6606PDA cells in WT mice compared to Ucp2 KO mice. In tumors in the knockout strain, higher levels of interferon-γ mRNA despite similar numbers of tumor-infiltrating T cells were observed. 6606PDA cells triggered a stronger stromal reaction in Ucp2 KO mice than in WT animals. Accordingly, pancreatic stellate cells from Ucp2 KO mice proliferated at a higher rate than cells of the WT strain when they were incubated with conditioned media from PDAC cells. CONCLUSIONS: Ucp2 modulates PDAC microenvironment in a way that favors tumor progression and implicates an altered stromal response as one of the underlying mechanisms.

Cell-Type-Specific Gene Regulatory Networks of Pro-Inflammatory and Pro-Resolving Lipid Mediator Biosynthesis in the Immune System.

Lipid mediators are important regulators in inflammatory responses, and their biosynthetic pathways are targeted by commonly used anti-inflammatory drugs. Switching from pro-inflammatory lipid mediators (PIMs) to specialized pro-resolving (SPMs) is a critical step toward acute inflammation resolution and preventing chronic inflammation. Although the biosynthetic pathways and enzymes for PIMs and SPMs have now been largely identified, the actual transcriptional profiles underlying the immune cell type-specific transcriptional profiles of these mediators are still unknown. Using the Atlas of Inflammation Resolution, we created a large network of gene regulatory interactions linked to the biosynthesis of SPMs and PIMs. By mapping single-cell sequencing data, we identified cell type-specific gene regulatory networks of the lipid mediator biosynthesis. Using machine learning approaches combined with network features, we identified cell clusters of similar transcriptional regulation and demonstrated how specific immune cell activation affects PIM and SPM profiles. We found substantial differences in regulatory networks in related cells, accounting for network-based preprocessing in functional single-cell analyses. Our results not only provide further insight into the gene regulation of lipid mediators in the immune response but also shed light on the contribution of selected cell types in their biosynthesis.

Exploring Thiazolopyridine AV25R: Unraveling of Biological Activities, Selective Anti-Cancer Properties and In Silico Target and Binding Prediction in Hematological Neoplasms.

Thiazolopyridines are a highly relevant class of small molecules, which have previously shown a wide range of biological activities. Besides their anti-tubercular, anti-microbial and anti-viral activities, they also show anti-cancerogenic properties, and play a role as inhibitors of cancer-related proteins. Herein, the biological effects of the thiazolopyridine AV25R, a novel small molecule with unknown biological effects, were characterized. Screening of a set of lymphoma (SUP-T1, SU-DHL-4) and B- acute leukemia cell lines (RS4;11, SEM) revealed highly selective effects of AV25R. The selective anti-proliferative and metabolism-modulating effects were observed in vitro for the B-ALL cell line RS4;11. Further, we were able to detect severe morphological changes and the induction of apoptosis. Gene expression analysis identified a large number of differentially expressed genes after AV25R exposure and significant differentially regulated cancer-related signaling pathways, such as VEGFA-VEGFR2 signaling and the EGF/EGFR pathway. Structure-based pharmacophore screening approaches using in silico modeling identified potential biological AV25R targets. Our results indicate that AV25R binds with several proteins known to regulate cell proliferation and tumor progression, such as FECH, MAP11, EGFR, TGFBR1 and MDM2. The molecular docking analyses indicates that AV25R has a higher binding affinity compared to many of the experimentally validated small molecule inhibitors of these targets. Thus, here we present in vitro and in silico analyses which characterize, for the first time, the molecular acting mechanism of AV25R, including cellular and molecular biologic effects. Additionally, this predicted the target binding of the molecule, revealing a high affinity to cancer-related proteins and, thus, classified AVR25 for targeted intervention approaches.

Structure-function relationships explain CTCF zinc finger mutation phenotypes in cancer.

CCCTC-binding factor (CTCF) plays fundamental roles in transcriptional regulation and chromatin architecture maintenance. CTCF is also a tumour suppressor frequently mutated in cancer, however, the structural and functional impact of mutations have not been examined. We performed molecular and structural characterisation of five cancer-specific CTCF missense zinc finger (ZF) mutations occurring within key intra- and inter-ZF residues. Functional characterisation of CTCF ZF mutations revealed a complete (L309P, R339W, R377H) or intermediate (R339Q) abrogation as well as an enhancement (G420D) of the anti-proliferative effects of CTCF. DNA binding at select sites was disrupted and transcriptional regulatory activities abrogated. Molecular docking and molecular dynamics confirmed that mutations in residues specifically contacting DNA bases or backbone exhibited loss of DNA binding. However, R339Q and G420D were stabilised by the formation of new primary DNA bonds, contributing to gain-of-function. Our data confirm that a spectrum of loss-, change- and gain-of-function impacts on CTCF zinc fingers are observed in cell growth regulation and gene regulatory activities. Hence, diverse cellular phenotypes of mutant CTCF are clearly explained by examining structure-function relationships.

The Effects of Hypoxic Preconditioned Murine Mesenchymal Stem Cells on Post-Infarct Arrhythmias in the Mouse Model.

Ventricular arrhythmias associated with myocardial infarction (MI) have a significant impact on mortality in patients following heart attack. Therefore, targeted reduction of arrhythmia represents a therapeutic approach for the prevention and treatment of severe events after infarction. Recent research transplanting mesenchymal stem cells (MSC) showed their potential in MI therapy. Our study aimed to investigate the effects of MSC injection on post-infarction arrhythmia. We used our murine double infarction model, which we previously established, to more closely mimic the clinical situation and intramyocardially injected hypoxic pre-conditioned murine MSC to the infarction border. Thereafter, various types of arrhythmias were recorded and analyzed. We observed a homogenous distribution of all types of arrhythmias after the first infarction, without any significant differences between the groups. Yet, MSC therapy after double infarction led to a highly significant reduction in simple and complex arrhythmias. Moreover, RNA-sequencing of samples from stem cell treated mice after re-infarction demonstrated a significant decline in most arrhythmias with reduced inflammatory pathways. Additionally, following stem-cell therapy we found numerous highly expressed genes to be either linked to lowering the risk of heart failure, cardiomyopathy or sudden cardiac death. Moreover, genes known to be associated with arrhythmogenesis and key mutations underlying arrhythmias were downregulated. In summary, our stem-cell therapy led to a reduction in cardiac arrhythmias after MI and showed a downregulation of already established inflammatory pathways. Furthermore, our study reveals gene regulation pathways that have a potentially direct influence on arrhythmogenesis after myocardial infarction.

Automated annotation of rare-cell types from single-cell RNA-sequencing data through synthetic oversampling.

BACKGROUND: The research landscape of single-cell and single-nuclei RNA-sequencing is evolving rapidly. In particular, the area for the detection of rare cells was highly facilitated by this technology. However, an automated, unbiased, and accurate annotation of rare subpopulations is challenging. Once rare cells are identified in one dataset, it is usually necessary to generate further specific datasets to enrich the analysis (e.g., with samples from other tissues). From a machine learning perspective, the challenge arises from the fact that rare-cell subpopulations constitute an imbalanced classification problem. We here introduce a Machine Learning (ML)-based oversampling method that uses gene expression counts of already identified rare cells as an input to generate synthetic cells to then identify similar (rare) cells in other publicly available experiments. We utilize single-cell synthetic oversampling (sc-SynO), which is based on the Localized Random Affine Shadowsampling (LoRAS) algorithm. The algorithm corrects for the overall imbalance ratio of the minority and majority class. RESULTS: We demonstrate the effectiveness of our method for three independent use cases, each consisting of already published datasets. The first use case identifies cardiac glial cells in snRNA-Seq data (17 nuclei out of 8635). This use case was designed to take a larger imbalance ratio (~1 to 500) into account and only uses single-nuclei data. The second use case was designed to jointly use snRNA-Seq data and scRNA-Seq on a lower imbalance ratio (~1 to 26) for the training step to likewise investigate the potential of the algorithm to consider both single-cell capture procedures and the impact of "less" rare-cell types. The third dataset refers to the murine data of the Allen Brain Atlas, including more than 1 million cells. For validation purposes only, all datasets have also been analyzed traditionally using common data analysis approaches, such as the Seurat workflow. CONCLUSIONS: In comparison to baseline testing without oversampling, our approach identifies rare-cells with a robust precision-recall balance, including a high accuracy and low false positive detection rate. A practical benefit of our algorithm is that it can be readily implemented in other and existing workflows. The code basis in R and Python is publicly available at FairdomHub, as well as GitHub, and can easily be transferred to identify other rare-cell types.

Harmonizing semantic annotations for computational models in biology.

Life science researchers use computational models to articulate and test hypotheses about the behavior of biological systems. Semantic annotation is a critical component for enhancing the interoperability and reusability of such models as well as for the integration of the data needed for model parameterization and validation. Encoded as machine-readable links to knowledge resource terms, semantic annotations describe the computational or biological meaning of what models and data represent. These annotations help researchers find and repurpose models, accelerate model composition and enable knowledge integration across model repositories and experimental data stores. However, realizing the potential benefits of semantic annotation requires the development of model annotation standards that adhere to a community-based annotation protocol. Without such standards, tool developers must account for a variety of annotation formats and approaches, a situation that can become prohibitively cumbersome and which can defeat the purpose of linking model elements to controlled knowledge resource terms. Currently, no consensus protocol for semantic annotation exists among the larger biological modeling community. Here, we report on the landscape of current annotation practices among the COmputational Modeling in BIology NEtwork community and provide a set of recommendations for building a consensus approach to semantic annotation.

GEMtractor: extracting views into genome-scale metabolic models.

SUMMARY: Computational metabolic models typically encode for graphs of species, reactions and enzymes. Comparing genome-scale models through topological analysis of multipartite graphs is challenging. However, in many practical cases it is not necessary to compare the full networks. The GEMtractor is a web-based tool to trim models encoded in SBML. It can be used to extract subnetworks, for example focusing on reaction- and enzyme-centric views into the model. AVAILABILITY AND IMPLEMENTATION: The GEMtractor is licensed under the terms of GPLv3 and developed at github.com/binfalse/GEMtractor-a public version is available at sbi.uni-rostock.de/gemtractor.

Protein-coding variants contribute to the risk of atopic dermatitis and skin-specific gene expression.

BACKGROUND: Fifteen percent of atopic dermatitis (AD) liability-scale heritability could be attributed to 31 susceptibility loci identified by using genome-wide association studies, with only 3 of them (IL13, IL-6 receptor [IL6R], and filaggrin [FLG]) resolved to protein-coding variants. OBJECTIVE: We examined whether a significant portion of unexplained AD heritability is further explained by low-frequency and rare variants in the gene-coding sequence. METHODS: We evaluated common, low-frequency, and rare protein-coding variants using exome chip and replication genotype data of 15,574 patients and 377,839 control subjects combined with whole-transcriptome data on lesional, nonlesional, and healthy skin samples of 27 patients and 38 control subjects. RESULTS: An additional 12.56% (SE, 0.74%) of AD heritability is explained by rare protein-coding variation. We identified docking protein 2 (DOK2) and CD200 receptor 1 (CD200R1) as novel genome-wide significant susceptibility genes. Rare coding variants associated with AD are further enriched in 5 genes (IL-4 receptor [IL4R], IL13, Janus kinase 1 [JAK1], JAK2, and tyrosine kinase 2 [TYK2]) of the IL13 pathway, all of which are targets for novel systemic AD therapeutics. Multiomics-based network and RNA sequencing analysis revealed DOK2 as a central hub interacting with, among others, CD200R1, IL6R, and signal transducer and activator of transcription 3 (STAT3). Multitissue gene expression profile analysis for 53 tissue types from the Genotype-Tissue Expression project showed that disease-associated protein-coding variants exert their greatest effect in skin tissues. CONCLUSION: Our discoveries highlight a major role of rare coding variants in AD acting independently of common variants. Further extensive functional studies are required to detect all potential causal variants and to specify the contribution of the novel susceptibility genes DOK2 and CD200R1 to overall disease susceptibility.

An integrative network-driven pipeline for systematic identification of lncRNA-associated regulatory network motifs in metastatic melanoma.

BACKGROUND: Melanoma phenotype and the dynamics underlying its progression are determined by a complex interplay between different types of regulatory molecules. In particular, transcription factors (TFs), microRNAs (miRNAs), and long non-coding RNAs (lncRNAs) interact in layers that coalesce into large molecular interaction networks. Our goal here is to study molecules associated with the cross-talk between various network layers, and their impact on tumor progression. RESULTS: To elucidate their contribution to disease, we developed an integrative computational pipeline to construct and analyze a melanoma network focusing on lncRNAs, their miRNA and protein targets, miRNA target genes, and TFs regulating miRNAs. In the network, we identified three-node regulatory loops each composed of lncRNA, miRNA, and TF. To prioritize these motifs for their role in melanoma progression, we integrated patient-derived RNAseq dataset from TCGA (SKCM) melanoma cohort, using a weighted multi-objective function. We investigated the expression profile of the top-ranked motifs and used them to classify patients into metastatic and non-metastatic phenotypes. CONCLUSIONS: The results of this study showed that network motif UCA1/AKT1/hsa-miR-125b-1 has the highest prediction accuracy (ACC = 0.88) for discriminating metastatic and non-metastatic melanoma phenotypes. The observation is also confirmed by the progression-free survival analysis where the patient group characterized by the metastatic-type expression profile of the motif suffers a significant reduction in survival. The finding suggests a prognostic value of network motifs for the classification and treatment of melanoma.

Label-Free Leukemia Monitoring by Computer Vision.

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. While there are a number of well-recognized prognostic biomarkers at diagnosis, the most powerful independent prognostic factor is the response of the leukemia to induction chemotherapy (Campana and Pui: Blood 129 (2017) 1913-1918). Given the potential for machine learning to improve precision medicine, we tested its capacity to monitor disease in children undergoing ALL treatment. Diagnostic and on-treatment bone marrow samples were labeled with an ALL-discriminating antibody combination and analyzed by imaging flow cytometry. Ignoring the fluorescent markers and using only features extracted from bright-field and dark-field cell images, a deep learning model was able to identify ALL cells at an accuracy of >88%. This antibody-free, single cell method is cheap, quick, and could be adapted to a simple, laser-free cytometer to allow automated, point-of-care testing to detect slow early responders. Adaptation to other types of leukemia is feasible, which would revolutionize residual disease monitoring. © 2020 The Authors. Cytometry Part A published by Wiley Periodicals, Inc. on behalf of International Society for Advancement of Cytometry.

The role of cooperativity in a p53-miR34 dynamical mathematical model.

The objective of this study is to evaluate the role of cooperativity, captured by the Hill coefficient, in a minimal mathematical model describing the interactions between p53 and miR-34a. The model equations are analyzed for negative, none and normal cooperativity using a specific version of bifurcation theory and they are solved numerically. Special attention is paid to the sign of so-called first Lyapunov value. Interpretations of the results are given, both according to dynamic theory and in biological terms. In terms of cell signaling, we propose the hypothesis that when the outgoing signal of a system spends a physiologically significant amount of time outside of its equilibrium state, then the value of that signal can be sampled at any point along the trajectory towards that equilibrium and indeed, at multiple points. Coupled with non-linear behavior, such as that caused by cooperativity, this feature can account for a complex and varied response, which p53 is known for. From dynamical point of view, we found that when cooperativity is negative, the system has only one stable equilibrium point. In the absence of cooperativity, there is a single unstable equilibrium point with a critical boundary of stability. In the case with normal cooperativity, the system can have one, two, or three steady states with both, bi-stability and bi-instability occurring.

Label-Free Identification of White Blood Cells Using Machine Learning.

White blood cell (WBC) differential counting is an established clinical routine to assess patient immune system status. Fluorescent markers and a flow cytometer are required for the current state-of-the-art method for determining WBC differential counts. However, this process requires several sample preparation steps and may adversely disturb the cells. We present a novel label-free approach using an imaging flow cytometer and machine learning algorithms, where live, unstained WBCs were classified. It achieved an average F1-score of 97% and two subtypes of WBCs, B and T lymphocytes, were distinguished from each other with an average F1-score of 78%, a task previously considered impossible for unlabeled samples. We provide an open-source workflow to carry out the procedure. We validated the WBC analysis with unstained samples from 85 donors. The presented method enables robust and highly accurate identification of WBCs, minimizing the disturbance to the cells and leaving marker channels free to answer other biological questions. It also opens the door to employing machine learning for liquid biopsy, here, using the rich information in cell morphology for a wide range of diagnostics of primary blood. © 2019 The Authors. Cytometry Part A published by Wiley Periodicals, Inc. on behalf of International Society for Advancement of Cytometry.

The RNA workbench: best practices for RNA and high-throughput sequencing bioinformatics in Galaxy.

RNA-based regulation has become a major research topic in molecular biology. The analysis of epigenetic and expression data is therefore incomplete if RNA-based regulation is not taken into account. Thus, it is increasingly important but not yet standard to combine RNA-centric data and analysis tools with other types of experimental data such as RNA-seq or ChIP-seq. Here, we present the RNA workbench, a comprehensive set of analysis tools and consolidated workflows that enable the researcher to combine these two worlds. Based on the Galaxy framework the workbench guarantees simple access, easy extension, flexible adaption to personal and security needs, and sophisticated analyses that are independent of command-line knowledge. Currently, it includes more than 50 bioinformatics tools that are dedicated to different research areas of RNA biology including RNA structure analysis, RNA alignment, RNA annotation, RNA-protein interaction, ribosome profiling, RNA-seq analysis and RNA target prediction. The workbench is developed and maintained by experts in RNA bioinformatics and the Galaxy framework. Together with the growing community evolving around this workbench, we are committed to keep the workbench up-to-date for future standards and needs, providing researchers with a reliable and robust framework for RNA data analysis. AVAILABILITY: The RNA workbench is available at https://github.com/bgruening/galaxy-rna-workbench.

Unveiling network-based functional features through integration of gene expression into protein networks.

Decoding health and disease phenotypes is one of the fundamental objectives in biomedicine. Whereas high-throughput omics approaches are available, it is evident that any single omics approach might not be adequate to capture the complexity of phenotypes. Therefore, integrated multi-omics approaches have been used to unravel genotype-phenotype relationships such as global regulatory mechanisms and complex metabolic networks in different eukaryotic organisms. Some of the progress and challenges associated with integrated omics studies have been reviewed previously in comprehensive studies. In this work, we highlight and review the progress, challenges and advantages associated with emerging approaches, integrating gene expression and protein-protein interaction networks to unravel network-based functional features. This includes identifying disease related genes, gene prioritization, clustering protein interactions, developing the modules, extract active subnetworks and static protein complexes or dynamic/temporal protein complexes. We also discuss how these approaches contribute to our understanding of the biology of complex traits and diseases. This article is part of a Special Issue entitled: Cardiac adaptations to obesity, diabetes and insulin resistance, edited by Professors Jan F.C. Glatz, Jason R.B. Dyck and Christine Des Rosiers.

A web platform for the network analysis of high-throughput data in melanoma and its use to investigate mechanisms of resistance to anti-PD1 immunotherapy.

Cellular phenotypes are established and controlled by complex and precisely orchestrated molecular networks. In cancer, mutations and dysregulations of multiple molecular factors perturb the regulation of these networks and lead to malignant transformation. High-throughput technologies are a valuable source of information to establish the complex molecular relationships behind the emergence of malignancy, but full exploitation of this massive amount of data requires bioinformatics tools that rely on network-based analyses. In this report we present the Virtual Melanoma Cell, an online tool developed to facilitate the mining and interpretation of high-throughput data on melanoma by biomedical researches. The platform is based on a comprehensive, manually generated and expert-validated regulatory map composed of signaling pathways important in malignant melanoma. The Virtual Melanoma Cell is a tool designed to accept, visualize and analyze user-generated datasets. It is available at: https://www.vcells.net/melanoma. To illustrate the utilization of the web platform and the regulatory map, we have analyzed a large publicly available dataset accounting for anti-PD1 immunotherapy treatment of malignant melanoma patients.

The RNA world in the 21st century-a systems approach to finding non-coding keys to clinical questions.

There was evidence that RNAs are a functionally rich class of molecules not only since the arrival of the next-generation sequencing technology. Non-coding RNAs (ncRNA) could be the key to accelerated diagnosis and enhanced prediction of disease and therapy outcomes as well as the design of advanced therapeutic strategies to overcome yet unsatisfactory approaches.In this review, we discuss the state of the art in RNA systems biology with focus on the application in the systems biomedicine field. We propose guidelines for analysing the role of microRNAs and long non-coding RNAs in human pathologies. We introduce RNA expression profiling and network approaches for the identification of stable and effective RNomics-based biomarkers, providing insights into the role of ncRNAs in disease regulation. Towards this, we discuss ways to model the dynamics of gene regulatory networks and signalling pathways that involve ncRNAs. We also describe data resources and computational methods for finding putative mechanisms of action of ncRNAs. Finally, we discuss avenues for the computer-aided design of novel RNA-based therapeutics.