ASU Scholarship Showcase

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.

Pay-for-performance (PFP) is a relatively new approach to agricultural conservation that attaches an incentive payment to quantified reductions in nutrient runoff from a participating farm. Similar to a payment for ecosystem services approach, PFP lends itself to providing incentives for the most beneficial practices at the field level. To date, PFP conservation in the U.S. has only been applied in small pilot programs. Because monitoring conservation performance for each field enrolled in a program would be cost-prohibitive, field-level modeling can provide cost-effective estimates of anticipated improvements in nutrient runoff. We developed a PFP system that uses a unique application of one of the leading agricultural models, the USDA's Soil and Water Assessment Tool, to evaluate the nutrient load reductions of potential farm practice changes based on field-level agronomic and management data. The initial phase of the project focused on simulating individual fields in the River Raisin watershed in southeastern Michigan. Here we present development of the modeling approach and results from the pilot year, 2015-2016. These results stress that (1) there is variability in practice effectiveness both within and between farms, and thus there is not one "best practice" for all farms, (2) conservation decisions are made most effectively at the scale of the farm field rather than the sub-watershed or watershed level, and (3) detailed, field-level management information is needed to accurately model and manage on-farm nutrient loadings.

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.