Spatial processes and evolutionary models: a critical review

Evolution is a fundamentally population level process in which variation, drift and selection produce both temporal and spatial patterns of change. Statistical model fitting is now commonly used to estimate which kind of evolutionary process best explains patterns of change through time using models like Brownian motion, stabilizing selection (Ornstein–Uhlenbeck) and directional selection on traits measured from stratigraphic sequences or on phylogenetic trees. But these models assume that the traits possessed by a species are homogeneous. Spatial processes such as dispersal, gene flow and geographical range changes can produce patterns of trait evolution that do not fit the expectations of standard models, even when evolution at the local‐population level is governed by drift or a typical OU model of selection. The basic properties of population level processes (variation, drift, selection and population size) are reviewed and the relationship between their spatial and temporal dynamics is discussed. Typical evolutionary models used in palaeontology incorporate the temporal component of these dynamics, but not the spatial. Range expansions and contractions introduce rate variability into drift processes, range expansion under a drift model can drive directional change in trait evolution, and spatial selection gradients can create spatial variation in traits that can produce long‐term directional trends and punctuation events depending on the balance between selection strength, gene flow, extirpation probability and model of speciation. Using computational modelling that spatial processes can create evolutionary outcomes that depart from basic population‐level notions from these standard macroevolutionary models.     

EVOLUTION OF VARIATION AND VARIABILITY UNDER FLUCTUATING,STABILIZING, AND DISRUPTIVE SELECTION

How variation and variability (the capacity to vary) may respond to selection remain open questions. Indeed, effects of different selection regimes on variational properties, such as canalization and developmental stability are under debate. We analyzed the patterns of among‐ and within‐individual variation in two wing‐shape characters in populations of Drosophila melanogaster maintained under fluctuating, disruptive, and stabilizing selection for more than 20 generations. Patterns of variation in wing size, which was not a direct target of selection, were also analyzed. Disruptive selection dramatically increased phenotypic variation in the two shape characters, but left phenotypic variation in wing size unaltered. Fluctuating and stabilizing selection consistently decreased phenotypic variation in all traits. In contrast, within‐individual variation, measured by the level of fluctuating asymmetry, increased for all traits under all selection regimes. These results suggest that canalization and developmental stability are evolvable and presumably controlled by different underlying genetic mechanisms, but the evolutionary responses are not consistent with an adaptive response to selection on variation. Selection also affected patterns of directional asymmetry, although inconsistently across traits and treatments.     

The form of sexual selection on male genitalia cannot be inferred from within-population variance and allometry - a case study in Aquarius remigis

Male genital morphology in insects and arachnids is characterized by static hypoallometry and low intrapopulational levels of phenotypic variation relative to other male traits. The one-size-fits-all model of genital evolution attributes these patterns to stabilizing sexual selection. This model relies on the assumption that the observed patterns of variation and allometry reflect the form of sexual selection acting these traits. We test this by examining the patterns of scaling and trait variation for a set of genitalic and somatic morphological traits in male water striders (Aquarius remigis). This suite of traits is of particular interest because previous work has shown that the genitalic traits are under strong directional selection whereas the somatic traits are under either weak directional or stabilizing selection. Because the selection regime for these traits is known, we can, for the first time, test the purported relationship between trait variation, scaling, and the form of sexual selection. We show that the patterns of variation and scaling of these traits differ sharply from those predicted for traits experiencing strong directional sexual selection. Specifically, the male genital structures show static hypoallometry and low intrapopulational levels of phenotypic variation relative to other male traits, in spite of consistent, strong, directional sexual selection. These scaling relationships and levels of variation are typical of genital traits in other insect species, where they have been presumed to reflect stabilizing sexual selection. Our data clearly refute the assumption of the one-size-fits-all hypothesis that hypoallometric scaling of genitalic traits implies stabilizing selection. We discuss the implications of this finding and propose future directions for improving our current understanding of genital evolution in arthropods.     

QUANTITATIVE GENETICS OF BRYOZOAN PHENOTYPIC EVOLUTION. I. RATE TESTS FOR RANDOM CHANGE VERSUS SELECTION IN DIFFERENTIATION OF LIVING SPECIES

The possible roles of random genetic change and natural selection in bryozoan speciation were analyzed using quantitative genetic methods on breeding data for traits of skeletal morphology in two closely related species of the cheilostome Stylopoma. The hypothesis that morphologic differences between the species are caused entirely by mutation and genetic drift could not be rejected for reasonable rates of mutation maintained for as few as 10³ to 10⁴ generations. Divergence times this short or shorter are consistent with the abrupt appearances of many invertebrate species in the fossil record, commonly followed by millions of years of morphologic stasis. To produce these differences over 10³ generations or fewer, directional selection acting alone would require unrealistically high levels of minimum selective mortality throughout divergence. Thus, selection is unnecessary to explain the divergence of these species, except as a means of accelerating the effects of random genetic change on shorter time scales (directional selection), or decelerating them over longer ones (stabilizing selection). These results are consistent with a variety of models of phenotypic evolution involving random shifts between multiple adaptive peaks. Similar results were obtained by substituting trait heritabilities and genetic covariances reconstructed by partitioning within- and among-colony phenotypic variance in place of the values based on breeding data. Quantitative genetic analysis of speciation in fossil bryozoan lineages is thus justified.     

The evolution of alternative developmental pathways: footprints of selection on life-history traits in a butterfly

Developmental pathways may evolve to optimize alternative phenotypes across environments. However, the maintenance of such adaptive plasticity under relaxed selection has received little study. We compare the expression of life-history traits across two developmental pathways in two populations of the butterfly Pararge aegeria where both populations express a diapause pathway but one never expresses direct development in nature. In the population with ongoing selection on both pathways, the difference between pathways in development time and growth rate was larger, whereas the difference in body size was smaller compared with the population experiencing relaxed selection on one pathway. This indicates that relaxed selection on the direct pathway has allowed life-history traits to drift towards values associated with lower fitness when following this pathway. Relaxed selection on direct development was also associated with a higher degree of genetic variation for protandry expressed as within-family sexual dimorphism in growth rate. Genetic correlations for larval growth rate across sexes and pathways were generally positive, with the notable exception of correlation estimates that involved directly developing males of the population that experienced relaxed selection on this pathway. We conclude that relaxed selection on one developmental pathway appears to have partly disrupted the developmental regulation of life-history trait expression. This in turn suggests that ongoing selection may be responsible for maintaining adaptive developmental regulation along alternative developmental pathways in these populations.     

Patterns of contemporary phenotypic selection and flower integration in the hummingbird‐pollinated Nicotiana glauca between populations with different flower–pollinator combinations

We studied six populations of the hummingbird‐pollinated Nicotiana glauca to determine if the marked differences in the degree of floral‐pollinator mismatch between populations promote divergences in the pattern of pollinator‐mediated phenotypic selection on single traits and on the evolution of complexes of many interacting floral traits. We found evidence that flower phenotype is being shaped by pollinator‐mediated phenotypic selection, since corolla length was consistently under contemporary directional or stabilizing selection. Weak directional selection for longer corollas was found in two populations with low flower–pollinator mismatch; much stronger directional selection was detected for shorter corollas in two populations with high flower–pollinator mismatch; finally, the remaining two populations with intermediate flower–pollinator mismatch showed stabilizing selection for corolla length. N. glauca populations differed in every flower character measured but variations in pollinator‐mediated selection among populations were only observed for corolla length. Multiple covariation among traits was favoured, as suggested by the predominately functional patterns of integration and selection of complexes of many interacting floral traits. This was consistent with the patterns of correlational selection exhibited by four of the six populations, where corolla length was under significant selection in combination with corolla width, style length or stamen length. Overall floral integration was relatively high in all populations but phenotypic integration patterns were not clearly accounted by the degree of flower–pollinator mismatch or type of phenotypic selection, suggesting that trait covariation at the entire flower level is not explained by the current scenario of pollinator‐mediated selection.     

VARIATION IN ACOUSTIC SIGNALLING TRAITS EXHIBITS FOOTPRINTS OF SEXUAL SELECTION

Phenotypic variation is ubiquitous in nature and a precondition for adaptive evolution. However, theory predicts that the extent of phenotypic variation should decrease with increasing strength of selection on a trait. Comparative analyses of trait variability have repeatedly used this expectation to infer the type or strength of selection. Yet, the suggested influence of selection on trait variability has rarely been tested empirically. In the present study, I compare estimates of sexual selection strength and trait variability from published data. I constricted the analysis to acoustic courtship traits in amphibians and insects with known variability and corresponding results of female binary choice experiments on these traits. Trait variability and strength of sexual selection were significantly correlated, and both were correlated with signal duration. Because traits under stronger selection had lower variation even after the effect of signal duration was eliminated, I conclude that traces of the strength of selection can be observed with respect to variation of acoustic signaling traits in insects and amphibians. The analysis also shows that traits under stabilizing selection have significantly lower phenotypic variability than traits under directional selection.     

Testing the molecular and evolutionary causes of a 'leapfrog' pattern of geographical variation in coloration

Understanding the mechanisms accounting for the evolution of phenotypic diversity is central to evolutionary biology. We use molecular and phenotypic data to test hypotheses for 'leapfrog' patterns of geographical variation, in which phenotypically similar, disjunct populations are separated by distinct populations of the same species. Phylogenetic reconstructions revealed independent evolution of melanic plumage characters in different populations in the Neotropical avian genus Arremon. Thus, phenotypic similarities between distant populations cannot be explained by close phylogenetic affinity. Nor can they be attributed to recurring mutations in the MC1R gene, a locus involved in melanic pigmentation. A coalescent analysis indicates that plumage traits have become fixed at a faster rate than expected under genetic drift, suggesting that selection underlies their repeated evolution. In contrast to views that genetic drift drives phenotypic differentiation in Neotropical montane birds, our results imply that geographical variation preceding speciation may reflect the action of deterministic selective processes.     

Natural selection on immune defense: A field experiment

Predicting the evolution of phenotypic traits requires an understanding of natural selection on them. Despite its indispensability in the fight against parasites, selection on host immune defense has remained understudied. Theory predicts immune traits to be under stabilizing selection due to associated trade‐offs with other fitness‐related traits. Empirical studies, however, report mainly positive directional selection. This discrepancy could be caused by low phenotypic variation in the examined individuals and/or variation in host resource level that confounds trade‐offs in empirical studies. In a field experiment where we maintained Lymnaea stagnalis snails individually in cages in a lake, we investigated phenotypic selection on two immune defense traits, phenoloxidase (PO)‐like activity and antibacterial activity, in hemolymph. We used a diverse laboratory population and manipulated snail resource level by limiting their food supply. For six weeks, we followed immune activity, growth, and two fitness components, survival and fecundity of snails. We found that PO‐like activity and growth were under stabilizing selection, while antibacterial activity was under positive directional selection. Selection on immune traits was mainly driven by variation in survival. The form of selection on immune defense apparently depends on the particular trait, possibly due to its importance for countering the present parasite community.     

REDUCED GENETIC VARIANCE AMONG HIGH FITNESS INDIVIDUALS: INFERRING STABILIZING SELECTION ON MALE SEXUAL DISPLAYS IN DROSOPHILA SERRATA

Directional selection is prevalent in nature, yet phenotypes tend to remain relatively constant, suggesting a limit to trait evolution. However, the genetic basis of this limit is unresolved. Given widespread pleiotropy, opposing selection on a trait may arise from the effects of the underlying alleles on other traits under selection, generating net stabilizing selection on trait genetic variance. These pleiotropic costs of trait exaggeration may arise through any number of other traits, making them hard to detect in phenotypic analyses. Stabilizing selection can be inferred, however, if genetic variance is greater among low‐ compared to high‐fitness individuals. We extend a recently suggested approach to provide a direct test of a difference in genetic variance for a suite of cuticular hydrocarbons (CHCs) in Drosophila serrata. Despite strong directional sexual selection on these traits, genetic variance differed between high‐ and low‐fitness individuals and was greater among the low‐fitness males for seven of eight CHCs, significantly more than expected by chance. Univariate tests of a difference in genetic variance were nonsignificant but likely have low power. Our results suggest that further CHC exaggeration in D. serrata in response to sexual selection is limited by pleiotropic costs mediated through other traits.     

The evolution of genetic architecture under frequency-dependent disruptive selection

We propose a model to analyze a quantitative trait under frequency-dependent disruptive selection. Selection on the trait is a combination of stabilizing selection and intraspecific competition, where competition is maximal between individuals with equal phenotypes. In addition, there is a density-dependent component induced by population regulation. The trait is determined additively by a number of biallelic loci, which can have different effects on the trait value. In contrast to most previous models, we assume that the allelic effects at the loci can evolve due to epistatic interactions with the genetic background. Using a modifier approach, we derive analytical results under the assumption of weak selection and constant population size, and we investigate the full model by numerical simulations. We find that frequency-dependent disruptive selection favors the evolution of a highly asymmetric genetic architecture, where most of the genetic variation is concentrated on a small number of loci. We show that the evolution of genetic architecture can be understood in terms of the ecological niches created by competition. The phenotypic distribution of a population with an adapted genetic architecture closely matches this niche structure. Thus, evolution of the genetic architecture seems to be a plausible way for populations to adapt to regimes of frequency-dependent disruptive selection. As such, it should be seen as a potential evolutionary pathway to discrete polymorphisms and as a potential alternative to other evolutionary responses, such as the evolution of sexual dimorphism or assortative mating.     

Sexual selection and genetic colour polymorphisms in animals

Genetic colour polymorphisms are widespread across animals and often subjected to complex selection regimes. Traditionally, colour morphs were used as simple visual markers to measure allele frequency changes in nature, selection, population divergence and speciation. With advances in sequencing technology and analysis methods, several model systems are emerging where the molecular targets of selection are being described. Here, we discuss recent studies on the genetics of sexually selected colour polymorphisms, aiming at (i) reviewing the evidence of sexual selection on colour polymorphisms, (ii) highlighting the genetic architecture, molecular and developmental basis underlying phenotypic colour diversification and (iii) discuss how the maintenance of such polymorphisms might be facilitated or constrained by these. Studies of the genetic architecture of colour polymorphism point towards the importance of tight clustering of colour loci with other trait loci, such as in the case of inversions and supergene structures. Other interesting findings include linkage between colour loci and mate preferences or sex determination, and the role of introgression and regulatory variation in fuelling polymorphisms. We highlight that more studies are needed that explicitly integrate fitness consequences of sexual selection on colour with the underlying molecular targets of colour to gain insights into the evolutionary consequences of sexual selection on polymorphism maintenance.     

MORPHOLOGICAL EVOLUTION OF THE DROSOPHILA VIRILIS SPECIES GROUP AS ASSESSED BY RATE TESTS FOR NATURAL SELECTION ON QUANTITATIVE CHARACTERS

Two rate tests for assessing natural selection on quantitative traits are discussed for their usefulness in macroevolutionary and adaptational studies. The underlying assumptions and parameter estimation for the constant-heritability (CH) and mutation-drift-equilibrium (MDE) models, which are the bases for these tests, are discussed. The purpose of these rate tests is to determine whether morphological change has occurred either too fast to be explained by neutral drift, which suggests directional selection, or too slow, which suggests stabilizing selection. Previous formulations of these rate tests have ignored the phylogenetic component. Several models of evolution are considered to help account for phylogeny in the context of rate tests. The MDE rate test for stabilizing selection was performed on nine morphological characters among several species of the Drosophila virilis species group. These tests can be interpreted to suggest that stabilizing selection has probably been a major factor in producing the observed similarity among the Drosophila species examined.     

A population genetics model for multiple quantitative traits exhibiting pleiotropy and epistasis

We study a population genetics model of an organism with a genome of L(tot)loci that determine the values of T quantitative traits. Each trait is controlled by a subset of L loci assigned randomly from the genome. There is an optimum value for each trait, and stabilizing selection acts on the phenotype as a whole to maintain actual trait values close to their optima. The model contains pleiotropic effects (loci can affect more than one trait) and epistasis in fitness. We use adaptive walk simulations to find high-fitness genotypes and to study the way these genotypes are distributed in sequence space. We then simulate the evolution of haploid and diploid populations on these fitness landscapes and show that the genotypes of populations are able to drift through sequence space despite stabilizing selection on the phenotype. We study the way the rate of drift and the extent of the accessible region of sequence space is affected by mutation rate, selection strength, population size, recombination rate, and the parameters L and T that control the landscape shape. There are three regimes of the model. If LTL(tot), there are many small peaks that can be spread over a wide region of sequence space. Compensatory neutral mutations are important in the population dynamics in this case.     

Speciation despite gene flow when developmental pathways evolve

Abstract.— Evolutionary biologists assume that species formation requires a drastic reduction in gene exchange between populations, but the rate sufficient to prevent speciation is unknown. To study speciation, we use a new class of population genetic models that incorporate simple developmental genetic rules, likely present in all organisms, to construct the phenotype. When we allow replicate populations to evolve in parallel to a new, shared optimal phenotype, often their hybrids acquire poorly regulated phenotypes: Dobzhansky-Muller incompatibilities arise and postzygotic reproductive isolation evolves. Here we show that, although gene exchange does inhibit this process, it is the proportion of migrants exchanged ( m ) rather than the number of migrants ( Nm ) that is critical, and rates as high as 16 individuals exchanged per generation still permit the evolution of postzygotic isolation. Stronger directional selection counters the inhibitory effect of gene flow, increasing the speciation probability. We see similar results when populations in a standard two-locus, two-allele Dobzhansky-Muller model are subject to simultaneous directional selection and gene flow. However, in developmental pathway models with more than two loci, gene flow is more able to impede speciation. Genetic incompatibilities arise as frequent by-products of adaptive evolution of traits determined by regulatory pathways, something that does not occur when phenotypes are modeled using the standard, additive genetic framework. Development therefore not only constrains the microevolutionary process, it also facilitates the interactions among genes and gene products that make speciation more likely–even in the face of strong gene flow.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

Adaptive evolution of reproductive and vegetative traits driven by breeding systems

The evolution of inflorescence size, a key trait in reproductive success, was studied in the genus Acer under a perspective of adaptive evolution. Breeding systems, hypothesized to indicate different levels of mating competition, were considered as the selective scenarios defining different optima of inflorescence size. Larger inflorescences, which increase male fitness by generating larger floral displays, were hypothesized to be selected under scenarios with higher competition with unisexuals. An identical approach was used to test if the same selective regimes could be driving the evolution of leaf size, a vegetative trait that was found to be correlated with inflorescence size. A Brownian motion model of inflorescence/leaf-size evolution (which cannot distinguish between changes caused by pure drift processes and changes caused by natural selection in rapidly and randomly changing environments) was compared with several adaptive Ornstein-Uhlenbeck (OU) models, which can quantify the effects of both stochasticity and natural selection. The best-fitting model for inflorescence/leaf-size evolution was an OU model with three optima that increased with the level of mating competition. Both traits evolved under the same selective regimes and in the same direction, confirming a pattern of correlated evolution. These results show that a selective regime hypothetically related to the evolution of a reproductive trait can also explain the evolution of a vegetative trait.     

Bimodal trait distributions with large variances question the reliability of trait-based aggregate models

Functionally diverse communities can adjust their species composition to altered environmental conditions, which may influence food web dynamics. Trait-based aggregate models cope with this complexity by ignoring details about species identities and focusing on their functional characteristics (traits). They describe the temporal changes of the aggregate properties of entire communities, including their total biomasses, mean trait values, and trait variances. The applicability of aggregate models depends on the validity of their underlying assumptions that trait distributions are normal and exhibit small variances. We investigated to what extent this can be expected to work by comparing an innovative model that accounts for the full trait distributions of predator and prey communities to a corresponding aggregate model. We used a food web structure with well-established trade-offs among traits promoting mutual adjustments between prey edibility and predator selectivity in response to selection. We altered the shape of the trade-offs to compare the outcome of the two models under different selection regimes, leading to trait distributions increasingly deviating from normality. Their biomass and trait dynamics agreed very well for stabilizing selection and reasonably well for directional selection, under which different trait values are favored at different times. However, for disruptive selection, the results of the aggregate model strongly deviated from the full trait distribution model that showed bimodal trait distributions with large variances. Hence, the outcome of aggregate models is reliable under ideal conditions but has to be questioned when confronted with more complex selection regimes and trait distributions, which are commonly observed in nature.     

Genomic Adaptations to Salinity Resist Gene Flow in the Evolution of Floridian Watersnakes

The migration-selection balance often governs the evolution of lineages, and speciation with gene flow is now considered common across the tree of life. Ecological speciation is a process that can facilitate divergence despite gene flow due to strong selective pressures caused by ecological differences; however, the exact traits under selection are often unknown. The transition from freshwater to saltwater habitats provides strong selection targeting traits with osmoregulatory function. Several lineages of North American watersnakes (Nerodia spp.) are known to occur in saltwater habitat and represent a useful system for studying speciation by providing an opportunity to investigate gene flow and evaluate how species boundaries are maintained or degraded. We use double digest restriction-site associated DNA sequencing to characterize the migration-selection balance and test for evidence of ecological divergence within the Nerodia fasciata-clarkii complex in Florida. We find evidence of high intraspecific gene flow with a pattern of isolation-by-distance underlying subspecific lineages. However, we identify genetic structure indicative of reduced gene flow between inland and coastal lineages suggesting divergence due to isolation-by-environment. This pattern is consistent with observed environmental differences where the amount of admixture decreases with increased salinity. Furthermore, we identify significantly enriched terms related to osmoregulatory function among a set of candidate loci, including several genes that have been previously implicated in adaptation to salinity stress. Collectively, our results demonstrate that ecological differences, likely driven by salinity, cause strong divergent selection which promotes divergence in the N. fasciata-clarkii complex despite significant gene flow.     

LARGE-SCALE CANDIDATE GENE SCAN REVEALS THE ROLE OF CHEMORECEPTOR GENES IN HOST PLANT SPECIALIZATION AND SPECIATION IN THE PEA APHID

Understanding the drivers of speciation is critical to interpreting patterns of biodiversity. The identification of the genetic changes underlying adaptation and reproductive isolation is necessary to link barriers to gene flow to the causal origins of divergence. Here, we present a novel approach to the genetics of speciation, which should complement the commonly used approaches of quantitative trait locus mapping and genome-wide scans for selection. We present a large-scale candidate gene approach by means of sequence capture, applied to identifying the genetic changes underlying reproductive isolation in the pea aphid, a model system for the study of ecological speciation. Targeted resequencing enabled us to scale up the candidate gene approach, specifically testing for the role of chemosensory gene families in host plant specialization. Screening for the signature of divergence under selection at 172 candidate and noncandidate loci, we revealed a handful of loci that show high levels of differentiation among host races, which almost all correspond to odorant and gustatory receptor genes. This study offers the first indication that some chemoreceptor genes, often tightly linked together in the genome, could play a key role in local adaptation and reproductive isolation in the pea aphid and potentially other phytophagous insects. Our approach opens a new route toward the functional genomics of ecological speciation.