TY - JOUR AB - Key innovations are fundamental to biological diversification, but their genetic basis is poorly understood. A recent transition from egg-laying to live-bearing in marine snails (Littorina spp.) provides the opportunity to study the genetic architecture of an innovation that has evolved repeatedly across animals. Individuals do not cluster by reproductive mode in a genome-wide phylogeny, but local genealogical analysis revealed numerous small genomic regions where all live-bearers carry the same core haplotype. Candidate regions show evidence for live-bearer–specific positive selection and are enriched for genes that are differentially expressed between egg-laying and live-bearing reproductive systems. Ages of selective sweeps suggest that live-bearer–specific alleles accumulated over more than 200,000 generations. Our results suggest that new functions evolve through the recruitment of many alleles rather than in a single evolutionary step. AU - Stankowski, Sean AU - Zagrodzka, Zuzanna B. AU - Garlovsky, Martin D. AU - Pal, Arka AU - Shipilina, Daria AU - Garcia Castillo, Diego Fernando AU - Lifchitz, Hila AU - Le Moan, Alan AU - Leder, Erica AU - Reeve, James AU - Johannesson, Kerstin AU - Westram, Anja M AU - Butlin, Roger K. ID - 14796 IS - 6678 JF - Science TI - The genetic basis of a recent transition to live-bearing in marine snails VL - 383 ER - TY - JOUR AB - Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems. AU - Johannesson, Kerstin AU - Faria, Rui AU - Le Moan, Alan AU - Rafajlović, Marina AU - Westram, Anja M AU - Butlin, Roger K. AU - Stankowski, Sean ID - 15099 JF - Trends in Genetics SN - 0168-9525 TI - Diverse pathways to speciation revealed by marine snails ER - TY - JOUR AB - Understanding population divergence that eventually leads to speciation is essential for evolutionary biology. High species diversity in the sea was regarded as a paradox when strict allopatry was considered necessary for most speciation events because geographical barriers seemed largely absent in the sea, and many marine species have high dispersal capacities. Combining genome-wide data with demographic modelling to infer the demographic history of divergence has introduced new ways to address this classical issue. These models assume an ancestral population that splits into two subpopulations diverging according to different scenarios that allow tests for periods of gene flow. Models can also test for heterogeneities in population sizes and migration rates along the genome to account, respectively, for background selection and selection against introgressed ancestry. To investigate how barriers to gene flow arise in the sea, we compiled studies modelling the demographic history of divergence in marine organisms and extracted preferred demographic scenarios together with estimates of demographic parameters. These studies show that geographical barriers to gene flow do exist in the sea but that divergence can also occur without strict isolation. Heterogeneity of gene flow was detected in most population pairs suggesting the predominance of semipermeable barriers during divergence. We found a weak positive relationship between the fraction of the genome experiencing reduced gene flow and levels of genome-wide differentiation. Furthermore, we found that the upper bound of the ‘grey zone of speciation’ for our dataset extended beyond that found before, implying that gene flow between diverging taxa is possible at higher levels of divergence than previously thought. Finally, we list recommendations for further strengthening the use of demographic modelling in speciation research. These include a more balanced representation of taxa, more consistent and comprehensive modelling, clear reporting of results and simulation studies to rule out nonbiological explanations for general results. AU - De Jode, Aurélien AU - Le Moan, Alan AU - Johannesson, Kerstin AU - Faria, Rui AU - Stankowski, Sean AU - Westram, Anja M AU - Butlin, Roger K. AU - Rafajlović, Marina AU - Fraisse, Christelle ID - 11479 IS - 2 JF - Evolutionary Applications TI - Ten years of demographic modelling of divergence and speciation in the sea VL - 16 ER - TY - JOUR AB - Inversions are structural mutations that reverse the sequence of a chromosome segment and reduce the effective rate of recombination in the heterozygous state. They play a major role in adaptation, as well as in other evolutionary processes such as speciation. Although inversions have been studied since the 1920s, they remain difficult to investigate because the reduced recombination conferred by them strengthens the effects of drift and hitchhiking, which in turn can obscure signatures of selection. Nonetheless, numerous inversions have been found to be under selection. Given recent advances in population genetic theory and empirical study, here we review how different mechanisms of selection affect the evolution of inversions. A key difference between inversions and other mutations, such as single nucleotide variants, is that the fitness of an inversion may be affected by a larger number of frequently interacting processes. This considerably complicates the analysis of the causes underlying the evolution of inversions. We discuss the extent to which these mechanisms can be disentangled, and by which approach. AU - Berdan, Emma L. AU - Barton, Nicholas H AU - Butlin, Roger AU - Charlesworth, Brian AU - Faria, Rui AU - Fragata, Inês AU - Gilbert, Kimberly J. AU - Jay, Paul AU - Kapun, Martin AU - Lotterhos, Katie E. AU - Mérot, Claire AU - Durmaz Mitchell, Esra AU - Pascual, Marta AU - Peichel, Catherine L. AU - Rafajlović, Marina AU - Westram, Anja M AU - Schaeffer, Stephen W. AU - Johannesson, Kerstin AU - Flatt, Thomas ID - 14556 JF - Journal of Evolutionary Biology SN - 1010-061X TI - How chromosomal inversions reorient the evolutionary process ER - TY - JOUR AB - Chromosomal rearrangements (CRs) have been known since almost the beginning of genetics. While an important role for CRs in speciation has been suggested, evidence primarily stems from theoretical and empirical studies focusing on the microevolutionary level (i.e., on taxon pairs where speciation is often incomplete). Although the role of CRs in eukaryotic speciation at a macroevolutionary level has been supported by associations between species diversity and rates of evolution of CRs across phylogenies, these findings are limited to a restricted range of CRs and taxa. Now that more broadly applicable and precise CR detection approaches have become available, we address the challenges in filling some of the conceptual and empirical gaps between micro- and macroevolutionary studies on the role of CRs in speciation. We synthesize what is known about the macroevolutionary impact of CRs and suggest new research avenues to overcome the pitfalls of previous studies to gain a more comprehensive understanding of the evolutionary significance of CRs in speciation across the tree of life. AU - Lucek, Kay AU - Giménez, Mabel D. AU - Joron, Mathieu AU - Rafajlović, Marina AU - Searle, Jeremy B. AU - Walden, Nora AU - Westram, Anja M AU - Faria, Rui ID - 14742 IS - 11 JF - Cold Spring Harbor Perspectives in Biology KW - General Biochemistry KW - Genetics and Molecular Biology SN - 1943-0264 TI - The impact of chromosomal rearrangements in speciation: From micro- to macroevolution VL - 15 ER - 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 - JOUR AB - Conflict over reproduction between females and males exists because of anisogamy and promiscuity. Together they generate differences in fitness optima between the sexes and result in antagonistic coevolution of female and male reproductive traits. Mounting duration is likely to be a compromise between male and female interests whose outcome depends on the intensity of sexual selection. The timing of sperm transfer during mounting is critical. For example, mountings may be interrupted before sperm is transferred as a consequence of female or male choice, or they may be prolonged to function as mate guarding. In the highly promiscuous intertidal snail Littorina saxatilis, mountings vary substantially in duration, from less than a minute to more than an hour, and it has been assumed that mountings of a few minutes do not result in any sperm being transferred. Here, we examined the timing of sperm transfer, a reproductive trait that is likely affected by sexual conflict. We performed time-controlled mounting trials using L. saxatilis males and virgin females, aiming to examine indirectly when the transfer of sperm starts. We observed the relationship between mounting duration and the proportion of developing embryos out of all eggs and embryos in the brood pouch. Developing embryos were observed in similar proportions in all treatments (i.e. 1, 5 and 10 or more minutes at which mountings were artificially interrupted), suggesting that sperm transfer begins rapidly (within 1 min) in L. saxatilis and very short matings do not result in sperm shortage in the females. We discuss how the observed pattern can be influenced by predation risk, population density, and female status and receptivity. AU - Perini, Samuel AU - Butlin, Rogerk AU - Westram, Anja M AU - Johannesson, Kerstin ID - 10926 IS - 1 JF - Journal of Molluscan Studies SN - 0260-1230 TI - Very short mountings are enough for sperm transfer in Littorina saxatilis VL - 88 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - Combining hybrid zone analysis with genomic data is a promising approach to understanding the genomic basis of adaptive divergence. It allows for the identification of genomic regions underlying barriers to gene flow. It also provides insights into spatial patterns of allele frequency change, informing about the interplay between environmental factors, dispersal and selection. However, when only a single hybrid zone is analysed, it is difficult to separate patterns generated by selection from those resulting from chance. Therefore, it is beneficial to look for repeatable patterns across replicate hybrid zones in the same system. We applied this approach to the marine snail Littorina saxatilis, which contains two ecotypes, adapted to wave-exposed rocks vs. high-predation boulder fields. The existence of numerous hybrid zones between ecotypes offered the opportunity to test for the repeatability of genomic architectures and spatial patterns of divergence. We sampled and phenotyped snails from seven replicate hybrid zones on the Swedish west coast and genotyped them for thousands of single nucleotide polymorphisms. Shell shape and size showed parallel clines across all zones. Many genomic regions showing steep clines and/or high differentiation were shared among hybrid zones, consistent with a common evolutionary history and extensive gene flow between zones, and supporting the importance of these regions for divergence. In particular, we found that several large putative inversions contribute to divergence in all locations. Additionally, we found evidence for consistent displacement of clines from the boulder–rock transition. Our results demonstrate patterns of spatial variation that would not be accessible without continuous spatial sampling, a large genomic data set and replicate hybrid zones. AU - Westram, Anja M AU - Faria, Rui AU - Johannesson, Kerstin AU - Butlin, Roger ID - 10838 IS - 15 JF - Molecular Ecology KW - Genetics KW - Ecology KW - Evolution KW - Behavior and Systematics SN - 0962-1083 TI - Using replicate hybrid zones to understand the genomic basis of adaptive divergence VL - 30 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 - 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 - JOUR AB - 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. AU - Perini, Samuel AU - Rafajlović, Marina AU - Westram, Anja M AU - Johannesson, Kerstin AU - Butlin, Roger K. ID - 7995 IS - 7 JF - Evolution SN - 00143820 TI - Assortative mating, sexual selection, and their consequences for gene flow in Littorina VL - 74 ER - TY - GEN AB - 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. AU - Perini, Samuel AU - Rafajlovic, Marina AU - Westram, Anja M AU - Johannesson, Kerstin AU - Butlin, Roger ID - 8809 TI - Data from: Assortative mating, sexual selection and their consequences for gene flow in Littorina ER - TY - JOUR AB - 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. AU - Kulmuni, Jonna AU - Butlin, Roger K. AU - Lucek, Kay AU - Savolainen, Vincent AU - Westram, Anja M ID - 8168 IS - 1806 JF - Philosophical Transactions of the Royal Society. Series B: Biological sciences SN - 0962-8436 TI - Towards the completion of speciation: The evolution of reproductive isolation beyond the first barriers VL - 375 ER - TY - JOUR AB - 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. AU - Stankowski, Sean AU - Westram, Anja M AU - Zagrodzka, Zuzanna B. AU - Eyres, Isobel AU - Broquet, Thomas AU - Johannesson, Kerstin AU - Butlin, Roger K. ID - 8167 IS - 1806 JF - Philosophical Transactions of the Royal Society. Series B: Biological Sciences TI - The evolution of strong reproductive isolation between sympatric intertidal snails VL - 375 ER - TY - JOUR AB - 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. AU - Johannesson, Kerstin AU - Zagrodzka, Zuzanna AU - Faria, Rui AU - Westram, Anja M AU - Butlin, Roger K. ID - 7205 IS - 3 JF - Journal of Evolutionary Biology SN - 1010061X TI - Is embryo abortion a post-zygotic barrier to gene flow between Littorina ecotypes? VL - 33 ER - TY - JOUR AB - 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. AU - Faria, Rui AU - Johannesson, Kerstin AU - Butlin, Roger K. AU - Westram, Anja M ID - 5911 IS - 3 JF - Trends in Ecology and Evolution SN - 01695347 TI - Evolving inversions VL - 34 ER - TY - JOUR AB - 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. AU - Faria, Rui AU - Chaube, Pragya AU - Morales, Hernán E. AU - Larsson, Tomas AU - Lemmon, Alan R. AU - Lemmon, Emily M. AU - Rafajlović, Marina AU - Panova, Marina AU - Ravinet, Mark AU - Johannesson, Kerstin AU - Westram, Anja M AU - Butlin, Roger K. ID - 6095 IS - 6 JF - Molecular Ecology SN - 0962-1083 TI - Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes VL - 28 ER - TY - GEN AB - 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 postzygotic barriers due to selection against hybrids remains unclear. High embryo abortion rates in this species could indicate the presence of such barriers. Postzygotic 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 1011 embryos (mean 130±123) and abortion rates varied between 0 and100% (mean 12%). We described female genotypes by using a hybrid index based on hundreds of SNPs differentiated between ecotypes with which we characterised 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 postzygotic barriers contributing to ecotype divergence and the high and variable abortion rates were marginally, if at all, explained by hybrid index of females. AU - Johannesson, Kerstin AU - Zagrodzka, Zuzanna AU - Faria, Rui AU - Westram, Anja M AU - Butlin, Roger ID - 13067 TI - Data from: Is embryo abortion a postzygotic barrier to gene flow between Littorina ecotypes? ER - TY - JOUR AB - The study of parallel ecological divergence provides important clues to the operation of natural selection. Parallel divergence often occurs in heterogeneous environments with different kinds of environmental gradients in different locations, but the genomic basis underlying this process is unknown. We investigated the genomics of rapid parallel adaptation in the marine snail Littorina saxatilis in response to two independent environmental axes (crab-predation versus wave-action and low-shore versus high-shore). Using pooled whole-genome resequencing, we show that sharing of genomic regions of high differentiation between environments is generally low but increases at smaller spatial scales. We identify different shared genomic regions of divergence for each environmental axis and show that most of these regions overlap with candidate chromosomal inversions. Several inversion regions are divergent and polymorphic across many localities. We argue that chromosomal inversions could store shared variation that fuels rapid parallel adaptation to heterogeneous environments, possibly as balanced polymorphism shared by adaptive gene flow. AU - Morales, Hernán E. AU - Faria, Rui AU - Johannesson, Kerstin AU - Larsson, Tomas AU - Panova, Marina AU - Westram, Anja M AU - Butlin, Roger K. ID - 7393 IS - 12 JF - Science Advances SN - 2375-2548 TI - Genomic architecture of parallel ecological divergence: Beyond a single environmental contrast VL - 5 ER - TY - GEN AB - 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. AU - Faria, Rui AU - Chaube, Pragya AU - Morales, Hernán E. AU - Larsson, Tomas AU - Lemmon, Alan R. AU - Lemmon, Emily M. AU - Rafajlović, Marina AU - Panova, Marina AU - Ravinet, Mark AU - Johannesson, Kerstin AU - Westram, Anja M AU - Butlin, Roger K. ID - 9837 TI - Data from: Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes ER - TY - JOUR AB - Adaptive divergence and speciation may happen despite opposition by gene flow. Identifying the genomic basis underlying divergence with gene flow is a major task in evolutionary genomics. Most approaches (e.g., outlier scans) focus on genomic regions of high differentiation. However, not all genomic architectures potentially underlying divergence are expected to show extreme differentiation. Here, we develop an approach that combines hybrid zone analysis (i.e., focuses on spatial patterns of allele frequency change) with system-specific simulations to identify loci inconsistent with neutral evolution. We apply this to a genome-wide SNP set from an ideally suited study organism, the intertidal snail Littorina saxatilis, which shows primary divergence between ecotypes associated with different shore habitats. We detect many SNPs with clinal patterns, most of which are consistent with neutrality. Among non-neutral SNPs, most are located within three large putative inversions differentiating ecotypes. Many non-neutral SNPs show relatively low levels of differentiation. We discuss potential reasons for this pattern, including loose linkage to selected variants, polygenic adaptation and a component of balancing selection within populations (which may be expected for inversions). Our work is in line with theory predicting a role for inversions in divergence, and emphasizes that genomic regions contributing to divergence may not always be accessible with methods purely based on allele frequency differences. These conclusions call for approaches that take spatial patterns of allele frequency change into account in other systems. AU - Westram, Anja M AU - Rafajlović, Marina AU - Chaube, Pragya AU - Faria, Rui AU - Larsson, Tomas AU - Panova, Marina AU - Ravinet, Mark AU - Blomberg, Anders AU - Mehlig, Bernhard AU - Johannesson, Kerstin AU - Butlin, Roger ID - 9917 IS - 4 JF - Evolution Letters SN - 2056-3744 TI - Clines on the seashore: The genomic architecture underlying rapid divergence in the face of gene flow VL - 2 ER - TY - JOUR AB - The evolution of assortative mating is a key part of the speciation process. Stronger assortment, or greater divergence in mating traits, between species pairs with overlapping ranges is commonly observed, but possible causes of this pattern of reproductive character displacement are difficult to distinguish. We use a multidisciplinary approach to provide a rare example where it is possible to distinguish among hypotheses concerning the evolution of reproductive character displacement. We build on an earlier comparative analysis that illustrated a strong pattern of greater divergence in penis form between pairs of sister species with overlapping ranges than between allopatric sister-species pairs, in a large clade of marine gastropods (Littorinidae). We investigate both assortative mating and divergence in male genitalia in one of the sister-species pairs, discriminating among three contrasting processes each of which can generate a pattern of reproductive character displacement: reinforcement, reproductive interference and the Templeton effect. We demonstrate reproductive character displacement in assortative mating, but not in genital form between this pair of sister species and use demographic models to distinguish among the different processes. Our results support a model with no gene flow since secondary contact and thus favor reproductive interference as the cause of reproductive character displacement for mate choice, rather than reinforcement. High gene flow within species argues against the Templeton effect. Secondary contact appears to have had little impact on genital divergence. AU - Hollander, Johan AU - Montaño-Rendón, Mauricio AU - Bianco, Giuseppe AU - Yang, Xi AU - Westram, Anja M AU - Duvaux, Ludovic AU - Reid, David G. AU - Butlin, Roger K. ID - 9915 IS - 6 JF - Evolution Letters SN - 2056-3744 TI - Are assortative mating and genital divergence driven by reinforcement? VL - 2 ER - TY - GEN AB - Adaptive divergence and speciation may happen despite opposition by gene flow. Identifying the genomic basis underlying divergence with gene flow is a major task in evolutionary genomics. Most approaches (e.g. outlier scans) focus on genomic regions of high differentiation. However, not all genomic architectures potentially underlying divergence are expected to show extreme differentiation. Here, we develop an approach that combines hybrid zone analysis (i.e. focuses on spatial patterns of allele frequency change) with system-specific simulations to identify loci inconsistent with neutral evolution. We apply this to a genome-wide SNP set from an ideally-suited study organism, the intertidal snail Littorina saxatilis, which shows primary divergence between ecotypes associated with different shore habitats. We detect many SNPs with clinal patterns, most of which are consistent with neutrality. Among non-neutral SNPs, most are located within three large putative inversions differentiating ecotypes. Many non-neutral SNPs show relatively low levels of differentiation. We discuss potential reasons for this pattern, including loose linkage to selected variants, polygenic adaptation and a component of balancing selection within populations (which may be expected for inversions). Our work is in line with theory predicting a role for inversions in divergence, and emphasises that genomic regions contributing to divergence may not always be accessible with methods purely based on allele frequency differences. These conclusions call for approaches that take spatial patterns of allele frequency change into account in other systems. AU - Westram, Anja M AU - Rafajlović, Marina AU - Chaube, Pragya AU - Faria, Rui AU - Larsson, Tomas AU - Panova, Marina AU - Ravinet, Mark AU - Blomberg, Anders AU - Mehlig, Bernhard AU - Johannesson, Kerstin AU - Butlin, Roger ID - 9930 TI - Data from: Clines on the seashore: the genomic architecture underlying rapid divergence in the face of gene flow ER - TY - GEN AB - The evolution of assortative mating is a key part of the speciation process. Stronger assortment, or greater divergence in mating traits, between species pairs with overlapping ranges is commonly observed, but possible causes of this pattern of reproductive character displacement are difficult to distinguish. We use a multidisciplinary approach to provide a rare example where it is possible to distinguish among hypotheses concerning the evolution of reproductive character displacement. We build on an earlier comparative analysis that illustrated a strong pattern of greater divergence in penis form between pairs of sister species with overlapping ranges than between allopatric sister-species pairs, in a large clade of marine gastropods (Littorinidae). We investigate both assortative mating and divergence in male genitalia in one of the sister-species pairs, discriminating among three contrasting processes each of which can generate a pattern of reproductive character displacement: reinforcement, reproductive interference and the Templeton effect. We demonstrate reproductive character displacement in assortative mating, but not in genital form between this pair of sister species and use demographic models to distinguish among the different processes. Our results support a model with no gene flow since secondary contact and thus favour reproductive interference as the cause of reproductive character displacement for mate choice, rather than reinforcement. High gene flow within species argues against the Templeton effect. Secondary contact appears to have had little impact on genital divergence. AU - Hollander, Johan AU - Montaño-Rendón, Mauricio AU - Bianco, Giuseppe AU - Yang, Xi AU - Westram, Anja M AU - Duvaux, Ludovic AU - Reid, David G. AU - Butlin, Roger K. ID - 9929 TI - Data from: Are assortative mating and genital divergence driven by reinforcement? ER -