@misc{9849, abstract = {This text provides additional information about the model, a derivation of the analytic results in Eq (4), and details about simulations of an additional parameter set.}, author = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago}, publisher = {Public Library of Science}, title = {{Modelling and simulation details}}, doi = {10.1371/journal.pcbi.1005609.s001}, year = {2017}, } @misc{9850, abstract = {In this text, we discuss how a cost of resistance and the possibility of lethal mutations impact our model.}, author = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago}, publisher = {Public Library of Science}, title = {{Extensions of the model}}, doi = {10.1371/journal.pcbi.1005609.s002}, year = {2017}, } @misc{9851, abstract = {Based on the intuitive derivation of the dynamics of SIM allele frequency pM in the main text, we present a heuristic prediction for the long-term SIM allele frequencies with χ > 1 stresses and compare it to numerical simulations.}, author = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago}, publisher = {Public Library of Science}, title = {{Heuristic prediction for multiple stresses}}, doi = {10.1371/journal.pcbi.1005609.s003}, year = {2017}, } @misc{9852, abstract = {We show how different combination strategies affect the fraction of individuals that are multi-resistant.}, author = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago}, publisher = {Public Library of Science}, title = {{Resistance frequencies for different combination strategies}}, doi = {10.1371/journal.pcbi.1005609.s004}, year = {2017}, } @misc{9856, author = {Schmidt, Tom and Barton, Nicholas H and Rasic, Gordana and Turley, Andrew and Montgomery, Brian and Iturbe Ormaetxe, Inaki and Cook, Peter and Ryan, Peter and Ritchie, Scott and Hoffmann, Ary and O’Neill, Scott and Turelli, Michael}, publisher = {Public Library of Science}, title = {{Supporting Information concerning additional likelihood analyses and results}}, doi = {10.1371/journal.pbio.2001894.s014}, year = {2017}, } @misc{9858, author = {Schmidt, Tom and Barton, Nicholas H and Rasic, Gordana and Turley, Andrew and Montgomery, Brian and Iturbe Ormaetxe, Inaki and Cook, Peter and Ryan, Peter and Ritchie, Scott and Hoffmann, Ary and O’Neill, Scott and Turelli, Michael}, publisher = {Public Library of Science}, title = {{Excel file with data on mosquito densities, Wolbachia infection status and housing characteristics}}, doi = {10.1371/journal.pbio.2001894.s016}, year = {2017}, } @article{990, abstract = {Assortative mating is an important driver of speciation in populations with gene flow and is predicted to evolve under certain conditions in few-locus models. However, the evolution of assortment is less understood for mating based on quantitative traits, which are often characterized by high genetic variability and extensive linkage disequilibrium between trait loci. We explore this scenario for a two-deme model with migration, by considering a single polygenic trait subject to divergent viability selection across demes, as well as assortative mating and sexual selection within demes, and investigate how trait divergence is shaped by various evolutionary forces. Our analysis reveals the existence of sharp thresholds of assortment strength, at which divergence increases dramatically. We also study the evolution of assortment via invasion of modifiers of mate discrimination and show that the ES assortment strength has an intermediate value under a range of migration-selection parameters, even in diverged populations, due to subtle effects which depend sensitively on the extent of phenotypic variation within these populations. The evolutionary dynamics of the polygenic trait is studied using the hypergeometric and infinitesimal models. We further investigate the sensitivity of our results to the assumptions of the hypergeometric model, using individual-based simulations.}, author = {Sachdeva, Himani and Barton, Nicholas H}, issn = {0014-3820}, journal = {Evolution; International Journal of Organic Evolution}, number = {6}, pages = {1478 -- 1493 }, publisher = {Wiley-Blackwell}, title = {{Divergence and evolution of assortative mating in a polygenic trait model of speciation with gene flow}}, doi = {10.1111/evo.13252}, volume = {71}, year = {2017}, } @article{951, abstract = {Dengue-suppressing Wolbachia strains are promising tools for arbovirus control, particularly as they have the potential to self-spread following local introductions. To test this, we followed the frequency of the transinfected Wolbachia strain wMel through Ae. aegypti in Cairns, Australia, following releases at 3 nonisolated locations within the city in early 2013. Spatial spread was analysed graphically using interpolation and by fitting a statistical model describing the position and width of the wave. For the larger 2 of the 3 releases (covering 0.97 km2 and 0.52 km2), we observed slow but steady spatial spread, at about 100–200 m per year, roughly consistent with theoretical predictions. In contrast, the smallest release (0.11 km2) produced erratic temporal and spatial dynamics, with little evidence of spread after 2 years. This is consistent with the prediction concerning fitness-decreasing Wolbachia transinfections that a minimum release area is needed to achieve stable local establishment and spread in continuous habitats. Our graphical and likelihood analyses produced broadly consistent estimates of wave speed and wave width. Spread at all sites was spatially heterogeneous, suggesting that environmental heterogeneity will affect large-scale Wolbachia transformations of urban mosquito populations. The persistence and spread of Wolbachia in release areas meeting minimum area requirements indicates the promise of successful large-scale population transfo}, author = {Schmidt, Tom and Barton, Nicholas H and Rasic, Gordana and Turley, Andrew and Montgomery, Brian and Iturbe Ormaetxe, Inaki and Cook, Peter and Ryan, Peter and Ritchie, Scott and Hoffmann, Ary and O’Neill, Scott and Turelli, Michael}, issn = {1544-9173}, journal = {PLoS Biology}, number = {5}, publisher = {Public Library of Science}, title = {{Local introduction and heterogeneous spatial spread of dengue-suppressing Wolbachia through an urban population of Aedes Aegypti}}, doi = {10.1371/journal.pbio.2001894}, volume = {15}, year = {2017}, } @misc{9857, author = {Schmidt, Tom and Barton, Nicholas H and Rasic, Gordana and Turley, Andrew and Montgomery, Brian and Iturbe Ormaetxe, Inaki and Cook, Peter and Ryan, Peter and Ritchie, Scott and Hoffmann, Ary and O’Neill, Scott and Turelli, Michael}, publisher = {Public Library of Science}, title = {{Supporting information concerning observed wMel frequencies and analyses of habitat variables}}, doi = {10.1371/journal.pbio.2001894.s015}, year = {2017}, } @article{696, abstract = {Mutator strains are expected to evolve when the availability and effect of beneficial mutations are high enough to counteract the disadvantage from deleterious mutations that will inevitably accumulate. As the population becomes more adapted to its environment, both availability and effect of beneficial mutations necessarily decrease and mutation rates are predicted to decrease. It has been shown that certain molecular mechanisms can lead to increased mutation rates when the organism finds itself in a stressful environment. While this may be a correlated response to other functions, it could also be an adaptive mechanism, raising mutation rates only when it is most advantageous. Here, we use a mathematical model to investigate the plausibility of the adaptive hypothesis. We show that such a mechanism can be mantained if the population is subjected to diverse stresses. By simulating various antibiotic treatment schemes, we find that combination treatments can reduce the effectiveness of second-order selection on stress-induced mutagenesis. We discuss the implications of our results to strategies of antibiotic therapy.}, author = {Lukacisinova, Marta and Novak, Sebastian and Paixao, Tiago}, issn = {1553-734X}, journal = {PLoS Computational Biology}, number = {7}, publisher = {Public Library of Science}, title = {{Stress induced mutagenesis: Stress diversity facilitates the persistence of mutator genes}}, doi = {10.1371/journal.pcbi.1005609}, volume = {13}, year = {2017}, } @article{1158, abstract = {Speciation results from the progressive accumulation of mutations that decrease the probability of mating between parental populations or reduce the fitness of hybrids—the so-called species barriers. The speciation genomic literature, however, is mainly a collection of case studies, each with its own approach and specificities, such that a global view of the gradual process of evolution from one to two species is currently lacking. Of primary importance is the prevalence of gene flow between diverging entities, which is central in most species concepts and has been widely discussed in recent years. Here, we explore the continuum of speciation thanks to a comparative analysis of genomic data from 61 pairs of populations/species of animals with variable levels of divergence. Gene flow between diverging gene pools is assessed under an approximate Bayesian computation (ABC) framework. We show that the intermediate "grey zone" of speciation, in which taxonomy is often controversial, spans from 0.5% to 2% of net synonymous divergence, irrespective of species life history traits or ecology. Thanks to appropriate modeling of among-locus variation in genetic drift and introgression rate, we clarify the status of the majority of ambiguous cases and uncover a number of cryptic species. Our analysis also reveals the high incidence in animals of semi-isolated species (when some but not all loci are affected by barriers to gene flow) and highlights the intrinsic difficulty, both statistical and conceptual, of delineating species in the grey zone of speciation.}, author = {Roux, Camille and Fraisse, Christelle and Romiguier, Jonathan and Anciaux, Youann and Galtier, Nicolas and Bierne, Nicolas}, journal = {PLoS Biology}, number = {12}, publisher = {Public Library of Science}, title = {{Shedding light on the grey zone of speciation along a continuum of genomic divergence}}, doi = {10.1371/journal.pbio.2000234}, volume = {14}, year = {2016}, } @article{1172, abstract = {A central issue in cell biology is the physico-chemical basis of organelle biogenesis in intracellular trafficking pathways, its most impressive manifestation being the biogenesis of Golgi cisternae. At a basic level, such morphologically and chemically distinct compartments should arise from an interplay between the molecular transport and chemical maturation. Here, we formulate analytically tractable, minimalist models, that incorporate this interplay between transport and chemical progression in physical space, and explore the conditions for de novo biogenesis of distinct cisternae. We propose new quantitative measures that can discriminate between the various models of transport in a qualitative manner-this includes measures of the dynamics in steady state and the dynamical response to perturbations of the kind amenable to live-cell imaging.}, author = {Sachdeva, Himani and Barma, Mustansir and Rao, Madan}, journal = {Scientific Reports}, publisher = {Nature Publishing Group}, title = {{Nonequilibrium description of de novo biogenesis and transport through Golgi-like cisternae}}, doi = {10.1038/srep38840}, volume = {6}, year = {2016}, } @article{1195, abstract = {The genetic analysis of experimentally evolving populations typically relies on short reads from pooled individuals (Pool-Seq). While this method provides reliable allele frequency estimates, the underlying haplotype structure remains poorly characterized. With small population sizes and adaptive variants that start from low frequencies, the interpretation of selection signatures in most Evolve and Resequencing studies remains challenging. To facilitate the characterization of selection targets, we propose a new approach that reconstructs selected haplotypes from replicated time series, using Pool-Seq data. We identify selected haplotypes through the correlated frequencies of alleles carried by them. Computer simulations indicate that selected haplotype-blocks of several Mb can be reconstructed with high confidence and low error rates, even when allele frequencies change only by 20% across three replicates. Applying this method to real data from D. melanogaster populations adapting to a hot environment, we identify a selected haplotype-block of 6.93 Mb. We confirm the presence of this haplotype-block in evolved populations by experimental haplotyping, demonstrating the power and accuracy of our haplotype reconstruction from Pool-Seq data. We propose that the combination of allele frequency estimates with haplotype information will provide the key to understanding the dynamics of adaptive alleles. }, author = {Franssen, Susan and Barton, Nicholas H and Schlötterer, Christian}, journal = {Molecular Biology and Evolution}, number = {1}, pages = {174 -- 184}, publisher = {Oxford University Press}, title = {{Reconstruction of haplotype-blocks selected during experimental evolution.}}, doi = {10.1093/molbev/msw210}, volume = {34}, year = {2016}, } @article{1224, abstract = {Sexual dimorphism in resource allocation is expected to change during the life cycle of dioecious plants because of temporal differences between the sexes in reproductive investment. Given the potential for sex-specific differences in reproductive costs, resource availability may contribute to variation in reproductive allocation in females and males. Here, we used Rumex hastatulus, a dioecious, wind-pollinated annual plant, to investigate whether sexual dimorphism varies with life-history stage and nutrient availability, and determine whether allocation patterns differ depending on reproductive commitment. To examine if the costs of reproduction varied between the sexes, reproduction was either allowed or prevented through bud removal, and biomass allocation was measured at maturity. In a second experiment to assess variation in sexual dimorphism across the life cycle, and whether this varied with resource availability, plants were grown in high and low nutrients and allocation to roots, aboveground vegetative growth and reproduction were measured at three developmental stages. Males prevented from reproducing compensated with increased above- and belowground allocation to a much larger degree than females, suggesting that male reproductive costs reduce vegetative growth. The proportional allocation to roots, reproductive structures and aboveground vegetative growth varied between the sexes and among life-cycle stages, but not with nutrient treatment. Females allocated proportionally more resources to roots than males at peak flowering, but this pattern was reversed at reproductive maturity under low-nutrient conditions. Our study illustrates the importance of temporal dynamics in sex-specific resource allocation and provides support for high male reproductive costs in wind-pollinated plants.}, author = {Teitel, Zachary and Pickup, Melinda and Field, David and Barrett, Spencer}, journal = {Plant Biology}, number = {1}, pages = {98 -- 103}, publisher = {Wiley-Blackwell}, title = {{The dynamics of resource allocation and costs of reproduction in a sexually dimorphic, wind-pollinated dioecious plant}}, doi = {10.1111/plb.12336}, volume = {18}, year = {2016}, } @article{1241, abstract = {How likely is it that a population escapes extinction through adaptive evolution? The answer to this question is of great relevance in conservation biology, where we aim at species’ rescue and the maintenance of biodiversity, and in agriculture and medicine, where we seek to hamper the emergence of pesticide or drug resistance. By reshuffling the genome, recombination has two antagonistic effects on the probability of evolutionary rescue: It generates and it breaks up favorable gene combinations. Which of the two effects prevails depends on the fitness effects of mutations and on the impact of stochasticity on the allele frequencies. In this article, we analyze a mathematical model for rescue after a sudden environmental change when adaptation is contingent on mutations at two loci. The analysis reveals a complex nonlinear dependence of population survival on recombination. We moreover find that, counterintuitively, a fast eradication of the wild type can promote rescue in the presence of recombination. The model also shows that two-step rescue is not unlikely to happen and can even be more likely than single-step rescue (where adaptation relies on a single mutation), depending on the circumstances.}, author = {Uecker, Hildegard and Hermisson, Joachim}, journal = {Genetics}, number = {2}, pages = {721 -- 732}, publisher = {Genetics Society of America}, title = {{The role of recombination in evolutionary rescue}}, doi = {10.1534/genetics.115.180299}, volume = {202}, year = {2016}, } @inproceedings{1349, abstract = {Crossing fitness valleys is one of the major obstacles to function optimization. In this paper we investigate how the structure of the fitness valley, namely its depth d and length ℓ, influence the runtime of different strategies for crossing these valleys. We present a runtime comparison between the (1+1) EA and two non-elitist nature-inspired algorithms, Strong Selection Weak Mutation (SSWM) and the Metropolis algorithm. While the (1+1) EA has to jump across the valley to a point of higher fitness because it does not accept decreasing moves, the non-elitist algorithms may cross the valley by accepting worsening moves. We show that while the runtime of the (1+1) EA algorithm depends critically on the length of the valley, the runtimes of the non-elitist algorithms depend crucially only on the depth of the valley. In particular, the expected runtime of both SSWM and Metropolis is polynomial in ℓ and exponential in d while the (1+1) EA is efficient only for valleys of small length. Moreover, we show that both SSWM and Metropolis can also efficiently optimize a rugged function consisting of consecutive valleys.}, author = {Oliveto, Pietro and Paixao, Tiago and Heredia, Jorge and Sudholt, Dirk and Trubenova, Barbora}, booktitle = {Proceedings of the Genetic and Evolutionary Computation Conference 2016 }, location = {Denver, CO, USA}, pages = {1163 -- 1170}, publisher = {ACM}, title = {{When non-elitism outperforms elitism for crossing fitness valleys}}, doi = {10.1145/2908812.2908909}, year = {2016}, } @article{1356, author = {Barton, Nicholas H}, journal = {Genetics}, number = {1}, pages = {3 -- 4}, publisher = {Genetics Society of America}, title = {{Sewall Wright on evolution in Mendelian populations and the “Shifting Balance”}}, doi = {10.1534/genetics.115.184796}, volume = {202}, year = {2016}, } @article{1357, author = {Barton, Nicholas H}, journal = {Genetics}, number = {3}, pages = {865 -- 866}, publisher = {Genetics Society of America}, title = {{Richard Hudson and Norman Kaplan on the coalescent process}}, doi = {10.1534/genetics.116.187542}, volume = {202}, year = {2016}, } @article{1358, abstract = {Gene regulation relies on the specificity of transcription factor (TF)–DNA interactions. Limited specificity may lead to crosstalk: a regulatory state in which a gene is either incorrectly activated due to noncognate TF–DNA interactions or remains erroneously inactive. As each TF can have numerous interactions with noncognate cis-regulatory elements, crosstalk is inherently a global problem, yet has previously not been studied as such. We construct a theoretical framework to analyse the effects of global crosstalk on gene regulation. We find that crosstalk presents a significant challenge for organisms with low-specificity TFs, such as metazoans. Crosstalk is not easily mitigated by known regulatory schemes acting at equilibrium, including variants of cooperativity and combinatorial regulation. Our results suggest that crosstalk imposes a previously unexplored global constraint on the functioning and evolution of regulatory networks, which is qualitatively distinct from the known constraints that act at the level of individual gene regulatory elements.}, author = {Friedlander, Tamar and Prizak, Roshan and Guet, Calin C and Barton, Nicholas H and Tkacik, Gasper}, journal = {Nature Communications}, publisher = {Nature Publishing Group}, title = {{Intrinsic limits to gene regulation by global crosstalk}}, doi = {10.1038/ncomms12307}, volume = {7}, year = {2016}, } @article{1359, abstract = {The role of gene interactions in the evolutionary process has long been controversial. Although some argue that they are not of importance, because most variation is additive, others claim that their effect in the long term can be substantial. Here, we focus on the long-term effects of genetic interactions under directional selection assuming no mutation or dominance, and that epistasis is symmetrical overall. We ask by how much the mean of a complex trait can be increased by selection and analyze two extreme regimes, in which either drift or selection dominate the dynamics of allele frequencies. In both scenarios, epistatic interactions affect the long-term response to selection by modulating the additive genetic variance. When drift dominates, we extend Robertson ’ s [Robertson A (1960) Proc R Soc Lond B Biol Sci 153(951):234 − 249] argument to show that, for any form of epistasis, the total response of a haploid population is proportional to the initial total genotypic variance. In contrast, the total response of a diploid population is increased by epistasis, for a given initial genotypic variance. When selection dominates, we show that the total selection response can only be increased by epistasis when s ome initially deleterious alleles become favored as the genetic background changes. We find a sim- ple approximation for this effect and show that, in this regime, it is the structure of the genotype - phenotype map that matters and not the variance components of the population.}, author = {Paixao, Tiago and Barton, Nicholas H}, journal = {PNAS}, number = {16}, pages = {4422 -- 4427}, publisher = {National Academy of Sciences}, title = {{The effect of gene interactions on the long-term response to selection}}, doi = {10.1073/pnas.1518830113}, volume = {113}, year = {2016}, }