Immune-related toxicity and soluble profile in patients affected by solid tumors: a network approach.

BACKGROUND: Immune checkpoint inhibitors (ICIs) have particular, immune-related adverse events (irAEs), as a consequence of interfering with self-tolerance mechanisms. The incidence of irAEs varies depending on ICI class, administered dose and treatment schedule. The aim of this study was to define a baseline (T0) immune profile (IP) predictive of irAE development. METHODS: A prospective, multicenter study evaluating the immune profile (IP) of 79 patients with advanced cancer and treated with anti-programmed cell death protein 1 (anti-PD-1) drugs as a first- or second-line setting was performed. The results were then correlated with irAEs onset. The IP was studied by means of multiplex assay, evaluating circulating concentration of 12 cytokines, 5 chemokines, 13 soluble immune checkpoints and 3 adhesion molecules. Indoleamine 2, 3-dioxygenase (IDO) activity was measured through a modified liquid chromatography-tandem mass spectrometry using the high-performance liquid chromatography-mass spectrometry (HPLC-MS/MS) method. A connectivity heatmap was obtained by calculating Spearman correlation coefficients. Two different networks of connectivity were constructed, based on the toxicity profile. RESULTS: Toxicity was predominantly of low/moderate grade. High-grade irAEs were relatively rare, while cumulative toxicity was high (35%). Positive and statistically significant correlations between the cumulative toxicity and IP10 and IL8, sLAG3, sPD-L2, sHVEM, sCD137, sCD27 and sICAM-1 serum concentration were found. Moreover, patients who experienced irAEs had a markedly different connectivity pattern, characterized by disruption of most of the paired connections between cytokines, chemokines and connections of sCD137, sCD27 and sCD28, while sPDL-2 pair-wise connectivity values seemed to be intensified. Network connectivity analysis identified a total of 187 statistically significant interactions in patients without toxicity and a total of 126 statistically significant interactions in patients with toxicity. Ninety-eight interactions were common to both networks, while 29 were specifically observed in patients who experienced toxicity. CONCLUSIONS: A particular, common pattern of immune dysregulation was defined in patients developing irAEs. This immune serological profile, if confirmed in a larger patient population, could lead to the design of a personalized therapeutic strategy in order to prevent, monitor and treat irAEs at an early stage.

SWIMmeR: an R-based software to unveiling crucial nodes in complex biological networks.

SUMMARY: We present SWIMmeR, an open-source version of its predecessor SWIM (SWitchMiner) that is a network-based tool for mining key (switch) genes that are associated with intriguing patterns of molecular co-abundance and may play a crucial role in phenotypic transitions in various biological settings. SWIM was originally written in MATLAB®, a proprietary programming language that requires the purchase of a license to install, manipulate, operate and run the software. Over the last years, SWIM has sparked a widespread interest within the scientific community thanks to the promising results obtained through its application in a broad range of phenotype-specific scenarios, spanning from complex diseases to grapevine berry maturation. This success has created the call for it to be distributed in a freely accessible, open-source, runtime environment, such as R, aimed at a general audience of non-expert users that cannot afford the leading proprietary solution. SWIMmeR is provided as a comprehensive collection of R functions and it also includes several additional features that make it less intensive in terms of computer time and more efficient in terms of usability and further implementation and extension. AVAILABILITY AND IMPLEMENTATION: The SWIMmeR source code is freely available at https://github.com/sportingCode/SWIMmeR.git, along with a practical user guide, including a usage example of its application on breast cancer dataset. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An ontology-based approach for modelling and querying Alzheimer's disease data.

BACKGROUND: The recent advances in biotechnology and computer science have led to an ever-increasing availability of public biomedical data distributed in large databases worldwide. However, these data collections are far from being "standardized" so to be harmonized or even integrated, making it impossible to fully exploit the latest machine learning technologies for the analysis of data themselves. Hence, facing this huge flow of biomedical data is a challenging task for researchers and clinicians due to their complexity and high heterogeneity. This is the case of neurodegenerative diseases and the Alzheimer's Disease (AD) in whose context specialized data collections such as the one by the Alzheimer's Disease Neuroimaging Initiative (ADNI) are maintained. METHODS: Ontologies are controlled vocabularies that allow the semantics of data and their relationships in a given domain to be represented. They are often exploited to aid knowledge and data management in healthcare research. Computational Ontologies are the result of the combination of data management systems and traditional ontologies. Our approach is i) to define a computational ontology representing a logic-based formal conceptual model of the ADNI data collection and ii) to provide a means for populating the ontology with the actual data in the Alzheimer Disease Neuroimaging Initiative (ADNI). These two components make it possible to semantically query the ADNI database in order to support data extraction in a more intuitive manner. RESULTS: We developed: i) a detailed computational ontology for clinical multimodal datasets from the ADNI repository in order to simplify the access to these data; ii) a means for populating this ontology with the actual ADNI data. Such computational ontology immediately makes it possible to facilitate complex queries to the ADNI files, obtaining new diagnostic knowledge about Alzheimer's disease. CONCLUSIONS: The proposed ontology will improve the access to the ADNI dataset, allowing queries to extract multivariate datasets to perform multidimensional and longitudinal statistical analyses. Moreover, the proposed ontology can be a candidate for supporting the design and implementation of new information systems for the collection and management of AD data and metadata, and for being a reference point for harmonizing or integrating data residing in different sources.

Hepatitis B Virus and microRNAs: A Bioinformatics Approach.

In recent decades, microRNAs (miRNAs) have emerged as key regulators of gene expression, and the identification of viral miRNAs (v-miRNAs) within some viruses, including hepatitis B virus (HBV), has attracted significant attention. HBV infections often progress to chronic states (CHB) and may induce fibrosis/cirrhosis and hepatocellular carcinoma (HCC). The presence of HBV can dysregulate host miRNA expression, influencing several biological pathways, such as apoptosis, innate and immune response, viral replication, and pathogenesis. Consequently, miRNAs are considered a promising biomarker for diagnostic, prognostic, and treatment response. The dynamics of miRNAs during HBV infection are multifaceted, influenced by host variability and miRNA interactions. Given the ability of miRNAs to target multiple messenger RNA (mRNA), understanding the viral-host (human) interplay is complex but essential to develop novel clinical applications. Therefore, bioinformatics can help to analyze, identify, and interpret a vast amount of miRNA data. This review explores the bioinformatics tools available for viral and host miRNA research. Moreover, we introduce a brief overview focusing on the role of miRNAs during HBV infection. In this way, this review aims to help the selection of the most appropriate bioinformatics tools based on requirements and research goals.

SPINNAKER: an R-based tool to highlight key RNA interactions in complex biological networks.

BACKGROUND: Recently, we developed a mathematical model for identifying putative competing endogenous RNA (ceRNA) interactions. This methodology has aroused a broad acknowledgment within the scientific community thanks to the encouraging results achieved when applied to breast invasive carcinoma, leading to the identification of PVT1, a long non-coding RNA functioning as ceRNA for the miR-200 family. The main shortcoming of the model is that it is no freely available and implemented in MATLAB®, a proprietary programming platform requiring a paid license for installing, operating, manipulating, and running the software. RESULTS: Breaking through these model limitations demands to distribute it in an open-source, freely accessible environment, such as R, designed for an ordinary audience of users that are not able to afford a proprietary solution. Here, we present SPINNAKER (SPongeINteractionNetworkmAKER), the open-source version of our widely established mathematical model for predicting ceRNAs crosstalk, that is released as an exhaustive collection of R functions. SPINNAKER has been even designed for providing many additional features that facilitate its usability, make it more efficient in terms of further implementation and extension, and less intense in terms of computational execution time. CONCLUSIONS: SPINNAKER source code is freely available at https://github.com/sportingCode/SPINNAKER.git together with a thoroughgoing PPT-based guideline. In order to help users get the key points more conveniently, also a practical R-styled plain-text guideline is provided. Finally, a short movie is available to help the user to set the own directory, properly.

SAveRUNNER: A network-based algorithm for drug repurposing and its application to COVID-19.

The novelty of new human coronavirus COVID-19/SARS-CoV-2 and the lack of effective drugs and vaccines gave rise to a wide variety of strategies employed to fight this worldwide pandemic. Many of these strategies rely on the repositioning of existing drugs that could shorten the time and reduce the cost compared to de novo drug discovery. In this study, we presented a new network-based algorithm for drug repositioning, called SAveRUNNER (Searching off-lAbel dRUg aNd NEtwoRk), which predicts drug-disease associations by quantifying the interplay between the drug targets and the disease-specific proteins in the human interactome via a novel network-based similarity measure that prioritizes associations between drugs and diseases locating in the same network neighborhoods. Specifically, we applied SAveRUNNER on a panel of 14 selected diseases with a consolidated knowledge about their disease-causing genes and that have been found to be related to COVID-19 for genetic similarity (i.e., SARS), comorbidity (e.g., cardiovascular diseases), or for their association to drugs tentatively repurposed to treat COVID-19 (e.g., malaria, HIV, rheumatoid arthritis). Focusing specifically on SARS subnetwork, we identified 282 repurposable drugs, including some the most rumored off-label drugs for COVID-19 treatments (e.g., chloroquine, hydroxychloroquine, tocilizumab, heparin), as well as a new combination therapy of 5 drugs (hydroxychloroquine, chloroquine, lopinavir, ritonavir, remdesivir), actually used in clinical practice. Furthermore, to maximize the efficiency of putative downstream validation experiments, we prioritized 24 potential anti-SARS-CoV repurposable drugs based on their network-based similarity values. These top-ranked drugs include ACE-inhibitors, monoclonal antibodies (e.g., anti-IFNγ, anti-TNFα, anti-IL12, anti-IL1ß, anti-IL6), and thrombin inhibitors. Finally, our findings were in-silico validated by performing a gene set enrichment analysis, which confirmed that most of the network-predicted repurposable drugs may have a potential treatment effect against human coronavirus infections.

A Comparison of Network-Based Methods for Drug Repurposing along with an Application to Human Complex Diseases.

Drug repurposing strategy, proposing a therapeutic switching of already approved drugs with known medical indications to new therapeutic purposes, has been considered as an efficient approach to unveil novel drug candidates with new pharmacological activities, significantly reducing the cost and shortening the time of de novo drug discovery. Meaningful computational approaches for drug repurposing exploit the principles of the emerging field of Network Medicine, according to which human diseases can be interpreted as local perturbations of the human interactome network, where the molecular determinants of each disease (disease genes) are not randomly scattered, but co-localized in highly interconnected subnetworks (disease modules), whose perturbation is linked to the pathophenotype manifestation. By interpreting drug effects as local perturbations of the interactome, for a drug to be on-target effective against a specific disease or to cause off-target adverse effects, its targets should be in the nearby of disease-associated genes. Here, we used the network-based proximity measure to compute the distance between the drug module and the disease module in the human interactome by exploiting five different metrics (minimum, maximum, mean, median, mode), with the aim to compare different frameworks for highlighting putative repurposable drugs to treat complex human diseases, including malignant breast and prostate neoplasms, schizophrenia, and liver cirrhosis. Whilst the standard metric (that is the minimum) for the network-based proximity remained a valid tool for efficiently screening off-label drugs, we observed that the other implemented metrics specifically predicted further interesting drug candidates worthy of investigation for yielding a potentially significant clinical benefit.

Repurposing Histaminergic Drugs in Multiple Sclerosis.

Multiple sclerosis is an autoimmune disease with a strong neuroinflammatory component that contributes to severe demyelination, neurodegeneration and lesions formation in white and grey matter of the spinal cord and brain. Increasing attention is being paid to the signaling of the biogenic amine histamine in the context of several pathological conditions. In multiple sclerosis, histamine regulates the differentiation of oligodendrocyte precursors, reduces demyelination, and improves the remyelination process. However, the concomitant activation of histamine H1-H4 receptors can sustain either damaging or favorable effects, depending on the specifically activated receptor subtype/s, the timing of receptor engagement, and the central versus peripheral target district. Conventional drug development has failed so far to identify curative drugs for multiple sclerosis, thus causing a severe delay in therapeutic options available to patients. In this perspective, drug repurposing offers an exciting and complementary alternative for rapidly approving some medicines already approved for other indications. In the present work, we have adopted a new network-medicine-based algorithm for drug repurposing called SAveRUNNER, for quantifying the interplay between multiple sclerosis-associated genes and drug targets in the human interactome. We have identified new histamine drug-disease associations and predicted off-label novel use of the histaminergic drugs amodiaquine, rupatadine, and diphenhydramine among others, for multiple sclerosis. Our work suggests that selected histamine-related molecules might get to the root causes of multiple sclerosis and emerge as new potential therapeutic strategies for the disease.

SAveRUNNER: an R-based tool for drug repurposing.

BACKGROUND: Currently, no proven effective drugs for the novel coronavirus disease COVID-19 exist and despite widespread vaccination campaigns, we are far short from herd immunity. The number of people who are still vulnerable to the virus is too high to hamper new outbreaks, leading a compelling need to find new therapeutic options devoted to combat SARS-CoV-2 infection. Drug repurposing represents an effective drug discovery strategy from existing drugs that could shorten the time and reduce the cost compared to de novo drug discovery. RESULTS: We developed a network-based tool for drug repurposing provided as a freely available R-code, called SAveRUNNER (Searching off-lAbel dRUg aNd NEtwoRk), with the aim to offer a promising framework to efficiently detect putative novel indications for currently marketed drugs against diseases of interest. SAveRUNNER predicts drug-disease associations by quantifying the interplay between the drug targets and the disease-associated proteins in the human interactome through the computation of a novel network-based similarity measure, which prioritizes associations between drugs and diseases located in the same network neighborhoods. CONCLUSIONS: The algorithm was successfully applied to predict off-label drugs to be repositioned against the new human coronavirus (2019-nCoV/SARS-CoV-2), and it achieved a high accuracy in the identification of well-known drug indications, thus revealing itself as a powerful tool to rapidly detect potential novel medical indications for various drugs that are worth of further investigation. SAveRUNNER source code is freely available at https://github.com/giuliafiscon/SAveRUNNER.git , along with a comprehensive user guide.

The New Paradigm of Network Medicine to Analyze Breast Cancer Phenotypes.

Breast cancer (BC) is a heterogeneous and complex disease as witnessed by the existence of different subtypes and clinical characteristics that poses significant challenges in disease management. The complexity of this tumor may rely on the highly interconnected nature of the various biological processes as stated by the new paradigm of Network Medicine. We explored The Cancer Genome Atlas (TCGA)-BRCA data set, by applying the network-based algorithm named SWItch Miner, and mapping the findings on the human interactome to capture the molecular interconnections associated with the disease modules. To characterize BC phenotypes, we constructed protein-protein interaction modules based on "hub genes", called switch genes, both common and specific to the four tumor subtypes. Transcriptomic profiles of patients were stratified according to both clinical (immunohistochemistry) and genetic (PAM50) classifications. 266 and 372 switch genes were identified from immunohistochemistry and PAM50 classifications, respectively. Moreover, the identified switch genes were functionally characterized to select an interconnected pathway of disease genes. By intersecting the common switch genes of the two classifications, we selected a unique signature of 28 disease genes that were BC subtype-independent and classification subtype-independent. Data were validated both in vitro (10 BC cell lines) and ex vivo (66 BC tissues) experiments. Results showed that four of these hub proteins (AURKA, CDC45, ESPL1, and RAD54L) were over-expressed in all tumor subtypes. Moreover, the inhibition of one of the identified switch genes (AURKA) similarly affected all BC subtypes. In conclusion, using a network-based approach, we identified a common BC disease module which might reflect its pathological signature, suggesting a new vision to face with the disease heterogeneity.

MIENTURNET: an interactive web tool for microRNA-target enrichment and network-based analysis.

BACKGROUND: miRNAs regulate the expression of several genes with one miRNA able to target multiple genes and with one gene able to be simultaneously targeted by more than one miRNA. Therefore, it has become indispensable to shorten the long list of miRNA-target interactions to put in the spotlight in order to gain insight into understanding the regulatory mechanism orchestrated by miRNAs in various cellular processes. A reasonable solution is certainly to prioritize miRNA-target interactions to maximize the effectiveness of the downstream analysis. RESULTS: We propose a new and easy-to-use web tool MIENTURNET (MicroRNA ENrichment TURned NETwork) that receives in input a list of miRNAs or mRNAs and tackles the problem of prioritizing miRNA-target interactions by performing a statistical analysis followed by a fully featured network-based visualization and analysis. The statistics is used to assess the significance of an over-representation of miRNA-target interactions and then MIENTURNET filters based on the statistical significance associated with each miRNA-target interaction. In addition, the holistic approach of the network theory is used to infer possible evidences of miRNA regulation by capturing emergent properties of the miRNA-target regulatory network that would be not evident through a pairwise analysis of the individual components. CONCLUSION: MIENTURNET offers the possibility to consistently perform both statistical and network-based analyses by using only a single tool leading to a more effective prioritization of the miRNA-target interactions. This has the potential to avoid researchers without computational and informatics skills to navigate multiple websites and thus to independently investigate miRNA activity in every cellular process of interest in an easy and at the same time exhaustive way thanks to the intuitive web interface. The web application along with a well-documented and comprehensive user guide are freely available at http://userver.bio.uniroma1.it/apps/mienturnet/ without any login requirement.

SWIM tool application to expression data of glioblastoma stem-like cell lines, corresponding primary tumors and conventional glioma cell lines.

BACKGROUND: It is well-known that glioblastoma contains self-renewing, stem-like subpopulation with the ability to sustain tumor growth. These cells - called cancer stem-like cells - share certain phenotypic characteristics with untransformed stem cells and are resistant to many conventional cancer therapies, which might explain the limitations in curing human malignancies. Thus, the identification of genes controlling the differentiation of these stem-like cells is becoming a successful therapeutic strategy, owing to the promise of novel targets for treating malignancies. METHODS: Recently, we developed SWIM, a software able to unveil a small pool of genes - called switch genes - critically associated with drastic changes in cell phenotype. Here, we applied SWIM to the expression profiling of glioblastoma stem-like cells and conventional glioma cell lines, in order to identify switch genes related to stem-like phenotype. RESULTS: SWIM identifies 171 switch genes that are all down-regulated in glioblastoma stem-like cells. This list encompasses genes like CAV1, COL5A1, COL6A3, FLNB, HMMR, ITGA3, ITGA5, MET, SDC1, THBS1, and VEGFC, involved in "ECM-receptor interaction" and "focal adhesion" pathways. The inhibition of switch genes highly correlates with the activation of genes related to neural development and differentiation, such as the 4-core OLIG2, POU3F2, SALL2, SOX2, whose induction has been shown to be sufficient to reprogram differentiated glioblastoma into stem-like cells. Among switch genes, the transcription factor FOSL1 appears as the brightest star since: it is down-regulated in stem-like cells; it highly negatively correlates with the 4-core genes that are all up-regulated in stem-like cells; the promoter regions of the 4-core genes harbor a consensus binding motif for FOSL1. CONCLUSIONS: We suggest that the inhibition of switch genes in stem-like cells could induce the deregulation of cell communication pathways, contributing to neoplastic progression and tumor invasiveness. Conversely, their activation could restore the physiological equilibrium between cell adhesion and migration, hampering the progression of cancer. Moreover, we posit FOSL1 as promising candidate to orchestrate the differentiation of cancer stem-like cells by repressing the 4-core genes' expression, which severely halts cancer growth and might affect the therapeutic outcome. We suggest FOSL1 as novel putative therapeutic and prognostic biomarker, worthy of further investigation.

Combining EEG signal processing with supervised methods for Alzheimer's patients classification.

BACKGROUND: Alzheimer's Disease (AD) is a neurodegenaritive disorder characterized by a progressive dementia, for which actually no cure is known. An early detection of patients affected by AD can be obtained by analyzing their electroencephalography (EEG) signals, which show a reduction of the complexity, a perturbation of the synchrony, and a slowing down of the rhythms. METHODS: In this work, we apply a procedure that exploits feature extraction and classification techniques to EEG signals, whose aim is to distinguish patient affected by AD from the ones affected by Mild Cognitive Impairment (MCI) and healthy control (HC) samples. Specifically, we perform a time-frequency analysis by applying both the Fourier and Wavelet Transforms on 109 samples belonging to AD, MCI, and HC classes. The classification procedure is designed with the following steps: (i) preprocessing of EEG signals; (ii) feature extraction by means of the Discrete Fourier and Wavelet Transforms; and (iii) classification with tree-based supervised methods. RESULTS: By applying our procedure, we are able to extract reliable human-interpretable classification models that allow to automatically assign the patients into their belonging class. In particular, by exploiting a Wavelet feature extraction we achieve 83%, 92%, and 79% of accuracy when dealing with HC vs AD, HC vs MCI, and MCI vs AD classification problems, respectively. CONCLUSIONS: Finally, by comparing the classification performances with both feature extraction methods, we find out that Wavelets analysis outperforms Fourier. Hence, we suggest it in combination with supervised methods for automatic patients classification based on their EEG signals for aiding the medical diagnosis of dementia.

TCGA2BED: extracting, extending, integrating, and querying The Cancer Genome Atlas.

BACKGROUND: Data extraction and integration methods are becoming essential to effectively access and take advantage of the huge amounts of heterogeneous genomics and clinical data increasingly available. In this work, we focus on The Cancer Genome Atlas, a comprehensive archive of tumoral data containing the results of high-throughout experiments, mainly Next Generation Sequencing, for more than 30 cancer types. RESULTS: We propose TCGA2BED a software tool to search and retrieve TCGA data, and convert them in the structured BED format for their seamless use and integration. Additionally, it supports the conversion in CSV, GTF, JSON, and XML standard formats. Furthermore, TCGA2BED extends TCGA data with information extracted from other genomic databases (i.e., NCBI Entrez Gene, HGNC, UCSC, and miRBase). We also provide and maintain an automatically updated data repository with publicly available Copy Number Variation, DNA-methylation, DNA-seq, miRNA-seq, and RNA-seq (V1,V2) experimental data of TCGA converted into the BED format, and their associated clinical and biospecimen meta data in attribute-value text format. CONCLUSIONS: The availability of the valuable TCGA data in BED format reduces the time spent in taking advantage of them: it is possible to efficiently and effectively deal with huge amounts of cancer genomic data integratively, and to search, retrieve and extend them with additional information. The BED format facilitates the investigators allowing several knowledge discovery analyses on all tumor types in TCGA with the final aim of understanding pathological mechanisms and aiding cancer treatments.

CAMUR: Knowledge extraction from RNA-seq cancer data through equivalent classification rules.

MOTIVATION: Nowadays, knowledge extraction methods from Next Generation Sequencing data are highly requested. In this work, we focus on RNA-seq gene expression analysis and specifically on case-control studies with rule-based supervised classification algorithms that build a model able to discriminate cases from controls. State of the art algorithms compute a single classification model that contains few features (genes). On the contrary, our goal is to elicit a higher amount of knowledge by computing many classification models, and therefore to identify most of the genes related to the predicted class. RESULTS: We propose CAMUR, a new method that extracts multiple and equivalent classification models. CAMUR iteratively computes a rule-based classification model, calculates the power set of the genes present in the rules, iteratively eliminates those combinations from the data set, and performs again the classification procedure until a stopping criterion is verified. CAMUR includes an ad-hoc knowledge repository (database) and a querying tool.We analyze three different types of RNA-seq data sets (Breast, Head and Neck, and Stomach Cancer) from The Cancer Genome Atlas (TCGA) and we validate CAMUR and its models also on non-TCGA data. Our experimental results show the efficacy of CAMUR: we obtain several reliable equivalent classification models, from which the most frequent genes, their relationships, and the relation with a particular cancer are deduced. AVAILABILITY AND IMPLEMENTATION: dmb.iasi.cnr.it/camur.php CONTACT: emanuel@iasi.cnr.it SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MONSTER v1.1: a tool to extract and search for RNA non-branching structures.

BACKGROUND: Detection of RNA structure similarities is still one of the major computational problems in the discovery of RNA functions. A case in point is the study of the new appreciated long non-coding RNAs (lncRNAs), emerging as new players involved in many cellular processes and molecular interactions. Among several mechanisms of action, some lncRNAs show specific substructures that are likely to be instrumental for their functioning. For instance, it has been reported in literature that some lncRNAs have a guiding or scaffolding role by binding chromatin-modifying protein complexes. Thus, a functionally characterized lncRNA (reference) can be used to infer the function of others that are functionally unknown (target), based on shared structural motifs. METHODS: In our previous work we presented a tool, MONSTER v1.0, able to identify structural motifs shared between two full-length RNAs. Our procedure is mainly composed of two ad-hoc developed algorithms: nbRSSP_extractor for characterizing the folding of an RNA sequence by means of a sequence-structure descriptor (i.e., an array of non-overlapping substructures located on the RNA sequence and coded by dot-bracket notation); and SSD_finder, to enable an effective search engine for groups of matches (i.e., chains) common to the reference and target RNA based on a dynamic programming approach with a new score function. Here, we present an updated version of the previous one (MONSTER v1.1) accounting for the peculiar feature of lncRNAs that are not expected to have a unique fold, but appear to fluctuate among a large number of equally-stable folds. In particular, we improved our SSD_finder algorithm in order to take into account all the alternative equally-stable structures. RESULTS: We present an application of MONSTER v1.1 on lincRNAs, which are a specific class of lncRNAs located in genomic regions which do not overlap protein-coding genes. In particular, we provide reliable predictions of the shared chains between HOTAIR, ANRIL and COLDAIR. The latter are lincRNAs which interact with the same protein complexes of the Polycomb group and hence they are expected to share structural motifs.

Circular RNA mediated gene regulation in human breast cancer: A bioinformatics analysis.

Circular RNAs (circRNAs) are a new acknowledged class of RNAs that has been shown to play a major role in several biological functions both in physiological and pathological conditions, operating as critical part of regulatory processes, like competing endogenous RNA (ceRNA) networks. The ceRNA hypothesis is a recently discovered molecular mechanism that adds a new key layer of post-transcriptional regulation, whereby various types of RNAs can reciprocally influence each other's expression competing for binding the same pool of microRNAs, even affecting disease development. In this study, we build a network of circRNA-miRNA-mRNA interactions in human breast cancer, called CERNOMA, that is a bipartite graph with one class of nodes corresponding to differentially expressed miRNAs (DEMs) and the other one corresponding to differentially expressed circRNAs (DEC) and mRNAs (DEGs). A link between a DEC (or DEG) and DEM is placed if it is predicted to be a target of the DEM and shows an opposite expression level trend with respect to the DEM. Within the CERNOMA, we highlighted an interesting deregulated circRNA-miRNA-mRNA triplet, including the up-regulated hsa_circRNA_102908 (BRCA1 associated RING domain 1), the down-regulated miR-410-3p, and the up-regulated ESM1, whose overexpression has been already shown to promote tumor dissemination and metastasis in breast cancer.

A Network of MicroRNAs and mRNAs Involved in Melanosome Maturation and Trafficking Defines the Lower Response of Pigmentable Melanoma Cells to Targeted Therapy.

BACKGROUND: The ability to increase their degree of pigmentation is an adaptive response that confers pigmentable melanoma cells higher resistance to BRAF inhibitors (BRAFi) compared to non-pigmentable melanoma cells. METHODS: Here, we compared the miRNome and the transcriptome profile of pigmentable 501Mel and SK-Mel-5 melanoma cells vs. non-pigmentable A375 melanoma cells, following treatment with the BRAFi vemurafenib (vem). In depth bioinformatic analyses (clusterProfiler, WGCNA and SWIMmeR) allowed us to identify the miRNAs, mRNAs and biological processes (BPs) that specifically characterize the response of pigmentable melanoma cells to the drug. Such BPs were studied using appropriate assays in vitro and in vivo (xenograft in zebrafish embryos). RESULTS: Upon vem treatment, miR-192-5p, miR-211-5p, miR-374a-5p, miR-486-5p, miR-582-5p, miR-1260a and miR-7977, as well as GPR143, OCA2, RAB27A, RAB32 and TYRP1 mRNAs, are differentially expressed only in pigmentable cells. These miRNAs and mRNAs belong to BPs related to pigmentation, specifically melanosome maturation and trafficking. In fact, an increase in the number of intracellular melanosomes-due to increased maturation and/or trafficking-confers resistance to vem. CONCLUSION: We demonstrated that the ability of pigmentable cells to increase the number of intracellular melanosomes fully accounts for their higher resistance to vem compared to non-pigmentable cells. In addition, we identified a network of miRNAs and mRNAs that are involved in melanosome maturation and/or trafficking. Finally, we provide the rationale for testing BRAFi in combination with inhibitors of these biological processes, so that pigmentable melanoma cells can be turned into more sensitive non-pigmentable cells.

A network approach to define the predictive role of immune profile on tumor response and toxicity of anti PD-1 single agent immunotherapy in patients with solid tumors.

Background: The immune profile of each patient could be considered as a portrait of the fitness of his/her own immune system. The predictive role of the immune profile in immune-related toxicities (irAEs) development and tumour response to treatment was investigated. Methods: A prospective, multicenter study evaluating, through a multiplex assay, the soluble immune profile at the baseline of 53 patients with advanced cancer, treated with immunotherapy as single agent was performed. Four connectivity heat maps and networks were obtained by calculating the Spearman correlation coefficients for each group: responder patients who developed cumulative toxicity (R-T), responders who did not develop cumulative toxicity (R-NT), non-responders who developed cumulative toxicity (NR-T), non-responders who did not develop cumulative toxicity (NR-NT). Results: A statistically significant up-regulation of IL-17A, sCTLA4, sCD80, I-CAM-1, sP-Selectin and sEselectin in NR-T was detected. A clear loss of connectivity of most of the soluble immune checkpoints and cytokines characterized the immune profile of patients with toxicity, while an inversion of the correlation for ICAM-1 and sP-selectin was observed in NR-T. Four connectivity networks were built for each group. The highest number of connections characterized the NR-T. Conclusions: A connectivity network of immune dysregulation was defined for each subgroup of patients, regardless of tumor type. In patients with the worst prognosis (NR-T) the peculiar connectivity model could facilitate their early and timely identification, as well as the design of a personalized treatment approach to improve outcomes or prevent irAEs.

The Impact of Drug-Drug Interactions on the Toxicity Profile of Combined Treatment with BRAF and MEK Inhibitors in Patients with BRAF-Mutated Metastatic Melanoma.

BACKGROUND: BRAF and MEK inhibition is a successful strategy in managing BRAF-mutant melanoma, even if the treatment-related toxicity is substantial. We analyzed the role of drug-drug interactions (DDI) on the toxicity profile of anti-BRAF/anti-MEK therapy. METHODS: In this multicenter, observational, and retrospective study, DDIs were assessed using Drug-PIN software (V 2/23). The association between the Drug-PIN continuous score or the Drug-PIN traffic light and the occurrence of treatment-related toxicities and oncological outcomes was evaluated. RESULTS: In total, 177 patients with advanced BRAF-mutated melanoma undergoing BRAF/MEK targeted therapy were included. All grade toxicity was registered in 79% of patients. Cardiovascular toxicities occurred in 31 patients (17.5%). Further, 94 (55.9%) patients had comorbidities requiring specific pharmacological treatments. The median Drug-PIN score significantly increased when the target combination was added to the patient's home therapy (p-value < 0.0001). Cardiovascular toxicity was significantly associated with the Drug-PIN score (p-value = 0.048). The Drug-PIN traffic light (p = 0.00821) and the Drug-PIN score (p = 0.0291) were seen to be significant predictors of cardiotoxicity. Patients with low-grade vs. high-grade interactions showed a better prognosis regarding overall survival (OS) (p = 0.0045) and progression-free survival (PFS) (p = 0.012). The survival analysis of the subgroup of patients with cardiological toxicity demonstrated that patients with low-grade vs. high-grade DDIs had better outcomes in terms of OS (p = 0.0012) and a trend toward significance in PFS (p = 0.068). CONCLUSIONS: DDIs emerged as a critical issue for the risk of treatment-related cardiovascular toxicity. Our findings support the utility of DDI assessment in melanoma patients treated with BRAF/MEK inhibitors.