@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}, } @article{910, abstract = {Frequency-independent selection is generally considered as a force that acts to reduce the genetic variation in evolving populations, yet rigorous arguments for this idea are scarce. When selection fluctuates in time, it is unclear whether frequency-independent selection may maintain genetic polymorphism without invoking additional mechanisms. We show that constant frequency-independent selection with arbitrary epistasis on a well-mixed haploid population eliminates genetic variation if we assume linkage equilibrium between alleles. To this end, we introduce the notion of frequency-independent selection at the level of alleles, which is sufficient to prove our claim and contains the notion of frequency-independent selection on haploids. When selection and recombination are weak but of the same order, there may be strong linkage disequilibrium; numerical calculations show that stable equilibria are highly unlikely. Using the example of a diallelic two-locus model, we then demonstrate that frequency-independent selection that fluctuates in time can maintain stable polymorphism if linkage disequilibrium changes its sign periodically. We put our findings in the context of results from the existing literature and point out those scenarios in which the possible role of frequency-independent selection in maintaining genetic variation remains unclear. }, author = {Novak, Sebastian and Barton, Nicholas H}, journal = {Genetics}, number = {2}, pages = {653 -- 668}, publisher = {Genetics Society of America}, title = {{When does frequency-independent selection maintain genetic variation?}}, doi = {10.1534/genetics.117.300129}, volume = {207}, year = {2017}, } @article{614, abstract = {Moths and butterflies (Lepidoptera) usually have a pair of differentiated WZ sex chromosomes. However, in most lineages outside of the division Ditrysia, as well as in the sister order Trichoptera, females lack a W chromosome. The W is therefore thought to have been acquired secondarily. Here we compare the genomes of three Lepidoptera species (one Dytrisia and two non-Dytrisia) to test three models accounting for the origin of the W: (1) a Z-autosome fusion; (2) a sex chromosome turnover; and (3) a non-canonical mechanism (e.g., through the recruitment of a B chromosome). We show that the gene content of the Z is highly conserved across Lepidoptera (rejecting a sex chromosome turnover) and that very few genes moved onto the Z in the common ancestor of the Ditrysia (arguing against a Z-autosome fusion). Our comparative genomics analysis therefore supports the secondary acquisition of the Lepidoptera W by a non-canonical mechanism, and it confirms the extreme stability of well-differentiated sex chromosomes.}, author = {Fraisse, Christelle and Picard, Marion A and Vicoso, Beatriz}, issn = {20411723}, journal = {Nature Communications}, number = {1}, publisher = {Nature Publishing Group}, title = {{The deep conservation of the Lepidoptera Z chromosome suggests a non canonical origin of the W}}, doi = {10.1038/s41467-017-01663-5}, volume = {8}, year = {2017}, } @misc{7163, abstract = {The de novo genome assemblies generated for this study, and the associated metadata.}, author = {Fraisse, Christelle}, publisher = {Institute of Science and Technology Austria}, title = {{Supplementary Files for "The deep conservation of the Lepidoptera Z chromosome suggests a non canonical origin of the W"}}, doi = {10.15479/AT:ISTA:7163}, 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 = {1553734X}, 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{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{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}, } @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{1409, author = {Abbott, Richard and Barton, Nicholas H and Good, Jeffrey}, journal = {Molecular Ecology}, number = {11}, pages = {2325 -- 2332}, publisher = {Wiley-Blackwell}, title = {{Genomics of hybridization and its evolutionary consequences}}, doi = {10.1111/mec.13685}, volume = {25}, year = {2016}, } @article{1420, abstract = {Selection, mutation, and random drift affect the dynamics of allele frequencies and consequently of quantitative traits. While the macroscopic dynamics of quantitative traits can be measured, the underlying allele frequencies are typically unobserved. Can we understand how the macroscopic observables evolve without following these microscopic processes? This problem has been studied previously by analogy with statistical mechanics: the allele frequency distribution at each time point is approximated by the stationary form, which maximizes entropy. We explore the limitations of this method when mutation is small (4Nμ < 1) so that populations are typically close to fixation, and we extend the theory in this regime to account for changes in mutation strength. We consider a single diallelic locus either under directional selection or with overdominance and then generalize to multiple unlinked biallelic loci with unequal effects. We find that the maximum-entropy approximation is remarkably accurate, even when mutation and selection change rapidly. }, author = {Bod'ová, Katarína and Tkacik, Gasper and Barton, Nicholas H}, journal = {Genetics}, number = {4}, pages = {1523 -- 1548}, publisher = {Genetics Society of America}, title = {{A general approximation for the dynamics of quantitative traits}}, doi = {10.1534/genetics.115.184127}, volume = {202}, year = {2016}, } @article{1518, abstract = {The inference of demographic history from genome data is hindered by a lack of efficient computational approaches. In particular, it has proved difficult to exploit the information contained in the distribution of genealogies across the genome. We have previously shown that the generating function (GF) of genealogies can be used to analytically compute likelihoods of demographic models from configurations of mutations in short sequence blocks (Lohse et al. 2011). Although the GF has a simple, recursive form, the size of such likelihood calculations explodes quickly with the number of individuals and applications of this framework have so far been mainly limited to small samples (pairs and triplets) for which the GF can be written by hand. Here we investigate several strategies for exploiting the inherent symmetries of the coalescent. In particular, we show that the GF of genealogies can be decomposed into a set of equivalence classes that allows likelihood calculations from nontrivial samples. Using this strategy, we automated blockwise likelihood calculations for a general set of demographic scenarios in Mathematica. These histories may involve population size changes, continuous migration, discrete divergence, and admixture between multiple populations. To give a concrete example, we calculate the likelihood for a model of isolation with migration (IM), assuming two diploid samples without phase and outgroup information. We demonstrate the new inference scheme with an analysis of two individual butterfly genomes from the sister species Heliconius melpomene rosina and H. cydno.}, author = {Lohse, Konrad and Chmelik, Martin and Martin, Simon and Barton, Nicholas H}, journal = {Genetics}, number = {2}, pages = {775 -- 786}, publisher = {Genetics Society of America}, title = {{Efficient strategies for calculating blockwise likelihoods under the coalescent}}, doi = {10.1534/genetics.115.183814}, volume = {202}, year = {2016}, } @article{1631, abstract = {Ancestral processes are fundamental to modern population genetics and spatial structure has been the subject of intense interest for many years. Despite this interest, almost nothing is known about the distribution of the locations of pedigree or genetic ancestors. Using both spatially continuous and stepping-stone models, we show that the distribution of pedigree ancestors approaches a travelling wave, for which we develop two alternative approximations. The speed and width of the wave are sensitive to the local details of the model. After a short time, genetic ancestors spread far more slowly than pedigree ancestors, ultimately diffusing out with radius ## rather than spreading at constant speed. In contrast to the wave of pedigree ancestors, the spread of genetic ancestry is insensitive to the local details of the models.}, author = {Kelleher, Jerome and Etheridge, Alison and Véber, Amandine and Barton, Nicholas H}, journal = {Theoretical Population Biology}, pages = {1 -- 12}, publisher = {Academic Press}, title = {{Spread of pedigree versus genetic ancestry in spatially distributed populations}}, doi = {10.1016/j.tpb.2015.10.008}, volume = {108}, year = {2016}, } @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}, } @misc{9862, author = {Roux, Camille and Fraisse, Christelle and Romiguier, Jonathan and Anciaux, Youann and Galtier, Nicolas and Bierne, Nicolas}, publisher = {Public Library of Science}, title = {{Simulation study to test the robustness of ABC in face of recent times of divergence}}, doi = {10.1371/journal.pbio.2000234.s016}, year = {2016}, } @misc{9863, author = {Roux, Camille and Fraisse, Christelle and Romiguier, Jonathan and Anciaux, Youann and Galtier, Nicolas and Bierne, Nicolas}, publisher = {Public Library of Science}, title = {{Accessions of surveyed individuals, geographic locations and summary statistics}}, doi = {10.1371/journal.pbio.2000234.s017}, year = {2016}, }