At the Biological Modeling and Simulation Frontier

We provide a rationale for and describe examples of synthetic modeling and simulation (M&S) of biological systems. We explain how synthetic methods are distinct from familiar inductive methods. Synthetic M&S is a means to better understand the mechanisms that generate normal and disease-related phenomena observed in research, and how compounds of interest interact with them to alter phenomena. An objective is to build better, working hypotheses of plausible mechanisms. A synthetic model is an extant hypothesis: execution produces an observable mechanism and phenomena. Mobile objects representing compounds carry information enabling components to distinguish between them and react accordingly when different compounds are studied simultaneously. We argue that the familiar inductive approaches contribute to the general inefficiencies being experienced by pharmaceutical R&D, and that use of synthetic approaches accelerates and improves R&D decision-making and thus the drug development process. A reason is that synthetic models encourage and facilitate abductive scientific reasoning, a primary means of knowledge creation and creative cognition. When synthetic models are executed, we observe different aspects of knowledge in action from different perspectives. These models can be tuned to reflect differences in experimental conditions and individuals, making translational research more concrete while moving us closer to personalized medicine

A Phylogeny and Timescale for the Evolution of Pseudocheiridae (Marsupialia: Diprotodontia) in Australia and New Guinea

Pseudocheiridae (Marsupialia: Diprotodontia) is a family of endemic Australasian arboreal folivores, more commonly known as ringtail possums. Seventeen extant species are grouped into six genera (Pseudocheirus, Pseudochirulus, Hemibelideus, Petauroides, Pseudochirops, Petropseudes). Pseudochirops and Pseudochirulus are the only genera with representatives on New Guinea and surrounding western islands. Here, we examine phylogenetic relationships among 13 of the 17 extant pseudocheirid species based on protein-coding portions of the ApoB, BRCA1, ENAM, IRBP, Rag1, and vWF genes. Maximum parsimony, maximum likelihood, and Bayesian methods were used to estimate phylogenetic relationships. Two different relaxed molecular clock methods were used to estimate divergence times. Bayesian and maximum parsimony methods were used to reconstruct ancestral character states for geographic provenance and maximum elevation occupied. We find robust support for the monophyly of Pseudocheirinae (Pseudochirulus + Pseudocheirus), Hemibelidinae (Hemibelideus + Petauroides), and Pseudochiropsinae (Pseudochirops + Petropseudes), respectively, and for an association of Pseudocheirinae and Hemibelidinae to the exclusion of Pseudochiropsinae. Within Pseudochiropsinae, Petropseudes grouped more closely with the New Guinean Pseudochirops spp. than with the Australian Pseudochirops archeri, rendering Pseudochirops paraphyletic. New Guinean species belonging to Pseudochirops are monophyletic, as are New Guinean species belonging to Pseudochirulus. Molecular dates and ancestral reconstructions of geographic provenance combine to suggest that the ancestors of extant New Guinean Pseudochirops spp. and Pseudochirulus spp. dispersed from Australia to New Guinea ∼12.1–6.5 Ma (Pseudochirops) and ∼6.0–2.4 Ma (Pseudochirulus). Ancestral state reconstructions support the hypothesis that occupation of high elevations (>3000 m) is a derived feature that evolved on the terminal branch leading to Pseudochirops cupreus, and either evolved in the ancestor of Pseudochirulus forbesi, Pseudochirulus mayeri, and Pseudochirulus caroli, with subsequent loss in P. caroli, or evolved independently in P. mayeri and P. forbesi. Divergence times within the New Guinean Pseudochirops clade are generally coincident with the uplift of the central cordillera and other highlands. Diversification within New Guinean Pseudochirulus occurred in the Plio-Pleistocene after the establishment of the Central Range and other highlands

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

Huntington's disease mice and human brain tissue exhibit increased G3BP1 granules and TDP43 mislocalization

Chronic cellular stress associated with neurodegenerative disease can result in the persistence of stress granule (SG) structures, membraneless organelles that form in response to cellular stress. In Huntington's disease (HD), chronic expression of mutant huntingtin generates various forms of cellular stress, including activation of the unfolded protein response and oxidative stress. However, it has yet to be determined whether SGs are a feature of HD neuropathology. We examined the miRNA composition of extracellular vesicles (EVs) present in the cerebrospinal fluid (CSF) of HD patients and show that a subset of their target mRNAs were differentially expressed in the prefrontal cortex of HD patients. Of these targets, SG components were enriched, including the SG nucleating Ras GTPase-activating protein-binding protein 1 (G3BP1). We investigated localization and levels of G3BP1 and found a significant increase in the density of G3BP1-positive granules in the cortex and hippocampus of R6/2 transgenic mice and in the superior frontal cortex of HD patient brains. Intriguingly, we also observed that the SG-associated TAR DNA-Binding Protein-43 (TDP43), a nuclear RNA/DNA binding protein, was mislocalized to the cytoplasm of G3BP1-granule positive HD cortical neurons. These findings suggest that G3BP1 SG dynamics may play a role in the pathophysiology of HD

Cryptic diversity of a widespread global pathogen reveals expanded threats to amphibian conservation

Biodiversity loss is one major outcome of human-mediated ecosystem disturbance. One way that humans have triggered wildlife declines is by transporting disease-causing agents to remote areas of the world. Amphibians have been hit particularly hard by disease due in part to a globally distributed pathogenic chytrid fungus (Batrachochytrium dendrobatidis [Bd]). Prior research has revealed important insights into the biology and distribution of Bd; however, there are still many outstanding questions in this system. Although we know that there are multiple divergent lineages of Bd that differ in pathogenicity, we know little about how these lineages are distributed around the world and where lineages may be coming into contact. Here, we implement a custom genotyping method for a global set of Bd samples. This method is optimized to amplify and sequence degraded DNA from noninvasive skin swab samples. We describe a divergent lineage of Bd, which we call BdASIA3, that appears to be widespread in Southeast Asia. This lineage co-occurs with the global panzootic lineage (BdGPL) in multiple localities. Additionally, we shed light on the global distribution of BdGPL and highlight the expanded range of another lineage, BdCAPE. Finally, we argue that more monitoring needs to take place where Bd lineages are coming into contact and where we know little about Bd lineage diversity. Monitoring need not use expensive or difficult field techniques but can use archived swab samples to further explore the history-and predict the future impacts-of this devastating pathogen