@article{11479, abstract = {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.}, author = {De Jode, Aurélien and Le Moan, Alan and Johannesson, Kerstin and Faria, Rui and Stankowski, Sean and Westram, Anja M and Butlin, Roger K. and Rafajlović, Marina and Fraisse, Christelle}, issn = {1752-4571}, journal = {Evolutionary Applications}, number = {2}, pages = {542--559}, publisher = {Wiley}, title = {{Ten years of demographic modelling of divergence and speciation in the sea}}, doi = {10.1111/eva.13428}, volume = {16}, year = {2023}, } @article{12514, abstract = {The concept of a “speciation continuum” has gained popularity in recent decades. It emphasizes speciation as a continuous process that may be studied by comparing contemporary population pairs that show differing levels of divergence. In their recent perspective article in Evolution, Stankowski and Ravinet provided a valuable service by formally defining the speciation continuum as a continuum of reproductive isolation, based on opinions gathered from a survey of speciation researchers. While we agree that the speciation continuum has been a useful concept to advance the understanding of the speciation process, some intrinsic limitations exist. Here, we advocate for a multivariate extension, the speciation hypercube, first proposed by Dieckmann et al. in 2004, but rarely used since. We extend the idea of the speciation cube and suggest it has strong conceptual and practical advantages over a one-dimensional model. We illustrate how the speciation hypercube can be used to visualize and compare different speciation trajectories, providing new insights into the processes and mechanisms of speciation. A key strength of the speciation hypercube is that it provides a unifying framework for speciation research, as it allows questions from apparently disparate subfields to be addressed in a single conceptual model.}, author = {Bolnick, Daniel I. and Hund, Amanda K. and Nosil, Patrik and Peng, Foen and Ravinet, Mark and Stankowski, Sean and Subramanian, Swapna and Wolf, Jochen B.W. and Yukilevich, Roman}, issn = {1558-5646}, journal = {Evolution: International journal of organic evolution}, number = {1}, pages = {318--328}, publisher = {Oxford University Press}, title = {{A multivariate view of the speciation continuum}}, doi = {10.1093/evolut/qpac004}, volume = {77}, 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{14452, abstract = {The classical infinitesimal model is a simple and robust model for the inheritance of quantitative traits. In this model, a quantitative trait is expressed as the sum of a genetic and an environmental component, and the genetic component of offspring traits within a family follows a normal distribution around the average of the parents’ trait values, and has a variance that is independent of the parental traits. In previous work, we showed that when trait values are determined by the sum of a large number of additive Mendelian factors, each of small effect, one can justify the infinitesimal model as a limit of Mendelian inheritance. In this paper, we show that this result extends to include dominance. We define the model in terms of classical quantities of quantitative genetics, before justifying it as a limit of Mendelian inheritance as the number, M, of underlying loci tends to infinity. As in the additive case, the multivariate normal distribution of trait values across the pedigree can be expressed in terms of variance components in an ancestral population and probabilities of identity by descent determined by the pedigree. Now, with just first-order dominance effects, we require two-, three-, and four-way identities. We also show that, even if we condition on parental trait values, the “shared” and “residual” components of trait values within each family will be asymptotically normally distributed as the number of loci tends to infinity, with an error of order 1/M−−√⁠. We illustrate our results with some numerical examples.}, author = {Barton, Nicholas H and Etheridge, Alison M. and Véber, Amandine}, issn = {1943-2631}, journal = {Genetics}, number = {2}, publisher = {Oxford Academic}, title = {{The infinitesimal model with dominance}}, doi = {10.1093/genetics/iyad133}, volume = {225}, year = {2023}, } @misc{12949, abstract = {The classical infinitesimal model is a simple and robust model for the inheritance of quantitative traits. In this model, a quantitative trait is expressed as the sum of a genetic and a non-genetic (environmental) component and the genetic component of offspring traits within a family follows a normal distribution around the average of the parents’ trait values, and has a variance that is independent of the trait values of the parents. Although the trait distribution across the whole population can be far from normal, the trait distributions within families are normally distributed with a variance-covariance matrix that is determined entirely by that in the ancestral population and the probabilities of identity determined by the pedigree. Moreover, conditioning on some of the trait values within the pedigree has predictable effects on the mean and variance within and between families. In previous work, Barton et al. (2017), we showed that when trait values are determined by the sum of a large number of Mendelian factors, each of small effect, one can justify the infinitesimal model as limit of Mendelian inheritance. It was also shown that under some forms of epistasis, trait values within a family are still normally distributed.}, author = {Barton, Nicholas H}, keywords = {Quantitative genetics, infinitesimal model}, publisher = {Institute of Science and Technology Austria}, title = {{The infinitesimal model with dominance}}, doi = {10.15479/AT:ISTA:12949}, year = {2023}, } @article{14556, abstract = {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.}, author = {Berdan, Emma L. and Barton, Nicholas H and Butlin, Roger and Charlesworth, Brian and Faria, Rui and Fragata, Inês and Gilbert, Kimberly J. and Jay, Paul and Kapun, Martin and Lotterhos, Katie E. and Mérot, Claire and Durmaz Mitchell, Esra and Pascual, Marta and Peichel, Catherine L. and Rafajlović, Marina and Westram, Anja M and Schaeffer, Stephen W. and Johannesson, Kerstin and Flatt, Thomas}, issn = {1420-9101}, journal = {Journal of Evolutionary Biology}, publisher = {Wiley}, title = {{How chromosomal inversions reorient the evolutionary process}}, doi = {10.1111/jeb.14242}, year = {2023}, } @article{14552, abstract = {Interactions between plants and herbivores are central in most ecosystems, but their strength is highly variable. The amount of variability within a system is thought to influence most aspects of plant-herbivore biology, from ecological stability to plant defense evolution. Our understanding of what influences variability, however, is limited by sparse data. We collected standardized surveys of herbivory for 503 plant species at 790 sites across 116° of latitude. With these data, we show that within-population variability in herbivory increases with latitude, decreases with plant size, and is phylogenetically structured. Differences in the magnitude of variability are thus central to how plant-herbivore biology varies across macroscale gradients. We argue that increased focus on interaction variability will advance understanding of patterns of life on Earth.}, author = {Robinson, M. L. and Hahn, P. G. and Inouye, B. D. and Underwood, N. and Whitehead, S. R. and Abbott, K. C. and Bruna, E. M. and Cacho, N. I. and Dyer, L. A. and Abdala-Roberts, L. and Allen, W. J. and Andrade, J. F. and Angulo, D. F. and Anjos, D. and Anstett, D. N. and Bagchi, R. and Bagchi, S. and Barbosa, M. and Barrett, S. and Baskett, Carina and Ben-Simchon, E. and Bloodworth, K. J. and Bronstein, J. L. and Buckley, Y. M. and Burghardt, K. T. and Bustos-Segura, C. and Calixto, E. S. and Carvalho, R. L. and Castagneyrol, B. and Chiuffo, M. C. and Cinoğlu, D. and Cinto Mejía, E. and Cock, M. C. and Cogni, R. and Cope, O. L. and Cornelissen, T. and Cortez, D. R. and Crowder, D. W. and Dallstream, C. and Dáttilo, W. and Davis, J. K. and Dimarco, R. D. and Dole, H. E. and Egbon, I. N. and Eisenring, M. and Ejomah, A. and Elderd, B. D. and Endara, M. J. and Eubanks, M. D. and Everingham, S. E. and Farah, K. N. and Farias, R. P. and Fernandes, A. P. and Fernandes, G. W. and Ferrante, M. and Finn, A. and Florjancic, G. A. and Forister, M. L. and Fox, Q. N. and Frago, E. and França, F. M. and Getman-Pickering, A. S. and Getman-Pickering, Z. and Gianoli, E. and Gooden, B. and Gossner, M. M. and Greig, K. A. and Gripenberg, S. and Groenteman, R. and Grof-Tisza, P. and Haack, N. and Hahn, L. and Haq, S. M. and Helms, A. M. and Hennecke, J. and Hermann, S. L. and Holeski, L. M. and Holm, S. and Hutchinson, M. C. and Jackson, E. E. and Kagiya, S. and Kalske, A. and Kalwajtys, M. and Karban, R. and Kariyat, R. and Keasar, T. and Kersch-Becker, M. F. and Kharouba, H. M. and Kim, T. N. and Kimuyu, D. M. and Kluse, J. and Koerner, S. E. and Komatsu, K. J. and Krishnan, S. and Laihonen, M. and Lamelas-López, L. and Lascaleia, M. C. and Lecomte, N. and Lehn, C. R. and Li, X. and Lindroth, R. L. and Lopresti, E. F. and Losada, M. and Louthan, A. M. and Luizzi, V. J. and Lynch, S. C. and Lynn, J. S. and Lyon, N. J. and Maia, L. F. and Maia, R. A. and Mannall, T. L. and Martin, B. S. and Massad, T. J. and Mccall, A. C. and Mcgurrin, K. and Merwin, A. C. and Mijango-Ramos, Z. and Mills, C. H. and Moles, A. T. and Moore, C. M. and Moreira, X. and Morrison, C. R. and Moshobane, M. C. and Muola, A. and Nakadai, R. and Nakajima, K. and Novais, S. and Ogbebor, C. O. and Ohsaki, H. and Pan, V. S. and Pardikes, N. A. and Pareja, M. and Parthasarathy, N. and Pawar, R. R. and Paynter, Q. and Pearse, I. S. and Penczykowski, R. M. and Pepi, A. A. and Pereira, C. C. and Phartyal, S. S. and Piper, F. I. and Poveda, K. and Pringle, E. G. and Puy, J. and Quijano, T. and Quintero, C. and Rasmann, S. and Rosche, C. and Rosenheim, L. Y. and Rosenheim, J. A. and Runyon, J. B. and Sadeh, A. and Sakata, Y. and Salcido, D. M. and Salgado-Luarte, C. and Santos, B. A. and Sapir, Y. and Sasal, Y. and Sato, Y. and Sawant, M. and Schroeder, H. and Schumann, I. and Segoli, M. and Segre, H. and Shelef, O. and Shinohara, N. and Singh, R. P. and Smith, D. S. and Sobral, M. and Stotz, G. C. and Tack, A. J.M. and Tayal, M. and Tooker, J. F. and Torrico-Bazoberry, D. and Tougeron, K. and Trowbridge, A. M. and Utsumi, S. and Uyi, O. and Vaca-Uribe, J. L. and Valtonen, A. and Van Dijk, L. J.A. and Vandvik, V. and Villellas, J. and Waller, L. P. and Weber, M. G. and Yamawo, A. and Yim, S. and Zarnetske, P. L. and Zehr, L. N. and Zhong, Z. and Wetzel, W. C.}, issn = {1095-9203}, journal = {Science}, number = {6671}, pages = {679--683}, publisher = {AAAS}, title = {{Plant size, latitude, and phylogeny explain within-population variability in herbivory}}, doi = {10.1126/science.adh8830}, volume = {382}, year = {2023}, } @misc{14579, abstract = {This is associated with our paper "Plant size, latitude, and phylogeny explain within-population variability in herbivory" published in Science. }, author = {Wetzel, William}, publisher = {Zenodo}, title = {{HerbVar-Network/HV-Large-Patterns-MS-public: v1.0.0}}, doi = {10.5281/ZENODO.8133117}, year = {2023}, } @phdthesis{14058, abstract = {Females and males across species are subject to divergent selective pressures arising from di↵erent reproductive interests and ecological niches. This often translates into a intricate array of sex-specific natural and sexual selection on traits that have a shared genetic basis between both sexes, causing a genetic sexual conflict. The resolution of this conflict mostly relies on the evolution of sex-specific expression of the shared genes, leading to phenotypic sexual dimorphism. Such sex-specific gene expression is thought to evolve via modifications of the genetic networks ultimately linked to sex-determining transcription factors. Although much empirical and theoretical evidence supports this standard picture of the molecular basis of sexual conflict resolution, there still are a few open questions regarding the complex array of selective forces driving phenotypic di↵erentiation between the sexes, as well as the molecular mechanisms underlying sexspecific adaptation. I address some of these open questions in my PhD thesis. First, how do patterns of phenotypic sexual dimorphism vary within populations, as a response to the temporal and spatial changes in sex-specific selective forces? To tackle this question, I analyze the patterns of sex-specific phenotypic variation along three life stages and across populations spanning the whole geographical range of Rumex hastatulus, a wind-pollinated angiosperm, in the first Chapter of the thesis. Second, how do gene expression patterns lead to phenotypic dimorphism, and what are the molecular mechanisms underlying the observed transcriptomic variation? I address this question by examining the sex- and tissue-specific expression variation in newly-generated datasets of sex-specific expression in heads and gonads of Drosophila melanogaster. I additionally used two complementary approaches for the study of the genetic basis of sex di↵erences in gene expression in the second and third Chapters of the thesis. Third, how does intersex correlation, thought to be one of the main aspects constraining the ability for the two sexes to decouple, interact with the evolution of sexual dimorphism? I develop models of sex-specific stabilizing selection, mutation and drift to formalize common intuition regarding the patterns of covariation between intersex correlation and sexual dimorphism in the fourth Chapter of the thesis. Alltogether, the work described in this PhD thesis provides useful insights into the links between genetic, transcriptomic and phenotypic layers of sex-specific variation, and contributes to our general understanding of the dynamics of sexual dimorphism evolution.}, author = {Puixeu Sala, Gemma}, isbn = {978-3-99078-035-0}, issn = {2663-337X}, pages = {230}, publisher = {Institute of Science and Technology Austria}, title = {{The molecular basis of sexual dimorphism: Experimental and theoretical characterization of phenotypic, transcriptomic and genetic patterns of sex-specific adaptation}}, doi = {10.15479/at:ista:14058}, year = {2023}, } @article{14077, abstract = {The regulatory architecture of gene expression is known to differ substantially between sexes in Drosophila, but most studies performed so far used whole-body data and only single crosses, which may have limited their scope to detect patterns that are robust across tissues and biological replicates. Here, we use allele-specific gene expression of parental and reciprocal hybrid crosses between 6 Drosophila melanogaster inbred lines to quantify cis- and trans-regulatory variation in heads and gonads of both sexes separately across 3 replicate crosses. Our results suggest that female and male heads, as well as ovaries, have a similar regulatory architecture. On the other hand, testes display more and substantially different cis-regulatory effects, suggesting that sex differences in the regulatory architecture that have been previously observed may largely derive from testis-specific effects. We also examine the difference in cis-regulatory variation of genes across different levels of sex bias in gonads and heads. Consistent with the idea that intersex correlations constrain expression and can lead to sexual antagonism, we find more cis variation in unbiased and moderately biased genes in heads. In ovaries, reduced cis variation is observed for male-biased genes, suggesting that cis variants acting on these genes in males do not lead to changes in ovary expression. Finally, we examine the dominance patterns of gene expression and find that sex- and tissue-specific patterns of inheritance as well as trans-regulatory variation are highly variable across biological crosses, although these were performed in highly controlled experimental conditions. This highlights the importance of using various genetic backgrounds to infer generalizable patterns.}, author = {Puixeu Sala, Gemma and Macon, Ariana and Vicoso, Beatriz}, issn = {2160-1836}, journal = {G3: Genes, Genomes, Genetics}, keywords = {Genetics (clinical), Genetics, Molecular Biology}, number = {8}, publisher = {Oxford University Press}, title = {{Sex-specific estimation of cis and trans regulation of gene expression in heads and gonads of Drosophila melanogaster}}, doi = {10.1093/g3journal/jkad121}, volume = {13}, year = {2023}, } @misc{12933, abstract = {Datasets of the publication "Sex-specific estimation of cis and trans regulation of gene expression in heads and gonads of Drosophila melanogaster".}, author = {Puixeu Sala, Gemma}, publisher = {Institute of Science and Technology Austria}, title = {{Data from: Sex-specific estimation of cis and trans regulation of gene expression in heads and gonads of Drosophila melanogaster}}, doi = {10.15479/AT:ISTA:12933}, year = {2023}, } @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}, 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}, year = {2023}, } @phdthesis{14651, abstract = {For self-incompatibility (SI) to be stable in a population, theory predicts that sufficient inbreeding depression (ID) is required: the fitness of offspring from self-mated individuals must be low enough to prevent the spread of self-compatibility (SC). Reviews of natural plant populations have supported this theory, with SI species generally showing high levels of ID. However, there is thought to be an under-sampling of self-incompatible taxa in the current literature. In this thesis, I study inbreeding depression in the SI plant species Antirrhinum majus using both greenhouse crosses and a large collected field dataset. Additionally, the gametophytic S-locus of A. majus is highly heterozygous and polymorphic, thus making assembly and discovery of S-alleles very difficult. Here, 206 new alleles of the male component SLFs are presented, along with a phylogeny showing the high conservation with alleles from another Antirrhinum species. Lastly, selected sites within the protein structure of SLFs are investigated, with one site in particular highlighted as potentially being involved in the SI recognition mechanism.}, author = {Arathoon, Louise S}, issn = {2663 - 337X}, pages = {96}, publisher = {Institute of Science and Technology Austria}, title = {{Investigating inbreeding depression and the self-incompatibility locus of Antirrhinum majus}}, doi = {10.15479/at:ista:14651}, year = {2023}, } @article{14742, abstract = {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.}, author = {Lucek, Kay and Giménez, Mabel D. and Joron, Mathieu and Rafajlović, Marina and Searle, Jeremy B. and Walden, Nora and Westram, Anja M and Faria, Rui}, issn = {1943-0264}, journal = {Cold Spring Harbor Perspectives in Biology}, keywords = {General Biochemistry, Genetics and Molecular Biology}, number = {11}, publisher = {Cold Spring Harbor Laboratory}, title = {{The impact of chromosomal rearrangements in speciation: From micro- to macroevolution}}, doi = {10.1101/cshperspect.a041447}, volume = {15}, 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{14833, abstract = {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.}, author = {Stankowski, Sean and Zagrodzka, Zuzanna B and Galindo, Juan and Montaño-Rendón, Mauricio and Faria, Rui and Mikhailova, Natalia and Blakeslee, April M H and Arnason, Einar and Broquet, Thomas and Morales, Hernán E and Grahame, John W and Westram, Anja M and Johannesson, Kerstin and Butlin, Roger K}, issn = {2752-938X}, journal = {Evolutionary Journal of the Linnean Society}, number = {1}, publisher = {Oxford University Press}, title = {{Whole-genome phylogeography of the intertidal snail Littorina saxatilis}}, doi = {10.1093/evolinnean/kzad002}, volume = {2}, year = {2023}, } @unpublished{14732, abstract = {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.}, author = {Olusanya, Oluwafunmilola O and Khudiakova, Kseniia and Sachdeva, Himani}, booktitle = {bioRxiv}, title = {{Genetic load, eco-evolutionary feedback and extinction in a metapopulation}}, doi = {10.1101/2023.12.02.569702}, year = {2023}, } @misc{14812, abstract = {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}, author = {Stankowski, Sean}, publisher = {Zenodo}, title = {{Data and code for: The genetic architecture of a recent transition to live-bearing in marine snails}}, doi = {10.5281/ZENODO.8318995}, year = {2023}, } @phdthesis{12800, abstract = {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. }, author = {Julseth, Mara}, issn = {2791-4585}, pages = {21}, publisher = {Institute of Science and Technology Austria}, title = {{The effect of local population structure on genetic variation at selected loci in the A. majus hybrid zone}}, doi = {10.15479/at:ista:12800}, year = {2023}, } @article{11702, abstract = {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.}, author = {Barton, Nicholas H}, issn = {1091-6490}, journal = {Proceedings of the National Academy of Sciences of the United States of America}, number = {30}, publisher = {Proceedings of the National Academy of Sciences}, title = {{The "New Synthesis"}}, doi = {10.1073/pnas.2122147119}, volume = {119}, year = {2022}, } @phdthesis{11128, abstract = {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.}, author = {Matejovicova, Lenka}, isbn = {978-3-99078-016-9}, issn = {2663-337X}, pages = {112}, publisher = {Institute of Science and Technology Austria}, title = {{Genetic basis of flower colour as a model for adaptive evolution}}, doi = {10.15479/at:ista:11128}, year = {2022}, } @article{10604, abstract = {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.}, author = {Turelli, Michael and Barton, Nicholas H}, issn = {2056-3744}, journal = {Evolution Letters}, keywords = {genetics, ecology, evolution, behavior and systematics}, number = {1}, pages = {92--105}, publisher = {Wiley}, title = {{Why did the Wolbachia transinfection cross the road? Drift, deterministic dynamics, and disease control}}, doi = {10.1002/evl3.270}, volume = {6}, year = {2022}, } @misc{11686, abstract = {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.}, author = {Turelli, Michael and Barton, Nicholas H}, keywords = {Biological sciences}, publisher = {Dryad}, title = {{Wolbachia frequency data from: Why did the Wolbachia transinfection cross the road? Drift, deterministic dynamics and disease control}}, doi = {10.25338/B81931}, year = {2022}, } @article{10736, abstract = {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.}, author = {Lagator, Mato and Sarikas, Srdjan and Steinrueck, Magdalena and Toledo-Aparicio, David and Bollback, Jonathan P and Guet, Calin C and Tkačik, Gašper}, issn = {2050-084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Predicting bacterial promoter function and evolution from random sequences}}, doi = {10.7554/eLife.64543}, volume = {11}, year = {2022}, } @article{11334, abstract = {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.}, author = {Freitas, Susana and Westram, Anja M and Schwander, Tanja and Arakelyan, Marine and Ilgaz, Çetin and Kumlutas, Yusuf and Harris, David James and Carretero, Miguel A. and Butlin, Roger K.}, issn = {1558-5646}, journal = {Evolution}, number = {5}, pages = {899--914}, publisher = {Wiley}, title = {{Parthenogenesis in Darevskia lizards: A rare outcome of common hybridization, not a common outcome of rare hybridization}}, doi = {10.1111/evo.14462}, volume = {76}, year = {2022}, } @article{11447, abstract = {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.}, author = {Saona Urmeneta, Raimundo J and Kondrashov, Fyodor and Khudiakova, Kseniia}, issn = {1522-9602}, journal = {Bulletin of Mathematical Biology}, keywords = {Computational Theory and Mathematics, General Agricultural and Biological Sciences, Pharmacology, General Environmental Science, General Biochemistry, Genetics and Molecular Biology, General Mathematics, Immunology, General Neuroscience}, number = {8}, publisher = {Springer Nature}, title = {{Relation between the number of peaks and the number of reciprocal sign epistatic interactions}}, doi = {10.1007/s11538-022-01029-z}, volume = {84}, year = {2022}, } @article{11546, abstract = {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.}, author = {Westram, Anja M and Faria, Rui and Johannesson, Kerstin and Butlin, Roger and Barton, Nicholas H}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society B: Biological Sciences}, keywords = {General Agricultural and Biological Sciences, General Biochemistry, Genetics and Molecular Biology}, number = {1856}, publisher = {Royal Society of London}, title = {{Inversions and parallel evolution}}, doi = {10.1098/rstb.2021.0203}, volume = {377}, year = {2022}, } @article{11640, abstract = {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.}, author = {Szep, Eniko and Trubenova, Barbora and Csilléry, Katalin}, issn = {1755-0998}, journal = {Molecular Ecology Resources}, number = {8}, pages = {2941--2955}, publisher = {Wiley}, title = {{Using gridCoal to assess whether standard population genetic theory holds in the presence of spatio-temporal heterogeneity in population size}}, doi = {10.1111/1755-0998.13676}, volume = {22}, year = {2022}, } @article{12001, abstract = {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.}, author = {Hearn, Katherine E. and Koch, Eva L. and Stankowski, Sean and Butlin, Roger K. and Faria, Rui and Johannesson, Kerstin and Westram, Anja M}, issn = {2056-3744}, journal = {Evolution Letters}, number = {5}, pages = {358--374}, publisher = {Oxford Academic}, title = {{Differing associations between sex determination and sex-linked inversions in two ecotypes of Littorina saxatilis}}, doi = {10.1002/evl3.295}, volume = {6}, year = {2022}, } @article{12157, abstract = {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.}, author = {Hayward, Laura and Sella, Guy}, issn = {2050-084X}, journal = {eLife}, keywords = {General Immunology and Microbiology, General Biochemistry, Genetics and Molecular Biology, General Medicine, General Neuroscience}, publisher = {eLife Sciences Publications}, title = {{Polygenic adaptation after a sudden change in environment}}, doi = {10.7554/elife.66697}, volume = {11}, year = {2022}, } @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 = {2022}, } @article{12234, abstract = {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.}, author = {Stankowski, Sean}, issn = {1558-5646}, journal = {Evolution}, keywords = {General Agricultural and Biological Sciences, Genetics, Ecology, Evolution, Behavior and Systematics}, number = {11}, pages = {2784--2785}, publisher = {Wiley}, title = {{Digest: On the origin of a possible hybrid species}}, doi = {10.1111/evo.14632}, volume = {76}, year = {2022}, } @article{12247, abstract = {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.}, author = {Koch, Eva L. and Ravinet, Mark and Westram, Anja M and Johannesson, Kerstin and Butlin, Roger K.}, issn = {1558-5646}, journal = {Evolution}, keywords = {General Agricultural and Biological Sciences, Genetics, Ecology, Evolution, Behavior and Systematics}, number = {10}, pages = {2332--2346}, publisher = {Wiley}, title = {{Genetic architecture of repeated phenotypic divergence in Littorina saxatilis evolution}}, doi = {10.1111/evo.14602}, volume = {76}, year = {2022}, } @misc{13066, abstract = {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.}, author = {Koch, Eva and Ravinet, Mark and Westram, Anja M and Jonannesson, Kerstin and Butlin, Roger}, publisher = {Dryad}, title = {{Data from: Genetic architecture of repeated phenotypic divergence in Littorina saxatilis ecotype evolution}}, doi = {10.5061/DRYAD.M905QFV4B}, year = {2022}, } @article{12264, abstract = {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.}, author = {Westram, Anja M and Stankowski, Sean and Surendranadh, Parvathy and Barton, Nicholas H}, issn = {1420-9101}, journal = {Journal of Evolutionary Biology}, keywords = {Ecology, Evolution, Behavior and Systematics}, number = {9}, pages = {1143--1164}, publisher = {Wiley}, title = {{What is reproductive isolation?}}, doi = {10.1111/jeb.14005}, volume = {35}, year = {2022}, } @article{12265, author = {Westram, Anja M and Stankowski, Sean and Surendranadh, Parvathy and Barton, Nicholas H}, issn = {1420-9101}, journal = {Journal of Evolutionary Biology}, keywords = {Ecology, Evolution, Behavior and Systematics}, number = {9}, pages = {1200--1205}, publisher = {Wiley}, title = {{Reproductive isolation, speciation, and the value of disagreement: A reply to the commentaries on ‘What is reproductive isolation?’}}, doi = {10.1111/jeb.14082}, volume = {35}, year = {2022}, } @article{10787, abstract = {A species distributed across diverse environments may adapt to local conditions. We ask how quickly such a species changes its range in response to changed conditions. Szép et al. (Szép E, Sachdeva H, Barton NH. 2021 Polygenic local adaptation in metapopulations: a stochastic eco-evolutionary model. Evolution75, 1030–1045 (doi:10.1111/evo.14210)) used the infinite island model to find the stationary distribution of allele frequencies and deme sizes. We extend this to find how a metapopulation responds to changes in carrying capacity, selection strength, or migration rate when deme sizes are fixed. We further develop a ‘fixed-state’ approximation. Under this approximation, polymorphism is only possible for a narrow range of habitat proportions when selection is weak compared to drift, but for a much wider range otherwise. When rates of selection or migration relative to drift change in a single deme of the metapopulation, the population takes a time of order m−1 to reach the new equilibrium. However, even with many loci, there can be substantial fluctuations in net adaptation, because at each locus, alleles randomly get lost or fixed. Thus, in a finite metapopulation, variation may gradually be lost by chance, even if it would persist in an infinite metapopulation. When conditions change across the whole metapopulation, there can be rapid change, which is predicted well by the fixed-state approximation. This work helps towards an understanding of how metapopulations extend their range across diverse environments. This article is part of the theme issue ‘Species’ ranges in the face of changing environments (Part II)’.}, author = {Barton, Nicholas H and Olusanya, Oluwafunmilola O}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society B: Biological Sciences}, keywords = {General Agricultural and Biological Sciences, General Biochemistry, Genetics and Molecular Biology}, number = {1848}, publisher = {The Royal Society}, title = {{The response of a metapopulation to a changing environment}}, doi = {10.1098/rstb.2021.0009}, volume = {377}, year = {2022}, } @article{10658, abstract = {We analyse how migration from a large mainland influences genetic load and population numbers on an island, in a scenario where fitness-affecting variants are unconditionally deleterious, and where numbers decline with increasing load. Our analysis shows that migration can have qualitatively different effects, depending on the total mutation target and fitness effects of deleterious variants. In particular, we find that populations exhibit a genetic Allee effect across a wide range of parameter combinations, when variants are partially recessive, cycling between low-load (large-population) and high-load (sink) states. Increased migration reduces load in the sink state (by increasing heterozygosity) but further inflates load in the large-population state (by hindering purging). We identify various critical parameter thresholds at which one or other stable state collapses, and discuss how these thresholds are influenced by the genetic versus demographic effects of migration. Our analysis is based on a ‘semi-deterministic’ analysis, which accounts for genetic drift but neglects demographic stochasticity. We also compare against simulations which account for both demographic stochasticity and drift. Our results clarify the importance of gene flow as a key determinant of extinction risk in peripheral populations, even in the absence of ecological gradients. This article is part of the theme issue ‘Species’ ranges in the face of changing environments (part I)’.}, author = {Sachdeva, Himani and Olusanya, Oluwafunmilola O and Barton, Nicholas H}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society B}, number = {1846}, publisher = {The Royal Society}, title = {{Genetic load and extinction in peripheral populations: The roles of migration, drift and demographic stochasticity}}, doi = {10.1098/rstb.2021.0010}, volume = {377}, year = {2022}, } @article{11411, abstract = {Many studies have quantified the distribution of heterozygosity and relatedness in natural populations, but few have examined the demographic processes driving these patterns. In this study, we take a novel approach by studying how population structure affects both pairwise identity and the distribution of heterozygosity in a natural population of the self-incompatible plant Antirrhinum majus. Excess variance in heterozygosity between individuals is due to identity disequilibrium, which reflects the variance in inbreeding between individuals; it is measured by the statistic g2. We calculated g2 together with FST and pairwise relatedness (Fij) using 91 SNPs in 22,353 individuals collected over 11 years. We find that pairwise Fij declines rapidly over short spatial scales, and the excess variance in heterozygosity between individuals reflects significant variation in inbreeding. Additionally, we detect an excess of individuals with around half the average heterozygosity, indicating either selfing or matings between close relatives. We use 2 types of simulation to ask whether variation in heterozygosity is consistent with fine-scale spatial population structure. First, by simulating offspring using parents drawn from a range of spatial scales, we show that the known pollen dispersal kernel explains g2. Second, we simulate a 1,000-generation pedigree using the known dispersal and spatial distribution and find that the resulting g2 is consistent with that observed from the field data. In contrast, a simulated population with uniform density underestimates g2, indicating that heterogeneous density promotes identity disequilibrium. Our study shows that heterogeneous density and leptokurtic dispersal can together explain the distribution of heterozygosity.}, author = {Surendranadh, Parvathy and Arathoon, Louise S and Baskett, Carina and Field, David and Pickup, Melinda and Barton, Nicholas H}, issn = {1943-2631}, journal = {Genetics}, number = {3}, publisher = {Oxford University Press}, title = {{Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus}}, doi = {10.1093/genetics/iyac083}, volume = {221}, year = {2022}, } @misc{11321, abstract = {Here are the research data underlying the publication "Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus" Further information are summed up in the README document. }, author = {Surendranadh, Parvathy and Arathoon, Louise S and Baskett, Carina and Field, David and Pickup, Melinda and Barton, Nicholas H}, publisher = {Institute of Science and Technology Austria}, title = {{Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus}}, doi = {10.15479/at:ista:11321}, year = {2022}, } @article{12081, abstract = {Selection accumulates information in the genome—it guides stochastically evolving populations toward states (genotype frequencies) that would be unlikely under neutrality. This can be quantified as the Kullback–Leibler (KL) divergence between the actual distribution of genotype frequencies and the corresponding neutral distribution. First, we show that this population-level information sets an upper bound on the information at the level of genotype and phenotype, limiting how precisely they can be specified by selection. Next, we study how the accumulation and maintenance of information is limited by the cost of selection, measured as the genetic load or the relative fitness variance, both of which we connect to the control-theoretic KL cost of control. The information accumulation rate is upper bounded by the population size times the cost of selection. This bound is very general, and applies across models (Wright–Fisher, Moran, diffusion) and to arbitrary forms of selection, mutation, and recombination. Finally, the cost of maintaining information depends on how it is encoded: Specifying a single allele out of two is expensive, but one bit encoded among many weakly specified loci (as in a polygenic trait) is cheap.}, author = {Hledik, Michal and Barton, Nicholas H and Tkačik, Gašper}, issn = {1091-6490}, journal = {Proceedings of the National Academy of Sciences}, number = {36}, publisher = {Proceedings of the National Academy of Sciences}, title = {{Accumulation and maintenance of information in evolution}}, doi = {10.1073/pnas.2123152119}, volume = {119}, year = {2022}, } @phdthesis{11388, abstract = {In evolve and resequence experiments, a population is sequenced, subjected to selection and then sequenced again, so that genetic changes before and after selection can be observed at the genetic level. Here, I use these studies to better understand the genetic basis of complex traits - traits which depend on more than a few genes. In the first chapter, I discuss the first evolve and resequence experiment, in which a population of mice, the so-called "Longshanks" mice, were selected for tibia length while their body mass was kept constant. The full pedigree is known. We observed a selection response on all chromosomes and used the infinitesimal model with linkage, a model which assumes an infinite number of genes with infinitesimally small effect sizes, as a null model. Results implied a very polygenic basis with a few loci of major effect standing out and changing in parallel. There was large variability between the different chromosomes in this study, probably due to LD. In chapter two, I go on to discuss the impact of LD, on the variability in an allele-frequency based summary statistic, giving an equation based on the initial allele frequencies, average pairwise LD, and the first four moments of the haplotype block copy number distribution. I describe this distribution by referring back to the founder generation. I then demonstrate how to infer selection via a maximum likelihood scheme on the example of a single locus and discuss how to extend this to more realistic scenarios. In chapter three, I discuss the second evolve and resequence experiment, in which a small population of Drosophila melanogaster was selected for increased pupal case size over 6 generations. The experiment was highly replicated with 27 lines selected within family and a known pedigree. We observed a phenotypic selection response of over one standard deviation. I describe the patterns in allele frequency data, including allele frequency changes and patterns of heterozygosity, and give ideas for future work.}, author = {Belohlavy, Stefanie}, isbn = {978-3-99078-018-3}, pages = {98}, publisher = {Institute of Science and Technology Austria}, title = {{The genetic basis of complex traits studied via analysis of evolve and resequence experiments}}, doi = {10.15479/at:ista:11388}, year = {2022}, } @article{10535, abstract = {Realistic models of biological processes typically involve interacting components on multiple scales, driven by changing environment and inherent stochasticity. Such models are often analytically and numerically intractable. We revisit a dynamic maximum entropy method that combines a static maximum entropy with a quasi-stationary approximation. This allows us to reduce stochastic non-equilibrium dynamics expressed by the Fokker-Planck equation to a simpler low-dimensional deterministic dynamics, without the need to track microscopic details. Although the method has been previously applied to a few (rather complicated) applications in population genetics, our main goal here is to explain and to better understand how the method works. We demonstrate the usefulness of the method for two widely studied stochastic problems, highlighting its accuracy in capturing important macroscopic quantities even in rapidly changing non-stationary conditions. For the Ornstein-Uhlenbeck process, the method recovers the exact dynamics whilst for a stochastic island model with migration from other habitats, the approximation retains high macroscopic accuracy under a wide range of scenarios in a dynamic environment.}, author = {Bod'ová, Katarína and Szep, Eniko and Barton, Nicholas H}, issn = {1553-7358}, journal = {PLoS Computational Biology}, number = {12}, publisher = {Public Library of Science}, title = {{Dynamic maximum entropy provides accurate approximation of structured population dynamics}}, doi = {10.1371/journal.pcbi.1009661}, volume = {17}, year = {2021}, } @article{8708, abstract = {The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study ‘replicated’ instances of secondary contact between closely related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry‐informative panel of such SNPs. We then compared their frequencies in newly sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi‐stable variants (Dobzhansky‐Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.}, author = {Simon, Alexis and Fraisse, Christelle and El Ayari, Tahani and Liautard‐Haag, Cathy and Strelkov, Petr and Welch, John J and Bierne, Nicolas}, issn = {14209101}, journal = {Journal of Evolutionary Biology}, number = {1}, pages = {208--223}, publisher = {Wiley}, title = {{How do species barriers decay? Concordance and local introgression in mosaic hybrid zones of mussels}}, doi = {10.1111/jeb.13709}, volume = {34}, year = {2021}, } @article{8743, abstract = {Montane cloud forests are areas of high endemism, and are one of the more vulnerable terrestrial ecosystems to climate change. Thus, understanding how they both contribute to the generation of biodiversity, and will respond to ongoing climate change, are important and related challenges. The widely accepted model for montane cloud forest dynamics involves upslope forcing of their range limits with global climate warming. However, limited climate data provides some support for an alternative model, where range limits are forced downslope with climate warming. Testing between these two models is challenging, due to the inherent limitations of climate and pollen records. We overcome this with an alternative source of historical information, testing between competing model predictions using genomic data and demographic analyses for a species of beetle tightly associated to an oceanic island cloud forest. Results unequivocally support the alternative model: populations that were isolated at higher elevation peaks during the Last Glacial Maximum are now in contact and hybridizing at lower elevations. Our results suggest that genomic data are a rich source of information to further understand how montane cloud forest biodiversity originates, and how it is likely to be impacted by ongoing climate change.}, author = {Salces-Castellano, Antonia and Stankowski, Sean and Arribas, Paula and Patino, Jairo and Karger, Dirk N. and Butlin, Roger and Emerson, Brent C.}, issn = {1558-5646}, journal = {Evolution}, number = {2}, pages = {231--244}, publisher = {Wiley}, title = {{Long-term cloud forest response to climate warming revealed by insect speciation history}}, doi = {10.1111/evo.14111}, volume = {75}, year = {2021}, } @article{8928, abstract = {Domestication is a human‐induced selection process that imprints the genomes of domesticated populations over a short evolutionary time scale and that occurs in a given demographic context. Reconstructing historical gene flow, effective population size changes and their timing is therefore of fundamental interest to understand how plant demography and human selection jointly shape genomic divergence during domestication. Yet, the comparison under a single statistical framework of independent domestication histories across different crop species has been little evaluated so far. Thus, it is unclear whether domestication leads to convergent demographic changes that similarly affect crop genomes. To address this question, we used existing and new transcriptome data on three crop species of Solanaceae (eggplant, pepper and tomato), together with their close wild relatives. We fitted twelve demographic models of increasing complexity on the unfolded joint allele frequency spectrum for each wild/crop pair, and we found evidence for both shared and species‐specific demographic processes between species. A convergent history of domestication with gene flow was inferred for all three species, along with evidence of strong reduction in the effective population size during the cultivation stage of tomato and pepper. The absence of any reduction in size of the crop in eggplant stands out from the classical view of the domestication process; as does the existence of a “protracted period” of management before cultivation. Our results also suggest divergent management strategies of modern cultivars among species as their current demography substantially differs. Finally, the timing of domestication is species‐specific and supported by the few historical records available.}, author = {Arnoux, Stéphanie and Fraisse, Christelle and Sauvage, Christopher}, issn = {14209101}, journal = {Journal of Evolutionary Biology}, number = {2}, pages = {270--283}, publisher = {Wiley}, title = {{Genomic inference of complex domestication histories in three Solanaceae species}}, doi = {10.1111/jeb.13723}, volume = {34}, year = {2021}, } @article{9100, abstract = {Marine environments are inhabited by a broad representation of the tree of life, yet our understanding of speciation in marine ecosystems is extremely limited compared with terrestrial and freshwater environments. Developing a more comprehensive picture of speciation in marine environments requires that we 'dive under the surface' by studying a wider range of taxa and ecosystems is necessary for a more comprehensive picture of speciation. Although studying marine evolutionary processes is often challenging, recent technological advances in different fields, from maritime engineering to genomics, are making it increasingly possible to study speciation of marine life forms across diverse ecosystems and taxa. Motivated by recent research in the field, including the 14 contributions in this issue, we highlight and discuss six axes of research that we think will deepen our understanding of speciation in the marine realm: (a) study a broader range of marine environments and organisms; (b) identify the reproductive barriers driving speciation between marine taxa; (c) understand the role of different genomic architectures underlying reproductive isolation; (d) infer the evolutionary history of divergence using model‐based approaches; (e) study patterns of hybridization and introgression between marine taxa; and (f) implement highly interdisciplinary, collaborative research programmes. In outlining these goals, we hope to inspire researchers to continue filling this critical knowledge gap surrounding the origins of marine biodiversity.}, author = {Faria, Rui and Johannesson, Kerstin and Stankowski, Sean}, issn = {14209101}, journal = {Journal of Evolutionary Biology}, number = {1}, pages = {4--15}, publisher = {Wiley}, title = {{Speciation in marine environments: Diving under the surface}}, doi = {10.1111/jeb.13756}, volume = {34}, year = {2021}, } @article{9168, abstract = {Interspecific crossing experiments have shown that sex chromosomes play a major role in reproductive isolation between many pairs of species. However, their ability to act as reproductive barriers, which hamper interspecific genetic exchange, has rarely been evaluated quantitatively compared to Autosomes. This genome-wide limitation of gene flow is essential for understanding the complete separation of species, and thus speciation. Here, we develop a mainland-island model of secondary contact between hybridizing species of an XY (or ZW) sexual system. We obtain theoretical predictions for the frequency of introgressed alleles, and the strength of the barrier to neutral gene flow for the two types of chromosomes carrying multiple interspecific barrier loci. Theoretical predictions are obtained for scenarios where introgressed alleles are rare. We show that the same analytical expressions apply for sex chromosomes and autosomes, but with different sex-averaged effective parameters. The specific features of sex chromosomes (hemizygosity and absence of recombination in the heterogametic sex) lead to reduced levels of introgression on the X (or Z) compared to autosomes. This effect can be enhanced by certain types of sex-biased forces, but it remains overall small (except when alleles causing incompatibilities are recessive). We discuss these predictions in the light of empirical data comprising model-based tests of introgression and cline surveys in various biological systems.}, author = {Fraisse, Christelle and Sachdeva, Himani}, issn = {1943-2631}, journal = {Genetics}, number = {2}, publisher = {Genetics Society of America}, title = {{The rates of introgression and barriers to genetic exchange between hybridizing species: Sex chromosomes vs autosomes}}, doi = {10.1093/genetics/iyaa025}, volume = {217}, year = {2021}, } @article{9119, abstract = {We present DILS, a deployable statistical analysis platform for conducting demographic inferences with linked selection from population genomic data using an Approximate Bayesian Computation framework. DILS takes as input single‐population or two‐population data sets (multilocus fasta sequences) and performs three types of analyses in a hierarchical manner, identifying: (a) the best demographic model to study the importance of gene flow and population size change on the genetic patterns of polymorphism and divergence, (b) the best genomic model to determine whether the effective size Ne and migration rate N, m are heterogeneously distributed along the genome (implying linked selection) and (c) loci in genomic regions most associated with barriers to gene flow. Also available via a Web interface, an objective of DILS is to facilitate collaborative research in speciation genomics. Here, we show the performance and limitations of DILS by using simulations and finally apply the method to published data on a divergence continuum composed by 28 pairs of Mytilus mussel populations/species.}, author = {Fraisse, Christelle and Popovic, Iva and Mazoyer, Clément and Spataro, Bruno and Delmotte, Stéphane and Romiguier, Jonathan and Loire, Étienne and Simon, Alexis and Galtier, Nicolas and Duret, Laurent and Bierne, Nicolas and Vekemans, Xavier and Roux, Camille}, issn = {17550998}, journal = {Molecular Ecology Resources}, pages = {2629--2644}, publisher = {Wiley}, title = {{DILS: Demographic inferences with linked selection by using ABC}}, doi = {10.1111/1755-0998.13323}, volume = {21}, year = {2021}, } @article{9375, abstract = {Genetic variation segregates as linked sets of variants, or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. And yet, genomic data often lack haplotype information, due to constraints in sequencing technologies. Here we present “haplotagging”, a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species the geographic clines for the major wing pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the centre of the hybrid zone. We propose that shared warning signalling (Müllerian mimicry) may couple the cline shifts seen in both species, and facilitate the parallel co-emergence of a novel hybrid morph in both co-mimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.}, author = {Meier, Joana I. and Salazar, Patricio A. and Kučka, Marek and Davies, Robert William and Dréau, Andreea and Aldás, Ismael and Power, Olivia Box and Nadeau, Nicola J. and Bridle, Jon R. and Rolian, Campbell and Barton, Nicholas H and McMillan, W. Owen and Jiggins, Chris D. and Chan, Yingguang Frank}, issn = {0027-8424}, journal = {PNAS}, number = {25}, publisher = {Proceedings of the National Academy of Sciences}, title = {{Haplotype tagging reveals parallel formation of hybrid races in two butterfly species}}, doi = {10.1073/pnas.2015005118}, volume = {118}, year = {2021}, } @article{9394, abstract = {Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.}, author = {Koch, Eva L. and Morales, Hernán E. and Larsson, Jenny and Westram, Anja M and Faria, Rui and Lemmon, Alan R. and Lemmon, E. Moriarty and Johannesson, Kerstin and Butlin, Roger K.}, issn = {2056-3744}, journal = {Evolution Letters}, number = {3}, pages = {196--213}, publisher = {Wiley}, title = {{Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis}}, doi = {10.1002/evl3.227}, volume = {5}, year = {2021}, } @article{9392, abstract = {Humans conceptualize the diversity of life by classifying individuals into types we call ‘species’1. The species we recognize influence political and financial decisions and guide our understanding of how units of diversity evolve and interact. Although the idea of species may seem intuitive, a debate about the best way to define them has raged even before Darwin2. So much energy has been devoted to the so-called ‘species problem’ that no amount of discourse will ever likely solve it2,3. Dozens of species concepts are currently recognized3, but we lack a concrete understanding of how much researchers actually disagree and the factors that cause them to think differently1,2. To address this, we used a survey to quantify the species problem for the first time. The results indicate that the disagreement is extensive: two randomly chosen respondents will most likely disagree on the nature of species. The probability of disagreement is not predicted by researcher experience or broad study system, but tended to be lower among researchers with similar focus, training and who study the same organism. Should we see this diversity of perspectives as a problem? We argue that we should not.}, author = {Stankowski, Sean and Ravinet, Mark}, issn = {18790445}, journal = {Current Biology}, number = {9}, pages = {R428--R429}, publisher = {Cell Press}, title = {{Quantifying the use of species concepts}}, doi = {10.1016/j.cub.2021.03.060}, volume = {31}, year = {2021}, } @misc{12987, abstract = {Chromosomal inversion polymorphisms, segments of chromosomes that are flipped in orientation and occur in reversed order in some individuals, have long been recognized to play an important role in local adaptation. They can reduce recombination in heterozygous individuals and thus help to maintain sets of locally adapted alleles. In a wide range of organisms, populations adapted to different habitats differ in frequency of inversion arrangements. However, getting a full understanding of the importance of inversions for adaptation requires confirmation of their influence on traits under divergent selection. Here, we studied a marine snail, Littorina saxatilis, that has evolved ecotypes adapted to wave exposure or crab predation. These two types occur in close proximity on different parts of the shore. Gene flow between them exists in contact zones. However, they exhibit strong phenotypic divergence in several traits under habitat-specific selection, including size, shape and behaviour. We used crosses between these ecotypes to identify genomic regions that explain variation in these traits by using QTL analysis and variance partitioning across linkage groups. We could show that previously detected inversion regions contribute to adaptive divergence. Some inversions influenced multiple traits suggesting that they contain sets of locally adaptive alleles. Our study also identified regions without known inversions that are important for phenotypic divergence. Thus, we provide a more complete overview of the importance of inversions in relation to the remaining genome.}, author = {Koch, Eva and Morales, Hernán E. and Larsson, Jenny and Westram, Anja M and Faria, Rui and Lemmon, Alan R. and Lemmon, E. Moriarty and Johannesson, Kerstin and Butlin, Roger K.}, publisher = {Dryad}, title = {{Data from: Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis}}, doi = {10.5061/DRYAD.ZGMSBCCB4}, year = {2021}, } @article{9410, abstract = {Antibiotic concentrations vary dramatically in the body and the environment. Hence, understanding the dynamics of resistance evolution along antibiotic concentration gradients is critical for predicting and slowing the emergence and spread of resistance. While it has been shown that increasing the concentration of an antibiotic slows resistance evolution, how adaptation to one antibiotic concentration correlates with fitness at other points along the gradient has not received much attention. Here, we selected populations of Escherichia coli at several points along a concentration gradient for three different antibiotics, asking how rapidly resistance evolved and whether populations became specialized to the antibiotic concentration they were selected on. Populations selected at higher concentrations evolved resistance more slowly but exhibited equal or higher fitness across the whole gradient. Populations selected at lower concentrations evolved resistance rapidly, but overall fitness in the presence of antibiotics was lower. However, these populations readily adapted to higher concentrations upon subsequent selection. Our results indicate that resistance management strategies must account not only for the rates of resistance evolution but also for the fitness of evolved strains.}, author = {Lagator, Mato and Uecker, Hildegard and Neve, Paul}, issn = {1744957X}, journal = {Biology letters}, number = {5}, publisher = {Royal Society of London}, title = {{Adaptation at different points along antibiotic concentration gradients}}, doi = {10.1098/rsbl.2020.0913}, volume = {17}, year = {2021}, } @article{9470, abstract = {A key step in understanding the genetic basis of different evolutionary outcomes (e.g., adaptation) is to determine the roles played by different mutation types (e.g., SNPs, translocations and inversions). To do this we must simultaneously consider different mutation types in an evolutionary framework. Here, we propose a research framework that directly utilizes the most important characteristics of mutations, their population genetic effects, to determine their relative evolutionary significance in a given scenario. We review known population genetic effects of different mutation types and show how these may be connected to different evolutionary outcomes. We provide examples of how to implement this framework and pinpoint areas where more data, theory and synthesis are needed. Linking experimental and theoretical approaches to examine different mutation types simultaneously is a critical step towards understanding their evolutionary significance.}, author = {Berdan, Emma L. and Blanckaert, Alexandre and Slotte, Tanja and Suh, Alexander and Westram, Anja M and Fragata, Inês}, issn = {1365294X}, journal = {Molecular Ecology}, number = {12}, pages = {2710--2723}, publisher = {Wiley}, title = {{Unboxing mutations: Connecting mutation types with evolutionary consequences}}, doi = {10.1111/mec.15936}, volume = {30}, year = {2021}, } @article{9816, abstract = {Aims: Mass antigen testing programs have been challenged because of an alleged insufficient specificity, leading to a large number of false positives. The objective of this study is to derive a lower bound of the specificity of the SD Biosensor Standard Q Ag-Test in large scale practical use. Methods: Based on county data from the nationwide tests for SARS-CoV-2 in Slovakia between 31.10.–1.11. 2020 we calculate a lower confidence bound for the specificity. As positive test results were not systematically verified by PCR tests, we base the lower bound on a worst case assumption, assuming all positives to be false positives. Results: 3,625,332 persons from 79 counties were tested. The lowest positivity rate was observed in the county of Rožňava where 100 out of 34307 (0.29%) tests were positive. This implies a test specificity of at least 99.6% (97.5% one-sided lower confidence bound, adjusted for multiplicity). Conclusion: The obtained lower bound suggests a higher specificity compared to earlier studies in spite of the underlying worst case assumption and the application in a mass testing setting. The actual specificity is expected to exceed 99.6% if the prevalence in the respective regions was non-negligible at the time of testing. To our knowledge, this estimate constitutes the first bound obtained from large scale practical use of an antigen test.}, author = {Hledik, Michal and Polechova, Jitka and Beiglböck, Mathias and Herdina, Anna Nele and Strassl, Robert and Posch, Martin}, issn = {1932-6203}, journal = {PLoS ONE}, number = {7}, publisher = {Public Library of Science}, title = {{Analysis of the specificity of a COVID-19 antigen test in the Slovak mass testing program}}, doi = {10.1371/journal.pone.0255267}, volume = {16}, year = {2021}, } @article{9252, abstract = {This paper analyses the conditions for local adaptation in a metapopulation with infinitely many islands under a model of hard selection, where population size depends on local fitness. Each island belongs to one of two distinct ecological niches or habitats. Fitness is influenced by an additive trait which is under habitat‐dependent directional selection. Our analysis is based on the diffusion approximation and accounts for both genetic drift and demographic stochasticity. By neglecting linkage disequilibria, it yields the joint distribution of allele frequencies and population size on each island. We find that under hard selection, the conditions for local adaptation in a rare habitat are more restrictive for more polygenic traits: even moderate migration load per locus at very many loci is sufficient for population sizes to decline. This further reduces the efficacy of selection at individual loci due to increased drift and because smaller populations are more prone to swamping due to migration, causing a positive feedback between increasing maladaptation and declining population sizes. Our analysis also highlights the importance of demographic stochasticity, which exacerbates the decline in numbers of maladapted populations, leading to population collapse in the rare habitat at significantly lower migration than predicted by deterministic arguments.}, author = {Szep, Eniko and Sachdeva, Himani and Barton, Nicholas H}, issn = {1558-5646}, journal = {Evolution}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics, General Agricultural and Biological Sciences}, number = {5}, pages = {1030--1045}, publisher = {Wiley}, title = {{Polygenic local adaptation in metapopulations: A stochastic eco‐evolutionary model}}, doi = {10.1111/evo.14210}, volume = {75}, year = {2021}, } @article{9374, abstract = {If there are no constraints on the process of speciation, then the number of species might be expected to match the number of available niches and this number might be indefinitely large. One possible constraint is the opportunity for allopatric divergence. In 1981, Felsenstein used a simple and elegant model to ask if there might also be genetic constraints. He showed that progress towards speciation could be described by the build‐up of linkage disequilibrium among divergently selected loci and between these loci and those contributing to other forms of reproductive isolation. Therefore, speciation is opposed by recombination, because it tends to break down linkage disequilibria. Felsenstein then introduced a crucial distinction between “two‐allele” models, which are subject to this effect, and “one‐allele” models, which are free from the recombination constraint. These fundamentally important insights have been the foundation for both empirical and theoretical studies of speciation ever since.}, author = {Butlin, Roger K. and Servedio, Maria R. and Smadja, Carole M. and Bank, Claudia and Barton, Nicholas H and Flaxman, Samuel M. and Giraud, Tatiana and Hopkins, Robin and Larson, Erica L. and Maan, Martine E. and Meier, Joana and Merrill, Richard and Noor, Mohamed A. F. and Ortiz‐Barrientos, Daniel and Qvarnström, Anna}, issn = {1558-5646}, journal = {Evolution}, keywords = {Genetics, Ecology, Evolution, Behavior and Systematics, General Agricultural and Biological Sciences}, number = {5}, pages = {978--988}, publisher = {Wiley}, title = {{Homage to Felsenstein 1981, or why are there so few/many species?}}, doi = {10.1111/evo.14235}, volume = {75}, year = {2021}, } @misc{13062, abstract = {This paper analyzes the conditions for local adaptation in a metapopulation with infinitely many islands under a model of hard selection, where population size depends on local fitness. Each island belongs to one of two distinct ecological niches or habitats. Fitness is influenced by an additive trait which is under habitat-dependent directional selection. Our analysis is based on the diffusion approximation and accounts for both genetic drift and demographic stochasticity. By neglecting linkage disequilibria, it yields the joint distribution of allele frequencies and population size on each island. We find that under hard selection, the conditions for local adaptation in a rare habitat are more restrictive for more polygenic traits: even moderate migration load per locus at very many loci is sufficient for population sizes to decline. This further reduces the efficacy of selection at individual loci due to increased drift and because smaller populations are more prone to swamping due to migration, causing a positive feedback between increasing maladaptation and declining population sizes. Our analysis also highlights the importance of demographic stochasticity, which exacerbates the decline in numbers of maladapted populations, leading to population collapse in the rare habitat at significantly lower migration than predicted by deterministic arguments.}, author = {Szep, Eniko and Sachdeva, Himani and Barton, Nicholas H}, publisher = {Dryad}, title = {{Supplementary code for: Polygenic local adaptation in metapopulations: A stochastic eco-evolutionary model}}, doi = {10.5061/DRYAD.8GTHT76P1}, year = {2021}, } @article{9383, abstract = {A primary roadblock to our understanding of speciation is that it usually occurs over a timeframe that is too long to study from start to finish. The idea of a speciation continuum provides something of a solution to this problem; rather than observing the entire process, we can simply reconstruct it from the multitude of speciation events that surround us. But what do we really mean when we talk about the speciation continuum, and can it really help us understand speciation? We explored these questions using a literature review and online survey of speciation researchers. Although most researchers were familiar with the concept and thought it was useful, our survey revealed extensive disagreement about what the speciation continuum actually tells us. This is due partly to the lack of a clear definition. Here, we provide an explicit definition that is compatible with the Biological Species Concept. That is, the speciation continuum is a continuum of reproductive isolation. After outlining the logic of the definition in light of alternatives, we explain why attempts to reconstruct the speciation process from present‐day populations will ultimately fail. We then outline how we think the speciation continuum concept can continue to act as a foundation for understanding the continuum of reproductive isolation that surrounds us.}, author = {Stankowski, Sean and Ravinet, Mark}, issn = {1558-5646}, journal = {Evolution}, number = {6}, pages = {1256--1273}, publisher = {Oxford University Press}, title = {{Defining the speciation continuum}}, doi = {10.1111/evo.14215}, volume = {75}, year = {2021}, } @inbook{14984, abstract = {Hybrid zones are narrow geographic regions where different populations, races or interbreeding species meet and mate, producing mixed ‘hybrid’ offspring. They are relatively common and can be found in a diverse range of organisms and environments. The study of hybrid zones has played an important role in our understanding of the origin of species, with hybrid zones having been described as ‘natural laboratories’. This is because they allow us to study,in situ, the conditions and evolutionary forces that enable divergent taxa to remain distinct despite some ongoing gene exchange between them.}, author = {Stankowski, Sean and Shipilina, Daria and Westram, Anja M}, booktitle = {Encyclopedia of Life Sciences}, isbn = {9780470016176}, publisher = {Wiley}, title = {{Hybrid Zones}}, doi = {10.1002/9780470015902.a0029355}, volume = {2}, year = {2021}, } @misc{9192, abstract = {Here are the research data underlying the publication " Effects of fine-scale population structure on inbreeding in a long-term study of snapdragons (Antirrhinum majus)." Further information are summed up in the README document.}, author = {Surendranadh, Parvathy and Arathoon, Louise S and Baskett, Carina and Field, David and Pickup, Melinda and Barton, Nicholas H}, publisher = {Institute of Science and Technology Austria}, title = {{Effects of fine-scale population structure on the distribution of heterozygosity in a long-term study of Antirrhinum majus}}, doi = {10.15479/AT:ISTA:9192}, year = {2021}, } @article{7651, abstract = {The growth of snail shells can be described by simple mathematical rules. Variation in a few parameters can explain much of the diversity of shell shapes seen in nature. However, empirical studies of gastropod shell shape variation typically use geometric morphometric approaches, which do not capture this growth pattern. We have developed a way to infer a set of developmentally descriptive shape parameters based on three-dimensional logarithmic helicospiral growth and using landmarks from two-dimensional shell images as input. We demonstrate the utility of this approach, and compare it to the geometric morphometric approach, using a large set of Littorina saxatilis shells in which locally adapted populations differ in shape. Our method can be modified easily to make it applicable to a wide range of shell forms, which would allow for investigations of the similarities and differences between and within many different species of gastropods.}, author = {Larsson, J. and Westram, Anja M and Bengmark, S. and Lundh, T. and Butlin, R. K.}, issn = {1742-5662}, journal = {Journal of The Royal Society Interface}, number = {163}, publisher = {The Royal Society}, title = {{A developmentally descriptive method for quantifying shape in gastropod shells}}, doi = {10.1098/rsif.2019.0721}, volume = {17}, year = {2020}, } @inbook{9123, abstract = {Inversions are chromosomal rearrangements where the order of genes is reversed. Inversions originate by mutation and can be under positive, negative or balancing selection. Selective effects result from potential disruptive effects on meiosis, gene disruption at inversion breakpoints and, importantly, the effects of inversions as modifiers of recombination rate: Recombination is strongly reduced in individuals heterozygous for an inversion, allowing for alleles at different loci to be inherited as a ‘block’. This may lead to a selective advantage whenever it is favourable to keep certain combinations of alleles associated, for example under local adaptation with gene flow. Inversions can cover a considerable part of a chromosome and contain numerous loci under different selection pressures, so that the resulting overall effects may be complex. Empirical data from various systems show that inversions may have a prominent role in local adaptation, speciation, parallel evolution, the maintenance of polymorphism and sex chromosome evolution.}, author = {Westram, Anja M and Faria, Rui and Butlin, Roger and Johannesson, Kerstin}, booktitle = {eLS}, isbn = {9780470016176}, publisher = {Wiley}, title = {{Inversions and Evolution}}, doi = {10.1002/9780470015902.a0029007}, year = {2020}, } @misc{13073, abstract = {The Mytilus complex of marine mussel species forms a mosaic of hybrid zones, found across temperate regions of the globe. This allows us to study "replicated" instances of secondary contact between closely-related species. Previous work on this complex has shown that local introgression is both widespread and highly heterogeneous, and has identified SNPs that are outliers of differentiation between lineages. Here, we developed an ancestry-informative panel of such SNPs. We then compared their frequencies in newly-sampled populations, including samples from within the hybrid zones, and parental populations at different distances from the contact. Results show that close to the hybrid zones, some outlier loci are near to fixation for the heterospecific allele, suggesting enhanced local introgression, or the local sweep of a shared ancestral allele. Conversely, genomic cline analyses, treating local parental populations as the reference, reveal a globally high concordance among loci, albeit with a few signals of asymmetric introgression. Enhanced local introgression at specific loci is consistent with the early transfer of adaptive variants after contact, possibly including asymmetric bi-stable variants (Dobzhansky-Muller incompatibilities), or haplotypes loaded with fewer deleterious mutations. Having escaped one barrier, however, these variants can be trapped or delayed at the next barrier, confining the introgression locally. These results shed light on the decay of species barriers during phases of contact.}, author = {Simon, Alexis and Fraisse, Christelle and El Ayari, Tahani and Liautard-Haag, Cathy and Strelkov, Petr and Welch, John and Bierne, Nicolas}, publisher = {Dryad}, title = {{How do species barriers decay? concordance and local introgression in mosaic hybrid zones of mussels}}, doi = {10.5061/DRYAD.R4XGXD29N}, year = {2020}, } @misc{13065, abstract = {Domestication is a human-induced selection process that imprints the genomes of domesticated populations over a short evolutionary time scale, and that occurs in a given demographic context. Reconstructing historical gene flow, effective population size changes and their timing is therefore of fundamental interest to understand how plant demography and human selection jointly shape genomic divergence during domestication. Yet, the comparison under a single statistical framework of independent domestication histories across different crop species has been little evaluated so far. Thus, it is unclear whether domestication leads to convergent demographic changes that similarly affect crop genomes. To address this question, we used existing and new transcriptome data on three crop species of Solanaceae (eggplant, pepper and tomato), together with their close wild relatives. We fitted twelve demographic models of increasing complexity on the unfolded joint allele frequency spectrum for each wild/crop pair, and we found evidence for both shared and species-specific demographic processes between species. A convergent history of domestication with gene-flow was inferred for all three species, along with evidence of strong reduction in the effective population size during the cultivation stage of tomato and pepper. The absence of any reduction in size of the crop in eggplant stands out from the classical view of the domestication process; as does the existence of a “protracted period” of management before cultivation. Our results also suggest divergent management strategies of modern cultivars among species as their current demography substantially differs. Finally, the timing of domestication is species-specific and supported by the few historical records available.}, author = {Arnoux, Stephanie and Fraisse, Christelle and Sauvage, Christopher}, publisher = {Dryad}, title = {{VCF files of synonymous SNPs related to: Genomic inference of complex domestication histories in three Solanaceae species}}, doi = {10.5061/DRYAD.Q2BVQ83HD}, year = {2020}, } @article{7995, abstract = {When divergent populations are connected by gene flow, the establishment of complete reproductive isolation usually requires the joint action of multiple barrier effects. One example where multiple barrier effects are coupled consists of a single trait that is under divergent natural selection and also mediates assortative mating. Such multiple‐effect traits can strongly reduce gene flow. However, there are few cases where patterns of assortative mating have been described quantitatively and their impact on gene flow has been determined. Two ecotypes of the coastal marine snail, Littorina saxatilis , occur in North Atlantic rocky‐shore habitats dominated by either crab predation or wave action. There is evidence for divergent natural selection acting on size, and size‐assortative mating has previously been documented. Here, we analyze the mating pattern in L. saxatilis with respect to size in intensively sampled transects across boundaries between the habitats. We show that the mating pattern is mostly conserved between ecotypes and that it generates both assortment and directional sexual selection for small male size. Using simulations, we show that the mating pattern can contribute to reproductive isolation between ecotypes but the barrier to gene flow is likely strengthened more by sexual selection than by assortment.}, author = {Perini, Samuel and Rafajlović, Marina and Westram, Anja M and Johannesson, Kerstin and Butlin, Roger K.}, issn = {15585646}, journal = {Evolution}, number = {7}, pages = {1482--1497}, publisher = {Wiley}, title = {{Assortative mating, sexual selection, and their consequences for gene flow in Littorina}}, doi = {10.1111/evo.14027}, volume = {74}, year = {2020}, } @misc{8809, abstract = {When divergent populations are connected by gene flow, the establishment of complete reproductive isolation usually requires the joint action of multiple barrier effects. One example where multiple barrier effects are coupled consists of a single trait that is under divergent natural selection and also mediates assortative mating. Such multiple-effect traits can strongly reduce gene flow. However, there are few cases where patterns of assortative mating have been described quantitatively and their impact on gene flow has been determined. Two ecotypes of the coastal marine snail, Littorina saxatilis, occur in North Atlantic rocky-shore habitats dominated by either crab predation or wave action. There is evidence for divergent natural selection acting on size, and size-assortative mating has previously been documented. Here, we analyze the mating pattern in L. saxatilis with respect to size in intensively-sampled transects across boundaries between the habitats. We show that the mating pattern is mostly conserved between ecotypes and that it generates both assortment and directional sexual selection for small male size. Using simulations, we show that the mating pattern can contribute to reproductive isolation between ecotypes but the barrier to gene flow is likely strengthened more by sexual selection than by assortment.}, author = {Perini, Samuel and Rafajlovic, Marina and Westram, Anja M and Johannesson, Kerstin and Butlin, Roger}, publisher = {Dryad}, title = {{Data from: Assortative mating, sexual selection and their consequences for gene flow in Littorina}}, doi = {10.5061/dryad.qrfj6q5cn}, year = {2020}, } @article{8112, author = {Barton, Nicholas H}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society. Series B: Biological Sciences}, number = {1806}, publisher = {The Royal Society}, title = {{On the completion of speciation}}, doi = {10.1098/rstb.2019.0530}, volume = {375}, year = {2020}, } @article{8168, abstract = {Speciation, that is, the evolution of reproductive barriers eventually leading to complete isolation, is a crucial process generating biodiversity. Recent work has contributed much to our understanding of how reproductive barriers begin to evolve, and how they are maintained in the face of gene flow. However, little is known about the transition from partial to strong reproductive isolation (RI) and the completion of speciation. We argue that the evolution of strong RI is likely to involve different processes, or new interactions among processes, compared with the evolution of the first reproductive barriers. Transition to strong RI may be brought about by changing external conditions, for example, following secondary contact. However, the increasing levels of RI themselves create opportunities for new barriers to evolve and, and interaction or coupling among barriers. These changing processes may depend on genomic architecture and leave detectable signals in the genome. We outline outstanding questions and suggest more theoretical and empirical work, considering both patterns and processes associated with strong RI, is needed to understand how speciation is completed.}, author = {Kulmuni, Jonna and Butlin, Roger K. and Lucek, Kay and Savolainen, Vincent and Westram, Anja M}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society. Series B: Biological sciences}, number = {1806}, publisher = {The Royal Society}, title = {{Towards the completion of speciation: The evolution of reproductive isolation beyond the first barriers}}, doi = {10.1098/rstb.2019.0528}, volume = {375}, year = {2020}, } @article{8167, abstract = {The evolution of strong reproductive isolation (RI) is fundamental to the origins and maintenance of biological diversity, especially in situations where geographical distributions of taxa broadly overlap. But what is the history behind strong barriers currently acting in sympatry? Using whole-genome sequencing and single nucleotide polymorphism genotyping, we inferred (i) the evolutionary relationships, (ii) the strength of RI, and (iii) the demographic history of divergence between two broadly sympatric taxa of intertidal snail. Despite being cryptic, based on external morphology, Littorina arcana and Littorina saxatilis differ in their mode of female reproduction (egg-laying versus brooding), which may generate a strong post-zygotic barrier. We show that egg-laying and brooding snails are closely related, but genetically distinct. Genotyping of 3092 snails from three locations failed to recover any recent hybrid or backcrossed individuals, confirming that RI is strong. There was, however, evidence for a very low level of asymmetrical introgression, suggesting that isolation remains incomplete. The presence of strong, asymmetrical RI was further supported by demographic analysis of these populations. Although the taxa are currently broadly sympatric, demographic modelling suggests that they initially diverged during a short period of geographical separation involving very low gene flow. Our study suggests that some geographical separation may kick-start the evolution of strong RI, facilitating subsequent coexistence of taxa in sympatry. The strength of RI needed to achieve sympatry and the subsequent effect of sympatry on RI remain open questions.}, author = {Stankowski, Sean and Westram, Anja M and Zagrodzka, Zuzanna B. and Eyres, Isobel and Broquet, Thomas and Johannesson, Kerstin and Butlin, Roger K.}, issn = {1471-2970}, journal = {Philosophical Transactions of the Royal Society. Series B: Biological Sciences}, number = {1806}, publisher = {The Royal Society}, title = {{The evolution of strong reproductive isolation between sympatric intertidal snails}}, doi = {10.1098/rstb.2019.0545}, volume = {375}, year = {2020}, } @article{8169, abstract = {Many recent studies have addressed the mechanisms operating during the early stages of speciation, but surprisingly few studies have tested theoretical predictions on the evolution of strong reproductive isolation (RI). To help address this gap, we first undertook a quantitative review of the hybrid zone literature for flowering plants in relation to reproductive barriers. Then, using Populus as an exemplary model group, we analysed genome-wide variation for phylogenetic tree topologies in both early- and late-stage speciation taxa to determine how these patterns may be related to the genomic architecture of RI. Our plant literature survey revealed variation in barrier complexity and an association between barrier number and introgressive gene flow. Focusing on Populus, our genome-wide analysis of tree topologies in speciating poplar taxa points to unusually complex genomic architectures of RI, consistent with earlier genome-wide association studies. These architectures appear to facilitate the ‘escape’ of introgressed genome segments from polygenic barriers even with strong RI, thus affecting their relationships with recombination rates. Placed within the context of the broader literature, our data illustrate how phylogenomic approaches hold great promise for addressing the evolution and temporary breakdown of RI during late stages of speciation.}, author = {Shang, Huiying and Hess, Jaqueline and Pickup, Melinda and Field, David and Ingvarsson, Pär K. and Liu, Jianquan and Lexer, Christian}, issn = {14712970}, journal = {Philosophical Transactions of the Royal Society. Series B: Biological Sciences}, number = {1806}, publisher = {The Royal Society}, title = {{Evolution of strong reproductive isolation in plants: Broad-scale patterns and lessons from a perennial model group}}, doi = {10.1098/rstb.2019.0544}, volume = {375}, year = {2020}, } @misc{9799, abstract = {Fitness interactions between mutations can influence a population’s evolution in many different ways. While epistatic effects are difficult to measure precisely, important information is captured by the mean and variance of log fitnesses for individuals carrying different numbers of mutations. We derive predictions for these quantities from a class of simple fitness landscapes, based on models of optimizing selection on quantitative traits. We also explore extensions to the models, including modular pleiotropy, variable effect sizes, mutational bias and maladaptation of the wild type. We illustrate our approach by reanalysing a large dataset of mutant effects in a yeast snoRNA. Though characterized by some large epistatic effects, these data give a good overall fit to the non-epistatic null model, suggesting that epistasis might have limited influence on the evolutionary dynamics in this system. We also show how the amount of epistasis depends on both the underlying fitness landscape and the distribution of mutations, and so is expected to vary in consistent ways between new mutations, standing variation and fixed mutations.}, author = {Fraisse, Christelle and Welch, John J.}, publisher = {Royal Society of London}, title = {{Simulation code for Fig S1 from the distribution of epistasis on simple fitness landscapes}}, doi = {10.6084/m9.figshare.7957469.v1}, year = {2020}, } @misc{9798, abstract = {Fitness interactions between mutations can influence a population’s evolution in many different ways. While epistatic effects are difficult to measure precisely, important information is captured by the mean and variance of log fitnesses for individuals carrying different numbers of mutations. We derive predictions for these quantities from a class of simple fitness landscapes, based on models of optimizing selection on quantitative traits. We also explore extensions to the models, including modular pleiotropy, variable effect sizes, mutational bias and maladaptation of the wild type. We illustrate our approach by reanalysing a large dataset of mutant effects in a yeast snoRNA. Though characterized by some large epistatic effects, these data give a good overall fit to the non-epistatic null model, suggesting that epistasis might have limited influence on the evolutionary dynamics in this system. We also show how the amount of epistasis depends on both the underlying fitness landscape and the distribution of mutations, and so is expected to vary in consistent ways between new mutations, standing variation and fixed mutations.}, author = {Fraisse, Christelle and Welch, John J.}, publisher = {Royal Society of London}, title = {{Simulation code for Fig S2 from the distribution of epistasis on simple fitness landscapes}}, doi = {10.6084/m9.figshare.7957472.v1}, year = {2020}, } @article{7236, abstract = {The biotic interactions hypothesis posits that biotic interactions are more important drivers of adaptation closer to the equator, evidenced by “stronger” contemporary interactions (e.g. greater interaction rates) and/or patterns of trait evolution consistent with a history of stronger interactions. Support for the hypothesis is mixed, but few studies span tropical and temperate regions while experimentally controlling for evolutionary history. Here, we integrate field observations and common garden experiments to quantify the relative importance of pollination and herbivory in a pair of tropical‐temperate congeneric perennial herbs. Phytolacca rivinoides and P. americana are pioneer species native to the Neotropics and the eastern USA, respectively. We compared plant‐pollinator and plant‐herbivore interactions between three tropical populations of P. rivinoides from Costa Rica and three temperate populations of P. americana from its northern range edge in Michigan and Ohio. For some metrics of interaction importance, we also included three subtropical populations of P. americana from its southern range edge in Florida. This approach confounds species and region but allows us, uniquely, to measure complementary proxies of interaction importance across a tropical‐temperate range in one system. To test the prediction that lower‐latitude plants are more reliant on insect pollinators, we quantified floral display and reward, insect visitation rates, and self‐pollination ability (autogamy). To test the prediction that lower‐latitude plants experience more herbivore pressure, we quantified herbivory rates, herbivore abundance, and leaf palatability. We found evidence supporting the biotic interactions hypothesis for most comparisons between P. rivinoides and north‐temperate P. americana (floral display, insect visitation, autogamy, herbivory, herbivore abundance, and young‐leaf palatability). Results for subtropical P. americana populations, however, were typically not intermediate between P. rivinoides and north‐temperate P. americana, as would be predicted by a linear latitudinal gradient in interaction importance. Subtropical young‐leaf palatability was intermediate, but subtropical mature leaves were the least palatable, and pollination‐related traits did not differ between temperate and subtropical regions. These nonlinear patterns of interaction importance suggest future work to relate interaction importance to climatic or biotic thresholds. In sum, we found that the biotic interactions hypothesis was more consistently supported at the larger spatial scale of our study.}, author = {Baskett, Carina and Schroeder, Lucy and Weber, Marjorie G. and Schemske, Douglas W.}, issn = {1557-7015}, journal = {Ecological Monographs}, number = {1}, publisher = {Wiley}, title = {{Multiple metrics of latitudinal patterns in insect pollination and herbivory for a tropical‐temperate congener pair}}, doi = {10.1002/ecm.1397}, volume = {90}, year = {2020}, } @article{7205, abstract = {Genetic incompatibilities contribute to reproductive isolation between many diverging populations, but it is still unclear to what extent they play a role if divergence happens with gene flow. In contact zones between the "Crab" and "Wave" ecotypes of the snail Littorina saxatilis, divergent selection forms strong barriers to gene flow, while the role of post‐zygotic barriers due to selection against hybrids remains unclear. High embryo abortion rates in this species could indicate the presence of such barriers. Post‐zygotic barriers might include genetic incompatibilities (e.g. Dobzhansky–Muller incompatibilities) but also maladaptation, both expected to be most pronounced in contact zones. In addition, embryo abortion might reflect physiological stress on females and embryos independent of any genetic stress. We examined all embryos of >500 females sampled outside and inside contact zones of three populations in Sweden. Females' clutch size ranged from 0 to 1,011 embryos (mean 130 ± 123), and abortion rates varied between 0% and 100% (mean 12%). We described female genotypes by using a hybrid index based on hundreds of SNPs differentiated between ecotypes with which we characterized female genotypes. We also calculated female SNP heterozygosity and inversion karyotype. Clutch size did not vary with female hybrid index, and abortion rates were only weakly related to hybrid index in two sites but not at all in a third site. No additional variation in abortion rate was explained by female SNP heterozygosity, but increased female inversion heterozygosity added slightly to increased abortion. Our results show only weak and probably biologically insignificant post‐zygotic barriers contributing to ecotype divergence, and the high and variable abortion rates were marginally, if at all, explained by hybrid index of females.}, author = {Johannesson, Kerstin and Zagrodzka, Zuzanna and Faria, Rui and Westram, Anja M and Butlin, Roger K.}, issn = {14209101}, journal = {Journal of Evolutionary Biology}, number = {3}, pages = {342--351}, publisher = {Wiley}, title = {{Is embryo abortion a post-zygotic barrier to gene flow between Littorina ecotypes?}}, doi = {10.1111/jeb.13570}, volume = {33}, year = {2020}, } @phdthesis{8574, abstract = {This thesis concerns itself with the interactions of evolutionary and ecological forces and the consequences on genetic diversity and the ultimate survival of populations. It is important to understand what signals processes leave on the genome and what we can infer from such data, which is usually abundant but noisy. Furthermore, understanding how and when populations adapt or go extinct is important for practical purposes, such as the genetic management of populations, as well as for theoretical questions, since local adaptation can be the first step toward speciation. In Chapter 2, we introduce the method of maximum entropy to approximate the demographic changes of a population in a simple setting, namely the logistic growth model with immigration. We show that this method is not only a powerful tool in physics but can be gainfully applied in an ecological framework. We investigate how well it approximates the real behavior of the system, and find that is does so, even in unexpected situations. Finally, we illustrate how it can model changing environments. In Chapter 3, we analyze the co-evolution of allele frequencies and population sizes in an infinite island model. We give conditions under which polygenic adaptation to a rare habitat is possible. The model we use is based on the diffusion approximation, considers eco-evolutionary feedback mechanisms (hard selection), and treats both drift and environmental fluctuations explicitly. We also look at limiting scenarios, for which we derive analytical expressions. In Chapter 4, we present a coalescent based simulation tool to obtain patterns of diversity in a spatially explicit subdivided population, in which the demographic history of each subpopulation can be specified. We compare the results to existing predictions, and explore the relative importance of time and space under a variety of spatial arrangements and demographic histories, such as expansion and extinction. In the last chapter, we give a brief outlook to further research. }, author = {Szep, Eniko}, issn = {2663-337X}, pages = {158}, publisher = {Institute of Science and Technology Austria}, title = {{Local adaptation in metapopulations}}, doi = {10.15479/AT:ISTA:8574}, year = {2020}, } @misc{8254, abstract = {Here are the research data underlying the publication "Estimating inbreeding and its effects in a long-term study of snapdragons (Antirrhinum majus)". Further information are summed up in the README document. The files for this record have been updated and are now found in the linked DOI https://doi.org/10.15479/AT:ISTA:9192.}, author = {Arathoon, Louise S}, publisher = {Institute of Science and Technology Austria}, title = {{Estimating inbreeding and its effects in a long-term study of snapdragons (Antirrhinum majus)}}, doi = {10.15479/AT:ISTA:8254}, year = {2020}, } @misc{9839, abstract = {More than 100 years after Grigg’s influential analysis of species’ borders, the causes of limits to species’ ranges still represent a puzzle that has never been understood with clarity. The topic has become especially important recently as many scientists have become interested in the potential for species’ ranges to shift in response to climate change—and yet nearly all of those studies fail to recognise or incorporate evolutionary genetics in a way that relates to theoretical developments. I show that range margins can be understood based on just two measurable parameters: (i) the fitness cost of dispersal—a measure of environmental heterogeneity—and (ii) the strength of genetic drift, which reduces genetic diversity. Together, these two parameters define an ‘expansion threshold’: adaptation fails when genetic drift reduces genetic diversity below that required for adaptation to a heterogeneous environment. When the key parameters drop below this expansion threshold locally, a sharp range margin forms. When they drop below this threshold throughout the species’ range, adaptation collapses everywhere, resulting in either extinction or formation of a fragmented metapopulation. Because the effects of dispersal differ fundamentally with dimension, the second parameter—the strength of genetic drift—is qualitatively different compared to a linear habitat. In two-dimensional habitats, genetic drift becomes effectively independent of selection. It decreases with ‘neighbourhood size’—the number of individuals accessible by dispersal within one generation. Moreover, in contrast to earlier predictions, which neglected evolution of genetic variance and/or stochasticity in two dimensions, dispersal into small marginal populations aids adaptation. This is because the reduction of both genetic and demographic stochasticity has a stronger effect than the cost of dispersal through increased maladaptation. The expansion threshold thus provides a novel, theoretically justified, and testable prediction for formation of the range margin and collapse of the species’ range.}, author = {Polechova, Jitka}, publisher = {Dryad}, title = {{Data from: Is the sky the limit? On the expansion threshold of a species' range}}, doi = {10.5061/dryad.5vv37}, year = {2019}, } @article{5911, abstract = {Empirical data suggest that inversions in many species contain genes important for intraspecific divergence and speciation, yet mechanisms of evolution remain unclear. While genes inside an inversion are tightly linked, inversions are not static but evolve separately from the rest of the genome by new mutations, recombination within arrangements, and gene flux between arrangements. Inversion polymorphisms are maintained by different processes, for example, divergent or balancing selection, or a mix of multiple processes. Moreover, the relative roles of selection, drift, mutation, and recombination will change over the lifetime of an inversion and within its area of distribution. We believe inversions are central to the evolution of many species, but we need many more data and new models to understand the complex mechanisms involved.}, author = {Faria, Rui and Johannesson, Kerstin and Butlin, Roger K. and Westram, Anja M}, issn = {01695347}, journal = {Trends in Ecology and Evolution}, number = {3}, pages = {239--248}, publisher = {Elsevier}, title = {{Evolving inversions}}, doi = {10.1016/j.tree.2018.12.005}, volume = {34}, year = {2019}, } @article{5680, abstract = {Pollinators display a remarkable diversity of foraging strategies with flowering plants, from primarily mutualistic interactions to cheating through nectar robbery. Despite numerous studies on the effect of nectar robbing on components of plant fitness, its contribution to reproductive isolation is unclear. We experimentally tested the impact of different pollinator strategies in a natural hybrid zone between two subspecies of Antirrhinum majus with alternate flower colour guides. On either side of a steep cline in flower colour between Antirrhinum majus pseudomajus (magenta) and A. m. striatum (yellow), we quantified the behaviour of all floral visitors at different time points during the flowering season. Using long-run camera surveys, we quantify the impact of nectar robbing on the number of flowers visited per inflorescence and the flower probing time. We further experimentally tested the effect of nectar robbing on female reproductive success by manipulating the intensity of robbing. While robbing increased over time the number of legitimate visitors tended to decrease concomitantly. We found that the number of flowers pollinated on a focal inflorescence decreased with the number of prior robbing events. However, in the manipulative experiment, fruit set and fruit volume did not vary significantly between low robbing and control treatments. Our findings challenge the idea that robbers have a negative impact on plant fitness through female function. This study also adds to our understanding of the components of pollinator-mediated reproductive isolation and the maintenance of Antirrhinum hybrid zones.}, author = {Andalo, Christophe and Burrus, Monique and Paute, Sandrine and Lauzeral, Christine and Field, David}, issn = {23818115}, journal = {Botany Letters}, number = {1}, pages = {80--92}, publisher = {Taylor and Francis}, title = {{Prevalence of legitimate pollinators and nectar robbers and the consequences for fruit set in an Antirrhinum majus hybrid zone}}, doi = {10.1080/23818107.2018.1545142}, volume = {166}, year = {2019}, } @article{6022, abstract = {The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.}, author = {Merrill, Richard M. and Rastas, Pasi and Martin, Simon H. and Melo Hurtado, Maria C and Barker, Sarah and Davey, John and Mcmillan, W. Owen and Jiggins, Chris D.}, journal = {PLoS Biology}, number = {2}, publisher = {Public Library of Science}, title = {{Genetic dissection of assortative mating behavior}}, doi = {10.1371/journal.pbio.2005902}, volume = {17}, year = {2019}, } @misc{9801, author = {Merrill, Richard M. and Rastas, Pasi and Martin, Simon H. and Melo Hurtado, Maria C and Barker, Sarah and Davey, John and Mcmillan, W. Owen and Jiggins, Chris D.}, publisher = {Public Library of Science}, title = {{Raw behavioral data}}, doi = {10.1371/journal.pbio.2005902.s006}, year = {2019}, } @article{6095, abstract = {Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients.}, author = {Faria, Rui and Chaube, Pragya and Morales, Hernán E. and Larsson, Tomas and Lemmon, Alan R. and Lemmon, Emily M. and Rafajlović, Marina and Panova, Marina and Ravinet, Mark and Johannesson, Kerstin and Westram, Anja M and Butlin, Roger K.}, issn = {1365-294X}, journal = {Molecular Ecology}, number = {6}, pages = {1375--1393}, publisher = {Wiley}, title = {{Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes}}, doi = {10.1111/mec.14972}, volume = {28}, year = {2019}, } @article{6230, abstract = {Great care is needed when interpreting claims about the genetic basis of human variation based on data from genome-wide association studies.}, author = {Barton, Nicholas H and Hermisson, Joachim and Nordborg, Magnus}, issn = {2050084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Why structure matters}}, doi = {10.7554/eLife.45380}, volume = {8}, year = {2019}, } @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 = {1365294X}, 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{6467, abstract = {Fitness interactions between mutations can influence a population’s evolution in many different ways. While epistatic effects are difficult to measure precisely, important information is captured by the mean and variance of log fitnesses for individuals carrying different numbers of mutations. We derive predictions for these quantities from a class of simple fitness landscapes, based on models of optimizing selection on quantitative traits. We also explore extensions to the models, including modular pleiotropy, variable effect sizes, mutational bias and maladaptation of the wild type. We illustrate our approach by reanalysing a large dataset of mutant effects in a yeast snoRNA (small nucleolar RNA). Though characterized by some large epistatic effects, these data give a good overall fit to the non-epistatic null model, suggesting that epistasis might have limited influence on the evolutionary dynamics in this system. We also show how the amount of epistasis depends on both the underlying fitness landscape and the distribution of mutations, and so is expected to vary in consistent ways between new mutations, standing variation and fixed mutations.}, author = {Fraisse, Christelle and Welch, John J.}, issn = {1744957X}, journal = {Biology Letters}, number = {4}, publisher = {Royal Society of London}, title = {{The distribution of epistasis on simple fitness landscapes}}, doi = {10.1098/rsbl.2018.0881}, volume = {15}, year = {2019}, } @article{6637, abstract = {The environment changes constantly at various time scales and, in order to survive, species need to keep adapting. Whether these species succeed in avoiding extinction is a major evolutionary question. Using a multilocus evolutionary model of a mutation‐limited population adapting under strong selection, we investigate the effects of the frequency of environmental fluctuations on adaptation. Our results rely on an “adaptive‐walk” approximation and use mathematical methods from evolutionary computation theory to investigate the interplay between fluctuation frequency, the similarity of environments, and the number of loci contributing to adaptation. First, we assume a linear additive fitness function, but later generalize our results to include several types of epistasis. We show that frequent environmental changes prevent populations from reaching a fitness peak, but they may also prevent the large fitness loss that occurs after a single environmental change. Thus, the population can survive, although not thrive, in a wide range of conditions. Furthermore, we show that in a frequently changing environment, the similarity of threats that a population faces affects the level of adaptation that it is able to achieve. We check and supplement our analytical results with simulations.}, author = {Trubenova, Barbora and Krejca, Martin and Lehre, Per Kristian and Kötzing, Timo}, journal = {Evolution}, number = {7}, pages = {1356--1374}, publisher = {Wiley}, title = {{Surfing on the seascape: Adaptation in a changing environment}}, doi = {10.1111/evo.13784}, volume = {73}, year = {2019}, } @article{6680, abstract = {This paper analyzes how partial selfing in a large source population influences its ability to colonize a new habitat via the introduction of a few founder individuals. Founders experience inbreeding depression due to partially recessive deleterious alleles as well as maladaptation to the new environment due to selection on a large number of additive loci. I first introduce a simplified version of the Inbreeding History Model (Kelly, 2007) in order to characterize mutation‐selection balance in a large, partially selfing source population under selection involving multiple non‐identical loci. I then use individual‐based simulations to study the eco‐evolutionary dynamics of founders establishing in the new habitat under a model of hard selection. The study explores how selfing rate shapes establishment probabilities of founders via effects on both inbreeding depression and adaptability to the new environment, and also distinguishes the effects of selfing on the initial fitness of founders from its effects on the long‐term adaptive response of the populations they found. A high rate of (but not complete) selfing is found to aid establishment over a wide range of parameters, even in the absence of mate limitation. The sensitivity of the results to assumptions about the nature of polygenic selection are discussed.}, author = {Sachdeva, Himani}, issn = {1558-5646}, journal = {Evolution}, number = {9}, pages = {1729--1745}, publisher = {Wiley}, title = {{Effect of partial selfing and polygenic selection on establishment in a new habitat}}, doi = {10.1111/evo.13812}, volume = {73}, year = {2019}, } @misc{9804, abstract = {Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response.}, author = {Castro, João Pl and Yancoskie, Michelle N. and Marchini, Marta and Belohlavy, Stefanie and Hiramatsu, Layla and Kučka, Marek and Beluch, William H. and Naumann, Ronald and Skuplik, Isabella and Cobb, John and Barton, Nicholas H and Rolian, Campbell and Chan, Yingguang Frank}, publisher = {Dryad}, title = {{Data from: An integrative genomic analysis of the Longshanks selection experiment for longer limbs in mice}}, doi = {10.5061/dryad.0q2h6tk}, year = {2019}, } @misc{9802, abstract = {This paper analyzes how partial selfing in a large source population influences its ability to colonize a new habitat via the introduction of a few founder individuals. Founders experience inbreeding depression due to partially recessive deleterious alleles as well as maladaptation to the new environment due to selection on a large number of additive loci. I first introduce a simplified version of the Inbreeding History Model (Kelly, 2007) in order to characterize mutation-selection balance in a large, partially selfing source population under selection involving multiple non-identical loci. I then use individual-based simulations to study the eco-evolutionary dynamics of founders establishing in the new habitat under a model of hard selection. The study explores how selfing rate shapes establishment probabilities of founders via effects on both inbreeding depression and adaptability to the new environment, and also distinguishes the effects of selfing on the initial fitness of founders from its effects on the long-term adaptive response of the populations they found. A high rate of (but not complete) selfing is found to aid establishment over a wide range of parameters, even in the absence of mate limitation. The sensitivity of the results to assumptions about the nature of polygenic selection are discussed.}, author = {Sachdeva, Himani}, publisher = {Dryad}, title = {{Data from: Effect of partial selfing and polygenic selection on establishment in a new habitat}}, doi = {10.5061/dryad.8tp0900}, year = {2019}, } @article{6795, abstract = {The green‐beard effect is one proposed mechanism predicted to underpin the evolu‐tion of altruistic behavior. It relies on the recognition and the selective help of altruists to each other in order to promote and sustain altruistic behavior. However, this mechanism has often been dismissed as unlikely or uncommon, as it is assumed that both the signaling trait and altruistic trait need to be encoded by the same gene or through tightly linked genes. Here, we use models of indirect genetic effects (IGEs) to find the minimum correlation between the signaling and altruistic trait required for the evolution of the latter. We show that this correlation threshold depends on the strength of the interaction (influence of the green beard on the expression of the altruistic trait), as well as the costs and benefits of the altruistic behavior. We further show that this correlation does not necessarily have to be high and support our analytical results by simulations.}, author = {Trubenova, Barbora and Hager, Reinmar}, issn = {20457758}, journal = {Ecology and Evolution}, number = {17}, pages = {9597--9608}, publisher = {Wiley}, title = {{Green beards in the light of indirect genetic effects}}, doi = {10.1002/ece3.5484}, volume = {9}, year = {2019}, } @article{6831, abstract = {* Understanding the mechanisms causing phenotypic differences between females and males has long fascinated evolutionary biologists. An extensive literature exists on animal sexual dimorphism but less information is known about sex differences in plants, particularly the extent of geographical variation in sexual dimorphism and its life‐cycle dynamics. * Here, we investigated patterns of genetically based sexual dimorphism in vegetative and reproductive traits of a wind‐pollinated dioecious plant, Rumex hastatulus, across three life‐cycle stages using open‐pollinated families from 30 populations spanning the geographic range and chromosomal variation (XY and XY1Y2) of the species. * The direction and degree of sexual dimorphism was highly variable among populations and life‐cycle stages. Sex‐specific differences in reproductive function explained a significant amount of temporal change in sexual dimorphism. For several traits, geographical variation in sexual dimorphism was associated with bioclimatic parameters, likely due to the differential responses of the sexes to climate. We found no systematic differences in sexual dimorphism between chromosome races. * Sex‐specific trait differences in dioecious plants largely result from a balance between sexual and natural selection on resource allocation. Our results indicate that abiotic factors associated with geographical context also play a role in modifying sexual dimorphism during the plant life‐cycle.}, author = {Puixeu Sala, Gemma and Pickup, Melinda and Field, David and Barrett, Spencer C.H.}, issn = {1469-8137}, journal = {New Phytologist}, number = {3}, pages = {1108--1120}, publisher = {Wiley}, title = {{Variation in sexual dimorphism in a wind-pollinated plant: The influence of geographical context and life-cycle dynamics}}, doi = {10.1111/nph.16050}, volume = {224}, year = {2019}, } @misc{9803, abstract = {Understanding the mechanisms causing phenotypic differences between females and males has long fascinated evolutionary biologists. An extensive literature exists on animal sexual dimorphism but less is known about sex differences in plants, particularly the extent of geographical variation in sexual dimorphism and its life-cycle dynamics. Here, we investigate patterns of genetically-based sexual dimorphism in vegetative and reproductive traits of a wind-pollinated dioecious plant, Rumex hastatulus, across three life-cycle stages using open-pollinated families from 30 populations spanning the geographic range and chromosomal variation (XY and XY1Y2) of the species. The direction and degree of sexual dimorphism was highly variable among populations and life-cycle stages. Sex-specific differences in reproductive function explained a significant amount of temporal change in sexual dimorphism. For several traits, geographical variation in sexual dimorphism was associated with bioclimatic parameters, likely due to the differential responses of the sexes to climate. We found no systematic differences in sexual dimorphism between chromosome races. Sex-specific trait differences in dioecious plants largely result from a balance between sexual and natural selection on resource allocation. Our results indicate that abiotic factors associated with geographical context also play a role in modifying sexual dimorphism during the plant life cycle.}, author = {Puixeu Sala, Gemma and Pickup, Melinda and Field, David and Barrett, Spencer C.H.}, publisher = {Dryad}, title = {{Data from: Variation in sexual dimorphism in a wind-pollinated plant: the influence of geographical context and life-cycle dynamics}}, doi = {10.5061/dryad.n1701c9}, year = {2019}, } @article{6855, abstract = {Many traits of interest are highly heritable and genetically complex, meaning that much of the variation they exhibit arises from differences at numerous loci in the genome. Complex traits and their evolution have been studied for more than a century, but only in the last decade have genome-wide association studies (GWASs) in humans begun to reveal their genetic basis. Here, we bring these threads of research together to ask how findings from GWASs can further our understanding of the processes that give rise to heritable variation in complex traits and of the genetic basis of complex trait evolution in response to changing selection pressures (i.e., of polygenic adaptation). Conversely, we ask how evolutionary thinking helps us to interpret findings from GWASs and informs related efforts of practical importance.}, author = {Sella, Guy and Barton, Nicholas H}, issn = {1545-293X}, journal = {Annual Review of Genomics and Human Genetics}, pages = {461--493}, publisher = {Annual Reviews}, title = {{Thinking about the evolution of complex traits in the era of genome-wide association studies}}, doi = {10.1146/annurev-genom-083115-022316}, volume = {20}, year = {2019}, } @article{6858, author = {Barton, Nicholas H}, issn = {2053-714X}, journal = {National Science Review}, number = {2}, pages = {291--292}, publisher = {Oxford University Press}, title = {{Is speciation driven by cycles of mixing and isolation?}}, doi = {10.1093/nsr/nwy113}, volume = {6}, year = {2019}, } @article{6857, abstract = {Gene Drives are regarded as future tools with a high potential for population control. Due to their inherent ability to overcome the rules of Mendelian inheritance, gene drives (GD) may spread genes rapidly through populations of sexually reproducing organisms. A release of organisms carrying a GD would constitute a paradigm shift in the handling of genetically modified organisms because gene drive organisms (GDO) are designed to drive their transgenes into wild populations and thereby increase the number of GDOs. The rapid development in this field and its focus on wild populations demand a prospective risk assessment with a focus on exposure related aspects. Presently, it is unclear how adequate risk management could be guaranteed to limit the spread of GDs in time and space, in order to avoid potential adverse effects in socio‐ecological systems. The recent workshop on the “Evaluation of Spatial and Temporal Control of Gene Drives” hosted by the Institute of Safety/Security and Risk Sciences (ISR) in Vienna aimed at gaining some insight into the potential population dynamic behavior of GDs and appropriate measures of control. Scientists from France, Germany, England, and the USA discussed both topics in this meeting on April 4–5, 2019. This article summarizes results of the workshop.}, author = {Giese, B and Friess, J L and Schetelig, M F and Barton, Nicholas H and Messer, Philip and Debarre, Florence and Meimberg, H and Windbichler, N and Boete, C}, issn = {1521-1878}, journal = {BioEssays}, number = {11}, publisher = {Wiley}, title = {{Gene Drives: Dynamics and regulatory matters – A report from the workshop “Evaluation of spatial and temporal control of Gene Drives”, 4 – 5 April 2019, Vienna}}, doi = {10.1002/bies.201900151}, volume = {41}, year = {2019}, } @misc{13067, abstract = {Genetic incompatibilities contribute to reproductive isolation between many diverging populations, but it is still unclear to what extent they play a role if divergence happens with gene flow. In contact zones between the "Crab" and "Wave" ecotypes of the snail Littorina saxatilis divergent selection forms strong barriers to gene flow, while the role of 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.}, author = {Johannesson, Kerstin and Zagrodzka, Zuzanna and Faria, Rui and Westram, Anja M and Butlin, Roger}, publisher = {Dryad}, title = {{Data from: Is embryo abortion a postzygotic barrier to gene flow between Littorina ecotypes?}}, doi = {10.5061/DRYAD.TB2RBNZWK}, year = {2019}, } @article{7393, abstract = {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.}, author = {Morales, Hernán E. and Faria, Rui and Johannesson, Kerstin and Larsson, Tomas and Panova, Marina and Westram, Anja M and Butlin, Roger K.}, issn = {2375-2548}, journal = {Science Advances}, number = {12}, publisher = {AAAS}, title = {{Genomic architecture of parallel ecological divergence: Beyond a single environmental contrast}}, doi = {10.1126/sciadv.aav9963}, volume = {5}, year = {2019}, } @inbook{8281, abstract = {We review the history of population genetics, starting with its origins a century ago from the synthesis between Mendel and Darwin's ideas, through to the recent development of sophisticated schemes of inference from sequence data, based on the coalescent. We explain the close relation between the coalescent and a diffusion process, which we illustrate by their application to understand spatial structure. We summarise the powerful methods available for analysis of multiple loci, when linkage equilibrium can be assumed, and then discuss approaches to the more challenging case, where associations between alleles require that we follow genotype, rather than allele, frequencies. Though we can hardly cover the whole of population genetics, we give an overview of the current state of the subject, and future challenges to it.}, author = {Barton, Nicholas H and Etheridge, Alison}, booktitle = {Handbook of statistical genomics}, editor = {Balding, David and Moltke, Ida and Marioni, John}, isbn = {9781119429142}, pages = {115--144}, publisher = {Wiley}, title = {{Mathematical models in population genetics}}, doi = {10.1002/9781119487845.ch4}, year = {2019}, }