Finding the shortest path with PesCa: A tool for network reconstruction

The growing dimension and complexity of the available experimental data generating biological networks have increased the need for tools that help in categorizing nodes by their topological relevance. Here we present CentiScaPe, a Cytoscape app specifically designed to calculate centrality indexes used for the identification of the most important nodes in a network. CentiScaPe is a comprehensive suite of algorithms dedicated to network nodes centrality analysis, computing several centralities for undirected, directed and weighted networks. The results of the topological analysis can be integrated with data set from lab experiments, like expression or phosphorylation levels for each protein represented in the network. Our app opens new perspectives in the analysis of biological networks, since the integration of topological analysis with lab experimental data enhance the predictive power of the bioinformatics analysis

Type II migration strikes back – an old paradigm for planet migration in discs

In this paper, we analyse giant gap-opening planet migration in proto-planetary discs, focusing on the type II migration regime. According to standard type II theory, planets migrate at the same rate as the gas in the disc, as they are coupled to the disc viscous evolution; however, recent studies questioned this paradigm, suggesting that planets migrate faster than the disc material. We study the problem through 2D long-time simulations of systems consistent with type II regime, using the hydrodynamical grid code FARGO3D. Even though our simulations confirm the presence of an initial phase characterized by fast migration, they also reveal that the migration velocity slows down and eventually reaches the theoretical prediction if we allow the system to evolve for enough time. We find the same tendency to evolve towards the theoretical predictions at later times when we analyse the mass flow through the gap and the torques acting on the planet. This transient is related to the initial conditions of our (and previous) simulations, and is due to the fact that the shape of the gap has to adjust to a new profile, once the planet is set into motion. Secondly, we test whether the type II theory expectation that giant planet migration is driven by viscosity is consistent with our simulation by comparing simulations with the same viscosity and different disc masses (or vice versa). We find a good agreement with the theory, since when the discs are characterized by the same viscosity, the migration properties are the same

Analysing omics data sets with weighted nodes networks (WNNets)

Current trends in biomedical research indicate data integration as a fundamental step towards precision medicine. In this context, network models allow representing and analysing complex biological processes. However, although effective in unveiling network properties, these models fail in considering the individual, biochemical variations occurring at molecular level. As a consequence, the analysis of these models partially loses its predictive power. To overcome these limitations, Weighted Nodes Networks (WNNets) were developed. WNNets allow to easily and effectively weigh nodes using experimental information from multiple conditions. In this study, the characteristics of WNNets were described and a proteomics data set was modelled and analysed. Results suggested that degree, an established centrality index, may offer a novel perspective about the functional role of nodes in WNNets. Indeed, degree allowed retrieving significant differences between experimental conditions, highlighting relevant proteins, and provided a novel interpretation for degree itself, opening new perspectives in experimental data modelling and analysis. Overall, WNNets may be used to model any high-throughput experimental data set requiring weighted nodes. Finally, improving the power of the analysis by using centralities such as betweenness may provide further biological insights and unveil novel, interesting characteristics of WNNets

Protoplanetary disc ‘isochrones’ and the evolution of discs in the M˙-Md plane

In this paper, we compare simple viscous diffusion models for the disc evolution with the results of recent surveys of the properties of young protoplanetary discs. We introduce the useful concept of \u2018disc isochrones\u2019 in the accretion rate\u2013disc mass plane and explore a set of Monte Carlo realization of disc initial conditions. We find that such simple viscous models can provide a remarkable agreement with the available data in the Lupus star forming region, with the key requirement that the average viscous evolutionary time-scale of the discs is comparable to the cluster age. Our models produce naturally a correlation between mass accretion rate and disc mass that is shallower than linear, contrary to previous results and in agreement with observations. We also predict that a linear correlation, with a tighter scatter, should be found for more evolved disc populations. Finally, we find that such viscous models can reproduce the observations in the Lupus region only in the assumption that the efficiency of angular momentum transport is a growing function of radius, thus putting interesting constraints on the nature of the microscopic processes that lead to disc accretion

Computational Identification of Phospho-Tyrosine Sub-Networks Related to Acanthocyte Generation in Neuroacanthocytosis

Acanthocytes, abnormal thorny red blood cells (RBC), are one of the biological hallmarks of neuroacanthocytosis syndromes (NA), a group of rare hereditary neurodegenerative disorders. Since RBCs are easily accessible, the study of acanthocytes in NA may provide insights into potential mechanisms of neurodegeneration. Previous studies have shown that changes in RBC membrane protein phosphorylation state affect RBC membrane mechanical stability and morphology. Here, we coupled tyrosine-phosphoproteomic analysis to topological network analysis. We aimed to predict signaling sub-networks possibly involved in the generation of acanthocytes in patients affected by the two core NA disorders, namely McLeod syndrome (MLS, XK-related, Xk protein) and chorea-acanthocytosis (ChAc, VPS13A-related, chorein protein). The experimentally determined phosphoproteomic data-sets allowed us to relate the subsequent network analysis to the pathogenetic background. To reduce the network complexity, we combined several algorithms of topological network analysis including cluster determination by shortest path analysis, protein categorization based on centrality indexes, along with annotation-based node filtering. We first identified XK- and VPS13A-related protein-protein interaction networks by identifying all the interactomic shortest paths linking Xk and chorein to the corresponding set of proteins whose tyrosine phosphorylation was altered in patients. These networks include the most likely paths of functional influence of Xk and chorein on phosphorylated proteins. We further refined the analysis by extracting restricted sets of highly interacting signaling proteins representing a common molecular background bridging the generation of acanthocytes in MLS and ChAc. The final analysis pointed to a novel, very restricted, signaling module of 14 highly interconnected kinases, whose alteration is possibly involved in generation of acanthocytes in MLS and ChAc

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale(1-3). Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter(4); identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation(5,6); analyses timings and patterns of tumour evolution(7); describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity(8,9); and evaluates a range of more-specialized features of cancer genomes(8,10-18).Peer reviewe

Identification of unique neoantigen qualities in long-term survivors of pancreatic cancer

Pancreatic ductal adenocarcinoma is a lethal cancer with fewer than 7% of patients surviving past 5 years. T-cell immunity has been linked to the exceptional outcome of the few long-term survivors1,2, yet the relevant antigens remain unknown. Here we use genetic, immunohistochemical and transcriptional immunoprofiling, computational biophysics, and functional assays to identify T-cell antigens in long-term survivors of pancreatic cancer. Using whole-exome sequencing and in silico neoantigen prediction, we found that tumours with both the highest neoantigen number and the most abundant CD8+ T-cell infiltrates, but neither alone, stratified patients with the longest survival. Investigating the specific neoantigen qualities promoting T-cell activation in long-term survivors, we discovered that these individuals were enriched in neoantigen qualities defined by a fitness model, and neoantigens in the tumour antigen MUC16 (also known as CA125). A neoantigen quality fitness model conferring greater immunogenicity to neoantigens with differential presentation and homology to infectious disease-derived peptides identified long-term survivors in two independent datasets, whereas a neoantigen quantity model ascribing greater immunogenicity to increasing neoantigen number alone did not. We detected intratumoural and lasting circulating T-cell reactivity to both high-quality and MUC16 neoantigens in long-term survivors of pancreatic cancer, including clones with specificity to both high-quality neoantigens and predicted cross-reactive microbial epitopes, consistent with neoantigen molecular mimicry. Notably, we observed selective loss of high-quality and MUC16 neoantigenic clones on metastatic progression, suggesting neoantigen immunoediting. Our results identify neoantigens with unique qualities as T-cell targets in pancreatic ductal adenocarcinoma. More broadly, we identify neoantigen quality as a biomarker for immunogenic tumours that may guide the application of immunotherapies

Targeting DNA Damage Response and Replication Stress in Pancreatic Cancer

Background and aims: Continuing recalcitrance to therapy cements pancreatic cancer (PC) as the most lethal malignancy, which is set to become the second leading cause of cancer death in our society. The study aim was to investigate the association between DNA damage response (DDR), replication stress and novel therapeutic response in PC to develop a biomarker driven therapeutic strategy targeting DDR and replication stress in PC. Methods: We interrogated the transcriptome, genome, proteome and functional characteristics of 61 novel PC patient-derived cell lines to define novel therapeutic strategies targeting DDR and replication stress. Validation was done in patient derived xenografts and human PC organoids. Results: Patient-derived cell lines faithfully recapitulate the epithelial component of pancreatic tumors including previously described molecular subtypes. Biomarkers of DDR deficiency, including a novel signature of homologous recombination deficiency, co-segregates with response to platinum (P < 0.001) and PARP inhibitor therapy (P < 0.001) in vitro and in vivo. We generated a novel signature of replication stress with which predicts response to ATR (P < 0.018) and WEE1 inhibitor (P < 0.029) treatment in both cell lines and human PC organoids. Replication stress was enriched in the squamous subtype of PC (P < 0.001) but not associated with DDR deficiency. Conclusions: Replication stress and DDR deficiency are independent of each other, creating opportunities for therapy in DDR proficient PC, and post-platinum therapy