ASU Scholarship Showcase

In this study, we present a novel methodology to infer indel parameters from multiple sequence alignments (MSAs) based on simulations. Our algorithm searches for the set of evolutionary parameters describing indel dynamics which best fits a given input MSA. In each step of the search, we use parametric bootstraps and the Mahalanobis distance to estimate how well a proposed set of parameters fits input data. Using simulations, we demonstrate that our methodology can accurately infer the indel parameters for a large variety of plausible settings. Moreover, using our methodology, we show that indel parameters substantially vary between three genomic data sets: Mammals, bacteria, and retroviruses. Finally, we demonstrate how our methodology can be used to simulate MSAs based on indel parameters inferred from real data sets.

Mathematical models of infectious diseases are a valuable tool in understanding the mechanisms and patterns of disease transmission. It is, however, a difficult subject to teach, requiring both mathematical expertise and extensive subject-matter knowledge of a variety of disease systems. In this article, we explore several uses of zombie epidemics to make mathematical modeling and infectious disease epidemiology more accessible to public health professionals, students, and the general public. We further introduce a web-based simulation, White Zed (http://cartwrig.ht/apps/whitezed/), that can be deployed in classrooms to allow students to explore models before implementing them. In our experience, zombie epidemics are familiar, approachable, flexible, and an ideal way to introduce basic concepts of infectious disease epidemiology.

Nonsense-mediated RNA decay (NMD) is a highly conserved pathway that selectively degrades specific subsets of RNA transcripts. Here, we provide evidence that NMD regulates early human developmental cell fate. We found that NMD factors tend to be expressed at higher levels in human pluripotent cells than in differentiated cells, raising the possibility that NMD must be downregulated to permit differentiation. Loss- and gain-of-function experiments in human embryonic stem cells (hESCs) demonstrated that, indeed, NMD downregulation is essential for efficient generation of definitive endoderm. RNA-seq analysis identified NMD target transcripts induced when NMD is suppressed in hESCs, including many encoding signaling components. This led us to test the role of TGF-β and BMP signaling, which we found NMD acts through to influence definitive endoderm versus mesoderm fate. Our results suggest that selective RNA decay is critical for specifying the developmental fate of specific human embryonic cell lineages.

Due to the limitation of current pharmacological therapeutic strategies, stem cell therapies have emerged as a viable option for treating many incurable neurological disorders. Specifically, human pluripotent stem cell (hPSC)-derived neural progenitor cells (hNPCs), a multipotent cell population that is capable of near indefinite expansion and subsequent differentiation into the various cell types that comprise the central nervous system (CNS), could provide an unlimited source of cells for such cell-based therapies. However the clinical application of these cells will require (i) defined, xeno-free conditions for their expansion and neuronal differentiation and (ii) scalable culture systems that enable their expansion and neuronal differentiation in numbers sufficient for regenerative medicine and drug screening purposes. Current extracellular matrix protein (ECMP)-based substrates for the culture of hNPCs are expensive, difficult to isolate, subject to batch-to-batch variations, and, therefore, unsuitable for clinical application of hNPCs. Using a high-throughput array-based screening approach, we identified a synthetic polymer, poly(4-vinyl phenol) (P4VP), that supported the long-term proliferation and self-renewal of hNPCs. The hNPCs cultured on P4VP maintained their characteristic morphology, expressed high levels of markers of multipotency, and retained their ability to differentiate into neurons. Such chemically defined substrates will eliminate critical roadblocks for the utilization of hNPCs for human neural regenerative repair, disease modeling, and drug discovery.

Mutation is the ultimate source of all genetic variation and is, therefore, central to evolutionary change. Previous work on Paramecium tetraurelia found an unusually low germline base-substitution mutation rate in this ciliate. Here, we tested the generality of this result among ciliates using Tetrahymena thermophila. We sequenced the genomes of 10 lines of T. thermophila that had each undergone approximately 1,000 generations of mutation accumulation (MA). We applied an existing mutation-calling pipeline and developed a new probabilistic mutation detection approach that directly models the design of an MA experiment and accommodates the noise introduced by mismapped reads. Our probabilistic mutation-calling method provides a straightforward way of estimating the number of sites at which a mutation could have been called if one was present, providing the denominator for our mutation rate calculations. From these methods, we find that T. thermophila has a germline base-substitution mutation rate of 7.61 × 10 ^-12 per-site, per cell division, which is consistent with the low base-substitution mutation rate in P. tetraurelia. Over the course of the evolution experiment, genomic exclusion lines derived from the MA lines experienced a fitness decline that cannot be accounted for by germline base-substitution mutations alone, suggesting that other genetic or epigenetic factors must be involved. Because selection can only operate to reduce mutation rates based upon the "visible" mutational load, asexual reproduction with a transcriptionally silent germline may allow ciliates to evolve extremely low germline mutation rates.

The most common evolutionary events at the molecular level are single-base substitutions, as well as insertions and deletions (indels) of short DNA segments. A large body of research has been devoted to develop probabilistic substitution models and to infer their parameters using likelihood and Bayesian approaches. In contrast, relatively little has been done to model indel dynamics, probably due to the difficulty in writing explicit likelihood functions. Here, we contribute to the effort of modeling indel dynamics by presenting SpartaABC, an approximate Bayesian computation (ABC) approach to infer indel parameters from sequence data (either aligned or unaligned). SpartaABC circumvents the need to use an explicit likelihood function by extracting summary statistics from simulated sequences. First, summary statistics are extracted from the input sequence data. Second, SpartaABC samples indel parameters from a prior distribution and uses them to simulate sequences. Third, it computes summary statistics from the simulated sets of sequences. By computing a distance between the summary statistics extracted from the input and each simulation, SpartaABC can provide an approximation to the posterior distribution of indel parameters as well as point estimates. We study the performance of our methodology and show that it provides accurate estimates of indel parameters in simulations. We next demonstrate the utility of SpartaABC by studying the impact of alignment errors on the inference of positive selection. A C ++ program implementing SpartaABC is freely available in http://spartaabc.tau.ac.il.