Evolution in heterogeneous environments: between soft and hard selection

Since their first formulations about half a century ago, the soft and hard selection models have become classical frameworks to study selection in subdivided populations. These models differ in the timing of density regulation and represent two extreme types of selection: density- and frequency-dependent selection (soft) and density- and frequency-independent selection (hard). Yet only few attempts have been made so far to model intermediate scenarios. Here, we design a model where migration may happen twice during the life cycle: before density regulation with probability d(J) (juvenile migration) and after density regulation with probability d(A) (adult migration). In the first step, we analyze the conditions for the coexistence of two specialists. We find that coexistence is possible under a large range of selection types, even when environmental heterogeneity is low. Then, we investigate the different possible outcomes obtained through gradual evolution. We show that polymorphism is more likely to evolve when the trade-off is weak, environmental heterogeneity is high, migration is low, and in particular when juvenile migration is low relative to adult migration, because the timing of migration affects the magnitude of frequency-dependent selection relative to gene flow. This model may provide a more general theoretical framework to experimentally study evolution in heterogeneous environments.     

Adaptive evolution of the Populus tremula photoperiod pathway

Perennial plants monitor seasonal changes through changes in environmental conditions such as the quantity and quality of light and genes in the photoperiodic pathway are known to be involved in controlling these processes. Here, we examine 25 of genes from the photoperiod pathway in Populus tremula (Salicaceae) for signatures of adaptive evolution. Overall, levels of synonymous polymorphism in the 25 genes are lower than at control loci selected randomly from the genome. This appears primarily to be caused by lower levels of synonymous polymorphism in genes associated with the circadian clock. Natural selection appears to play an important role in shaping protein evolution at several of the genes in the photoperiod pathways, which is highlighted by the fact that approximately 40% of the genes from the photoperiod pathway have estimates of selection on nonsynonymous polymorphisms that are significantly different from zero. A surprising observation we make is that circadian clock-associated genes appear to be over-represented among the genes showing elevated rates of protein evolution; seven genes are evolving under positive selection and all but one of these genes are involved in the circadian clock of Populus.     

Solvent accessibility and purifying selection within proteins of Escherichia coli and Salmonella enterica

The neutral theory of molecular evolution predicts that variation within species is inversely related to the strength of purifying selection, but the strength of purifying selection itself must be related to physical constraints imposed by protein folding and function. In this paper, we analyzed five enzymes for which polymorphic sequence variation within Escherichia coli and/or Salmonella enterica was available, along with a protein structure. Single and multivariate logistic regression models are presented that evaluate amino acid size, physicochemical properties, solvent accessibility, and secondary structure as predictors of polymorphism. A model that contains a positive coefficient of association between polymorphism and solvent accessibility and separate intercepts for each secondary-structure element is sufficient to explain the observed variation in polymorphism between sites. The model predicts an increase in the probability of amino acid polymorphism with increasing solvent accessibility for each protein regardless of physicochemical properties, secondary-structure element, or size of the amino acid. This result, when compared with the distribution of synonymous polymorphism, which shows no association with solvent accessibility, suggests a strong decrease in purifying selection with increasing solvent accessibility.     

ADAPTIVE DYNAMICS IN ALLELE SPACE: EVOLUTION OF GENETIC POLYMORPHISM BY SMALL MUTATIONS IN A HETEROGENEOUS ENVIRONMENT

We demonstrate how a genetic polymorphism of distinctly different alleles can develop during long-term frequency-dependent evolution in an initially monomorphic diploid population, if mutations have only small phenotypic effect. As a specific example, we use a version of Levene's (1953) soft selection model, where stabilizing selection acts on a continuous trait within each of two habitats. If the optimal phenotypes within the habitats are sufficiently different, then two distinctly different alleles evolve gradually from a single ancestral allele. In a wide range of parameter values, the two locally optimal phenotypes will be realized by one of the homozygotes and the heterozygote, rather than by the two homozygotes. Unlike in the haploid analogue of the model, there can be multiple polymorphic evolutionary attractors with different probabilities of convergence. Our results differ from the population genetic models of short-term evolution in two aspects: (1) a polymorphism that is population genetically stable may be invaded by a new mutant allele and, as a consequence, the population may fall back to monomorphism, (2) long-term evolution by allele substitutions may lead from a population where polymorphism is not possible into one where polymorphism is possible.     

Testing the role of ecological selection on colour pattern variation in the butterfly Parnassius clodius

Colour pattern has served as an important phenotype in understanding the process of natural selection, particularly in brightly coloured and variable species like butterflies. However, different selective forces operate on aspects of colour pattern, for example by favouring warning colours in eyespots or alternatively favoring investment in thermoregulatory properties of melanin. Additionally, genetic drift influences colour phenotypes, especially in populations undergoing population size change. Here, we investigated the relative roles of genetic drift and ecological selection in generating the phenotypic diversity of the butterfly Parnassius clodius. Genome‐wide patterns of single nucleotide polymorphism data show that P. clodius forms three population clusters, which experienced a period of population expansion following the last glacial maximum and have since remained relatively stable in size. After correcting for relatedness, morphological variation is best explained by climatic predictor variables, suggesting ecological selection generates trait variability. Solar radiation and precipitation are both negatively correlated with increasing total melanin in both sexes, supporting a thermoregulatory function of melanin. Similarly, wing size traits are significantly larger in warmer habitats for both sexes, supporting a Converse Bergmann Rule pattern. Bright red coloration is negatively correlated with temperature seasonality and solar radiation in males, and weakly associated with insectivorous avian predators in univariate models, providing mixed evidence that selection is linked to warning coloration and predator avoidance. Together, these results suggest that elements of butterfly wing phenotypes respond independently to different sources of selection and that thermoregulation is an important driver of phenotypic differentiation in Parnassian butterflies.     

Frequency-dependent selection and the maintenance of genetic variation: exploring the parameter space of the multiallelic pairwise interaction model

When individuals' fitnesses depend on the genetic composition of the population in which they are found, selection is then frequency dependent. Frequency-dependent selection (FDS) is often invoked as a heuristic explanation for the maintenance of large numbers of alleles at a locus. The pairwise interaction model is a general model of FDS via intraspecific competition at the genotypic level. Here we use a parameter-space approach to investigate the full potential for the maintenance of multiallelic equilibria under the pairwise interaction model. We find that FDS maintains full polymorphism more often than classic constant-selection models and produces more skewed equilibrium allele frequencies. Fitness sets with some degree of rare advantage maintained full polymorphism most often, but a wide variety of nonobvious fitness patterns were also found to have positive potential for polymorphism. An example is put forth suggesting possible explanations for multiallelic polymorphisms maintained despite positive FDS on individual alleles.     

Polymorphism in heterogeneous environments,evolution of habitat selection and sympatric speciation: Soft and hard selection models

Summary The adaptation to a variable environment has been studied within soft and hard selection frameworks. It is shown that an epistatically determined habitat preference, following a Markovian process, always leads to the maintenance of an adaptive polymorphism, in a soft selection context. Although local mating does not alter the conditions for polymorphism maintenance, it is shown that, in that case, habitat selection also leads to the evolution of isolated reproductive units within each available habitat. Habitat selection, however, cannot evolve in the total absence of adaptive polymorphism. This represents a theoretical problem for all models assuming habitat selection to be an initially fixed trait, and means that within a soft selection framework, all the available habitats will be exploited, even the less favourable ones.On the other hand, polymorphism cannot be maintained when selection is hard, even when all individuals select their habitat. Here, the evolution of habitat selection does not need any prerequisite polymorphism, and always leads to the exploitation of only one habitat by the most specialized genotype. It appears then that hard selection can account for the existence of empty habitat and for an easier evolution of habitat specialization.     

Ancient and recent adaptive evolution of primate non-homologous end joining genes

In human cells, DNA double-strand breaks are repaired primarily by the non-homologous end joining (NHEJ) pathway. Given their critical nature, we expected NHEJ proteins to be evolutionarily conserved, with relatively little sequence change over time. Here, we report that while critical domains of these proteins are conserved as expected, the sequence of NHEJ proteins has also been shaped by recurrent positive selection, leading to rapid sequence evolution in other protein domains. In order to characterize the molecular evolution of the human NHEJ pathway, we generated large simian primate sequence datasets for NHEJ genes. Codon-based models of gene evolution yielded statistical support for the recurrent positive selection of five NHEJ genes during primate evolution: XRCC4, NBS1, Artemis, POLλ, and CtIP. Analysis of human polymorphism data using the composite of multiple signals (CMS) test revealed that XRCC4 has also been subjected to positive selection in modern humans. Crystal structures are available for XRCC4, Nbs1, and Polλ; and residues under positive selection fall exclusively on the surfaces of these proteins. Despite the positive selection of such residues, biochemical experiments with variants of one positively selected site in Nbs1 confirm that functions necessary for DNA repair and checkpoint signaling have been conserved. However, many viruses interact with the proteins of the NHEJ pathway as part of their infectious lifecycle. We propose that an ongoing evolutionary arms race between viruses and NHEJ genes may be driving the surprisingly rapid evolution of these critical genes.     

The butterfly Danaus chrysippus (L.) in East Africa: polymorphism and morph-ratio clines within a complex, extensive and dynamic hybrid zone

Samples of the polymorphic butterfly Danaus chrysippus are analysed from six well separated sites in East Africa. Morph-ratio clines are described for four diallelic genes A, B, C and L, each of which influences the visual phenotype. Each of the four clines has a different orientation, consistent with an hypothesis that the polymorphism originated from hybridization between a number of polytypic denies which have at various times undergone range expansion. Allopatric subspeciation in isolated Pleistocene refugia is postulated. The phenotype of each geographical race is shared with one of the morphs within the hybrid zone; other sympatrically maintained polymorphic forms are normally confined to the hybrid zone. Wright's isolation-by-distance model best explains the present distribution of gene frequencies. Morph-ratios differ significantly between the sexes and are sometimes associated with heterozygote excess; gametic and genotypic disequilibria are general throughout the region and suggest the clines are maintained by strong natural selection. Seasonal cycling of phenotype frequency is believed to result from extensive migratory movements rather than natural selection. Female-biased sex-ratio, which is also seasonal, and Haldanc rule effects, result from hybrid breakdown when genetically distinct demes meet and interbreed. Oscillating sex-ratios and frequency of colour genes are functionally linked by negative feedback. The polymorphism owes its origin to allopatrie evolution but is now maintained sympatrically.     

Multiregional comparison of the ecological and phylogenetic structure of butterfly species richness gradients

Aim To examine butterfly species richness gradients in seven regions/countries and to quantify geographic mean root distance (MRD) patterns. My primary goal is to determine the extent to which an explanation for butterfly richness patterns based on tropical niche conservatism and the evolution of cold tolerance, proposed for the fauna of Canada and the USA, applies to other parts of the world. Location USA/Canada, Mexico, Europe/NW Africa, Transbaikal Siberia, Chile, South Africa and Australia. Methods Digitized range maps for butterfly species in each region were used to map richness patterns in summer (for all areas) and winter (for USA/Canada, Europe/NW Africa and Australia). A phylogeny resolved to subfamily was used to map the geographic MRD patterns. Regression trees and general linear models examined climatic and vegetation correlates of species richness and MRD within and among regions. Results Various combinations of climate and vegetation were strong predictors of species richness gradients within regions, but unresolved ‘regional’ factors contributed to the multiregional pattern. Regionally based differences in phylogenetic structure also exist, but MRD is negatively correlated with temperature both within and across areas. MRD patterns consistent with tropical niche conservatism occur in most areas. With a possible partial exception of Mexico, faunas in cold climates and in mountains are more derived than faunas in lowlands and tropical/subtropical climates. In USA/Canada, Europe and Australia, winter faunas are more derived than summer faunas. Main conclusions The phylogenetic pattern previously found in the USA and Canada is widespread in both the Northern and Southern Hemispheres, and niche conservatism and the evolution of cold tolerance is the likely explanation for the development of the global butterfly species richness gradient over evolutionary time. Contemporary climate also influences species richness patterns but is unlikely to be a complete explanation globally. The importance of climate is also manifested in the seasonal loss of more basal butterfly elements outside the tropics in winter.     

A model for colour pattern formation in the butterfly wing of Papilio dardanus

The butterfly Papilio dardanus is well known for the spectacular phenotypic polymorphism in the female of the species. We show that numerical simulations of a reaction diffusion model on a geometrically accurate wing domain produce spatial patterns that are consistent with many of those observed on the butterfly. Our results suggest that the wing coloration is due to a simple underlying stripe-like pattern of some pigment-inducing morphogen. We focus on the effect of key factors such as parameter values for mode selection, threshold values which determine colour, wing shape and boundary conditions. The generality of our approach should allow us to investigate other butterfly species. The relationship between these key factors and gene activities is discussed in the context of recent biological advances.     

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast‐Z effect in Heliconius

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex‐ and tissue‐specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male‐biased than female‐biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic‐specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast‐Z effect where hemizygosity influences the rate of adaptive substitutions.     

Rates and patterns of chromosomal evolution in Drosophila pseudoobscura and D. miranda

Comparisons of gene orders between species permit estimation of the rate of chromosomal evolution since their divergence from a common ancestor. We have compared gene orders on three chromosomes of Drosophila pseudoobscura with its close relative, D. miranda, and the distant outgroup species, D. melanogaster, by using the public genome sequences of D. pseudoobscura and D. melanogaster and approximately 50 in situ hybridizations of gene probes in D. miranda. We find no evidence for extensive transfer of genes among chromosomes in D. miranda. The rates of chromosomal rearrangements between D. miranda and D. pseudoobscura are far higher than those found before in Drosophila and approach those for nematodes, the fastest rates among higher eukaryotes. In addition, we find that the D. pseudoobscura chromosome with the highest level of inversion polymorphism (Muller's element C) does not show an unusually fast rate of evolution with respect to chromosome structure, suggesting that this classic case of inversion polymorphism reflects selection rather than mutational processes. On the basis of our results, we propose possible ancestral arrangements for the D. pseudoobscura C chromosome, which are different from those in the current literature. We also describe a new method for correcting for rearrangements that are not detected with a limited set of markers.     

Evidence for frequency‐dependent selection maintaining polymorphism in the Batesian mimic Papilio polytes in multiple islands in the Ryukyus,Japan

Batesian mimicry is a well‐studied adaptation for predation avoidance, in which a mimetic species resembles an unpalatable model species. Batesian mimicry can be under positive selection because of the protection gained against predators, due to resemblance to unpalatable model species. However, in some mimetic species, nonmimetic individuals are present in populations, despite the benefits of mimicry. The mechanism for evolution of such mimetic polymorphism remains an open question. Here, we address the hypothesis that the abundance of mimics is limited by that of the models, leading to mimetic polymorphism. In addition, other forces such as the effects of common ancestry and/or isolation by distance may explain this phenomenon. To investigate this question, we focused on the butterfly, Papilio polytes, that exhibits mimetic polymorphism on multiple islands of the Ryukyus, Japan, and performed field surveys and genetic analysis. We found that the mimic ratio of P. polytes was strongly correlated with the model abundance observed on each of the five islands, suggesting negative frequency‐dependent selection is driving the evolution of polymorphism in P. polytes populations. Molecular phylogenetic analysis indicated that the southern island populations are the major source of genetic diversity, and the middle and northern island populations arose by relatively recent migration. This view was also supported by mismatch distribution and Tajima's D analyses, suggesting a recent population expansion on the middle and northern islands, and stable population persistence on the southern islands. The frequency of the mimetic forms within P. polytes populations is thus explained by variations in the model abundance rather than by population structure. Thus, we propose that predation pressure, rather than neutral forces, have shaped the Batesian mimicry polymorphism in P. polytes observed in the Ryukyus.     

鸟类羽色多态现象：概念及其进化机制假说

多态现象对了解物种的遗传、变异和进化有着重要意义,鸟类羽色多态现象的进化机制,是多态现象研究中非常重要的内容。鸟类羽色多态现象的进化机制假说主要有：分化选择假说、异类选择假说和非随机性交配假说等。本文对鸟类羽色多态现象的概念和以上进化机制假说进行了综述,并针对鸟类羽色多态现象进化机制研究中存在争议的问题及研究的不足提出了自己的看法。     

Natural selection in mimicry

Biological mimicry has served as a salient example of natural selection for over a century, providing us with a dazzling array of very different examples across many unrelated taxa. We provide a conceptual framework that brings together apparently disparate examples of mimicry in a single model for the purpose of comparing how natural selection affects models, mimics and signal receivers across different interactions. We first analyse how model–mimic resemblance likely affects the fitness of models, mimics and receivers across diverse examples. These include classic Batesian and Müllerian butterfly systems, nectarless orchids that mimic Hymenoptera or nectar‐producing plants, caterpillars that mimic inert objects unlikely to be perceived as food, plants that mimic abiotic objects like carrion or dung and aggressive mimicry where predators mimic food items of their own prey. From this, we construct a conceptual framework of the selective forces that form the basis of all mimetic interactions. These interactions between models, mimics and receivers may follow four possible evolutionary pathways in terms of the direction of selection resulting from model–mimic resemblance. Two of these pathways correspond to the selective pressures associated with what is widely regarded as Batesian and Müllerian mimicry. The other two pathways suggest mimetic interactions underpinned by distinct selective pressures that have largely remained unrecognized. Each pathway is characterized by theoretical differences in how model–mimic resemblance influences the direction of selection acting on mimics, models and signal receivers, and the potential for consequent (co)evolutionary relationships between these three protagonists. The final part of this review describes how selective forces generated through model–mimic resemblance can be opposed by the basic ecology of interacting organisms and how those forces may affect the symmetry, strength and likelihood of (co)evolution between the three protagonists within the confines of the four broad evolutionary possibilities. We provide a clear and pragmatic visualization of selection pressures that portrays how different mimicry types may evolve. This conceptual framework provides clarity on how different selective forces acting on mimics, models and receivers are likely to interact and ultimately shape the evolutionary pathways taken by mimetic interactions, as well as the constraints inherent within these interactions.     

The long-term evolution of multilocus traits under frequency-dependent disruptive selection

Frequency-dependent disruptive selection is widely recognized as an important source of genetic variation. Its evolutionary consequences have been extensively studied using phenotypic evolutionary models, based on quantitative genetics, game theory, or adaptive dynamics. However, the genetic assumptions underlying these approaches are highly idealized and, even worse, predict different consequences of frequency-dependent disruptive selection. Population genetic models, by contrast, enable genotypic evolutionary models, but traditionally assume constant fitness values. Only a minority of these models thus addresses frequency-dependent selection, and only a few of these do so in a multilocus context. An inherent limitation of these remaining studies is that they only investigate the short-term maintenance of genetic variation. Consequently, the long-term evolution of multilocus characters under frequency-dependent disruptive selection remains poorly understood. We aim to bridge this gap between phenotypic and genotypic models by studying a multilocus version of Levene's soft-selection model. Individual-based simulations and deterministic approximations based on adaptive dynamics theory provide insights into the underlying evolutionary dynamics. Our analysis uncovers a general pattern of polymorphism formation and collapse, likely to apply to a wide variety of genetic systems: after convergence to a fitness minimum and the subsequent establishment of genetic polymorphism at multiple loci, genetic variation becomes increasingly concentrated on a few loci, until eventually only a single polymorphic locus remains. This evolutionary process combines features observed in quantitative genetics and adaptive dynamics models, and it can be explained as a consequence of changes in the selection regime that are inherent to frequency-dependent disruptive selection. Our findings demonstrate that the potential of frequency-dependent disruptive selection to maintain polygenic variation is considerably smaller than previously expected.     

Optimality models of phage life history and parallels in disease evolution

Optimality models constitute one of the simplest approaches to understanding phenotypic evolution. Yet they have shortcomings that are not easily evaluated in most organisms. Most importantly, the genetic basis of phenotype evolution is almost never understood, and phenotypic selection experiments are rarely possible. Both limitations can be overcome with bacteriophages. However, phages have such elementary life histories that few phenotypes seem appropriate for optimality approaches. Here we develop optimality models of two phage life history traits, lysis time and host range. The lysis time models show that the optimum is less sensitive to differences in host density than suggested by earlier analytical work. Host range evolution is approached from the perspective of whether the virus should avoid particular hosts, and the results match optimal foraging theory: there is an optimal "diet" in which host types are either strictly included or excluded, depending on their infection qualities. Experimental tests of both models are feasible, and phages provide concrete illustrations of many ways that optimality models can guide understanding and explanation. Phage genetic systems already support the perspective that lysis time and host range can evolve readily and evolve without greatly affecting other traits, one of the main tenets of optimality theory. The models can be extended to more general properties of infection, such as the evolution of virulence and tissue tropism.     

Studying Patterns of Recent Evolution at Synonymous Sites and Intronic Sites in Drosophila melanogaster

Most previous studies of the evolution of codon usage bias (CUB) and intronic GC content (iGC) in Drosophila melanogaster were based on between-species comparisons, reflecting long-term evolutionary events. However, a complete picture of the evolution of CUB and iGC cannot be drawn without knowledge of their more recent evolutionary history. Here, we used a polymorphism dataset collected from Zimbabwe to study patterns of the recent evolution of CUB and iGC. Analyzing coding and intronic data jointly with a model which can simultaneously estimate selection, mutational, and demographic parameters, we have found that: (1) natural selection is probably acting on synonymous codons; (2) a constant population size model seems to be sufficient to explain most of the observed synonymous polymorphism patterns; (3) GC is favored over AT in introns. In agreement with the long-term evolutionary patterns, ongoing selection acting on X-linked synonymous codons is stronger than that acting on autosomal codons. The selective differences between preferred and unpreferred codons tend to be greater than the differences between GC and AT in introns, suggesting that natural selection, not just biased gene conversion, may have influenced the evolution of CUB. Interestingly, evidence for non-equilibrium evolution comes exclusively from the intronic data. However, three different models, an equilibrium model with two classes of selected sites and two non-equilibrium models with changes in either population size or mutational parameters, fit the intronic data equally well. These results show that using inadequate selection (or demographic) models can result in incorrect estimates of demographic (or selection) parameters.     

Human hemoglobin: polymorphism, neutrality of variants, evolutionary aspect

Integrated analysis of the polymorphism of human hemoglobin has been made using populational genetics, hematological, physiological, protein chemistry and molecular biology data. The known variants of human hemoglobin are conventionally classified as "widely common", "less common" and "rare", depending on their contribution to polymorphism. The importance of homeostasis and compensatory reactions for maintaining the resistance of the human body against mutant hemoglobins is emphasized. Hb D Punjab and Hb O Arab being relatively neutral, the genetic structure of populations may restrict their distribution. A hypothesis is put forward concerning the possible role of an increased local conformational mobility of protein in creating neutral protein variants. It is proposed to discriminate between truly neutral and pseudoneutral protein variants. In case of possible changes in the genetic and environmental factors, the former are not subject to selection, while the latter may be. Contribution to neutral evolution can be made only by truly neutral variants. In a compensated heterozygotic state the truly neutral and pseudoneutral variants may give rise to new functions and adaptively valuable properties in protein. The evolution of proteins is believed to proceed from a stage which is consistent with M. Kimura's concept of neutrality of protein polymorphism toward a stage which is consistent with the concept of selectionism. It is concluded, that the currently observed degree of polymorphism of human hemoglobin corresponds to the present stage of molecular evolution of the protein.