TY - THES AB - Although we often see studies focusing on simple or even discrete traits in studies of colouration, the variation of “appearance” phenotypes found in nature is often more complex, continuous and high-dimensional. Therefore, we developed automated methods suitable for large datasets of genomes and images, striving to account for their complex nature, while minimising human bias. We used these methods on a dataset of more than 20, 000 plant SNP genomes and corresponding fower images from a hybrid zone of two subspecies of Antirrhinum majus with distinctly coloured fowers to improve our understanding of the genetic nature of the fower colour in our study system. Firstly, we use the advantage of large numbers of genotyped plants to estimate the haplotypes in the main fower colour regulating region. We study colour- and geography-related characteristics of the estimated haplotypes and how they connect to their relatedness. We show discrepancies from the expected fower colour distributions given the genotype and identify particular haplotypes leading to unexpected phenotypes. We also confrm a signifcant defcit of the double recessive recombinant and quite surprisingly, we show that haplotypes of the most frequent parental type are much less variable than others. Secondly, we introduce our pipeline capable of processing tens of thousands of full fower images without human interaction and summarising each image into a set of informative scores. We show the compatibility of these machine-measured fower colour scores with the previously used manual scores and study impact of external efect on the resulting scores. Finally, we use the machine-measured fower colour scores to ft and examine a phenotype cline across the hybrid zone in Planoles using full fower images as opposed to discrete, manual scores and compare it with the genotypic cline. AU - Matejovicova, Lenka ID - 11128 SN - 2663-337X TI - Genetic basis of flower colour as a model for adaptive evolution ER - TY - JOUR AB - Maternally inherited Wolbachia transinfections are being introduced into natural mosquito populations to reduce the transmission of dengue, Zika, and other arboviruses. Wolbachia-induced cytoplasmic incompatibility provides a frequency-dependent reproductive advantage to infected females that can spread transinfections within and among populations. However, because transinfections generally reduce host fitness, they tend to spread within populations only after their frequency exceeds a critical threshold. This produces bistability with stable equilibrium frequencies at both 0 and 1, analogous to the bistability produced by underdominance between alleles or karyotypes and by population dynamics under Allee effects. Here, we analyze how stochastic frequency variation produced by finite population size can facilitate the local spread of variants with bistable dynamics into areas where invasion is unexpected from deterministic models. Our exemplar is the establishment of wMel Wolbachia in the Aedes aegypti population of Pyramid Estates (PE), a small community in far north Queensland, Australia. In 2011, wMel was stably introduced into Gordonvale, separated from PE by barriers to A. aegypti dispersal. After nearly 6 years during which wMel was observed only at low frequencies in PE, corresponding to an apparent equilibrium between immigration and selection, wMel rose to fixation by 2018. Using analytic approximations and statistical analyses, we demonstrate that the observed fixation of wMel at PE is consistent with both stochastic transition past an unstable threshold frequency and deterministic transformation produced by steady immigration at a rate just above the threshold required for deterministic invasion. The indeterminacy results from a delicate balance of parameters needed to produce the delayed transition observed. Our analyses suggest that once Wolbachia transinfections are established locally through systematic introductions, stochastic “threshold crossing” is likely to only minimally enhance spatial spread, providing a local ratchet that slightly—but systematically—aids area-wide transformation of disease-vector populations in heterogeneous landscapes. AU - Turelli, Michael AU - Barton, Nicholas H ID - 10604 IS - 1 JF - Evolution Letters KW - genetics KW - ecology KW - evolution KW - behavior and systematics TI - Why did the Wolbachia transinfection cross the road? Drift, deterministic dynamics, and disease control VL - 6 ER - TY - GEN AB - Maternally inherited Wolbachia transinfections are being introduced into natural mosquito populations to reduce the transmission of dengue, Zika and other arboviruses. Wolbachia-induced cytoplasmic incompatibility provides a frequency-dependent reproductive advantage to infected females that can spread transinfections within and among populations. However, because transinfections generally reduce host fitness, they tend to spread within populations only after their frequency exceeds a critical threshold. This produces bistability with stable equilibrium frequencies at both 0 and 1, analogous to the bistability produced by underdominance between alleles or karyotypes and by population dynamics under Allee effects. Here, we analyze how stochastic frequency variation produced by finite population size can facilitate the local spread of variants with bistable dynamics into areas where invasion is unexpected from deterministic models. Our exemplar is the establishment of wMel Wolbachia in the Aedes aegypti population of Pyramid Estates (PE), a small community in far north Queensland, Australia. In 2011, wMel was stably introduced into Gordonvale, separated from PE by barriers to Ae. aegypti dispersal. After nearly six years during which wMel was observed only at low frequencies in PE, corresponding to an apparent equilibrium between immigration and selection, wMel rose to fixation by 2018. Using analytic approximations and statistical analyses, we demonstrate that the observed fixation of wMel at PE is consistent with both stochastic transition past an unstable threshold frequency and deterministic transformation produced by steady immigration at a rate just above the threshold required for deterministic invasion. The indeterminacy results from a delicate balance of parameters needed to produce the delayed transition observed. Our analyses suggest that once Wolbachia transinfections are established locally through systematic introductions, stochastic “threshold crossing” is likely to only minimally enhance spatial spread, providing a local ratchet that slightly – but systematically – aids area-wide transformation of disease-vector populations in heterogeneous landscapes. AU - Turelli, Michael AU - Barton, Nicholas H ID - 11686 KW - Biological sciences TI - Wolbachia frequency data from: Why did the Wolbachia transinfection cross the road? Drift, deterministic dynamics and disease control ER - TY - JOUR AB - Predicting function from sequence is a central problem of biology. Currently, this is possible only locally in a narrow mutational neighborhood around a wildtype sequence rather than globally from any sequence. Using random mutant libraries, we developed a biophysical model that accounts for multiple features of σ70 binding bacterial promoters to predict constitutive gene expression levels from any sequence. We experimentally and theoretically estimated that 10–20% of random sequences lead to expression and ~80% of non-expressing sequences are one mutation away from a functional promoter. The potential for generating expression from random sequences is so pervasive that selection acts against σ70-RNA polymerase binding sites even within inter-genic, promoter-containing regions. This pervasiveness of σ70-binding sites implies that emergence of promoters is not the limiting step in gene regulatory evolution. Ultimately, the inclusion of novel features of promoter function into a mechanistic model enabled not only more accurate predictions of gene expression levels, but also identified that promoters evolve more rapidly than previously thought. AU - Lagator, Mato AU - Sarikas, Srdjan AU - Steinrueck, Magdalena AU - Toledo-Aparicio, David AU - Bollback, Jonathan P AU - Guet, Calin C AU - Tkačik, Gašper ID - 10736 JF - eLife TI - Predicting bacterial promoter function and evolution from random sequences VL - 11 ER - TY - JOUR AB - Hybridization is a common evolutionary process with multiple possible outcomes. In vertebrates, interspecific hybridization has repeatedly generated parthenogenetic hybrid species. However, it is unknown whether the generation of parthenogenetic hybrids is a rare outcome of frequent hybridization between sexual species within a genus or the typical outcome of rare hybridization events. Darevskia is a genus of rock lizards with both hybrid parthenogenetic and sexual species. Using capture sequencing, we estimate phylogenetic relationships and gene flow among the sexual species, to determine how introgressive hybridization relates to the origins of parthenogenetic hybrids. We find evidence for widespread hybridization with gene flow, both between recently diverged species and deep branches. Surprisingly, we find no signal of gene flow between parental species of the parthenogenetic hybrids, suggesting that the parental pairs were either reproductively or geographically isolated early in their divergence. The generation of parthenogenetic hybrids in Darevskia is, then, a rare outcome of the total occurrence of hybridization within the genus, but the typical outcome when specific species pairs hybridize. Our results question the conventional view that parthenogenetic lineages are generated by hybridization in a window of divergence. Instead, they suggest that some lineages possess specific properties that underpin successful parthenogenetic reproduction. AU - Freitas, Susana AU - Westram, Anja M AU - Schwander, Tanja AU - Arakelyan, Marine AU - Ilgaz, Çetin AU - Kumlutas, Yusuf AU - Harris, David James AU - Carretero, Miguel A. AU - Butlin, Roger K. ID - 11334 IS - 5 JF - Evolution SN - 0014-3820 TI - Parthenogenesis in Darevskia lizards: A rare outcome of common hybridization, not a common outcome of rare hybridization VL - 76 ER - TY - JOUR AB - Empirical essays of fitness landscapes suggest that they may be rugged, that is having multiple fitness peaks. Such fitness landscapes, those that have multiple peaks, necessarily have special local structures, called reciprocal sign epistasis (Poelwijk et al. in J Theor Biol 272:141–144, 2011). Here, we investigate the quantitative relationship between the number of fitness peaks and the number of reciprocal sign epistatic interactions. Previously, it has been shown (Poelwijk et al. in J Theor Biol 272:141–144, 2011) that pairwise reciprocal sign epistasis is a necessary but not sufficient condition for the existence of multiple peaks. Applying discrete Morse theory, which to our knowledge has never been used in this context, we extend this result by giving the minimal number of reciprocal sign epistatic interactions required to create a given number of peaks. AU - Saona Urmeneta, Raimundo J AU - Kondrashov, Fyodor AU - Khudiakova, Kseniia ID - 11447 IS - 8 JF - Bulletin of Mathematical Biology KW - Computational Theory and Mathematics KW - General Agricultural and Biological Sciences KW - Pharmacology KW - General Environmental Science KW - General Biochemistry KW - Genetics and Molecular Biology KW - General Mathematics KW - Immunology KW - General Neuroscience SN - 0092-8240 TI - Relation between the number of peaks and the number of reciprocal sign epistatic interactions VL - 84 ER - TY - JOUR AB - Local adaptation leads to differences between populations within a species. In many systems, similar environmental contrasts occur repeatedly, sometimes driving parallel phenotypic evolution. Understanding the genomic basis of local adaptation and parallel evolution is a major goal of evolutionary genomics. It is now known that by preventing the break-up of favourable combinations of alleles across multiple loci, genetic architectures that reduce recombination, like chromosomal inversions, can make an important contribution to local adaptation. However, little is known about whether inversions also contribute disproportionately to parallel evolution. Our aim here is to highlight this knowledge gap, to showcase existing studies, and to illustrate the differences between genomic architectures with and without inversions using simple models. We predict that by generating stronger effective selection, inversions can sometimes speed up the parallel adaptive process or enable parallel adaptation where it would be impossible otherwise, but this is highly dependent on the spatial setting. We highlight that further empirical work is needed, in particular to cover a broader taxonomic range and to understand the relative importance of inversions compared to genomic regions without inversions. AU - Westram, Anja M AU - Faria, Rui AU - Johannesson, Kerstin AU - Butlin, Roger AU - Barton, Nicholas H ID - 11546 IS - 1856 JF - Philosophical Transactions of the Royal Society B: Biological Sciences KW - General Agricultural and Biological Sciences KW - General Biochemistry KW - Genetics and Molecular Biology SN - 0962-8436 TI - Inversions and parallel evolution VL - 377 ER - TY - JOUR AB - Spatially explicit population genetic models have long been developed, yet have rarely been used to test hypotheses about the spatial distribution of genetic diversity or the genetic divergence between populations. Here, we use spatially explicit coalescence simulations to explore the properties of the island and the two-dimensional stepping stone models under a wide range of scenarios with spatio-temporal variation in deme size. We avoid the simulation of genetic data, using the fact that under the studied models, summary statistics of genetic diversity and divergence can be approximated from coalescence times. We perform the simulations using gridCoal, a flexible spatial wrapper for the software msprime (Kelleher et al., 2016, Theoretical Population Biology, 95, 13) developed herein. In gridCoal, deme sizes can change arbitrarily across space and time, as well as migration rates between individual demes. We identify different factors that can cause a deviation from theoretical expectations, such as the simulation time in comparison to the effective deme size and the spatio-temporal autocorrelation across the grid. Our results highlight that FST, a measure of the strength of population structure, principally depends on recent demography, which makes it robust to temporal variation in deme size. In contrast, the amount of genetic diversity is dependent on the distant past when Ne is large, therefore longer run times are needed to estimate Ne than FST. Finally, we illustrate the use of gridCoal on a real-world example, the range expansion of silver fir (Abies alba Mill.) since the last glacial maximum, using different degrees of spatio-temporal variation in deme size. AU - Szep, Eniko AU - Trubenova, Barbora AU - Csilléry, Katalin ID - 11640 IS - 8 JF - Molecular Ecology Resources SN - 1755-098X TI - Using gridCoal to assess whether standard population genetic theory holds in the presence of spatio-temporal heterogeneity in population size VL - 22 ER - TY - JOUR AB - Sexual antagonism is a common hypothesis for driving the evolution of sex chromosomes, whereby recombination suppression is favored between sexually antagonistic loci and the sex-determining locus to maintain beneficial combinations of alleles. This results in the formation of a sex-determining region. Chromosomal inversions may contribute to recombination suppression but their precise role in sex chromosome evolution remains unclear. Because local adaptation is frequently facilitated through the suppression of recombination between adaptive loci by chromosomal inversions, there is potential for inversions that cover sex-determining regions to be involved in local adaptation as well, particularly if habitat variation creates environment-dependent sexual antagonism. With these processes in mind, we investigated sex determination in a well-studied example of local adaptation within a species: the intertidal snail, Littorina saxatilis. Using SNP data from a Swedish hybrid zone, we find novel evidence for a female-heterogametic sex determination system that is restricted to one ecotype. Our results suggest that four putative chromosomal inversions, two previously described and two newly discovered, span the putative sex chromosome pair. We determine their differing associations with sex, which suggest distinct strata of differing ages. The same inversions are found in the second ecotype but do not show any sex association. The striking disparity in inversion-sex associations between ecotypes that are connected by gene flow across a habitat transition that is just a few meters wide indicates a difference in selective regime that has produced a distinct barrier to the spread of the newly discovered sex-determining region between ecotypes. Such sex chromosome-environment interactions have not previously been uncovered in L. saxatilis and are known in few other organisms. A combination of both sex-specific selection and divergent natural selection is required to explain these highly unusual patterns. AU - Hearn, Katherine E. AU - Koch, Eva L. AU - Stankowski, Sean AU - Butlin, Roger K. AU - Faria, Rui AU - Johannesson, Kerstin AU - Westram, Anja M ID - 12001 IS - 5 JF - Evolution Letters TI - Differing associations between sex determination and sex-linked inversions in two ecotypes of Littorina saxatilis VL - 6 ER - TY - JOUR AB - Polygenic adaptation is thought to be ubiquitous, yet remains poorly understood. Here, we model this process analytically, in the plausible setting of a highly polygenic, quantitative trait that experiences a sudden shift in the fitness optimum. We show how the mean phenotype changes over time, depending on the effect sizes of loci that contribute to variance in the trait, and characterize the allele dynamics at these loci. Notably, we describe the two phases of the allele dynamics: The first is a rapid phase, in which directional selection introduces small frequency differences between alleles whose effects are aligned with or opposed to the shift, ultimately leading to small differences in their probability of fixation during a second, longer phase, governed by stabilizing selection. As we discuss, key results should hold in more general settings and have important implications for efforts to identify the genetic basis of adaptation in humans and other species. AU - Hayward, Laura AU - Sella, Guy ID - 12157 JF - eLife KW - General Immunology and Microbiology KW - General Biochemistry KW - Genetics and Molecular Biology KW - General Medicine KW - General Neuroscience TI - Polygenic adaptation after a sudden change in environment VL - 11 ER -