@article{20102, abstract = {Speciation is rarely observable directly. A way forward is to compare pairs of ecotypes that evolved in parallel in similar contexts but have reached different degrees of reproductive isolation. Such comparisons are possible in the marine snail Littorina saxatilis by contrasting barriers to gene flow between parallel ecotypes in Spain and Sweden. In both countries, divergent ecotypes have evolved to withstand either crab predation or wave action. Here, we explore transects spanning contact zones between the Crab and the Wave ecotypes using low-coverage whole-genome sequencing, morphological and behavioural traits. Despite parallel phenotypic divergence, distinct patterns of differentiation between the ecotypes emerged: a continuous cline in Sweden indicating a weak barrier to gene flow, but two highly genetically and phenotypically divergent, and partly spatially overlapping clusters in Spain suggesting a much stronger barrier to gene flow. The absence of Spanish early-generation hybrids supported strong isolation, but a low level of gene flow is evident from molecular data. In both countries, highly differentiated loci were located in both shared and country-specific chromosomal inversions but were also present in collinear regions. Despite being considered the same species and showing similar levels of phenotypic divergence, the Spanish ecotypes are much closer to full reproductive isolation than the Swedish ones. Barriers to gene flow of very different strengths between ecotypes within the same species might be explained by dissimilarities in the spatial arrangement of habitats, the selection gradients or the ages of the systems.}, author = {Raffini, Francesca and De Jode, Aurélien and Johannesson, Kerstin and Faria, Rui and Zagrodzka, Zuzanna B. and Westram, Anja M and Galindo, Juan and Rolán-Alvarez, Emilio and Butlin, Roger K.}, issn = {1365-294X}, journal = {Molecular Ecology}, number = {21}, publisher = {Wiley}, title = {{Phenotypic divergence and genomic architecture between parallel ecotypes at two different points on the speciation continuum in a marine snail}}, doi = {10.1111/mec.70025}, volume = {34}, year = {2025}, } @article{20325, abstract = {Inferring genealogical relationships of wild populations is useful because it gives direct estimates of mating patterns and variance in reproductive success. Inference can be improved by including information about parentage shared between siblings, or by modelling phenotypes or population data related to mating. However, we currently lack a framework to infer parent–offspring relationships, sibships and population parameters in a single analysis. To address this, we here extend a previous method, Fractional Analysis of Paternity and Sibships, to include population data for the case where one parent is known. We illustrate this with the example of pollen dispersal in a natural hybrid zone population of the snapdragon Antirrhinum majus. Pollen dispersal is leptokurtic, with half of mating events occurring within 30 m, but with a long tail of mating events up to 859 m. Using simulations, we find that both sibship and population information substantially improve pedigree reconstruction, and that we can expect to resolve median dispersal distances with high accuracy.}, author = {Ellis, Thomas and Field, David and Barton, Nicholas H}, issn = {1365-294X}, journal = {Molecular Ecology}, number = {15}, publisher = {Wiley}, title = {{Joint estimation of paternity, sibships and pollen dispersal in a snapdragon hybrid zone}}, doi = {10.1111/mec.70051}, volume = {34}, year = {2025}, } @article{20190, abstract = {A major goal of speciation research is identifying loci that underpin barriers to gene flow. Population genomics takes a ‘bottom-up’ approach, scanning the genome for molecular signatures of processes that drive or maintain divergence. However, interpreting the ‘genomic landscape’ of speciation is complicated, because genome scans conflate multiple processes, most of which are not informative about gene flow. However, studying replicated population contrasts, including multiple incidences of secondary contact, can strengthen inferences. In this paper, we use linked-read sequencing (haplotagging), FST scans and genealogical methods to characterise the genomic landscape associated with replicate hybrid zone formation. We studied two flower colour varieties of the common snapdragon, Antirrhinum majus subspecies majus, that form secondary hybrid zones in multiple independent valleys in the Pyrenees. Consistent with past work, we found very low differentiation at one well-studied zone (Planoles). However, at a second zone (Avellanet), we found stronger differentiation and greater heterogeneity, which we argue is due to differences in the amount of introgression following secondary contact. Topology weighting of genealogical trees identified loci where haplotype diversity was associated with the two snapdragon varieties. Two of the strongest associations were at previously identified flower colour loci: Flavia, that affects yellow pigmentation, and Rosea/Eluta, two linked loci that affect magenta pigmentation. Preliminary analysis of coalescence times provides additional evidence for selective sweeps at these loci and barriers to gene flow. Our study highlights the impact of demographic history on the differentiation landscape, emphasising the need to distinguish between historical divergence and recent introgression.}, author = {Pal, Arka and Shipilina, Daria and Le Moan, Alan and Mcnairn, Adrian J. and Grenier, Jennifer K. and Kucka, Marek and Coop, Graham and Chan, Yingguang Frank and Barton, Nicholas H and Field, David and Stankowski, Sean}, issn = {1365-294X}, journal = {Molecular Ecology}, number = {22}, publisher = {Wiley}, title = {{Genealogical analysis of replicate flower colour hybrid zones in Antirrhinum}}, doi = {10.1111/mec.70067}, volume = {34}, year = {2025}, } @article{14463, abstract = {Inversions are thought to play a key role in adaptation and speciation, suppressing recombination between diverging populations. Genes influencing adaptive traits cluster in inversions, and changes in inversion frequencies are associated with environmental differences. However, in many organisms, it is unclear if inversions are geographically and taxonomically widespread. The intertidal snail, Littorina saxatilis, is one such example. Strong associations between putative polymorphic inversions and phenotypic differences have been demonstrated between two ecotypes of L. saxatilis in Sweden and inferred elsewhere, but no direct evidence for inversion polymorphism currently exists across the species range. Using whole genome data from 107 snails, most inversion polymorphisms were found to be widespread across the species range. The frequencies of some inversion arrangements were significantly different among ecotypes, suggesting a parallel adaptive role. Many inversions were also polymorphic in the sister species, L. arcana, hinting at an ancient origin.}, author = {Reeve, James and Butlin, Roger K. and Koch, Eva L. and Stankowski, Sean and Faria, Rui}, issn = {1365-294X}, journal = {Molecular Ecology}, number = {24}, publisher = {Wiley}, title = {{Chromosomal inversion polymorphisms are widespread across the species ranges of rough periwinkles (Littorina saxatilis and L. arcana)}}, doi = {10.1111/mec.17160}, volume = {33}, year = {2024}, } @article{12166, abstract = {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.}, author = {Westram, Anja M and Butlin, Roger}, issn = {1365-294X}, journal = {Molecular Ecology}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics}, number = {1}, pages = {26--29}, publisher = {Wiley}, title = {{Professor Kerstin Johannesson–winner of the 2022 Molecular Ecology Prize}}, doi = {10.1111/mec.16779}, volume = {32}, year = {2023}, } @article{14787, abstract = {Understanding the phenotypic and genetic architecture of reproductive isolation is a long‐standing goal of speciation research. In several systems, large‐effect loci contributing to barrier phenotypes have been characterized, but such causal connections are rarely known for more complex genetic architectures. In this study, we combine “top‐down” and “bottom‐up” approaches with demographic modelling toward an integrated understanding of speciation across a monkeyflower hybrid zone. Previous work suggests that pollinator visitation acts as a primary barrier to gene flow between two divergent red‐ and yellow‐flowered ecotypes ofMimulus aurantiacus. Several candidate isolating traits and anonymous single nucleotide polymorphism loci under divergent selection have been identified, but their genomic positions remain unknown. Here, we report findings from demographic analyses that indicate this hybrid zone formed by secondary contact, but that subsequent gene flow was restricted by widespread barrier loci across the genome. Using a novel, geographic cline‐based genome scan, we demonstrate that candidate barrier loci are broadly distributed across the genome, rather than mapping to one or a few “islands of speciation.” Quantitative trait locus (QTL) mapping reveals that most floral traits are highly polygenic, with little evidence that QTL colocalize, indicating that most traits are genetically independent. Finally, we find little evidence that QTL and candidate barrier loci overlap, suggesting that some loci contribute to other forms of reproductive isolation. Our findings highlight the challenges of understanding the genetic architecture of reproductive isolation and reveal that barriers to gene flow other than pollinator isolation may play an important role in this system.}, author = {Stankowski, Sean and Chase, Madeline A. and McIntosh, Hanna and Streisfeld, Matthew A.}, issn = {1365-294X}, journal = {Molecular Ecology}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics}, number = {8}, pages = {2041--2054}, publisher = {Wiley}, title = {{Integrating top‐down and bottom‐up approaches to understand the genetic architecture of speciation across a monkeyflower hybrid zone}}, doi = {10.1111/mec.16849}, volume = {32}, year = {2023}, } @article{12159, abstract = {The term “haplotype block” is commonly used in the developing field of haplotype-based inference methods. We argue that the term should be defined based on the structure of the Ancestral Recombination Graph (ARG), which contains complete information on the ancestry of a sample. We use simulated examples to demonstrate key features of the relationship between haplotype blocks and ancestral structure, emphasizing the stochasticity of the processes that generate them. Even the simplest cases of neutrality or of a “hard” selective sweep produce a rich structure, often missed by commonly used statistics. We highlight a number of novel methods for inferring haplotype structure, based on the full ARG, or on a sequence of trees, and illustrate how they can be used to define haplotype blocks using an empirical data set. While the advent of new, computationally efficient methods makes it possible to apply these concepts broadly, they (and additional new methods) could benefit from adding features to explore haplotype blocks, as we define them. Understanding and applying the concept of the haplotype block will be essential to fully exploit long and linked-read sequencing technologies.}, author = {Shipilina, Daria and Pal, Arka and Stankowski, Sean and Chan, Yingguang Frank and Barton, Nicholas H}, issn = {1365-294X}, journal = {Molecular Ecology}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics}, number = {6}, pages = {1441--1457}, publisher = {Wiley}, title = {{On the origin and structure of haplotype blocks}}, doi = {10.1111/mec.16793}, volume = {32}, year = {2023}, } @article{10838, abstract = {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.}, author = {Westram, Anja M and Faria, Rui and Johannesson, Kerstin and Butlin, Roger}, issn = {1365-294X}, journal = {Molecular Ecology}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics}, number = {15}, pages = {3797--3814}, publisher = {Wiley}, title = {{Using replicate hybrid zones to understand the genomic basis of adaptive divergence}}, doi = {10.1111/mec.15861}, volume = {30}, year = {2021}, } @article{9470, abstract = {A key step in understanding the genetic basis of different evolutionary outcomes (e.g., adaptation) is to determine the roles played by different mutation types (e.g., SNPs, translocations and inversions). To do this we must simultaneously consider different mutation types in an evolutionary framework. Here, we propose a research framework that directly utilizes the most important characteristics of mutations, their population genetic effects, to determine their relative evolutionary significance in a given scenario. We review known population genetic effects of different mutation types and show how these may be connected to different evolutionary outcomes. We provide examples of how to implement this framework and pinpoint areas where more data, theory and synthesis are needed. Linking experimental and theoretical approaches to examine different mutation types simultaneously is a critical step towards understanding their evolutionary significance.}, author = {Berdan, Emma L. and Blanckaert, Alexandre and Slotte, Tanja and Suh, Alexander and Westram, Anja M and Fragata, Inês}, issn = {1365-294X}, journal = {Molecular Ecology}, number = {12}, pages = {2710--2723}, publisher = {Wiley}, title = {{Unboxing mutations: Connecting mutation types with evolutionary consequences}}, doi = {10.1111/mec.15936}, volume = {30}, year = {2021}, } @article{6095, abstract = {Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients.}, author = {Faria, Rui and Chaube, Pragya and Morales, Hernán E. and Larsson, Tomas and Lemmon, Alan R. and Lemmon, Emily M. and Rafajlović, Marina and Panova, Marina and Ravinet, Mark and Johannesson, Kerstin and Westram, Anja M and Butlin, Roger K.}, issn = {1365-294X}, journal = {Molecular Ecology}, number = {6}, pages = {1375--1393}, publisher = {Wiley}, title = {{Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes}}, doi = {10.1111/mec.14972}, volume = {28}, year = {2019}, } @article{7421, abstract = {X and Y chromosomes can diverge when rearrangements block recombination between them. Here we present the first genomic view of a reciprocal translocation that causes two physically unconnected pairs of chromosomes to be coinherited as sex chromosomes. In a population of the common frog (Rana temporaria), both pairs of X and Y chromosomes show extensive sequence differentiation, but not degeneration of the Y chromosomes. A new method based on gene trees shows both chromosomes are sex‐linked. Furthermore, the gene trees from the two Y chromosomes have identical topologies, showing they have been coinherited since the reciprocal translocation occurred. Reciprocal translocations can thus reshape sex linkage on a much greater scale compared with inversions, the type of rearrangement that is much better known in sex chromosome evolution, and they can greatly amplify the power of sexually antagonistic selection to drive genomic rearrangement. Two more populations show evidence of other rearrangements, suggesting that this species has unprecedented structural polymorphism in its sex chromosomes.}, author = {Toups, Melissa A and Rodrigues, Nicolas and Perrin, Nicolas and Kirkpatrick, Mark}, issn = {1365-294X}, journal = {Molecular Ecology}, number = {8}, pages = {1877--1889}, publisher = {Wiley}, title = {{A reciprocal translocation radically reshapes sex‐linked inheritance in the common frog}}, doi = {10.1111/mec.14990}, volume = {28}, year = {2019}, } @article{6466, abstract = {One of the most striking and consistent results in speciation genomics is the heterogeneous divergence observed across the genomes of closely related species. This pattern was initially attributed to different levels of gene exchange—with divergence preserved at loci generating a barrier to gene flow but homogenized at unlinked neutral loci. Although there is evidence to support this model, it is now recognized that interpreting patterns of divergence across genomes is not so straightforward. One problem is that heterogenous divergence between populations can also be generated by other processes (e.g. recurrent selective sweeps or background selection) without any involvement of differential gene flow. Thus, integrated studies that identify which loci are likely subject to divergent selection are required to shed light on the interplay between selection and gene flow during the early phases of speciation. In this issue of Molecular Ecology, Rifkin et al. (2019) confront this challenge using a pair of sister morning glory species. They wisely design their sampling to take the geographic context of individuals into account, including geographically isolated (allopatric) and co‐occurring (sympatric) populations. This enabled them to show that individuals are phenotypically less differentiated in sympatry. They also found that the loci that resist introgression are enriched for those most differentiated in allopatry and loci that exhibit signals of divergent selection. One great strength of the study is the combination of methods from population genetics and molecular evolution, including the development of a model to simultaneously infer admixture proportions and selfing rates.}, author = {Field, David and Fraisse, Christelle}, issn = {1365-294X}, journal = {Molecular ecology}, number = {7}, pages = {1579--1581}, publisher = {Wiley}, title = {{Breaking down barriers in morning glories}}, doi = {10.1111/mec.15048}, volume = {28}, year = {2019}, } @article{7739, abstract = {Currently, there is much debate on the genetic architecture of quantitative traits in wild populations. Is trait variation influenced by many genes of small effect or by a few genes of major effect? Where is additive genetic variation located in the genome? Do the same loci cause similar phenotypic variation in different populations? Great tits (Parus major) have been studied extensively in long‐term studies across Europe and consequently are considered an ecological ‘model organism’. Recently, genomic resources have been developed for the great tit, including a custom SNP chip and genetic linkage map. In this study, we used a suite of approaches to investigate the genetic architecture of eight quantitative traits in two long‐term study populations of great tits—one in the Netherlands and the other in the United Kingdom. Overall, we found little evidence for the presence of genes of large effects in either population. Instead, traits appeared to be influenced by many genes of small effect, with conservative estimates of the number of contributing loci ranging from 31 to 310. Despite concordance between population‐specific heritabilities, we found no evidence for the presence of loci having similar effects in both populations. While population‐specific genetic architectures are possible, an undetected shared architecture cannot be rejected because of limited power to map loci of small and moderate effects. This study is one of few examples of genetic architecture analysis in replicated wild populations and highlights some of the challenges and limitations researchers will face when attempting similar molecular quantitative genetic studies in free‐living populations.}, author = {Santure, Anna W. and Poissant, Jocelyn and De Cauwer, Isabelle and van Oers, Kees and Robinson, Matthew Richard and Quinn, John L. and Groenen, Martien A. M. and Visser, Marcel E. and Sheldon, Ben C. and Slate, Jon}, issn = {0962-1083}, journal = {Molecular Ecology}, pages = {6148--6162}, publisher = {Wiley}, title = {{Replicated analysis of the genetic architecture of quantitative traits in two wild great tit populations}}, doi = {10.1111/mec.13452}, volume = {24}, year = {2015}, } @article{7745, abstract = {The underlying basis of genetic variation in quantitative traits, in terms of the number of causal variants and the size of their effects, is largely unknown in natural populations. The expectation is that complex quantitative trait variation is attributable to many, possibly interacting, causal variants, whose effects may depend upon the sex, age and the environment in which they are expressed. A recently developed methodology in animal breeding derives a value of relatedness among individuals from high‐density genomic marker data, to estimate additive genetic variance within livestock populations. Here, we adapt and test the effectiveness of these methods to partition genetic variation for complex traits across genomic regions within ecological study populations where individuals have varying degrees of relatedness. We then apply this approach for the first time to a natural population and demonstrate that genetic variation in wing length in the great tit (Parus major) reflects contributions from multiple genomic regions. We show that a polygenic additive mode of gene action best describes the patterns observed, and we find no evidence of dosage compensation for the sex chromosome. Our results suggest that most of the genomic regions that influence wing length have the same effects in both sexes. We found a limited amount of genetic variance in males that is attributed to regions that have no effects in females, which could facilitate the sexual dimorphism observed for this trait. Although this exploratory work focuses on one complex trait, the methodology is generally applicable to any trait for any laboratory or wild population, paving the way for investigating sex‐, age‐ and environment‐specific genetic effects and thus the underlying genetic architecture of phenotype in biological study systems.}, author = {Robinson, Matthew Richard and Santure, Anna W. and DeCauwer, Isabelle and Sheldon, Ben C. and Slate, Jon}, issn = {0962-1083}, journal = {Molecular Ecology}, number = {15}, pages = {3963--3980}, publisher = {Wiley}, title = {{Partitioning of genetic variation across the genome using multimarker methods in a wild bird population}}, doi = {10.1111/mec.12375}, volume = {22}, year = {2013}, } @article{7746, abstract = {Clutch size and egg mass are life history traits that have been extensively studied in wild bird populations, as life history theory predicts a negative trade‐off between them, either at the phenotypic or at the genetic level. Here, we analyse the genomic architecture of these heritable traits in a wild great tit (Parus major) population, using three marker‐based approaches – chromosome partitioning, quantitative trait locus (QTL) mapping and a genome‐wide association study (GWAS). The variance explained by each great tit chromosome scales with predicted chromosome size, no location in the genome contains genome‐wide significant QTL, and no individual SNPs are associated with a large proportion of phenotypic variation, all of which may suggest that variation in both traits is due to many loci of small effect, located across the genome. There is no evidence that any regions of the genome contribute significantly to both traits, which combined with a small, nonsignificant negative genetic covariance between the traits, suggests the absence of genetic constraints on the independent evolution of these traits. Our findings support the hypothesis that variation in life history traits in natural populations is likely to be determined by many loci of small effect spread throughout the genome, which are subject to continued input of variation by mutation and migration, although we cannot exclude the possibility of an additional input of major effect genes influencing either trait.}, author = {Santure, Anna W. and De Cauwer, Isabelle and Robinson, Matthew Richard and Poissant, Jocelyn and Sheldon, Ben C. and Slate, Jon}, issn = {0962-1083}, journal = {Molecular Ecology}, number = {15}, pages = {3949--3962}, publisher = {Wiley}, title = {{Genomic dissection of variation in clutch size and egg mass in a wild great tit (Parus major) population}}, doi = {10.1111/mec.12376}, volume = {22}, year = {2013}, }