@article{12265, author = {Westram, Anja M and Stankowski, Sean and Surendranadh, Parvathy and Barton, Nicholas H}, issn = {1420-9101}, journal = {Journal of Evolutionary Biology}, keywords = {Ecology, Evolution, Behavior and Systematics}, number = {9}, pages = {1200--1205}, publisher = {Wiley}, title = {{Reproductive isolation, speciation, and the value of disagreement: A reply to the commentaries on ‘What is reproductive isolation?’}}, doi = {10.1111/jeb.14082}, volume = {35}, year = {2022}, } @article{10787, abstract = {A species distributed across diverse environments may adapt to local conditions. We ask how quickly such a species changes its range in response to changed conditions. Szép et al. (Szép E, Sachdeva H, Barton NH. 2021 Polygenic local adaptation in metapopulations: a stochastic eco-evolutionary model. Evolution75, 1030–1045 (doi:10.1111/evo.14210)) used the infinite island model to find the stationary distribution of allele frequencies and deme sizes. We extend this to find how a metapopulation responds to changes in carrying capacity, selection strength, or migration rate when deme sizes are fixed. We further develop a ‘fixed-state’ approximation. Under this approximation, polymorphism is only possible for a narrow range of habitat proportions when selection is weak compared to drift, but for a much wider range otherwise. When rates of selection or migration relative to drift change in a single deme of the metapopulation, the population takes a time of order m−1 to reach the new equilibrium. However, even with many loci, there can be substantial fluctuations in net adaptation, because at each locus, alleles randomly get lost or fixed. Thus, in a finite metapopulation, variation may gradually be lost by chance, even if it would persist in an infinite metapopulation. When conditions change across the whole metapopulation, there can be rapid change, which is predicted well by the fixed-state approximation. This work helps towards an understanding of how metapopulations extend their range across diverse environments. This article is part of the theme issue ‘Species’ ranges in the face of changing environments (Part II)’.}, author = {Barton, Nicholas H and Olusanya, Oluwafunmilola O}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society B: Biological Sciences}, keywords = {General Agricultural and Biological Sciences, General Biochemistry, Genetics and Molecular Biology}, number = {1848}, publisher = {The Royal Society}, title = {{The response of a metapopulation to a changing environment}}, doi = {10.1098/rstb.2021.0009}, volume = {377}, year = {2022}, } @article{10658, abstract = {We analyse how migration from a large mainland influences genetic load and population numbers on an island, in a scenario where fitness-affecting variants are unconditionally deleterious, and where numbers decline with increasing load. Our analysis shows that migration can have qualitatively different effects, depending on the total mutation target and fitness effects of deleterious variants. In particular, we find that populations exhibit a genetic Allee effect across a wide range of parameter combinations, when variants are partially recessive, cycling between low-load (large-population) and high-load (sink) states. Increased migration reduces load in the sink state (by increasing heterozygosity) but further inflates load in the large-population state (by hindering purging). We identify various critical parameter thresholds at which one or other stable state collapses, and discuss how these thresholds are influenced by the genetic versus demographic effects of migration. Our analysis is based on a ‘semi-deterministic’ analysis, which accounts for genetic drift but neglects demographic stochasticity. We also compare against simulations which account for both demographic stochasticity and drift. Our results clarify the importance of gene flow as a key determinant of extinction risk in peripheral populations, even in the absence of ecological gradients. This article is part of the theme issue ‘Species’ ranges in the face of changing environments (part I)’.}, author = {Sachdeva, Himani and Olusanya, Oluwafunmilola O and Barton, Nicholas H}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society B}, number = {1846}, publisher = {The Royal Society}, title = {{Genetic load and extinction in peripheral populations: The roles of migration, drift and demographic stochasticity}}, doi = {10.1098/rstb.2021.0010}, volume = {377}, year = {2022}, } @article{11411, abstract = {Many studies have quantified the distribution of heterozygosity and relatedness in natural populations, but few have examined the demographic processes driving these patterns. In this study, we take a novel approach by studying how population structure affects both pairwise identity and the distribution of heterozygosity in a natural population of the self-incompatible plant Antirrhinum majus. Excess variance in heterozygosity between individuals is due to identity disequilibrium, which reflects the variance in inbreeding between individuals; it is measured by the statistic g2. We calculated g2 together with FST and pairwise relatedness (Fij) using 91 SNPs in 22,353 individuals collected over 11 years. We find that pairwise Fij declines rapidly over short spatial scales, and the excess variance in heterozygosity between individuals reflects significant variation in inbreeding. Additionally, we detect an excess of individuals with around half the average heterozygosity, indicating either selfing or matings between close relatives. We use 2 types of simulation to ask whether variation in heterozygosity is consistent with fine-scale spatial population structure. First, by simulating offspring using parents drawn from a range of spatial scales, we show that the known pollen dispersal kernel explains g2. Second, we simulate a 1,000-generation pedigree using the known dispersal and spatial distribution and find that the resulting g2 is consistent with that observed from the field data. In contrast, a simulated population with uniform density underestimates g2, indicating that heterogeneous density promotes identity disequilibrium. Our study shows that heterogeneous density and leptokurtic dispersal can together explain the distribution of heterozygosity.}, author = {Surendranadh, Parvathy and Arathoon, Louise S and Baskett, Carina and Field, David and Pickup, Melinda and Barton, Nicholas H}, issn = {1943-2631}, journal = {Genetics}, number = {3}, publisher = {Oxford University Press}, title = {{Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus}}, doi = {10.1093/genetics/iyac083}, volume = {221}, year = {2022}, } @misc{11321, abstract = {Here are the research data underlying the publication "Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus" Further information are summed up in the README document. }, author = {Surendranadh, Parvathy and Arathoon, Louise S and Baskett, Carina and Field, David and Pickup, Melinda and Barton, Nicholas H}, publisher = {Institute of Science and Technology Austria}, title = {{Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus}}, doi = {10.15479/at:ista:11321}, year = {2022}, } @article{12081, abstract = {Selection accumulates information in the genome—it guides stochastically evolving populations toward states (genotype frequencies) that would be unlikely under neutrality. This can be quantified as the Kullback–Leibler (KL) divergence between the actual distribution of genotype frequencies and the corresponding neutral distribution. First, we show that this population-level information sets an upper bound on the information at the level of genotype and phenotype, limiting how precisely they can be specified by selection. Next, we study how the accumulation and maintenance of information is limited by the cost of selection, measured as the genetic load or the relative fitness variance, both of which we connect to the control-theoretic KL cost of control. The information accumulation rate is upper bounded by the population size times the cost of selection. This bound is very general, and applies across models (Wright–Fisher, Moran, diffusion) and to arbitrary forms of selection, mutation, and recombination. Finally, the cost of maintaining information depends on how it is encoded: Specifying a single allele out of two is expensive, but one bit encoded among many weakly specified loci (as in a polygenic trait) is cheap.}, author = {Hledik, Michal and Barton, Nicholas H and Tkačik, Gašper}, issn = {1091-6490}, journal = {Proceedings of the National Academy of Sciences}, number = {36}, publisher = {Proceedings of the National Academy of Sciences}, title = {{Accumulation and maintenance of information in evolution}}, doi = {10.1073/pnas.2123152119}, volume = {119}, year = {2022}, } @phdthesis{11388, abstract = {In evolve and resequence experiments, a population is sequenced, subjected to selection and then sequenced again, so that genetic changes before and after selection can be observed at the genetic level. Here, I use these studies to better understand the genetic basis of complex traits - traits which depend on more than a few genes. In the first chapter, I discuss the first evolve and resequence experiment, in which a population of mice, the so-called "Longshanks" mice, were selected for tibia length while their body mass was kept constant. The full pedigree is known. We observed a selection response on all chromosomes and used the infinitesimal model with linkage, a model which assumes an infinite number of genes with infinitesimally small effect sizes, as a null model. Results implied a very polygenic basis with a few loci of major effect standing out and changing in parallel. There was large variability between the different chromosomes in this study, probably due to LD. In chapter two, I go on to discuss the impact of LD, on the variability in an allele-frequency based summary statistic, giving an equation based on the initial allele frequencies, average pairwise LD, and the first four moments of the haplotype block copy number distribution. I describe this distribution by referring back to the founder generation. I then demonstrate how to infer selection via a maximum likelihood scheme on the example of a single locus and discuss how to extend this to more realistic scenarios. In chapter three, I discuss the second evolve and resequence experiment, in which a small population of Drosophila melanogaster was selected for increased pupal case size over 6 generations. The experiment was highly replicated with 27 lines selected within family and a known pedigree. We observed a phenotypic selection response of over one standard deviation. I describe the patterns in allele frequency data, including allele frequency changes and patterns of heterozygosity, and give ideas for future work.}, author = {Belohlavy, Stefanie}, isbn = {978-3-99078-018-3}, pages = {98}, publisher = {Institute of Science and Technology Austria}, title = {{The genetic basis of complex traits studied via analysis of evolve and resequence experiments}}, doi = {10.15479/at:ista:11388}, year = {2022}, } @article{10535, abstract = {Realistic models of biological processes typically involve interacting components on multiple scales, driven by changing environment and inherent stochasticity. Such models are often analytically and numerically intractable. We revisit a dynamic maximum entropy method that combines a static maximum entropy with a quasi-stationary approximation. This allows us to reduce stochastic non-equilibrium dynamics expressed by the Fokker-Planck equation to a simpler low-dimensional deterministic dynamics, without the need to track microscopic details. Although the method has been previously applied to a few (rather complicated) applications in population genetics, our main goal here is to explain and to better understand how the method works. We demonstrate the usefulness of the method for two widely studied stochastic problems, highlighting its accuracy in capturing important macroscopic quantities even in rapidly changing non-stationary conditions. For the Ornstein-Uhlenbeck process, the method recovers the exact dynamics whilst for a stochastic island model with migration from other habitats, the approximation retains high macroscopic accuracy under a wide range of scenarios in a dynamic environment.}, author = {Bod'ová, Katarína and Szep, Eniko and Barton, Nicholas H}, issn = {1553-7358}, journal = {PLoS Computational Biology}, number = {12}, publisher = {Public Library of Science}, title = {{Dynamic maximum entropy provides accurate approximation of structured population dynamics}}, doi = {10.1371/journal.pcbi.1009661}, volume = {17}, year = {2021}, } @article{8708, abstract = {The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study ‘replicated’ instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry‐informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi‐stable variants (Dobzhansky‐Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.}, author = {Simon, Alexis and Fraisse, Christelle and El Ayari, Tahani and Liautard‐Haag, Cathy and Strelkov, Petr and Welch, John J and Bierne, Nicolas}, issn = {14209101}, journal = {Journal of Evolutionary Biology}, number = {1}, pages = {208--223}, publisher = {Wiley}, title = {{How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels}}, doi = {10.1111/jeb.13709}, volume = {34}, year = {2021}, } @article{8743, abstract = {Montane cloud forests are areas of high endemism, and are one of the more vulnerable terrestrial ecosystems to climate change. Thus, understanding how they both contribute to the generation of biodiversity, and will respond to ongoing climate change, are important and related challenges. The widely accepted model for montane cloud forest dynamics involves upslope forcing of their range limits with global climate warming. However, limited climate data provides some support for an alternative model, where range limits are forced downslope with climate warming. Testing between these two models is challenging, due to the inherent limitations of climate and pollen records. We overcome this with an alternative source of historical information, testing between competing model predictions using genomic data and demographic analyses for a species of beetle tightly associated to an oceanic island cloud forest. Results unequivocally support the alternative model: populations that were isolated at higher elevation peaks during the Last Glacial Maximum are now in contact and hybridizing at lower elevations. Our results suggest that genomic data are a rich source of information to further understand how montane cloud forest biodiversity originates, and how it is likely to be impacted by ongoing climate change.}, author = {Salces-Castellano, Antonia and Stankowski, Sean and Arribas, Paula and Patino, Jairo and Karger, Dirk N. and Butlin, Roger and Emerson, Brent C.}, issn = {1558-5646}, journal = {Evolution}, number = {2}, pages = {231--244}, publisher = {Wiley}, title = {{Long-term cloud forest response to climate warming revealed by insect speciation history}}, doi = {10.1111/evo.14111}, volume = {75}, year = {2021}, } @article{8928, abstract = {Domestication is a human‐induced selection process that imprints the genomes of domesticated populations over a short evolutionary time scale and that occurs in a given demographic context. Reconstructing historical gene flow, effective population size changes and their timing is therefore of fundamental interest to understand how plant demography and human selection jointly shape genomic divergence during domestication. Yet, the comparison under a single statistical framework of independent domestication histories across different crop species has been little evaluated so far. Thus, it is unclear whether domestication leads to convergent demographic changes that similarly affect crop genomes. To address this question, we used existing and new transcriptome data on three crop species of Solanaceae (eggplant, pepper and tomato), together with their close wild relatives. We fitted twelve demographic models of increasing complexity on the unfolded joint allele frequency spectrum for each wild/crop pair, and we found evidence for both shared and species‐specific demographic processes between species. A convergent history of domestication with gene flow was inferred for all three species, along with evidence of strong reduction in the effective population size during the cultivation stage of tomato and pepper. The absence of any reduction in size of the crop in eggplant stands out from the classical view of the domestication process; as does the existence of a “protracted period” of management before cultivation. Our results also suggest divergent management strategies of modern cultivars among species as their current demography substantially differs. Finally, the timing of domestication is species‐specific and supported by the few historical records available.}, author = {Arnoux, Stéphanie and Fraisse, Christelle and Sauvage, Christopher}, issn = {14209101}, journal = {Journal of Evolutionary Biology}, number = {2}, pages = {270--283}, publisher = {Wiley}, title = {{Genomic inference of complex domestication histories in three Solanaceae species}}, doi = {10.1111/jeb.13723}, volume = {34}, year = {2021}, } @article{9100, abstract = {Marine environments are inhabited by a broad representation of the tree of life, yet our understanding of speciation in marine ecosystems is extremely limited compared with terrestrial and freshwater environments. Developing a more comprehensive picture of speciation in marine environments requires that we 'dive under the surface' by studying a wider range of taxa and ecosystems is necessary for a more comprehensive picture of speciation. Although studying marine evolutionary processes is often challenging, recent technological advances in different fields, from maritime engineering to genomics, are making it increasingly possible to study speciation of marine life forms across diverse ecosystems and taxa. Motivated by recent research in the field, including the 14 contributions in this issue, we highlight and discuss six axes of research that we think will deepen our understanding of speciation in the marine realm: (a) study a broader range of marine environments and organisms; (b) identify the reproductive barriers driving speciation between marine taxa; (c) understand the role of different genomic architectures underlying reproductive isolation; (d) infer the evolutionary history of divergence using model‐based approaches; (e) study patterns of hybridization and introgression between marine taxa; and (f) implement highly interdisciplinary, collaborative research programmes. In outlining these goals, we hope to inspire researchers to continue filling this critical knowledge gap surrounding the origins of marine biodiversity.}, author = {Faria, Rui and Johannesson, Kerstin and Stankowski, Sean}, issn = {14209101}, journal = {Journal of Evolutionary Biology}, number = {1}, pages = {4--15}, publisher = {Wiley}, title = {{Speciation in marine environments: Diving under the surface}}, doi = {10.1111/jeb.13756}, volume = {34}, year = {2021}, } @article{9168, abstract = {Interspecific crossing experiments have shown that sex chromosomes play a major role in reproductive isolation between many pairs of species. However, their ability to act as reproductive barriers, which hamper interspecific genetic exchange, has rarely been evaluated quantitatively compared to Autosomes. This genome-wide limitation of gene flow is essential for understanding the complete separation of species, and thus speciation. Here, we develop a mainland-island model of secondary contact between hybridizing species of an XY (or ZW) sexual system. We obtain theoretical predictions for the frequency of introgressed alleles, and the strength of the barrier to neutral gene flow for the two types of chromosomes carrying multiple interspecific barrier loci. Theoretical predictions are obtained for scenarios where introgressed alleles are rare. We show that the same analytical expressions apply for sex chromosomes and autosomes, but with different sex-averaged effective parameters. The specific features of sex chromosomes (hemizygosity and absence of recombination in the heterogametic sex) lead to reduced levels of introgression on the X (or Z) compared to autosomes. This effect can be enhanced by certain types of sex-biased forces, but it remains overall small (except when alleles causing incompatibilities are recessive). We discuss these predictions in the light of empirical data comprising model-based tests of introgression and cline surveys in various biological systems.}, author = {Fraisse, Christelle and Sachdeva, Himani}, issn = {1943-2631}, journal = {Genetics}, number = {2}, publisher = {Genetics Society of America}, title = {{The rates of introgression and barriers to genetic exchange between hybridizing species: Sex chromosomes vs autosomes}}, doi = {10.1093/genetics/iyaa025}, volume = {217}, year = {2021}, } @article{9119, abstract = {We present DILS, a deployable statistical analysis platform for conducting demographic inferences with linked selection from population genomic data using an Approximate Bayesian Computation framework. DILS takes as input single‐population or two‐population data sets (multilocus fasta sequences) and performs three types of analyses in a hierarchical manner, identifying: (a) the best demographic model to study the importance of gene flow and population size change on the genetic patterns of polymorphism and divergence, (b) the best genomic model to determine whether the effective size Ne and migration rate N, m are heterogeneously distributed along the genome (implying linked selection) and (c) loci in genomic regions most associated with barriers to gene flow. Also available via a Web interface, an objective of DILS is to facilitate collaborative research in speciation genomics. Here, we show the performance and limitations of DILS by using simulations and finally apply the method to published data on a divergence continuum composed by 28 pairs of Mytilus mussel populations/species.}, author = {Fraisse, Christelle and Popovic, Iva and Mazoyer, Clément and Spataro, Bruno and Delmotte, Stéphane and Romiguier, Jonathan and Loire, Étienne and Simon, Alexis and Galtier, Nicolas and Duret, Laurent and Bierne, Nicolas and Vekemans, Xavier and Roux, Camille}, issn = {17550998}, journal = {Molecular Ecology Resources}, pages = {2629--2644}, publisher = {Wiley}, title = {{DILS: Demographic inferences with linked selection by using ABC}}, doi = {10.1111/1755-0998.13323}, volume = {21}, year = {2021}, } @article{9375, abstract = {Genetic variation segregates as linked sets of variants, or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. And yet, genomic data often lack haplotype information, due to constraints in sequencing technologies. Here we present “haplotagging”, a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species the geographic clines for the major wing pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the centre of the hybrid zone. We propose that shared warning signalling (Müllerian mimicry) may couple the cline shifts seen in both species, and facilitate the parallel co-emergence of a novel hybrid morph in both co-mimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.}, author = {Meier, Joana I. and Salazar, Patricio A. and Kučka, Marek and Davies, Robert William and Dréau, Andreea and Aldás, Ismael and Power, Olivia Box and Nadeau, Nicola J. and Bridle, Jon R. and Rolian, Campbell and Barton, Nicholas H and McMillan, W. Owen and Jiggins, Chris D. and Chan, Yingguang Frank}, issn = {0027-8424}, journal = {PNAS}, number = {25}, publisher = {Proceedings of the National Academy of Sciences}, title = {{Haplotype tagging reveals parallel formation of hybrid races in two butterfly species}}, doi = {10.1073/pnas.2015005118}, volume = {118}, year = {2021}, } @article{9394, abstract = {Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.}, author = {Koch, Eva L. and Morales, Hernán E. and Larsson, Jenny and Westram, Anja M and Faria, Rui and Lemmon, Alan R. and Lemmon, E. Moriarty and Johannesson, Kerstin and Butlin, Roger K.}, issn = {2056-3744}, journal = {Evolution Letters}, number = {3}, pages = {196--213}, publisher = {Wiley}, title = {{Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis}}, doi = {10.1002/evl3.227}, volume = {5}, year = {2021}, } @article{9392, abstract = {Humans conceptualize the diversity of life by classifying individuals into types we call ‘species’1. The species we recognize influence political and financial decisions and guide our understanding of how units of diversity evolve and interact. Although the idea of species may seem intuitive, a debate about the best way to define them has raged even before Darwin2. So much energy has been devoted to the so-called ‘species problem’ that no amount of discourse will ever likely solve it2,3. Dozens of species concepts are currently recognized3, but we lack a concrete understanding of how much researchers actually disagree and the factors that cause them to think differently1,2. To address this, we used a survey to quantify the species problem for the first time. The results indicate that the disagreement is extensive: two randomly chosen respondents will most likely disagree on the nature of species. The probability of disagreement is not predicted by researcher experience or broad study system, but tended to be lower among researchers with similar focus, training and who study the same organism. Should we see this diversity of perspectives as a problem? We argue that we should not.}, author = {Stankowski, Sean and Ravinet, Mark}, issn = {18790445}, journal = {Current Biology}, number = {9}, pages = {R428--R429}, publisher = {Cell Press}, title = {{Quantifying the use of species concepts}}, doi = {10.1016/j.cub.2021.03.060}, volume = {31}, year = {2021}, } @misc{12987, abstract = {Chromosomal inversion polymorphisms, segments of chromosomes that are flipped in orientation and occur in reversed order in some individuals, have long been recognized to play an important role in local adaptation. They can reduce recombination in heterozygous individuals and thus help to maintain sets of locally adapted alleles. In a wide range of organisms, populations adapted to different habitats differ in frequency of inversion arrangements. However, getting a full understanding of the importance of inversions for adaptation requires confirmation of their influence on traits under divergent selection. Here, we studied a marine snail, Littorina saxatilis, that has evolved ecotypes adapted to wave exposure or crab predation. These two types occur in close proximity on different parts of the shore. Gene flow between them exists in contact zones. However, they exhibit strong phenotypic divergence in several traits under habitat-specific selection, including size, shape and behaviour. We used crosses between these ecotypes to identify genomic regions that explain variation in these traits by using QTL analysis and variance partitioning across linkage groups. We could show that previously detected inversion regions contribute to adaptive divergence. Some inversions influenced multiple traits suggesting that they contain sets of locally adaptive alleles. Our study also identified regions without known inversions that are important for phenotypic divergence. Thus, we provide a more complete overview of the importance of inversions in relation to the remaining genome.}, author = {Koch, Eva and Morales, Hernán E. and Larsson, Jenny and Westram, Anja M and Faria, Rui and Lemmon, Alan R. and Lemmon, E. Moriarty and Johannesson, Kerstin and Butlin, Roger K.}, publisher = {Dryad}, title = {{Data from: Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis}}, doi = {10.5061/DRYAD.ZGMSBCCB4}, year = {2021}, } @article{9410, abstract = {Antibiotic concentrations vary dramatically in the body and the environment. Hence, understanding the dynamics of resistance evolution along antibiotic concentration gradients is critical for predicting and slowing the emergence and spread of resistance. While it has been shown that increasing the concentration of an antibiotic slows resistance evolution, how adaptation to one antibiotic concentration correlates with fitness at other points along the gradient has not received much attention. Here, we selected populations of Escherichia coli at several points along a concentration gradient for three different antibiotics, asking how rapidly resistance evolved and whether populations became specialized to the antibiotic concentration they were selected on. Populations selected at higher concentrations evolved resistance more slowly but exhibited equal or higher fitness across the whole gradient. Populations selected at lower concentrations evolved resistance rapidly, but overall fitness in the presence of antibiotics was lower. However, these populations readily adapted to higher concentrations upon subsequent selection. Our results indicate that resistance management strategies must account not only for the rates of resistance evolution but also for the fitness of evolved strains.}, author = {Lagator, Mato and Uecker, Hildegard and Neve, Paul}, issn = {1744957X}, journal = {Biology letters}, number = {5}, publisher = {Royal Society of London}, title = {{Adaptation at different points along antibiotic concentration gradients}}, doi = {10.1098/rsbl.2020.0913}, volume = {17}, year = {2021}, } @article{9470, abstract = {A key step in understanding the genetic basis of different evolutionary outcomes (e.g., adaptation) is to determine the roles played by different mutation types (e.g., SNPs, translocations and inversions). To do this we must simultaneously consider different mutation types in an evolutionary framework. Here, we propose a research framework that directly utilizes the most important characteristics of mutations, their population genetic effects, to determine their relative evolutionary significance in a given scenario. We review known population genetic effects of different mutation types and show how these may be connected to different evolutionary outcomes. We provide examples of how to implement this framework and pinpoint areas where more data, theory and synthesis are needed. Linking experimental and theoretical approaches to examine different mutation types simultaneously is a critical step towards understanding their evolutionary significance.}, author = {Berdan, Emma L. and Blanckaert, Alexandre and Slotte, Tanja and Suh, Alexander and Westram, Anja M and Fragata, Inês}, issn = {1365294X}, journal = {Molecular Ecology}, number = {12}, pages = {2710--2723}, publisher = {Wiley}, title = {{Unboxing mutations: Connecting mutation types with evolutionary consequences}}, doi = {10.1111/mec.15936}, volume = {30}, year = {2021}, } @article{9816, abstract = {Aims: Mass antigen testing programs have been challenged because of an alleged insufficient specificity, leading to a large number of false positives. The objective of this study is to derive a lower bound of the specificity of the SD Biosensor Standard Q Ag-Test in large scale practical use. Methods: Based on county data from the nationwide tests for SARS-CoV-2 in Slovakia between 31.10.–1.11. 2020 we calculate a lower confidence bound for the specificity. As positive test results were not systematically verified by PCR tests, we base the lower bound on a worst case assumption, assuming all positives to be false positives. Results: 3,625,332 persons from 79 counties were tested. The lowest positivity rate was observed in the county of Rožňava where 100 out of 34307 (0.29%) tests were positive. This implies a test specificity of at least 99.6% (97.5% one-sided lower confidence bound, adjusted for multiplicity). Conclusion: The obtained lower bound suggests a higher specificity compared to earlier studies in spite of the underlying worst case assumption and the application in a mass testing setting. The actual specificity is expected to exceed 99.6% if the prevalence in the respective regions was non-negligible at the time of testing. To our knowledge, this estimate constitutes the first bound obtained from large scale practical use of an antigen test.}, author = {Hledik, Michal and Polechova, Jitka and Beiglböck, Mathias and Herdina, Anna Nele and Strassl, Robert and Posch, Martin}, issn = {1932-6203}, journal = {PLoS ONE}, number = {7}, publisher = {Public Library of Science}, title = {{Analysis of the specificity of a COVID-19 antigen test in the Slovak mass testing program}}, doi = {10.1371/journal.pone.0255267}, volume = {16}, year = {2021}, } @article{9252, abstract = {This paper analyses the conditions for local adaptation in a metapopulation with infinitely many islands under a model of hard selection, where population size depends on local fitness. Each island belongs to one of two distinct ecological niches or habitats. Fitness is influenced by an additive trait which is under habitat‐dependent directional selection. Our analysis is based on the diffusion approximation and accounts for both genetic drift and demographic stochasticity. By neglecting linkage disequilibria, it yields the joint distribution of allele frequencies and population size on each island. We find that under hard selection, the conditions for local adaptation in a rare habitat are more restrictive for more polygenic traits: even moderate migration load per locus at very many loci is sufficient for population sizes to decline. This further reduces the efficacy of selection at individual loci due to increased drift and because smaller populations are more prone to swamping due to migration, causing a positive feedback between increasing maladaptation and declining population sizes. Our analysis also highlights the importance of demographic stochasticity, which exacerbates the decline in numbers of maladapted populations, leading to population collapse in the rare habitat at significantly lower migration than predicted by deterministic arguments.}, author = {Szep, Eniko and Sachdeva, Himani and Barton, Nicholas H}, issn = {1558-5646}, journal = {Evolution}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics, General Agricultural and Biological Sciences}, number = {5}, pages = {1030--1045}, publisher = {Wiley}, title = {{Polygenic local adaptation in metapopulations: A stochastic eco‐evolutionary model}}, doi = {10.1111/evo.14210}, volume = {75}, year = {2021}, } @article{9374, abstract = {If there are no constraints on the process of speciation, then the number of species might be expected to match the number of available niches and this number might be indefinitely large. One possible constraint is the opportunity for allopatric divergence. In 1981, Felsenstein used a simple and elegant model to ask if there might also be genetic constraints. He showed that progress towards speciation could be described by the build‐up of linkage disequilibrium among divergently selected loci and between these loci and those contributing to other forms of reproductive isolation. Therefore, speciation is opposed by recombination, because it tends to break down linkage disequilibria. Felsenstein then introduced a crucial distinction between “two‐allele” models, which are subject to this effect, and “one‐allele” models, which are free from the recombination constraint. These fundamentally important insights have been the foundation for both empirical and theoretical studies of speciation ever since.}, author = {Butlin, Roger K. and Servedio, Maria R. and Smadja, Carole M. and Bank, Claudia and Barton, Nicholas H and Flaxman, Samuel M. and Giraud, Tatiana and Hopkins, Robin and Larson, Erica L. and Maan, Martine E. and Meier, Joana and Merrill, Richard and Noor, Mohamed A. F. and Ortiz‐Barrientos, Daniel and Qvarnström, Anna}, issn = {1558-5646}, journal = {Evolution}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics, General Agricultural and Biological Sciences}, number = {5}, pages = {978--988}, publisher = {Wiley}, title = {{Homage to Felsenstein 1981, or why are there so few/many species?}}, doi = {10.1111/evo.14235}, volume = {75}, year = {2021}, } @misc{13062, abstract = {This paper analyzes the conditions for local adaptation in a metapopulation with infinitely many islands under a model of hard selection, where population size depends on local fitness. Each island belongs to one of two distinct ecological niches or habitats. Fitness is influenced by an additive trait which is under habitat-dependent directional selection. Our analysis is based on the diffusion approximation and accounts for both genetic drift and demographic stochasticity. By neglecting linkage disequilibria, it yields the joint distribution of allele frequencies and population size on each island. We find that under hard selection, the conditions for local adaptation in a rare habitat are more restrictive for more polygenic traits: even moderate migration load per locus at very many loci is sufficient for population sizes to decline. This further reduces the efficacy of selection at individual loci due to increased drift and because smaller populations are more prone to swamping due to migration, causing a positive feedback between increasing maladaptation and declining population sizes. Our analysis also highlights the importance of demographic stochasticity, which exacerbates the decline in numbers of maladapted populations, leading to population collapse in the rare habitat at significantly lower migration than predicted by deterministic arguments.}, author = {Szep, Eniko and Sachdeva, Himani and Barton, Nicholas H}, publisher = {Dryad}, title = {{Supplementary code for: Polygenic local adaptation in metapopulations: A stochastic eco-evolutionary model}}, doi = {10.5061/DRYAD.8GTHT76P1}, year = {2021}, } @article{9383, abstract = {A primary roadblock to our understanding of speciation is that it usually occurs over a timeframe that is too long to study from start to finish. The idea of a speciation continuum provides something of a solution to this problem; rather than observing the entire process, we can simply reconstruct it from the multitude of speciation events that surround us. But what do we really mean when we talk about the speciation continuum, and can it really help us understand speciation? We explored these questions using a literature review and online survey of speciation researchers. Although most researchers were familiar with the concept and thought it was useful, our survey revealed extensive disagreement about what the speciation continuum actually tells us. This is due partly to the lack of a clear definition. Here, we provide an explicit definition that is compatible with the Biological Species Concept. That is, the speciation continuum is a continuum of reproductive isolation. After outlining the logic of the definition in light of alternatives, we explain why attempts to reconstruct the speciation process from present‐day populations will ultimately fail. We then outline how we think the speciation continuum concept can continue to act as a foundation for understanding the continuum of reproductive isolation that surrounds us.}, author = {Stankowski, Sean and Ravinet, Mark}, issn = {1558-5646}, journal = {Evolution}, number = {6}, pages = {1256--1273}, publisher = {Oxford University Press}, title = {{Defining the speciation continuum}}, doi = {10.1111/evo.14215}, volume = {75}, year = {2021}, } @inbook{14984, abstract = {Hybrid zones are narrow geographic regions where different populations, races or interbreeding species meet and mate, producing mixed ‘hybrid’ offspring. They are relatively common and can be found in a diverse range of organisms and environments. The study of hybrid zones has played an important role in our understanding of the origin of species, with hybrid zones having been described as ‘natural laboratories’. This is because they allow us to study,in situ, the conditions and evolutionary forces that enable divergent taxa to remain distinct despite some ongoing gene exchange between them.}, author = {Stankowski, Sean and Shipilina, Daria and Westram, Anja M}, booktitle = {Encyclopedia of Life Sciences}, isbn = {9780470016176}, publisher = {Wiley}, title = {{Hybrid Zones}}, doi = {10.1002/9780470015902.a0029355}, volume = {2}, year = {2021}, } @misc{9192, abstract = {Here are the research data underlying the publication " Effects of fine-scale population structure on inbreeding in a long-term study of snapdragons (Antirrhinum majus)." Further information are summed up in the README document.}, author = {Surendranadh, Parvathy and Arathoon, Louise S and Baskett, Carina and Field, David and Pickup, Melinda and Barton, Nicholas H}, publisher = {Institute of Science and Technology Austria}, title = {{Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus}}, doi = {10.15479/AT:ISTA:9192}, year = {2021}, } @article{7651, abstract = {The growth of snail shells can be described by simple mathematical rules. Variation in a few parameters can explain much of the diversity of shell shapes seen in nature. However, empirical studies of gastropod shell shape variation typically use geometric morphometric approaches, which do not capture this growth pattern. We have developed a way to infer a set of developmentally descriptive shape parameters based on three-dimensional logarithmic helicospiral growth and using landmarks from two-dimensional shell images as input. We demonstrate the utility of this approach, and compare it to the geometric morphometric approach, using a large set of Littorina saxatilis shells in which locally adapted populations differ in shape. Our method can be modified easily to make it applicable to a wide range of shell forms, which would allow for investigations of the similarities and differences between and within many different species of gastropods.}, author = {Larsson, J. and Westram, Anja M and Bengmark, S. and Lundh, T. and Butlin, R. K.}, issn = {1742-5662}, journal = {Journal of The Royal Society Interface}, number = {163}, publisher = {The Royal Society}, title = {{A developmentally descriptive method for quantifying shape in gastropod shells}}, doi = {10.1098/rsif.2019.0721}, volume = {17}, year = {2020}, } @inbook{9123, abstract = {Inversions are chromosomal rearrangements where the order of genes is reversed. Inversions originate by mutation and can be under positive, negative or balancing selection. Selective effects result from potential disruptive effects on meiosis, gene disruption at inversion breakpoints and, importantly, the effects of inversions as modifiers of recombination rate: Recombination is strongly reduced in individuals heterozygous for an inversion, allowing for alleles at different loci to be inherited as a ‘block’. This may lead to a selective advantage whenever it is favourable to keep certain combinations of alleles associated, for example under local adaptation with gene flow. Inversions can cover a considerable part of a chromosome and contain numerous loci under different selection pressures, so that the resulting overall effects may be complex. Empirical data from various systems show that inversions may have a prominent role in local adaptation, speciation, parallel evolution, the maintenance of polymorphism and sex chromosome evolution.}, author = {Westram, Anja M and Faria, Rui and Butlin, Roger and Johannesson, Kerstin}, booktitle = {eLS}, isbn = {9780470016176}, publisher = {Wiley}, title = {{Inversions and Evolution}}, doi = {10.1002/9780470015902.a0029007}, year = {2020}, } @misc{13073, abstract = {The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study "replicated" instances of secondary contact between closely-related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly-sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi-stable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.}, author = {Simon, Alexis and Fraisse, Christelle and El Ayari, Tahani and Liautard-Haag, Cathy and Strelkov, Petr and Welch, John and Bierne, Nicolas}, publisher = {Dryad}, title = {{How do species barriers decay? concordance and local introgression in mosaic hybrid zones of mussels}}, doi = {10.5061/DRYAD.R4XGXD29N}, year = {2020}, } @misc{13065, abstract = {Domestication is a human-induced selection process that imprints the genomes of domesticated populations over a short evolutionary time scale, and that occurs in a given demographic context. Reconstructing historical gene flow, effective population size changes and their timing is therefore of fundamental interest to understand how plant demography and human selection jointly shape genomic divergence during domestication. Yet, the comparison under a single statistical framework of independent domestication histories across different crop species has been little evaluated so far. Thus, it is unclear whether domestication leads to convergent demographic changes that similarly affect crop genomes. To address this question, we used existing and new transcriptome data on three crop species of Solanaceae (eggplant, pepper and tomato), together with their close wild relatives. We fitted twelve demographic models of increasing complexity on the unfolded joint allele frequency spectrum for each wild/crop pair, and we found evidence for both shared and species-specific demographic processes between species. A convergent history of domestication with gene-flow was inferred for all three species, along with evidence of strong reduction in the effective population size during the cultivation stage of tomato and pepper. The absence of any reduction in size of the crop in eggplant stands out from the classical view of the domestication process; as does the existence of a “protracted period” of management before cultivation. Our results also suggest divergent management strategies of modern cultivars among species as their current demography substantially differs. Finally, the timing of domestication is species-specific and supported by the few historical records available.}, author = {Arnoux, Stephanie and Fraisse, Christelle and Sauvage, Christopher}, publisher = {Dryad}, title = {{VCF files of synonymous SNPs related to: Genomic inference of complex domestication histories in three Solanaceae species}}, doi = {10.5061/DRYAD.Q2BVQ83HD}, year = {2020}, } @article{7995, abstract = {When divergent populations are connected by gene flow, the establishment of complete reproductive isolation usually requires the joint action of multiple barrier effects. One example where multiple barrier effects are coupled consists of a single trait that is under divergent natural selection and also mediates assortative mating. Such multiple‐effect traits can strongly reduce gene flow. However, there are few cases where patterns of assortative mating have been described quantitatively and their impact on gene flow has been determined. Two ecotypes of the coastal marine snail, Littorina saxatilis , occur in North Atlantic rocky‐shore habitats dominated by either crab predation or wave action. There is evidence for divergent natural selection acting on size, and size‐assortative mating has previously been documented. Here, we analyze the mating pattern in L. saxatilis with respect to size in intensively sampled transects across boundaries between the habitats. We show that the mating pattern is mostly conserved between ecotypes and that it generates both assortment and directional sexual selection for small male size. Using simulations, we show that the mating pattern can contribute to reproductive isolation between ecotypes but the barrier to gene flow is likely strengthened more by sexual selection than by assortment.}, author = {Perini, Samuel and Rafajlović, Marina and Westram, Anja M and Johannesson, Kerstin and Butlin, Roger K.}, issn = {15585646}, journal = {Evolution}, number = {7}, pages = {1482--1497}, publisher = {Wiley}, title = {{Assortative mating, sexual selection, and their consequences for gene flow in Littorina}}, doi = {10.1111/evo.14027}, volume = {74}, year = {2020}, } @misc{8809, abstract = {When divergent populations are connected by gene flow, the establishment of complete reproductive isolation usually requires the joint action of multiple barrier effects. One example where multiple barrier effects are coupled consists of a single trait that is under divergent natural selection and also mediates assortative mating. Such multiple-effect traits can strongly reduce gene flow. However, there are few cases where patterns of assortative mating have been described quantitatively and their impact on gene flow has been determined. Two ecotypes of the coastal marine snail, Littorina saxatilis, occur in North Atlantic rocky-shore habitats dominated by either crab predation or wave action. There is evidence for divergent natural selection acting on size, and size-assortative mating has previously been documented. Here, we analyze the mating pattern in L. saxatilis with respect to size in intensively-sampled transects across boundaries between the habitats. We show that the mating pattern is mostly conserved between ecotypes and that it generates both assortment and directional sexual selection for small male size. Using simulations, we show that the mating pattern can contribute to reproductive isolation between ecotypes but the barrier to gene flow is likely strengthened more by sexual selection than by assortment.}, author = {Perini, Samuel and Rafajlovic, Marina and Westram, Anja M and Johannesson, Kerstin and Butlin, Roger}, publisher = {Dryad}, title = {{Data from: Assortative mating, sexual selection and their consequences for gene flow in Littorina}}, doi = {10.5061/dryad.qrfj6q5cn}, year = {2020}, } @article{8112, author = {Barton, Nicholas H}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society. Series B: Biological Sciences}, number = {1806}, publisher = {The Royal Society}, title = {{On the completion of speciation}}, doi = {10.1098/rstb.2019.0530}, volume = {375}, year = {2020}, } @article{8168, abstract = {Speciation, that is, the evolution of reproductive barriers eventually leading to complete isolation, is a crucial process generating biodiversity. Recent work has contributed much to our understanding of how reproductive barriers begin to evolve, and how they are maintained in the face of gene flow. However, little is known about the transition from partial to strong reproductive isolation (RI) and the completion of speciation. We argue that the evolution of strong RI is likely to involve different processes, or new interactions among processes, compared with the evolution of the first reproductive barriers. Transition to strong RI may be brought about by changing external conditions, for example, following secondary contact. However, the increasing levels of RI themselves create opportunities for new barriers to evolve and, and interaction or coupling among barriers. These changing processes may depend on genomic architecture and leave detectable signals in the genome. We outline outstanding questions and suggest more theoretical and empirical work, considering both patterns and processes associated with strong RI, is needed to understand how speciation is completed.}, author = {Kulmuni, Jonna and Butlin, Roger K. and Lucek, Kay and Savolainen, Vincent and Westram, Anja M}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society. Series B: Biological sciences}, number = {1806}, publisher = {The Royal Society}, title = {{Towards the completion of speciation: The evolution of reproductive isolation beyond the first barriers}}, doi = {10.1098/rstb.2019.0528}, volume = {375}, year = {2020}, } @article{8167, abstract = {The evolution of strong reproductive isolation (RI) is fundamental to the origins and maintenance of biological diversity, especially in situations where geographical distributions of taxa broadly overlap. But what is the history behind strong barriers currently acting in sympatry? Using whole-genome sequencing and single nucleotide polymorphism genotyping, we inferred (i) the evolutionary relationships, (ii) the strength of RI, and (iii) the demographic history of divergence between two broadly sympatric taxa of intertidal snail. Despite being cryptic, based on external morphology, Littorina arcana and Littorina saxatilis differ in their mode of female reproduction (egg-laying versus brooding), which may generate a strong post-zygotic barrier. We show that egg-laying and brooding snails are closely related, but genetically distinct. Genotyping of 3092 snails from three locations failed to recover any recent hybrid or backcrossed individuals, confirming that RI is strong. There was, however, evidence for a very low level of asymmetrical introgression, suggesting that isolation remains incomplete. The presence of strong, asymmetrical RI was further supported by demographic analysis of these populations. Although the taxa are currently broadly sympatric, demographic modelling suggests that they initially diverged during a short period of geographical separation involving very low gene flow. Our study suggests that some geographical separation may kick-start the evolution of strong RI, facilitating subsequent coexistence of taxa in sympatry. The strength of RI needed to achieve sympatry and the subsequent effect of sympatry on RI remain open questions.}, author = {Stankowski, Sean and Westram, Anja M and Zagrodzka, Zuzanna B. and Eyres, Isobel and Broquet, Thomas and Johannesson, Kerstin and Butlin, Roger K.}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society. Series B: Biological Sciences}, number = {1806}, publisher = {The Royal Society}, title = {{The evolution of strong reproductive isolation between sympatric intertidal snails}}, doi = {10.1098/rstb.2019.0545}, volume = {375}, year = {2020}, } @article{8169, abstract = {Many recent studies have addressed the mechanisms operating during the early stages of speciation, but surprisingly few studies have tested theoretical predictions on the evolution of strong reproductive isolation (RI). To help address this gap, we first undertook a quantitative review of the hybrid zone literature for flowering plants in relation to reproductive barriers. Then, using Populus as an exemplary model group, we analysed genome-wide variation for phylogenetic tree topologies in both early- and late-stage speciation taxa to determine how these patterns may be related to the genomic architecture of RI. Our plant literature survey revealed variation in barrier complexity and an association between barrier number and introgressive gene flow. Focusing on Populus, our genome-wide analysis of tree topologies in speciating poplar taxa points to unusually complex genomic architectures of RI, consistent with earlier genome-wide association studies. These architectures appear to facilitate the ‘escape’ of introgressed genome segments from polygenic barriers even with strong RI, thus affecting their relationships with recombination rates. Placed within the context of the broader literature, our data illustrate how phylogenomic approaches hold great promise for addressing the evolution and temporary breakdown of RI during late stages of speciation.}, author = {Shang, Huiying and Hess, Jaqueline and Pickup, Melinda and Field, David and Ingvarsson, Pär K. and Liu, Jianquan and Lexer, Christian}, issn = {14712970}, journal = {Philosophical Transactions of the Royal Society. Series B: Biological Sciences}, number = {1806}, publisher = {The Royal Society}, title = {{Evolution of strong reproductive isolation in plants: Broad-scale patterns and lessons from a perennial model group}}, doi = {10.1098/rstb.2019.0544}, volume = {375}, year = {2020}, } @misc{9799, abstract = {Fitness interactions between mutations can influence a population’s evolution in many different ways. While epistatic effects are difficult to measure precisely, important information is captured by the mean and variance of log fitnesses for individuals carrying different numbers of mutations. We derive predictions for these quantities from a class of simple fitness landscapes, based on models of optimizing selection on quantitative traits. We also explore extensions to the models, including modular pleiotropy, variable effect sizes, mutational bias and maladaptation of the wild type. We illustrate our approach by reanalysing a large dataset of mutant effects in a yeast snoRNA. Though characterized by some large epistatic effects, these data give a good overall fit to the non-epistatic null model, suggesting that epistasis might have limited influence on the evolutionary dynamics in this system. We also show how the amount of epistasis depends on both the underlying fitness landscape and the distribution of mutations, and so is expected to vary in consistent ways between new mutations, standing variation and fixed mutations.}, author = {Fraisse, Christelle and Welch, John J.}, publisher = {Royal Society of London}, title = {{Simulation code for Fig S1 from the distribution of epistasis on simple fitness landscapes}}, doi = {10.6084/m9.figshare.7957469.v1}, year = {2020}, } @misc{9798, abstract = {Fitness interactions between mutations can influence a population’s evolution in many different ways. While epistatic effects are difficult to measure precisely, important information is captured by the mean and variance of log fitnesses for individuals carrying different numbers of mutations. We derive predictions for these quantities from a class of simple fitness landscapes, based on models of optimizing selection on quantitative traits. We also explore extensions to the models, including modular pleiotropy, variable effect sizes, mutational bias and maladaptation of the wild type. We illustrate our approach by reanalysing a large dataset of mutant effects in a yeast snoRNA. Though characterized by some large epistatic effects, these data give a good overall fit to the non-epistatic null model, suggesting that epistasis might have limited influence on the evolutionary dynamics in this system. We also show how the amount of epistasis depends on both the underlying fitness landscape and the distribution of mutations, and so is expected to vary in consistent ways between new mutations, standing variation and fixed mutations.}, author = {Fraisse, Christelle and Welch, John J.}, publisher = {Royal Society of London}, title = {{Simulation code for Fig S2 from the distribution of epistasis on simple fitness landscapes}}, doi = {10.6084/m9.figshare.7957472.v1}, year = {2020}, } @article{7236, abstract = {The biotic interactions hypothesis posits that biotic interactions are more important drivers of adaptation closer to the equator, evidenced by “stronger” contemporary interactions (e.g. greater interaction rates) and/or patterns of trait evolution consistent with a history of stronger interactions. Support for the hypothesis is mixed, but few studies span tropical and temperate regions while experimentally controlling for evolutionary history. Here, we integrate field observations and common garden experiments to quantify the relative importance of pollination and herbivory in a pair of tropical‐temperate congeneric perennial herbs. Phytolacca rivinoides and P. americana are pioneer species native to the Neotropics and the eastern USA, respectively. We compared plant‐pollinator and plant‐herbivore interactions between three tropical populations of P. rivinoides from Costa Rica and three temperate populations of P. americana from its northern range edge in Michigan and Ohio. For some metrics of interaction importance, we also included three subtropical populations of P. americana from its southern range edge in Florida. This approach confounds species and region but allows us, uniquely, to measure complementary proxies of interaction importance across a tropical‐temperate range in one system. To test the prediction that lower‐latitude plants are more reliant on insect pollinators, we quantified floral display and reward, insect visitation rates, and self‐pollination ability (autogamy). To test the prediction that lower‐latitude plants experience more herbivore pressure, we quantified herbivory rates, herbivore abundance, and leaf palatability. We found evidence supporting the biotic interactions hypothesis for most comparisons between P. rivinoides and north‐temperate P. americana (floral display, insect visitation, autogamy, herbivory, herbivore abundance, and young‐leaf palatability). Results for subtropical P. americana populations, however, were typically not intermediate between P. rivinoides and north‐temperate P. americana, as would be predicted by a linear latitudinal gradient in interaction importance. Subtropical young‐leaf palatability was intermediate, but subtropical mature leaves were the least palatable, and pollination‐related traits did not differ between temperate and subtropical regions. These nonlinear patterns of interaction importance suggest future work to relate interaction importance to climatic or biotic thresholds. In sum, we found that the biotic interactions hypothesis was more consistently supported at the larger spatial scale of our study.}, author = {Baskett, Carina and Schroeder, Lucy and Weber, Marjorie G. and Schemske, Douglas W.}, issn = {1557-7015}, journal = {Ecological Monographs}, number = {1}, publisher = {Wiley}, title = {{Multiple metrics of latitudinal patterns in insect pollination and herbivory for a tropical‐temperate congener pair}}, doi = {10.1002/ecm.1397}, volume = {90}, year = {2020}, } @article{7205, abstract = {Genetic incompatibilities contribute to reproductive isolation between many diverging populations, but it is still unclear to what extent they play a role if divergence happens with gene flow. In contact zones between the "Crab" and "Wave" ecotypes of the snail Littorina saxatilis, divergent selection forms strong barriers to gene flow, while the role of post‐zygotic barriers due to selection against hybrids remains unclear. High embryo abortion rates in this species could indicate the presence of such barriers. Post‐zygotic barriers might include genetic incompatibilities (e.g. Dobzhansky–Muller incompatibilities) but also maladaptation, both expected to be most pronounced in contact zones. In addition, embryo abortion might reflect physiological stress on females and embryos independent of any genetic stress. We examined all embryos of >500 females sampled outside and inside contact zones of three populations in Sweden. Females' clutch size ranged from 0 to 1,011 embryos (mean 130 ± 123), and abortion rates varied between 0% and 100% (mean 12%). We described female genotypes by using a hybrid index based on hundreds of SNPs differentiated between ecotypes with which we characterized female genotypes. We also calculated female SNP heterozygosity and inversion karyotype. Clutch size did not vary with female hybrid index, and abortion rates were only weakly related to hybrid index in two sites but not at all in a third site. No additional variation in abortion rate was explained by female SNP heterozygosity, but increased female inversion heterozygosity added slightly to increased abortion. Our results show only weak and probably biologically insignificant post‐zygotic barriers contributing to ecotype divergence, and the high and variable abortion rates were marginally, if at all, explained by hybrid index of females.}, author = {Johannesson, Kerstin and Zagrodzka, Zuzanna and Faria, Rui and Westram, Anja M and Butlin, Roger K.}, issn = {14209101}, journal = {Journal of Evolutionary Biology}, number = {3}, pages = {342--351}, publisher = {Wiley}, title = {{Is embryo abortion a post-zygotic barrier to gene flow between Littorina ecotypes?}}, doi = {10.1111/jeb.13570}, volume = {33}, year = {2020}, } @phdthesis{8574, abstract = {This thesis concerns itself with the interactions of evolutionary and ecological forces and the consequences on genetic diversity and the ultimate survival of populations. It is important to understand what signals processes leave on the genome and what we can infer from such data, which is usually abundant but noisy. Furthermore, understanding how and when populations adapt or go extinct is important for practical purposes, such as the genetic management of populations, as well as for theoretical questions, since local adaptation can be the first step toward speciation. In Chapter 2, we introduce the method of maximum entropy to approximate the demographic changes of a population in a simple setting, namely the logistic growth model with immigration. We show that this method is not only a powerful tool in physics but can be gainfully applied in an ecological framework. We investigate how well it approximates the real behavior of the system, and find that is does so, even in unexpected situations. Finally, we illustrate how it can model changing environments. In Chapter 3, we analyze the co-evolution of allele frequencies and population sizes in an infinite island model. We give conditions under which polygenic adaptation to a rare habitat is possible. The model we use is based on the diffusion approximation, considers eco-evolutionary feedback mechanisms (hard selection), and treats both drift and environmental fluctuations explicitly. We also look at limiting scenarios, for which we derive analytical expressions. In Chapter 4, we present a coalescent based simulation tool to obtain patterns of diversity in a spatially explicit subdivided population, in which the demographic history of each subpopulation can be specified. We compare the results to existing predictions, and explore the relative importance of time and space under a variety of spatial arrangements and demographic histories, such as expansion and extinction. In the last chapter, we give a brief outlook to further research. }, author = {Szep, Eniko}, issn = {2663-337X}, pages = {158}, publisher = {Institute of Science and Technology Austria}, title = {{Local adaptation in metapopulations}}, doi = {10.15479/AT:ISTA:8574}, year = {2020}, } @misc{8254, abstract = {Here are the research data underlying the publication "Estimating inbreeding and its effects in a long-term study of snapdragons (Antirrhinum majus)". Further information are summed up in the README document. The files for this record have been updated and are now found in the linked DOI https://doi.org/10.15479/AT:ISTA:9192.}, author = {Arathoon, Louise S}, publisher = {Institute of Science and Technology Austria}, title = {{Estimating inbreeding and its effects in a long-term study of snapdragons (Antirrhinum majus)}}, doi = {10.15479/AT:ISTA:8254}, year = {2020}, } @misc{9839, abstract = {More than 100 years after Grigg’s influential analysis of species’ borders, the causes of limits to species’ ranges still represent a puzzle that has never been understood with clarity. The topic has become especially important recently as many scientists have become interested in the potential for species’ ranges to shift in response to climate change—and yet nearly all of those studies fail to recognise or incorporate evolutionary genetics in a way that relates to theoretical developments. I show that range margins can be understood based on just two measurable parameters: (i) the fitness cost of dispersal—a measure of environmental heterogeneity—and (ii) the strength of genetic drift, which reduces genetic diversity. Together, these two parameters define an ‘expansion threshold’: adaptation fails when genetic drift reduces genetic diversity below that required for adaptation to a heterogeneous environment. When the key parameters drop below this expansion threshold locally, a sharp range margin forms. When they drop below this threshold throughout the species’ range, adaptation collapses everywhere, resulting in either extinction or formation of a fragmented metapopulation. Because the effects of dispersal differ fundamentally with dimension, the second parameter—the strength of genetic drift—is qualitatively different compared to a linear habitat. In two-dimensional habitats, genetic drift becomes effectively independent of selection. It decreases with ‘neighbourhood size’—the number of individuals accessible by dispersal within one generation. Moreover, in contrast to earlier predictions, which neglected evolution of genetic variance and/or stochasticity in two dimensions, dispersal into small marginal populations aids adaptation. This is because the reduction of both genetic and demographic stochasticity has a stronger effect than the cost of dispersal through increased maladaptation. The expansion threshold thus provides a novel, theoretically justified, and testable prediction for formation of the range margin and collapse of the species’ range.}, author = {Polechova, Jitka}, publisher = {Dryad}, title = {{Data from: Is the sky the limit? On the expansion threshold of a species' range}}, doi = {10.5061/dryad.5vv37}, year = {2019}, } @article{5911, abstract = {Empirical data suggest that inversions in many species contain genes important for intraspecific divergence and speciation, yet mechanisms of evolution remain unclear. While genes inside an inversion are tightly linked, inversions are not static but evolve separately from the rest of the genome by new mutations, recombination within arrangements, and gene flux between arrangements. Inversion polymorphisms are maintained by different processes, for example, divergent or balancing selection, or a mix of multiple processes. Moreover, the relative roles of selection, drift, mutation, and recombination will change over the lifetime of an inversion and within its area of distribution. We believe inversions are central to the evolution of many species, but we need many more data and new models to understand the complex mechanisms involved.}, author = {Faria, Rui and Johannesson, Kerstin and Butlin, Roger K. and Westram, Anja M}, issn = {01695347}, journal = {Trends in Ecology and Evolution}, number = {3}, pages = {239--248}, publisher = {Elsevier}, title = {{Evolving inversions}}, doi = {10.1016/j.tree.2018.12.005}, volume = {34}, year = {2019}, } @article{5680, abstract = {Pollinators display a remarkable diversity of foraging strategies with flowering plants, from primarily mutualistic interactions to cheating through nectar robbery. Despite numerous studies on the effect of nectar robbing on components of plant fitness, its contribution to reproductive isolation is unclear. We experimentally tested the impact of different pollinator strategies in a natural hybrid zone between two subspecies of Antirrhinum majus with alternate flower colour guides. On either side of a steep cline in flower colour between Antirrhinum majus pseudomajus (magenta) and A. m. striatum (yellow), we quantified the behaviour of all floral visitors at different time points during the flowering season. Using long-run camera surveys, we quantify the impact of nectar robbing on the number of flowers visited per inflorescence and the flower probing time. We further experimentally tested the effect of nectar robbing on female reproductive success by manipulating the intensity of robbing. While robbing increased over time the number of legitimate visitors tended to decrease concomitantly. We found that the number of flowers pollinated on a focal inflorescence decreased with the number of prior robbing events. However, in the manipulative experiment, fruit set and fruit volume did not vary significantly between low robbing and control treatments. Our findings challenge the idea that robbers have a negative impact on plant fitness through female function. This study also adds to our understanding of the components of pollinator-mediated reproductive isolation and the maintenance of Antirrhinum hybrid zones.}, author = {Andalo, Christophe and Burrus, Monique and Paute, Sandrine and Lauzeral, Christine and Field, David}, issn = {23818115}, journal = {Botany Letters}, number = {1}, pages = {80--92}, publisher = {Taylor and Francis}, title = {{Prevalence of legitimate pollinators and nectar robbers and the consequences for fruit set in an Antirrhinum majus hybrid zone}}, doi = {10.1080/23818107.2018.1545142}, volume = {166}, year = {2019}, } @article{6022, abstract = {The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.}, author = {Merrill, Richard M. and Rastas, Pasi and Martin, Simon H. and Melo Hurtado, Maria C and Barker, Sarah and Davey, John and Mcmillan, W. Owen and Jiggins, Chris D.}, journal = {PLoS Biology}, number = {2}, publisher = {Public Library of Science}, title = {{Genetic dissection of assortative mating behavior}}, doi = {10.1371/journal.pbio.2005902}, volume = {17}, year = {2019}, } @misc{9801, author = {Merrill, Richard M. and Rastas, Pasi and Martin, Simon H. and Melo Hurtado, Maria C and Barker, Sarah and Davey, John and Mcmillan, W. Owen and Jiggins, Chris D.}, publisher = {Public Library of Science}, title = {{Raw behavioral data}}, doi = {10.1371/journal.pbio.2005902.s006}, year = {2019}, } @article{6095, abstract = {Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients.}, author = {Faria, Rui and Chaube, Pragya and Morales, Hernán E. and Larsson, Tomas and Lemmon, Alan R. and Lemmon, Emily M. and Rafajlović, Marina and Panova, Marina and Ravinet, Mark and Johannesson, Kerstin and Westram, Anja M and Butlin, Roger K.}, issn = {1365-294X}, journal = {Molecular Ecology}, number = {6}, pages = {1375--1393}, publisher = {Wiley}, title = {{Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes}}, doi = {10.1111/mec.14972}, volume = {28}, year = {2019}, } @article{6230, abstract = {Great care is needed when interpreting claims about the genetic basis of human variation based on data from genome-wide association studies.}, author = {Barton, Nicholas H and Hermisson, Joachim and Nordborg, Magnus}, issn = {2050084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Why structure matters}}, doi = {10.7554/eLife.45380}, volume = {8}, year = {2019}, }