@article{14850, abstract = {Elaborate sexual signals are thought to have evolved and be maintained to serve as honest indicators of signaller quality. One measure of quality is health, which can be affected by parasite infection. Cnemaspis mysoriensis is a diurnal gecko that is often infested with ectoparasites in the wild, and males of this species express visual (coloured gular patches) and chemical (femoral gland secretions) traits that receivers could assess during social interactions. In this paper, we tested whether ectoparasites affect individual health, and whether signal quality is an indicator of ectoparasite levels. In wild lizards, we found that ectoparasite level was negatively correlated with body condition in both sexes. Moreover, some characteristics of both visual and chemical traits in males were strongly associated with ectoparasite levels. Specifically, males with higher ectoparasite levels had yellow gular patches with lower brightness and chroma, and chemical secretions with a lower proportion of aromatic compounds. We then determined whether ectoparasite levels in males influence female behaviour. Using sequential choice trials, wherein females were provided with either the visual or the chemical signals of wild-caught males that varied in ectoparasite level, we found that only chemical secretions evoked an elevated female response towards less parasitised males. Simultaneous choice trials in which females were exposed to the chemical secretions from males that varied in parasite level further confirmed a preference for males with lower parasites loads. Overall, we find that although health (body condition) or ectoparasite load can be honestly advertised through multiple modalities, the parasite-mediated female response is exclusively driven by chemical signals.}, author = {Pal, Arka and Joshi, Mihir and Thaker, Maria}, issn = {0022-0949}, journal = {Journal of Experimental Biology}, keywords = {Insect Science, Molecular Biology, Animal Science and Zoology, Aquatic Science, Physiology, Ecology, Evolution, Behavior and Systematics}, number = {1}, publisher = {The Company of Biologists}, title = {{Too much information? Males convey parasite levels using more signal modalities than females utilise}}, doi = {10.1242/jeb.246217}, volume = {227}, year = {2024}, } @phdthesis{14711, abstract = {In nature, different species find their niche in a range of environments, each with its unique characteristics. While some thrive in uniform (homogeneous) landscapes where environmental conditions stay relatively consistent across space, others traverse the complexities of spatially heterogeneous terrains. Comprehending how species are distributed and how they interact within these landscapes holds the key to gaining insights into their evolutionary dynamics while also informing conservation and management strategies. For species inhabiting heterogeneous landscapes, when the rate of dispersal is low compared to spatial fluctuations in selection pressure, localized adaptations may emerge. Such adaptation in response to varying selection strengths plays an important role in the persistence of populations in our rapidly changing world. Hence, species in nature are continuously in a struggle to adapt to local environmental conditions, to ensure their continued survival. Natural populations can often adapt in time scales short enough for evolutionary changes to influence ecological dynamics and vice versa, thereby creating a feedback between evolution and demography. The analysis of this feedback and the relative contributions of gene flow, demography, drift, and natural selection to genetic variation and differentiation has remained a recurring theme in evolutionary biology. Nevertheless, the effective role of these forces in maintaining variation and shaping patterns of diversity is not fully understood. Even in homogeneous environments devoid of local adaptations, such understanding remains elusive. Understanding this feedback is crucial, for example in determining the conditions under which extinction risk can be mitigated in peripheral populations subject to deleterious mutation accumulation at the edges of species’ ranges as well as in highly fragmented populations. In this thesis we explore both uniform and spatially heterogeneous metapopulations, investigating and providing theoretical insights into the dynamics of local adaptation in the latter and examining the dynamics of load and extinction as well as the impact of joint ecological and evolutionary (eco-evolutionary) dynamics in the former. The thesis is divided into 5 chapters. Chapter 1 provides a general introduction into the subject matter, clarifying concepts and ideas used throughout the thesis. In chapter 2, we explore how fast a species distributed across a heterogeneous landscape adapts to changing conditions marked by alterations in carrying capacity, selection pressure, and migration rate. In chapter 3, we investigate how migration selection and drift influences adaptation and the maintenance of variation in a metapopulation with three habitats, an extension of previous models of adaptation in two habitats. We further develop analytical approximations for the critical threshold required for polymorphism to persist. The focus of chapter 4 of the thesis is on understanding the interplay between ecology and evolution as coupled processes. We investigate how eco-evolutionary feedback between migration, selection, drift, and demography influences eco-evolutionary outcomes in marginal populations subject to deleterious mutation accumulation. Using simulations as well as theoretical approximations of the coupled dynamics of population size and allele frequency, we analyze how gene flow from a large mainland source influences genetic load and population size on an island (i.e., in a marginal population) under genetically realistic assumptions. Analyses of this sort are important because small isolated populations, are repeatedly affected by complex interactions between ecological and evolutionary processes, which can lead to their death. Understanding these interactions can therefore provide an insight into the conditions under which extinction risk can be mitigated in peripheral populations thus, contributing to conservation and restoration efforts. Chapter 5 extends the analysis in chapter 4 to consider the dynamics of load (due to deleterious mutation accumulation) and extinction risk in a metapopulation. We explore the role of gene flow, selection, and dominance on load and extinction risk and further pinpoint critical thresholds required for metapopulation persistence. Overall this research contributes to our understanding of ecological and evolutionary mechanisms that shape species’ persistence in fragmented landscapes, a crucial foundation for successful conservation efforts and biodiversity management.}, author = {Olusanya, Oluwafunmilola O}, issn = {2663 - 337X}, pages = {183}, publisher = {Institute of Science and Technology Austria}, title = {{Local adaptation, genetic load and extinction in metapopulations}}, doi = {10.15479/at:ista:14711}, year = {2024}, } @article{14796, abstract = {Key innovations are fundamental to biological diversification, but their genetic basis is poorly understood. A recent transition from egg-laying to live-bearing in marine snails (Littorina spp.) provides the opportunity to study the genetic architecture of an innovation that has evolved repeatedly across animals. Individuals do not cluster by reproductive mode in a genome-wide phylogeny, but local genealogical analysis revealed numerous small genomic regions where all live-bearers carry the same core haplotype. Candidate regions show evidence for live-bearer–specific positive selection and are enriched for genes that are differentially expressed between egg-laying and live-bearing reproductive systems. Ages of selective sweeps suggest that live-bearer–specific alleles accumulated over more than 200,000 generations. Our results suggest that new functions evolve through the recruitment of many alleles rather than in a single evolutionary step.}, author = {Stankowski, Sean and Zagrodzka, Zuzanna B. and Garlovsky, Martin D. and Pal, Arka and Shipilina, Daria and Garcia Castillo, Diego Fernando and Lifchitz, Hila and Le Moan, Alan and Leder, Erica and Reeve, James and Johannesson, Kerstin and Westram, Anja M and Butlin, Roger K.}, issn = {1095-9203}, journal = {Science}, number = {6678}, pages = {114--119}, publisher = {American Association for the Advancement of Science}, title = {{The genetic basis of a recent transition to live-bearing in marine snails}}, doi = {10.1126/science.adi2982}, volume = {383}, year = {2024}, } @phdthesis{15020, abstract = {This thesis consists of four distinct pieces of work within theoretical biology, with two themes in common: the concept of optimization in biological systems, and the use of information-theoretic tools to quantify biological stochasticity and statistical uncertainty. Chapter 2 develops a statistical framework for studying biological systems which we believe to be optimized for a particular utility function, such as retinal neurons conveying information about visual stimuli. We formalize such beliefs as maximum-entropy Bayesian priors, constrained by the expected utility. We explore how such priors aid inference of system parameters with limited data and enable optimality hypothesis testing: is the utility higher than by chance? Chapter 3 examines the ultimate biological optimization process: evolution by natural selection. As some individuals survive and reproduce more successfully than others, populations evolve towards fitter genotypes and phenotypes. We formalize this as accumulation of genetic information, and use population genetics theory to study how much such information can be accumulated per generation and maintained in the face of random mutation and genetic drift. We identify the population size and fitness variance as the key quantities that control information accumulation and maintenance. Chapter 4 reuses the concept of genetic information from Chapter 3, but from a different perspective: we ask how much genetic information organisms actually need, in particular in the context of gene regulation. For example, how much information is needed to bind transcription factors at correct locations within the genome? Population genetics provides us with a refined answer: with an increasing population size, populations achieve higher fitness by maintaining more genetic information. Moreover, regulatory parameters experience selection pressure to optimize the fitness-information trade-off, i.e. minimize the information needed for a given fitness. This provides an evolutionary derivation of the optimization priors introduced in Chapter 2. Chapter 5 proves an upper bound on mutual information between a signal and a communication channel output (such as neural activity). Mutual information is an important utility measure for biological systems, but its practical use can be difficult due to the large dimensionality of many biological channels. Sometimes, a lower bound on mutual information is computed by replacing the high-dimensional channel outputs with decodes (signal estimates). Our result provides a corresponding upper bound, provided that the decodes are the maximum posterior estimates of the signal.}, author = {Hledik, Michal}, issn = {2663 - 337X}, keywords = {Theoretical biology, Optimality, Evolution, Information}, pages = {158}, publisher = {Institute of Science and Technology Austria}, title = {{Genetic information and biological optimization}}, doi = {10.15479/at:ista:15020}, year = {2024}, } @article{15099, abstract = {Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.}, author = {Johannesson, Kerstin and Faria, Rui and Le Moan, Alan and Rafajlović, Marina and Westram, Anja M and Butlin, Roger K. and Stankowski, Sean}, issn = {1362-4555}, journal = {Trends in Genetics}, publisher = {Cell Press}, title = {{Diverse pathways to speciation revealed by marine snails}}, doi = {10.1016/j.tig.2024.01.002}, year = {2024}, } @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}, } @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}, } @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}, } @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}, }