A seven-step guide to spatial, agent-based modelling of tumour evolution.

Spatial agent-based models are frequently used to investigate the evolution of solid tumours subject to localized cell-cell interactions and microenvironmental heterogeneity. As spatial genomic, transcriptomic and proteomic technologies gain traction, spatial computational models are predicted to become ever more necessary for making sense of complex clinical and experimental data sets, for predicting clinical outcomes, and for optimizing treatment strategies. Here we present a non-technical step by step guide to developing such a model from first principles. Stressing the importance of tailoring the model structure to that of the biological system, we describe methods of increasing complexity, from the basic Eden growth model up to off-lattice simulations with diffusible factors. We examine choices that unavoidably arise in model design, such as implementation, parameterization, visualization and reproducibility. Each topic is illustrated with examples drawn from recent research studies and state of the art modelling platforms. We emphasize the benefits of simpler models that aim to match the complexity of the phenomena of interest, rather than that of the entire biological system. Our guide is aimed at both aspiring modellers and other biologists and oncologists who wish to understand the assumptions and limitations of the models on which major cancer studies now so often depend.

MAPP unravels frequent co-regulation of splicing and polyadenylation by RNA-binding proteins and their dysregulation in cancer.

Maturation of eukaryotic pre-mRNAs via splicing and polyadenylation is modulated across cell types and conditions by a variety of RNA-binding proteins (RBPs). Although there exist over 1,500 RBPs in human cells, their binding motifs and functions still remain to be elucidated, especially in the complex environment of tissues and in the context of diseases. To overcome the lack of methods for the systematic and automated detection of sequence motif-guided pre-mRNA processing regulation from RNA sequencing (RNA-Seq) data we have developed MAPP (Motif Activity on Pre-mRNA Processing). Applying MAPP to RBP knock-down experiments reveals that many RBPs regulate both splicing and polyadenylation of nascent transcripts by acting on similar sequence motifs. MAPP not only infers these sequence motifs, but also unravels the position-dependent impact of the RBPs on pre-mRNA processing. Interestingly, all investigated RBPs that act on both splicing and 3' end processing exhibit a consistently repressive or activating effect on both processes, providing a first glimpse on the underlying mechanism. Applying MAPP to normal and malignant brain tissue samples unveils that the motifs bound by the PTBP1 and RBFOX RBPs coordinately drive the oncogenic splicing program active in glioblastomas demonstrating that MAPP paves the way for characterizing pre-mRNA processing regulators under physiological and pathological conditions.

The ubiquity of clonal interference in cancer and other range expansions.

Evolution during range expansions shapes the fate of many biological systems including tumours, microbial communities, and invasive species. A fundamental process of interest is the selective sweep, in which an advantageous mutation evades clonal interference and spreads through the population to fixation. However, most theoretical investigations of selective sweeps have assumed constant population size or have ignored spatial structure. Here we use mathematical modelling and analysis to investigate selective sweep probabilities in populations that spread outwards as they evolve. In the case of constant radial expansion speed, we derive probability distributions for the arrival time and location of the first surviving mutant and hence find surprisingly simple approximate and exact expressions for selective sweep probabilities in one, two and three dimensions, which are independent of mutation rate. Namely, the selective sweep probability is approximately 1-cwt/cmd, where cwt and cm are the wildtype and mutant radial expansion speeds, and d the spatial dimension. Using agent-based simulations, we show that our analytical results accurately predict selective sweep frequencies in the two-dimensional spatial Moran process. We further compare our results with those obtained for alternative growth models. Parameterizing our model for human tumours, we find that selective sweeps are predicted to be rare except during very early solid tumour growth, thus providing a general, pan-cancer explanation for findings from recent sequencing studies.

Alternative splicing is not a key source of chemerin isoforms diversity.

BACKGROUND: Chemerin is a chemoattractant protein with adipokine and antimicrobial properties encoded by the retinoic acid receptor responder 2 (RARRES2) gene. Chemerin bioactivity largely depends on carboxyl-terminal proteolytic processing that generates chemerin isoforms with different chemotactic, regulatory, and antimicrobial potentials. While these mechanisms are relatively well known, the role of alternative splicing in generating isoform diversity remains obscure. METHODS AND RESULTS: Using rapid amplification of cDNA ends (RACE) PCR, we determined RARRES2 transcript variants present in mouse and human tissues and identified novel transcript variant 4 of mouse Rarres2 encoding mChem153K. Moreover, analyses of real-time quantitative PCR (RT-qPCR) and publicly-available next-generation RNA sequencing (RNA-seq) datasets showed that different alternatively spliced variants of mouse Rarres2 are present in mouse tissues and their expression patterns were unaffected by inflammatory and infectious stimuli except brown adipose tissue. However, only one transcript variant of human RARRES2 was present in liver and adipose tissue. CONCLUSION: Our findings indicate a limited role for alternative splicing in generating chemerin isoform diversity under all tested conditions.

CFIm-mediated alternative polyadenylation remodels cellular signaling and miRNA biogenesis.

The mammalian cleavage factor I (CFIm) has been implicated in alternative polyadenylation (APA) in a broad range of contexts, from cancers to learning deficits and parasite infections. To determine how the CFIm expression levels are translated into these diverse phenotypes, we carried out a multi-omics analysis of cell lines in which the CFIm25 (NUDT21) or CFIm68 (CPSF6) subunits were either repressed by siRNA-mediated knockdown or over-expressed from stably integrated constructs. We established that >800 genes undergo coherent APA in response to changes in CFIm levels, and they cluster in distinct functional classes related to protein metabolism. The activity of the ERK pathway traces the CFIm concentration, and explains some of the fluctuations in cell growth and metabolism that are observed upon CFIm perturbations. Furthermore, multiple transcripts encoding proteins from the miRNA pathway are targets of CFIm-dependent APA. This leads to an increased biogenesis and repressive activity of miRNAs at the same time as some 3' UTRs become shorter and presumably less sensitive to miRNA-mediated repression. Our study provides a first systematic assessment of a core set of APA targets that respond coherently to changes in CFIm protein subunit levels (CFIm25/CFIm68). We describe the elicited signaling pathways downstream of CFIm, which improve our understanding of the key role of CFIm in integrating RNA processing with other cellular activities.

The methylation status of the chemerin promoter region located from - 252 to + 258 bp regulates constitutive but not acute-phase cytokine-inducible chemerin expression levels.

Chemerin is a chemoattractant protein with adipokine properties encoded by the retinoic acid receptor responder 2 (RARRES2) gene. It has gained more attention in the past few years due to its multilevel impact on metabolism and immune responses. However, mechanisms controlling the constitutive and regulated expression of RARRES2 in a variety of cell types remain obscure. To our knowledge, this report is the first to show that DNA methylation plays an important role in the cell-specific expression of RARRES2 in adipocytes, hepatocytes, and B lymphocytes. Using luciferase reporter assays, we determined the proximal fragment of the RARRES2 gene promoter, located from - 252 to + 258 bp, to be a key regulator of transcription. Moreover, we showed that chemerin expression is regulated in murine adipocytes by acute-phase cytokines, interleukin 1ß and oncostatin M. In contrast with adipocytes, these cytokines exerted a weak, if any, response in mouse hepatocytes, suggesting that the effects of IL-1ß and OSM on chemerin expression is specific to fat tissue. Together, our findings highlight previously uncharacterized mediators and mechanisms that control chemerin expression.

The Alazami Syndrome-Associated Protein LARP7 Guides U6 Small Nuclear RNA Modification and Contributes to Splicing Robustness.

The La-related protein 7 (LARP7) forms a complex with the nuclear 7SK RNA to regulate RNA polymerase II transcription. It has been implicated in cancer and the Alazami syndrome, a severe developmental disorder. Here, we report a so far unknown role of this protein in RNA modification. We show that LARP7 physically connects the spliceosomal U6 small nuclear RNA (snRNA) with a distinct subset of box C/D small nucleolar RNAs (snoRNAs) guiding U6 2'-O-methylation. Consistently, these modifications are severely compromised in the absence of LARP7. Although general splicing remains largely unaffected, transcriptome-wide analysis revealed perturbations in alternative splicing in LARP7-depleted cells. Importantly, we identified defects in 2'-O-methylation of the U6 snRNA in Alazami syndrome siblings carrying a LARP7 mutation. Our data identify LARP7 as a bridging factor for snoRNA-guided modification of the U6 snRNA and suggest that alterations in splicing fidelity contribute to the etiology of the Alazami syndrome.

Host Taxon Predictor - A Tool for Predicting Taxon of the Host of a Newly Discovered Virus.

Recent advances in metagenomics provided a valuable alternative to culture-based approaches for better sampling viral diversity. However, some of newly identified viruses lack sequence similarity to any of previously sequenced ones, and cannot be easily assigned to their hosts. Here we present a bioinformatic approach to this problem. We developed classifiers capable of distinguishing eukaryotic viruses from the phages achieving almost 95% prediction accuracy. The classifiers are wrapped in Host Taxon Predictor (HTP) software written in Python which is freely available at https://github.com/wojciech-galan/viruses_classifier . HTP's performance was later demonstrated on a collection of newly identified viral genomes and genome fragments. In summary, HTP is a culture- and alignment-free approach for distinction between phages and eukaryotic viruses. We have also shown that it is possible to further extend our method to go up the evolutionary tree and predict whether a virus can infect narrower taxa.

Differential inflammatory microRNA and cytokine expression in pulmonary sarcoidosis.

Sarcoidosis is a granulomatous disease of unknown etiology. The disease has an important inflammatory and immune component; however, its immunopathogenesis is not completely understood. Recently, the role of microRNAs (miRNAs), the small non-coding RNAs, has attracted attention as both being involved in pathogenesis and serving as disease markers. Accordingly, changes in the expression of some miRNAs have been also associated with different autoimmune pathologies. However, not much is known about the role of miRNAs in sarcoidosis. Therefore, the aim of this study was to compare the level of expression of selected miRNAs in healthy individuals and patients with sarcoidosis. We detected significantly increased level of miR-34a in peripheral blood mononuclear cells isolated from sarcoidosis patients. Moreover, significantly up-regulated levels of interferon (IFN)-γ, IFN-γ inducible protein (IP-10) and vascular endothelial growth factor were detected in sera of patients when compared to healthy subjects. Our results add to a known inflammatory component in sarcoidosis. Changes in the levels of miR-34a may suggest its involvement in the pathology of this disease.