@article{1169,
  abstract     = {Dispersal is a crucial factor in natural evolution, since it determines the habitat experienced by any population and defines the spatial scale of interactions between individuals. There is compelling evidence for systematic differences in dispersal characteristics within the same population, i.e., genotype-dependent dispersal. The consequences of genotype-dependent dispersal on other evolutionary phenomena, however, are poorly understood. In this article we investigate the effect of genotype-dependent dispersal on spatial gene frequency patterns, using a generalization of the classical diffusion model of selection and dispersal. Dispersal is characterized by the variance of dispersal (diffusion coefficient) and the mean displacement (directional advection term). We demonstrate that genotype-dependent dispersal may change the qualitative behavior of Fisher waves, which change from being “pulled” to being “pushed” wave fronts as the discrepancy in dispersal between genotypes increases. The speed of any wave is partitioned into components due to selection, genotype-dependent variance of dispersal, and genotype-dependent mean displacement. We apply our findings to wave fronts maintained by selection against heterozygotes. Furthermore, we identify a benefit of increased variance of dispersal, quantify its effect on the speed of the wave, and discuss the implications for the evolution of dispersal strategies.},
  author       = {Novak, Sebastian and Kollár, Richard},
  issn         = {00166731},
  journal      = {Genetics},
  number       = {1},
  pages        = {367 -- 374},
  publisher    = {Genetics Society of America},
  title        = {{Spatial gene frequency waves under genotype dependent dispersal}},
  doi          = {10.1534/genetics.116.193946},
  volume       = {205},
  year         = {2017},
}

@article{1191,
  abstract     = {Variation in genotypes may be responsible for differences in dispersal rates, directional biases, and growth rates of individuals. These traits may favor certain genotypes and enhance their spatiotemporal spreading into areas occupied by the less advantageous genotypes. We study how these factors influence the speed of spreading in the case of two competing genotypes under the assumption that spatial variation of the total population is small compared to the spatial variation of the frequencies of the genotypes in the population. In that case, the dynamics of the frequency of one of the genotypes is approximately described by a generalized Fisher–Kolmogorov–Petrovskii–Piskunov (F–KPP) equation. This generalized F–KPP equation with (nonlinear) frequency-dependent diffusion and advection terms admits traveling wave solutions that characterize the invasion of the dominant genotype. Our existence results generalize the classical theory for traveling waves for the F–KPP with constant coefficients. Moreover, in the particular case of the quadratic (monostable) nonlinear growth–decay rate in the generalized F–KPP we study in detail the influence of the variance in diffusion and mean displacement rates of the two genotypes on the minimal wave propagation speed.},
  author       = {Kollár, Richard and Novak, Sebastian},
  journal      = {Bulletin of Mathematical Biology},
  number       = {3},
  pages        = {525--559},
  publisher    = {Springer},
  title        = {{Existence of traveling waves for the generalized F–KPP equation}},
  doi          = {10.1007/s11538-016-0244-3},
  volume       = {79},
  year         = {2017},
}

@article{1199,
  abstract     = {Much of quantitative genetics is based on the ‘infinitesimal model’, under which selection has a negligible effect on the genetic variance. This is typically justified by assuming a very large number of loci with additive effects. However, it applies even when genes interact, provided that the number of loci is large enough that selection on each of them is weak relative to random drift. In the long term, directional selection will change allele frequencies, but even then, the effects of epistasis on the ultimate change in trait mean due to selection may be modest. Stabilising selection can maintain many traits close to their optima, even when the underlying alleles are weakly selected. However, the number of traits that can be optimised is apparently limited to ~4Ne by the ‘drift load’, and this is hard to reconcile with the apparent complexity of many organisms. Just as for the mutation load, this limit can be evaded by a particular form of negative epistasis. A more robust limit is set by the variance in reproductive success. This suggests that selection accumulates information most efficiently in the infinitesimal regime, when selection on individual alleles is weak, and comparable with random drift. A review of evidence on selection strength suggests that although most variance in fitness may be because of alleles with large Nes, substantial amounts of adaptation may be because of alleles in the infinitesimal regime, in which epistasis has modest effects.},
  author       = {Barton, Nicholas H},
  journal      = {Heredity},
  pages        = {96 -- 109},
  publisher    = {Nature Publishing Group},
  title        = {{How does epistasis influence the response to selection?}},
  doi          = {10.1038/hdy.2016.109},
  volume       = {118},
  year         = {2017},
}

@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{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},
}

@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},
}

@article{1382,
  abstract     = {Background and aims Angiosperms display remarkable diversity in flower colour, implying that transitions between pigmentation phenotypes must have been common. Despite progress in understanding transitions between anthocyanin (blue, purple, pink or red) and unpigmented (white) flowers, little is known about the evolutionary patterns of flower-colour transitions in lineages with both yellow and anthocyanin-pigmented flowers. This study investigates the relative rates of evolutionary transitions between different combinations of yellow- and anthocyanin-pigmentation phenotypes in the tribe Antirrhineae. Methods We surveyed taxonomic literature for data on anthocyanin and yellow floral pigmentation for 369 species across the tribe. We then reconstructed the phylogeny of 169 taxa and used phylogenetic comparative methods to estimate transition rates among pigmentation phenotypes across the phylogeny. Key Results In contrast to previous studies we found a bias towards transitions involving a gain in pigmentation, although transitions to phenotypes with both anthocyanin and yellow taxa are nevertheless extremely rare. Despite the dominance of yellow and anthocyanin-pigmented taxa, transitions between these phenotypes are constrained to move through a white intermediate stage, whereas transitions to double-pigmentation are very rare. The most abundant transitions are between anthocyanin-pigmented and unpigmented flowers, and similarly the most abundant polymorphic taxa were those with anthocyanin-pigmented and unpigmented flowers. Conclusions Our findings show that pigment evolution is limited by the presence of other floral pigments. This interaction between anthocyanin and yellow pigments constrains the breadth of potential floral diversity observed in nature. In particular, they suggest that selection has repeatedly acted to promote the spread of single-pigmented phenotypes across the Antirrhineae phylogeny. Furthermore, the correlation between transition rates and polymorphism suggests that the forces causing and maintaining variance in the short term reflect evolutionary processes on longer time scales.},
  author       = {Ellis, Thomas and Field, David},
  journal      = {Annals of Botany},
  number       = {7},
  pages        = {1133 -- 1140},
  publisher    = {Oxford University Press},
  title        = {{Repeated gains in yellow and anthocyanin pigmentation in flower colour transitions in the Antirrhineae}},
  doi          = {10.1093/aob/mcw043},
  volume       = {117},
  year         = {2016},
}

@phdthesis{1398,
  abstract     = {Hybrid zones represent evolutionary laboratories, where recombination brings together alleles in combinations which have not previously been tested by selection. This provides an excellent opportunity to test the effect of molecular variation on fitness, and how this variation is able to spread through populations in a natural context. The snapdragon Antirrhinum majus is polymorphic in the wild for two loci controlling the distribution of yellow and magenta floral pigments. Where the yellow A. m. striatum and the magenta A. m. pseudomajus meet along a valley in the Spanish Pyrenees they form a stable hybrid zone Alleles at these loci recombine to give striking transgressive variation for flower colour. The sharp transition in phenotype over ~1km implies strong selection maintaining the hybrid zone. An indirect assay of pollinator visitation in the field found that pollinators forage in a positive-frequency dependent manner on Antirrhinum, matching previous data on fruit set. Experimental arrays and paternity analysis of wild-pollinated seeds demonstrated assortative mating for pigmentation alleles, and that pollinator behaviour alone is sufficient to explain this pattern. Selection by pollinators should be sufficiently strong to maintain the hybrid zone, although other mechanisms may be at work. At a broader scale I examined evolutionary transitions between yellow and anthocyanin pigmentation in the tribe Antirrhinae, and found that selection has acted strate that pollinators are a major determinant of reproductive success and mating patterns in wild Antirrhinum.},
  author       = {Ellis, Thomas},
  issn         = {2663-337X},
  pages        = {130},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{The role of pollinator-mediated selection in the maintenance of a flower color polymorphism in an Antirrhinum majus hybrid zone}},
  doi          = {10.15479/AT:ISTA:TH_526 },
  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μ &lt; 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},
}

@misc{5550,
  abstract     = {We collected flower colour information on species in the tribe Antirrhineae from taxonomic literature. We also retreived molecular data from GenBank for as many of these species as possible to estimate phylogenetic relationships among these taxa. We then used the R package 'diversitree' to examine patterns of evolutionary transitions between anthocyanin and yellow pigmentation across the phylogeny.

For full details of the methods see:
Ellis TJ and Field DL "Repeated gains in yellow and anthocyanin pigmentation in flower colour transitions in the Antirrhineae”, Annals of Botany (in press)},
  author       = {Ellis, Thomas and Field, David},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Flower colour data and phylogeny (NEXUS) files}},
  doi          = {10.15479/AT:ISTA:34},
  year         = {2016},
}

@misc{5551,
  abstract     = {Data from array experiments investigating pollinator behaviour on snapdragons in controlled conditions, and their effect on plant mating. Data were collected as part of Tom Ellis' PhD thesis , submitted February 2016.

We placed a total of 36 plants in a grid inside a closed organza tent, with a single hive of commercially bred bumblebees (Bombus hortorum). We used only the yellow-flowered Antirrhinum majus striatum and the magenta-flowered Antirrhinum majus pseudomajus, at ratios of 6:36, 12:24, 18:18, 24:12 and 30:6.

After 24 hours to learn how to deal with snapdragons, I observed pollinators foraging on plants, and recorded the transitions between plants. Thereafter seeds on plants were allowed to develops. A sample of these were grown to maturity when their flower colour could be determined, and they were scored as yellow, magenta, or hybrid.},
  author       = {Ellis, Thomas},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Data on pollinator observations and offpsring phenotypes}},
  doi          = {10.15479/AT:ISTA:35},
  year         = {2016},
}

@misc{5552,
  abstract     = {Data on pollinator visitation to wild snapdragons in a natural hybrid zone, collected as part of Tom Ellis' PhD thesis (submitted February 2016).

Snapdragon flowers have a mouth-like structure which pollinators must open to access nectar. We placed 5mm cellophane tags in these mouths, which are held in place by the pressure of the flower until a pollinator visits. When she opens the flower, the tag drops out, and one can infer a visit. We surveyed plants over multiple days in 2010, 2011 and 2012.

Also included are data on phenotypic and demographic variables which may be explanatory variables for pollinator visitation.},
  author       = {Ellis, Thomas},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Pollinator visitation data for wild Antirrhinum majus plants, with phenotypic and frequency data.}},
  doi          = {10.15479/AT:ISTA:36},
  year         = {2016},
}

@misc{5553,
  abstract     = {Genotypic, phenotypic and demographic data for 2128 wild snapdragons and 1127 open-pollinated progeny from a natural hybrid zone, collected as part of Tom Ellis' PhD thesis (submitted) February 2016).

Tissue samples were sent to LGC Genomics in Berlin for DNA extraction, and genotyping at 70 SNP markers by KASPR genotyping. 29 of these SNPs failed to amplify reliably, and have been removed from this dataset.

Other data were retreived from an online database of this population at www.antspec.org.},
  author       = {Field, David and Ellis, Thomas},
  keywords     = {paternity assignment, pedigree, matting patterns, assortative mating, Antirrhinum majus, frequency-dependent selection, plant-pollinator interaction},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Inference of mating patterns among wild snapdragons in a natural hybrid zone in 2012}},
  doi          = {10.15479/AT:ISTA:37},
  year         = {2016},
}

@misc{5554,
  abstract     = {The data stored here is used in Murat Tugrul's PhD thesis (Chapter 3), which is related to the evolution of bacterial RNA polymerase binding.
Magdalena Steinrueck (PhD Student in Calin Guet's group at IST Austria) performed the experiments and created the data on de novo promoter evolution. Fabienne Jesse (PhD Student in Jon Bollback's group at IST Austria) performed the experiments and created the data on lac promoter evolution.},
  author       = {Tugrul, Murat},
  keywords     = {RNAP binding, de novo promoter evolution, lac promoter},
  publisher    = {Institute of Science and Technology Austria},
  title        = {{Experimental Data for Binding Site Evolution of Bacterial RNA Polymerase}},
  doi          = {10.15479/AT:ISTA:43},
  year         = {2016},
}

