TY - JOUR AB - Understanding the factors that have shaped the current distributions and diversity of species is a central and longstanding aim of evolutionary biology. The recent inclusion of genomic data into phylogeographic studies has dramatically improved our understanding in organisms where evolutionary relationships have been challenging to infer. We used whole-genome sequences to study the phylogeography of the intertidal snail Littorina saxatilis, which has successfully colonized and diversified across a broad range of coastal environments in the Northern Hemisphere amid repeated cycles of glaciation. Building on past studies based on short DNA sequences, we used genome-wide data to provide a clearer picture of the relationships among samples spanning most of the species natural range. Our results confirm the trans-Atlantic colonization of North America from Europe, and have allowed us to identify rough locations of glacial refugia and to infer likely routes of colonization within Europe. We also investigated the signals in different datasets to account for the effects of genomic architecture and non-neutral evolution, which provides new insights about diversification of four ecotypes of L. saxatilis (the crab, wave, barnacle, and brackish ecotypes) at different spatial scales. Overall, we provide a much clearer picture of the biogeography of L. saxatilis, providing a foundation for more detailed phylogenomic and demographic studies. AU - Stankowski, Sean AU - Zagrodzka, Zuzanna B AU - Galindo, Juan AU - Montaño-Rendón, Mauricio AU - Faria, Rui AU - Mikhailova, Natalia AU - Blakeslee, April M H AU - Arnason, Einar AU - Broquet, Thomas AU - Morales, Hernán E AU - Grahame, John W AU - Westram, Anja M AU - Johannesson, Kerstin AU - Butlin, Roger K ID - 14833 IS - 1 JF - Evolutionary Journal of the Linnean Society TI - Whole-genome phylogeography of the intertidal snail Littorina saxatilis VL - 2 ER - TY - GEN AB - Fragmented landscapes pose a significant threat to the persistence of species as they are highly susceptible to heightened risk of extinction due to the combined effects of genetic and demographic factors such as genetic drift and demographic stochasticity. This paper explores the intricate interplay between genetic load and extinction risk within metapopulations with a focus on understanding the impact of eco-evolutionary feedback mechanisms. We distinguish between two models of selection: soft selection, characterised by subpopulations maintaining carrying capacity despite load, and hard selection, where load can significantly affect population size. Within the soft selection framework, we investigate the impact of gene flow on genetic load at a single locus, while also considering the effect of selection strength and dominance coefficient. We subsequently build on this to examine how gene flow influences both population size and load under hard selection as well as identify critical thresholds for metapopulation persistence. Our analysis employs the diffusion, semi-deterministic and effective migration approximations. Our findings reveal that under soft selection, even modest levels of migration can significantly alleviate the burden of load. In sharp contrast, with hard selection, a much higher degree of gene flow is required to mitigate load and prevent the collapse of the metapopulation. Overall, this study sheds light into the crucial role migration plays in shaping the dynamics of genetic load and extinction risk in fragmented landscapes, offering valuable insights for conservation strategies and the preservation of diversity in a changing world. AU - Olusanya, Oluwafunmilola O AU - Khudiakova, Kseniia AU - Sachdeva, Himani ID - 14732 T2 - bioRxiv TI - Genetic load, eco-evolutionary feedback and extinction in a metapopulation ER - TY - GEN AB - This repository contains the code and VCF files needed to conduct the analyses in our MS. Each folder contains a readMe document explaining the nature of each file and dataset and the results and analyses that they relate to. The same anlaysis code (but not VCF files) is also available at https://github.com/seanstankowski/Littorina_reproductive_mode AU - Stankowski, Sean ID - 14812 TI - Data and code for: The genetic architecture of a recent transition to live-bearing in marine snails ER - TY - THES AB - The evolutionary processes that brought about today’s plethora of living species and the many billions more ancient ones all underlie biology. Evolutionary pathways are neither directed nor deterministic, but rather an interplay between selection, migration, mutation, genetic drift and other environmental factors. Hybrid zones, as natural crossing experiments, offer a great opportunity to use cline analysis to deduce different evolutionary processes - for example, selection strength. Theoretical cline models, largely assuming uniform distribution of individuals, often lack the capability of incorporating population structure. Since in reality organisms mostly live in patchy distributions and their dispersal is hardly ever Gaussian, it is necessary to unravel the effect of these different elements of population structure on cline parameters and shape. In this thesis, I develop a simulation inspired by the A. majus hybrid zone of a single selected locus under frequency dependent selection. This simulation enables us to untangle the effects of different elements of population structure as for example a low-density center and long-range dispersal. This thesis is therefore a first step towards theoretically untangling the effects of different elements of population structure on cline parameters and shape. AU - Julseth, Mara ID - 12800 SN - 2791-4585 TI - The effect of local population structure on genetic variation at selected loci in the A. majus hybrid zone ER - TY - JOUR AB - When Mendel’s work was rediscovered in 1900, and extended to establish classical genetics, it was initially seen in opposition to Darwin’s theory of evolution by natural selection on continuous variation, as represented by the biometric research program that was the foundation of quantitative genetics. As Fisher, Haldane, and Wright established a century ago, Mendelian inheritance is exactly what is needed for natural selection to work efficiently. Yet, the synthesis remains unfinished. We do not understand why sexual reproduction and a fair meiosis predominate in eukaryotes, or how far these are responsible for their diversity and complexity. Moreover, although quantitative geneticists have long known that adaptive variation is highly polygenic, and that this is essential for efficient selection, this is only now becoming appreciated by molecular biologists—and we still do not have a good framework for understanding polygenic variation or diffuse function. AU - Barton, Nicholas H ID - 11702 IS - 30 JF - Proceedings of the National Academy of Sciences of the United States of America SN - 0027-8424 TI - The "New Synthesis" VL - 119 ER - TY - THES AB - Although we often see studies focusing on simple or even discrete traits in studies of colouration, the variation of “appearance” phenotypes found in nature is often more complex, continuous and high-dimensional. Therefore, we developed automated methods suitable for large datasets of genomes and images, striving to account for their complex nature, while minimising human bias. We used these methods on a dataset of more than 20, 000 plant SNP genomes and corresponding fower images from a hybrid zone of two subspecies of Antirrhinum majus with distinctly coloured fowers to improve our understanding of the genetic nature of the fower colour in our study system. Firstly, we use the advantage of large numbers of genotyped plants to estimate the haplotypes in the main fower colour regulating region. We study colour- and geography-related characteristics of the estimated haplotypes and how they connect to their relatedness. We show discrepancies from the expected fower colour distributions given the genotype and identify particular haplotypes leading to unexpected phenotypes. We also confrm a signifcant defcit of the double recessive recombinant and quite surprisingly, we show that haplotypes of the most frequent parental type are much less variable than others. Secondly, we introduce our pipeline capable of processing tens of thousands of full fower images without human interaction and summarising each image into a set of informative scores. We show the compatibility of these machine-measured fower colour scores with the previously used manual scores and study impact of external efect on the resulting scores. Finally, we use the machine-measured fower colour scores to ft and examine a phenotype cline across the hybrid zone in Planoles using full fower images as opposed to discrete, manual scores and compare it with the genotypic cline. AU - Matejovicova, Lenka ID - 11128 SN - 2663-337X TI - Genetic basis of flower colour as a model for adaptive evolution ER - TY - JOUR AB - Maternally inherited Wolbachia transinfections are being introduced into natural mosquito populations to reduce the transmission of dengue, Zika, and other arboviruses. Wolbachia-induced cytoplasmic incompatibility provides a frequency-dependent reproductive advantage to infected females that can spread transinfections within and among populations. However, because transinfections generally reduce host fitness, they tend to spread within populations only after their frequency exceeds a critical threshold. This produces bistability with stable equilibrium frequencies at both 0 and 1, analogous to the bistability produced by underdominance between alleles or karyotypes and by population dynamics under Allee effects. Here, we analyze how stochastic frequency variation produced by finite population size can facilitate the local spread of variants with bistable dynamics into areas where invasion is unexpected from deterministic models. Our exemplar is the establishment of wMel Wolbachia in the Aedes aegypti population of Pyramid Estates (PE), a small community in far north Queensland, Australia. In 2011, wMel was stably introduced into Gordonvale, separated from PE by barriers to A. aegypti dispersal. After nearly 6 years during which wMel was observed only at low frequencies in PE, corresponding to an apparent equilibrium between immigration and selection, wMel rose to fixation by 2018. Using analytic approximations and statistical analyses, we demonstrate that the observed fixation of wMel at PE is consistent with both stochastic transition past an unstable threshold frequency and deterministic transformation produced by steady immigration at a rate just above the threshold required for deterministic invasion. The indeterminacy results from a delicate balance of parameters needed to produce the delayed transition observed. Our analyses suggest that once Wolbachia transinfections are established locally through systematic introductions, stochastic “threshold crossing” is likely to only minimally enhance spatial spread, providing a local ratchet that slightly—but systematically—aids area-wide transformation of disease-vector populations in heterogeneous landscapes. AU - Turelli, Michael AU - Barton, Nicholas H ID - 10604 IS - 1 JF - Evolution Letters KW - genetics KW - ecology KW - evolution KW - behavior and systematics TI - Why did the Wolbachia transinfection cross the road? Drift, deterministic dynamics, and disease control VL - 6 ER - TY - GEN AB - Maternally inherited Wolbachia transinfections are being introduced into natural mosquito populations to reduce the transmission of dengue, Zika and other arboviruses. Wolbachia-induced cytoplasmic incompatibility provides a frequency-dependent reproductive advantage to infected females that can spread transinfections within and among populations. However, because transinfections generally reduce host fitness, they tend to spread within populations only after their frequency exceeds a critical threshold. This produces bistability with stable equilibrium frequencies at both 0 and 1, analogous to the bistability produced by underdominance between alleles or karyotypes and by population dynamics under Allee effects. Here, we analyze how stochastic frequency variation produced by finite population size can facilitate the local spread of variants with bistable dynamics into areas where invasion is unexpected from deterministic models. Our exemplar is the establishment of wMel Wolbachia in the Aedes aegypti population of Pyramid Estates (PE), a small community in far north Queensland, Australia. In 2011, wMel was stably introduced into Gordonvale, separated from PE by barriers to Ae. aegypti dispersal. After nearly six years during which wMel was observed only at low frequencies in PE, corresponding to an apparent equilibrium between immigration and selection, wMel rose to fixation by 2018. Using analytic approximations and statistical analyses, we demonstrate that the observed fixation of wMel at PE is consistent with both stochastic transition past an unstable threshold frequency and deterministic transformation produced by steady immigration at a rate just above the threshold required for deterministic invasion. The indeterminacy results from a delicate balance of parameters needed to produce the delayed transition observed. Our analyses suggest that once Wolbachia transinfections are established locally through systematic introductions, stochastic “threshold crossing” is likely to only minimally enhance spatial spread, providing a local ratchet that slightly – but systematically – aids area-wide transformation of disease-vector populations in heterogeneous landscapes. AU - Turelli, Michael AU - Barton, Nicholas H ID - 11686 KW - Biological sciences TI - Wolbachia frequency data from: Why did the Wolbachia transinfection cross the road? Drift, deterministic dynamics and disease control ER - TY - JOUR AB - Predicting function from sequence is a central problem of biology. Currently, this is possible only locally in a narrow mutational neighborhood around a wildtype sequence rather than globally from any sequence. Using random mutant libraries, we developed a biophysical model that accounts for multiple features of σ70 binding bacterial promoters to predict constitutive gene expression levels from any sequence. We experimentally and theoretically estimated that 10–20% of random sequences lead to expression and ~80% of non-expressing sequences are one mutation away from a functional promoter. The potential for generating expression from random sequences is so pervasive that selection acts against σ70-RNA polymerase binding sites even within inter-genic, promoter-containing regions. This pervasiveness of σ70-binding sites implies that emergence of promoters is not the limiting step in gene regulatory evolution. Ultimately, the inclusion of novel features of promoter function into a mechanistic model enabled not only more accurate predictions of gene expression levels, but also identified that promoters evolve more rapidly than previously thought. AU - Lagator, Mato AU - Sarikas, Srdjan AU - Steinrueck, Magdalena AU - Toledo-Aparicio, David AU - Bollback, Jonathan P AU - Guet, Calin C AU - Tkačik, Gašper ID - 10736 JF - eLife TI - Predicting bacterial promoter function and evolution from random sequences VL - 11 ER - TY - JOUR AB - Hybridization is a common evolutionary process with multiple possible outcomes. In vertebrates, interspecific hybridization has repeatedly generated parthenogenetic hybrid species. However, it is unknown whether the generation of parthenogenetic hybrids is a rare outcome of frequent hybridization between sexual species within a genus or the typical outcome of rare hybridization events. Darevskia is a genus of rock lizards with both hybrid parthenogenetic and sexual species. Using capture sequencing, we estimate phylogenetic relationships and gene flow among the sexual species, to determine how introgressive hybridization relates to the origins of parthenogenetic hybrids. We find evidence for widespread hybridization with gene flow, both between recently diverged species and deep branches. Surprisingly, we find no signal of gene flow between parental species of the parthenogenetic hybrids, suggesting that the parental pairs were either reproductively or geographically isolated early in their divergence. The generation of parthenogenetic hybrids in Darevskia is, then, a rare outcome of the total occurrence of hybridization within the genus, but the typical outcome when specific species pairs hybridize. Our results question the conventional view that parthenogenetic lineages are generated by hybridization in a window of divergence. Instead, they suggest that some lineages possess specific properties that underpin successful parthenogenetic reproduction. AU - Freitas, Susana AU - Westram, Anja M AU - Schwander, Tanja AU - Arakelyan, Marine AU - Ilgaz, Çetin AU - Kumlutas, Yusuf AU - Harris, David James AU - Carretero, Miguel A. AU - Butlin, Roger K. ID - 11334 IS - 5 JF - Evolution SN - 0014-3820 TI - Parthenogenesis in Darevskia lizards: A rare outcome of common hybridization, not a common outcome of rare hybridization VL - 76 ER - TY - JOUR AB - Empirical essays of fitness landscapes suggest that they may be rugged, that is having multiple fitness peaks. Such fitness landscapes, those that have multiple peaks, necessarily have special local structures, called reciprocal sign epistasis (Poelwijk et al. in J Theor Biol 272:141–144, 2011). Here, we investigate the quantitative relationship between the number of fitness peaks and the number of reciprocal sign epistatic interactions. Previously, it has been shown (Poelwijk et al. in J Theor Biol 272:141–144, 2011) that pairwise reciprocal sign epistasis is a necessary but not sufficient condition for the existence of multiple peaks. Applying discrete Morse theory, which to our knowledge has never been used in this context, we extend this result by giving the minimal number of reciprocal sign epistatic interactions required to create a given number of peaks. AU - Saona Urmeneta, Raimundo J AU - Kondrashov, Fyodor AU - Khudiakova, Kseniia ID - 11447 IS - 8 JF - Bulletin of Mathematical Biology KW - Computational Theory and Mathematics KW - General Agricultural and Biological Sciences KW - Pharmacology KW - General Environmental Science KW - General Biochemistry KW - Genetics and Molecular Biology KW - General Mathematics KW - Immunology KW - General Neuroscience SN - 0092-8240 TI - Relation between the number of peaks and the number of reciprocal sign epistatic interactions VL - 84 ER - TY - JOUR AB - Local adaptation leads to differences between populations within a species. In many systems, similar environmental contrasts occur repeatedly, sometimes driving parallel phenotypic evolution. Understanding the genomic basis of local adaptation and parallel evolution is a major goal of evolutionary genomics. It is now known that by preventing the break-up of favourable combinations of alleles across multiple loci, genetic architectures that reduce recombination, like chromosomal inversions, can make an important contribution to local adaptation. However, little is known about whether inversions also contribute disproportionately to parallel evolution. Our aim here is to highlight this knowledge gap, to showcase existing studies, and to illustrate the differences between genomic architectures with and without inversions using simple models. We predict that by generating stronger effective selection, inversions can sometimes speed up the parallel adaptive process or enable parallel adaptation where it would be impossible otherwise, but this is highly dependent on the spatial setting. We highlight that further empirical work is needed, in particular to cover a broader taxonomic range and to understand the relative importance of inversions compared to genomic regions without inversions. AU - Westram, Anja M AU - Faria, Rui AU - Johannesson, Kerstin AU - Butlin, Roger AU - Barton, Nicholas H ID - 11546 IS - 1856 JF - Philosophical Transactions of the Royal Society B: Biological Sciences KW - General Agricultural and Biological Sciences KW - General Biochemistry KW - Genetics and Molecular Biology SN - 0962-8436 TI - Inversions and parallel evolution VL - 377 ER - TY - JOUR AB - Spatially explicit population genetic models have long been developed, yet have rarely been used to test hypotheses about the spatial distribution of genetic diversity or the genetic divergence between populations. Here, we use spatially explicit coalescence simulations to explore the properties of the island and the two-dimensional stepping stone models under a wide range of scenarios with spatio-temporal variation in deme size. We avoid the simulation of genetic data, using the fact that under the studied models, summary statistics of genetic diversity and divergence can be approximated from coalescence times. We perform the simulations using gridCoal, a flexible spatial wrapper for the software msprime (Kelleher et al., 2016, Theoretical Population Biology, 95, 13) developed herein. In gridCoal, deme sizes can change arbitrarily across space and time, as well as migration rates between individual demes. We identify different factors that can cause a deviation from theoretical expectations, such as the simulation time in comparison to the effective deme size and the spatio-temporal autocorrelation across the grid. Our results highlight that FST, a measure of the strength of population structure, principally depends on recent demography, which makes it robust to temporal variation in deme size. In contrast, the amount of genetic diversity is dependent on the distant past when Ne is large, therefore longer run times are needed to estimate Ne than FST. Finally, we illustrate the use of gridCoal on a real-world example, the range expansion of silver fir (Abies alba Mill.) since the last glacial maximum, using different degrees of spatio-temporal variation in deme size. AU - Szep, Eniko AU - Trubenova, Barbora AU - Csilléry, Katalin ID - 11640 IS - 8 JF - Molecular Ecology Resources SN - 1755-098X TI - Using gridCoal to assess whether standard population genetic theory holds in the presence of spatio-temporal heterogeneity in population size VL - 22 ER - TY - JOUR AB - Sexual antagonism is a common hypothesis for driving the evolution of sex chromosomes, whereby recombination suppression is favored between sexually antagonistic loci and the sex-determining locus to maintain beneficial combinations of alleles. This results in the formation of a sex-determining region. Chromosomal inversions may contribute to recombination suppression but their precise role in sex chromosome evolution remains unclear. Because local adaptation is frequently facilitated through the suppression of recombination between adaptive loci by chromosomal inversions, there is potential for inversions that cover sex-determining regions to be involved in local adaptation as well, particularly if habitat variation creates environment-dependent sexual antagonism. With these processes in mind, we investigated sex determination in a well-studied example of local adaptation within a species: the intertidal snail, Littorina saxatilis. Using SNP data from a Swedish hybrid zone, we find novel evidence for a female-heterogametic sex determination system that is restricted to one ecotype. Our results suggest that four putative chromosomal inversions, two previously described and two newly discovered, span the putative sex chromosome pair. We determine their differing associations with sex, which suggest distinct strata of differing ages. The same inversions are found in the second ecotype but do not show any sex association. The striking disparity in inversion-sex associations between ecotypes that are connected by gene flow across a habitat transition that is just a few meters wide indicates a difference in selective regime that has produced a distinct barrier to the spread of the newly discovered sex-determining region between ecotypes. Such sex chromosome-environment interactions have not previously been uncovered in L. saxatilis and are known in few other organisms. A combination of both sex-specific selection and divergent natural selection is required to explain these highly unusual patterns. AU - Hearn, Katherine E. AU - Koch, Eva L. AU - Stankowski, Sean AU - Butlin, Roger K. AU - Faria, Rui AU - Johannesson, Kerstin AU - Westram, Anja M ID - 12001 IS - 5 JF - Evolution Letters TI - Differing associations between sex determination and sex-linked inversions in two ecotypes of Littorina saxatilis VL - 6 ER - TY - JOUR AB - Polygenic adaptation is thought to be ubiquitous, yet remains poorly understood. Here, we model this process analytically, in the plausible setting of a highly polygenic, quantitative trait that experiences a sudden shift in the fitness optimum. We show how the mean phenotype changes over time, depending on the effect sizes of loci that contribute to variance in the trait, and characterize the allele dynamics at these loci. Notably, we describe the two phases of the allele dynamics: The first is a rapid phase, in which directional selection introduces small frequency differences between alleles whose effects are aligned with or opposed to the shift, ultimately leading to small differences in their probability of fixation during a second, longer phase, governed by stabilizing selection. As we discuss, key results should hold in more general settings and have important implications for efforts to identify the genetic basis of adaptation in humans and other species. AU - Hayward, Laura AU - Sella, Guy ID - 12157 JF - eLife KW - General Immunology and Microbiology KW - General Biochemistry KW - Genetics and Molecular Biology KW - General Medicine KW - General Neuroscience TI - Polygenic adaptation after a sudden change in environment VL - 11 ER - TY - JOUR AB - Kerstin Johannesson is a marine ecologist and evolutionary biologist based at the Tjärnö Marine Laboratory of the University of Gothenburg, which is situated in the beautiful Kosterhavet National Park on the Swedish west coast. Her work, using marine periwinkles (especially Littorina saxatilis and L. fabalis) as main model systems, has made a remarkable contribution to marine evolutionary biology and our understanding of local adaptation and its genetic underpinnings. AU - Westram, Anja M AU - Butlin, Roger ID - 12166 IS - 1 JF - Molecular Ecology KW - Genetics KW - Ecology KW - Evolution KW - Behavior and Systematics SN - 0962-1083 TI - Professor Kerstin Johannesson–winner of the 2022 Molecular Ecology Prize VL - 32 ER - TY - JOUR AB - Hybrid speciation—the origin of new species resulting from the hybridization of genetically divergent lineages—was once considered rare, but genomic data suggest that it may occur more often than once thought. In this study, Noguerales and Ortego found genomic evidence supporting the hybrid origin of a grasshopper that is able to exploit a broader range of host plants than either of its putative parents. AU - Stankowski, Sean ID - 12234 IS - 11 JF - Evolution KW - General Agricultural and Biological Sciences KW - Genetics KW - Ecology KW - Evolution KW - Behavior and Systematics SN - 0014-3820 TI - Digest: On the origin of a possible hybrid species VL - 76 ER - TY - JOUR AB - Chromosomal inversions have been shown to play a major role in a local adaptation by suppressing recombination between alternative arrangements and maintaining beneficial allele combinations. However, so far, their importance relative to the remaining genome remains largely unknown. Understanding the genetic architecture of adaptation requires better estimates of how loci of different effect sizes contribute to phenotypic variation. Here, we used three Swedish islands where the marine snail Littorina saxatilis has repeatedly evolved into two distinct ecotypes along a habitat transition. We estimated the contribution of inversion polymorphisms to phenotypic divergence while controlling for polygenic effects in the remaining genome using a quantitative genetics framework. We confirmed the importance of inversions but showed that contributions of loci outside inversions are of similar magnitude, with variable proportions dependent on the trait and the population. Some inversions showed consistent effects across all sites, whereas others exhibited site-specific effects, indicating that the genomic basis for replicated phenotypic divergence is only partly shared. The contributions of sexual dimorphism as well as environmental factors to phenotypic variation were significant but minor compared to inversions and polygenic background. Overall, this integrated approach provides insight into the multiple mechanisms contributing to parallel phenotypic divergence. AU - Koch, Eva L. AU - Ravinet, Mark AU - Westram, Anja M AU - Johannesson, Kerstin AU - Butlin, Roger K. ID - 12247 IS - 10 JF - Evolution KW - General Agricultural and Biological Sciences KW - Genetics KW - Ecology KW - Evolution KW - Behavior and Systematics SN - 0014-3820 TI - Genetic architecture of repeated phenotypic divergence in Littorina saxatilis evolution VL - 76 ER - TY - GEN AB - Chromosomal inversions have been shown to play a major role in local adaptation by suppressing recombination between alternative arrangements and maintaining beneficial allele combinations. However, so far, their importance relative to the remaining genome remains largely unknown. Understanding the genetic architecture of adaptation requires better estimates of how loci of different effect sizes contribute to phenotypic variation. Here, we used three Swedish islands where the marine snail Littorina saxatilis has repeatedly evolved into two distinct ecotypes along a habitat transition. We estimated the contribution of inversion polymorphisms to phenotypic divergence while controlling for polygenic effects in the remaining genome using a quantitative genetics framework. We confirmed the importance of inversions but showed that contributions of loci outside inversions are of similar magnitude, with variable proportions dependent on the trait and the population. Some inversions showed consistent effects across all sites, whereas others exhibited site-specific effects, indicating that the genomic basis for replicated phenotypic divergence is only partly shared. The contributions of sexual dimorphism as well as environmental factors to phenotypic variation were significant but minor compared to inversions and polygenic background. Overall, this integrated approach provides insight into the multiple mechanisms contributing to parallel phenotypic divergence. AU - Koch, Eva AU - Ravinet, Mark AU - Westram, Anja M AU - Jonannesson, Kerstin AU - Butlin, Roger ID - 13066 TI - Data from: Genetic architecture of repeated phenotypic divergence in Littorina saxatilis ecotype evolution ER - TY - JOUR AB - Reproductive isolation (RI) is a core concept in evolutionary biology. It has been the central focus of speciation research since the modern synthesis and is the basis by which biological species are defined. Despite this, the term is used in seemingly different ways, and attempts to quantify RI have used very different approaches. After showing that the field lacks a clear definition of the term, we attempt to clarify key issues, including what RI is, how it can be quantified in principle, and how it can be measured in practice. Following other definitions with a genetic focus, we propose that RI is a quantitative measure of the effect that genetic differences between populations have on gene flow. Specifically, RI compares the flow of neutral alleles in the presence of these genetic differences to the flow without any such differences. RI is thus greater than zero when genetic differences between populations reduce the flow of neutral alleles between populations. We show how RI can be quantified in a range of scenarios. A key conclusion is that RI depends strongly on circumstances—including the spatial, temporal and genomic context—making it difficult to compare across systems. After reviewing methods for estimating RI from data, we conclude that it is difficult to measure in practice. We discuss our findings in light of the goals of speciation research and encourage the use of methods for estimating RI that integrate organismal and genetic approaches. AU - Westram, Anja M AU - Stankowski, Sean AU - Surendranadh, Parvathy AU - Barton, Nicholas H ID - 12264 IS - 9 JF - Journal of Evolutionary Biology KW - Ecology KW - Evolution KW - Behavior and Systematics SN - 1010-061X TI - What is reproductive isolation? VL - 35 ER - TY - JOUR AU - Westram, Anja M AU - Stankowski, Sean AU - Surendranadh, Parvathy AU - Barton, Nicholas H ID - 12265 IS - 9 JF - Journal of Evolutionary Biology KW - Ecology KW - Evolution KW - Behavior and Systematics SN - 1010-061X TI - Reproductive isolation, speciation, and the value of disagreement: A reply to the commentaries on ‘What is reproductive isolation?’ VL - 35 ER - TY - JOUR AB - 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)’. AU - Barton, Nicholas H AU - Olusanya, Oluwafunmilola O ID - 10787 IS - 1848 JF - Philosophical Transactions of the Royal Society B: Biological Sciences KW - General Agricultural and Biological Sciences KW - General Biochemistry KW - Genetics and Molecular Biology SN - 0962-8436 TI - The response of a metapopulation to a changing environment VL - 377 ER - TY - JOUR AB - 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)’. AU - Sachdeva, Himani AU - Olusanya, Oluwafunmilola O AU - Barton, Nicholas H ID - 10658 IS - 1846 JF - Philosophical Transactions of the Royal Society B SN - 0962-8436 TI - Genetic load and extinction in peripheral populations: The roles of migration, drift and demographic stochasticity VL - 377 ER - TY - JOUR AB - 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. AU - Surendranadh, Parvathy AU - Arathoon, Louise S AU - Baskett, Carina AU - Field, David AU - Pickup, Melinda AU - Barton, Nicholas H ID - 11411 IS - 3 JF - Genetics TI - Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus VL - 221 ER - TY - DATA AB - 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. AU - Surendranadh, Parvathy AU - Arathoon, Louise S AU - Baskett, Carina AU - Field, David AU - Pickup, Melinda AU - Barton, Nicholas H ID - 11321 TI - Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus ER - TY - JOUR AB - 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. AU - Hledik, Michal AU - Barton, Nicholas H AU - Tkačik, Gašper ID - 12081 IS - 36 JF - Proceedings of the National Academy of Sciences SN - 0027-8424 TI - Accumulation and maintenance of information in evolution VL - 119 ER - TY - THES AB - 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. AU - Belohlavy, Stefanie ID - 11388 SN - 978-3-99078-018-3 TI - The genetic basis of complex traits studied via analysis of evolve and resequence experiments ER - TY - JOUR AB - 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. AU - Bod'ová, Katarína AU - Szep, Eniko AU - Barton, Nicholas H ID - 10535 IS - 12 JF - PLoS Computational Biology SN - 1553-734X TI - Dynamic maximum entropy provides accurate approximation of structured population dynamics VL - 17 ER - TY - JOUR AB - 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. AU - Simon, Alexis AU - Fraisse, Christelle AU - El Ayari, Tahani AU - Liautard‐Haag, Cathy AU - Strelkov, Petr AU - Welch, John J AU - Bierne, Nicolas ID - 8708 IS - 1 JF - Journal of Evolutionary Biology SN - 1010061X TI - How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels VL - 34 ER - TY - JOUR AB - 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. AU - Salces-Castellano, Antonia AU - Stankowski, Sean AU - Arribas, Paula AU - Patino, Jairo AU - Karger, Dirk N. AU - Butlin, Roger AU - Emerson, Brent C. ID - 8743 IS - 2 JF - Evolution SN - 0014-3820 TI - Long-term cloud forest response to climate warming revealed by insect speciation history VL - 75 ER - TY - JOUR AB - 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. AU - Arnoux, Stéphanie AU - Fraisse, Christelle AU - Sauvage, Christopher ID - 8928 IS - 2 JF - Journal of Evolutionary Biology SN - 1010061X TI - Genomic inference of complex domestication histories in three Solanaceae species VL - 34 ER - TY - JOUR AB - 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. AU - Faria, Rui AU - Johannesson, Kerstin AU - Stankowski, Sean ID - 9100 IS - 1 JF - Journal of Evolutionary Biology SN - 1010061X TI - Speciation in marine environments: Diving under the surface VL - 34 ER - TY - JOUR AB - 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. AU - Fraisse, Christelle AU - Sachdeva, Himani ID - 9168 IS - 2 JF - Genetics SN - 1943-2631 TI - The rates of introgression and barriers to genetic exchange between hybridizing species: Sex chromosomes vs autosomes VL - 217 ER - TY - JOUR AB - 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. AU - Fraisse, Christelle AU - Popovic, Iva AU - Mazoyer, Clément AU - Spataro, Bruno AU - Delmotte, Stéphane AU - Romiguier, Jonathan AU - Loire, Étienne AU - Simon, Alexis AU - Galtier, Nicolas AU - Duret, Laurent AU - Bierne, Nicolas AU - Vekemans, Xavier AU - Roux, Camille ID - 9119 JF - Molecular Ecology Resources SN - 1755098X TI - DILS: Demographic inferences with linked selection by using ABC VL - 21 ER - TY - JOUR AB - 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. AU - Meier, Joana I. AU - Salazar, Patricio A. AU - Kučka, Marek AU - Davies, Robert William AU - Dréau, Andreea AU - Aldás, Ismael AU - Power, Olivia Box AU - Nadeau, Nicola J. AU - Bridle, Jon R. AU - Rolian, Campbell AU - Barton, Nicholas H AU - McMillan, W. Owen AU - Jiggins, Chris D. AU - Chan, Yingguang Frank ID - 9375 IS - 25 JF - PNAS TI - Haplotype tagging reveals parallel formation of hybrid races in two butterfly species VL - 118 ER - TY - JOUR AB - 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. AU - Koch, Eva L. AU - Morales, Hernán E. AU - Larsson, Jenny AU - Westram, Anja M AU - Faria, Rui AU - Lemmon, Alan R. AU - Lemmon, E. Moriarty AU - Johannesson, Kerstin AU - Butlin, Roger K. ID - 9394 IS - 3 JF - Evolution Letters TI - Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis VL - 5 ER - TY - JOUR AB - 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. AU - Stankowski, Sean AU - Ravinet, Mark ID - 9392 IS - 9 JF - Current Biology SN - 09609822 TI - Quantifying the use of species concepts VL - 31 ER - TY - GEN AB - 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. AU - Koch, Eva AU - Morales, Hernán E. AU - Larsson, Jenny AU - Westram, Anja M AU - Faria, Rui AU - Lemmon, Alan R. AU - Lemmon, E. Moriarty AU - Johannesson, Kerstin AU - Butlin, Roger K. ID - 12987 TI - Data from: Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis ER - TY - JOUR AB - 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. AU - Lagator, Mato AU - Uecker, Hildegard AU - Neve, Paul ID - 9410 IS - 5 JF - Biology letters TI - Adaptation at different points along antibiotic concentration gradients VL - 17 ER - TY - JOUR AB - 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. AU - Berdan, Emma L. AU - Blanckaert, Alexandre AU - Slotte, Tanja AU - Suh, Alexander AU - Westram, Anja M AU - Fragata, Inês ID - 9470 IS - 12 JF - Molecular Ecology SN - 09621083 TI - Unboxing mutations: Connecting mutation types with evolutionary consequences VL - 30 ER - TY - JOUR AB - 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. AU - Hledik, Michal AU - Polechova, Jitka AU - Beiglböck, Mathias AU - Herdina, Anna Nele AU - Strassl, Robert AU - Posch, Martin ID - 9816 IS - 7 JF - PLoS ONE TI - Analysis of the specificity of a COVID-19 antigen test in the Slovak mass testing program VL - 16 ER - TY - JOUR AB - 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. AU - Szep, Eniko AU - Sachdeva, Himani AU - Barton, Nicholas H ID - 9252 IS - 5 JF - Evolution KW - Genetics KW - Ecology KW - Evolution KW - Behavior and Systematics KW - General Agricultural and Biological Sciences SN - 0014-3820 TI - Polygenic local adaptation in metapopulations: A stochastic eco‐evolutionary model VL - 75 ER - TY - JOUR AB - 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. AU - Butlin, Roger K. AU - Servedio, Maria R. AU - Smadja, Carole M. AU - Bank, Claudia AU - Barton, Nicholas H AU - Flaxman, Samuel M. AU - Giraud, Tatiana AU - Hopkins, Robin AU - Larson, Erica L. AU - Maan, Martine E. AU - Meier, Joana AU - Merrill, Richard AU - Noor, Mohamed A. F. AU - Ortiz‐Barrientos, Daniel AU - Qvarnström, Anna ID - 9374 IS - 5 JF - Evolution KW - Genetics KW - Ecology KW - Evolution KW - Behavior and Systematics KW - General Agricultural and Biological Sciences SN - 0014-3820 TI - Homage to Felsenstein 1981, or why are there so few/many species? VL - 75 ER - TY - GEN AB - 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. AU - Szep, Eniko AU - Sachdeva, Himani AU - Barton, Nicholas H ID - 13062 TI - Supplementary code for: Polygenic local adaptation in metapopulations: A stochastic eco-evolutionary model ER - TY - JOUR AB - 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. AU - Stankowski, Sean AU - Ravinet, Mark ID - 9383 IS - 6 JF - Evolution SN - 0014-3820 TI - Defining the speciation continuum VL - 75 ER - TY - CHAP AB - 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. AU - Stankowski, Sean AU - Shipilina, Daria AU - Westram, Anja M ID - 14984 SN - 9780470016176 T2 - Encyclopedia of Life Sciences TI - Hybrid Zones VL - 2 ER - TY - DATA AB - 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. AU - Surendranadh, Parvathy AU - Arathoon, Louise S AU - Baskett, Carina AU - Field, David AU - Pickup, Melinda AU - Barton, Nicholas H ID - 9192 TI - Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus ER - TY - JOUR AB - 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. AU - Larsson, J. AU - Westram, Anja M AU - Bengmark, S. AU - Lundh, T. AU - Butlin, R. K. ID - 7651 IS - 163 JF - Journal of The Royal Society Interface SN - 1742-5689 TI - A developmentally descriptive method for quantifying shape in gastropod shells VL - 17 ER - TY - CHAP AB - 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. AU - Westram, Anja M AU - Faria, Rui AU - Butlin, Roger AU - Johannesson, Kerstin ID - 9123 SN - 9780470016176 T2 - eLS TI - Inversions and Evolution ER - TY - GEN AB - 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. AU - Simon, Alexis AU - Fraisse, Christelle AU - El Ayari, Tahani AU - Liautard-Haag, Cathy AU - Strelkov, Petr AU - Welch, John AU - Bierne, Nicolas ID - 13073 TI - How do species barriers decay? concordance and local introgression in mosaic hybrid zones of mussels ER -