Evolutionary history of a dispersal‐associated locus across sympatric and allopatric divergent populations of a wing‐polymorphic beetle across Atlantic Europe

Studying the evolutionary history of trait divergence, in particular those related to dispersal capacity, is of major interest for the process of local adaptation and metapopulation dynamics. Here, we reconstruct the evolution of different alleles at the nuclear‐encoded mitochondrial NADP⁺‐dependent isocitrate dehydrogenase (mtIdh) locus of the ground beetle Pogonus chalceus that are differentially and repeatedly selected in short‐ and long‐winged populations in response to different hydrological regimes at both allopatric and sympatric scales along the Atlantic European coasts. We sequenced 2788 bp of the mtIdh locus spanning a ~7‐kb genome region and compared its variation with that of two supposedly neutral genes. mtIdh sequences show (i) monophyletic clustering of the short‐winged associated mtIDH‐DE haplotypes within the long‐winged associated mtIDH‐AB haplotypes, (ii) a more than tenfold lower haplotype diversity associated with the mtIDH‐DE alleles compared to the mtIDH‐AB alleles and (iii) a high number of fixed nucleotide differences between both mtIDH haplotype clusters. Coalescent simulations suggest that this observed sequence variation in the mtIdh locus is most consistent with a singular origin in a partially isolated subpopulation, followed by a relatively recent spread of the mtIDH‐DE allele in short‐winged populations along the Atlantic coast. These results demonstrate that even traits associated with decreased dispersal capacity can rapidly spread and that reuse of adaptive alleles plays an important role in the adaptive potential within this sympatric mosaic of P. chalceus populations.     

Parallel adaptation to climate above the 35th parallel

Independent or parallel evolution of similar traits is key to understanding the genetics and limitations of adaptation. Adaptation from the same genetic changes in different populations defines parallel evolution. Such genetic changes can derive from standing ancestral variation or de novo mutations and excludes instances of adaptive introgression. In this issue of Molecular Ecology, Walden et al.(2020) investigate the scale of parallel climate adaptation from standing genetic variation between two North American Arabidopsis lyrata lineages, each formed by a distinct evolutionary history during the last glacial cycle. By identifying adaptive variants correlated with three ecologically significant climatic gradients, they show that instead of the same genetic variants or even genes, parallel evolution is only observed at the level of biological processes. The evolution of independent adaptive variants to climate in two genetically close lineages is explained by their different post‐glacial demographic histories. Separate glacial refugia and strong population bottlenecks were probably sufficient to change the landscape of shared allele frequencies, hindering the possibility of parallel evolution.     

PATHOGEN LIFE‐HISTORY TRADE‐OFFS REVEALED IN ALLOPATRY

Trade‐offs in life‐history traits is a central tenet in evolutionary biology, yet their ubiquity and relevance to realized fitness in natural populations remains questioned. Trade‐offs in pathogens are of particular interest because they may constrain the evolution and epidemiology of diseases. Here, we studied life‐history traits determining transmission in the obligate fungal pathogen, Podosphaera plantaginis, infecting Plantago lanceolata. We find that although traits are positively associated on sympatric host genotypes, on allopatric host genotypes relationships between infectivity and subsequent transmission traits change shape, becoming even negative. The epidemiological prediction of this change in life‐history relationships in allopatry is lower disease prevalence in newly established pathogen populations. An analysis of the natural pathogen metapopulation confirms that disease prevalence is lower in newly established pathogen populations and they are more prone to go extinct during winter than older pathogen populations. Hence, life‐history trade‐offs mediated by pathogen local adaptation may influence epidemiological dynamics at both population and metapopulation levels.     

Calls,colours, shape,and genes: a multi‐trait approach to the study of geographic variation in the Amazonian frog Allobates femoralis

Evolutionary divergence in behavioural traits related to mating may represent the initial stage of speciation. Direct selective forces are usually invoked to explain divergence in mate‐recognition traits, often neglecting a role for neutral processes or concomitant differentiation in ecological traits. We adopted a multi‐trait approach to obtain a deeper understanding of the mechanisms behind allopatric divergence in the Amazonian frog, Allobates femoralis. We tested the null hypothesis that geographic distance between populations correlates with genetic and phenotypic divergence, and compared divergence between mate‐recognition (acoustic) and ecological (coloration, body‐shape) traits. We quantified geographic variation in 39 phenotypic traits and a mitochondrial DNA marker among 125 individuals representing eight populations. Geographic variation in acoustic traits was pronounced and tracked the spatial genetic variation, which appeared to be neutral. Thus, the evolution of acoustic traits tracked the shared history of the populations, which is unexpected for pan‐Amazonian taxa or for mate‐recognition traits. Divergence in coloration appeared uncorrelated with genetic distance, and might be partly attributed to local selective pressures, and perhaps to Batesian mimicry. Divergence in body‐shape traits was low. The results obtained depict a complex evolutionary scenario and emphasize the importance of considering multiple traits when disentangling the forces behind allopatric divergence. ©2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 98 , 826–838.     

How phenotypic convergence arises in experimental evolution

Evolutionary convergence is a core issue in the study of adaptive evolution, as well as a highly debated topic at present. Few studies have analyzed this issue using a “real‐time” or evolutionary trajectory approach. Do populations that are initially differentiated converge to a similar adaptive state when experiencing a common novel environment? Drosophila subobscura populations founded from different locations and years showed initial differences and variation in evolutionary rates in several traits during short‐term (～20 generations) laboratory adaptation. Here, we extend that analysis to 40 more generations to analyze (1) how differences in evolutionary dynamics among populations change between shorter and longer time spans, and (2) whether evolutionary convergence occurs after 60 generations of evolution in a common environment. We found substantial variation in longer term evolutionary trajectories and differences between short‐ and longer term evolutionary dynamics. Although we observed pervasive patterns of convergence toward the character values of long‐established populations, populations still remain differentiated for several traits at the final generations analyzed. This pattern might involve transient divergence, as we report in some cases, indicating that more generations should lead to final convergence. These findings highlight the importance of longer term studies for understanding convergent evolution.     

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ～ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.     

A highly pleiotropic amino acid polymorphism in the Drosophila insulin receptor contributes to life‐history adaptation

Finding the specific nucleotides that underlie adaptive variation is a major goal in evolutionary biology, but polygenic traits pose a challenge because the complex genotype–phenotype relationship can obscure the effects of individual alleles. However, natural selection working in large wild populations can shift allele frequencies and indicate functional regions of the genome. Previously, we showed that the two most common alleles of a complex amino acid insertion–deletion polymorphism in the Drosophila insulin receptor show independent, parallel clines in frequency across the North American and Australian continents. Here, we report that the cline is stable over at least a five‐year period and that the polymorphism also demonstrates temporal shifts in allele frequency concurrent with seasonal change. We tested the alleles for effects on levels of insulin signaling, fecundity, development time, body size, stress tolerance, and life span. We find that the alleles are associated with predictable differences in these traits, consistent with patterns of Drosophila life‐history variation across geography that likely reflect adaptation to the heterogeneous climatic environment. These results implicate insulin signaling as a major mediator of life‐history adaptation in Drosophila, and suggest that life‐history trade‐offs can be explained by extensive pleiotropy at a single locus.     

Similar patterns of local barn owl adaptation in the Middle East and Europe with respect to melanic coloration

The maintenance of phenotypic variation is a central question in evolutionary biology. A commonly suggested mechanism is that of local adaptation, whereby different phenotypes are adapted to alternative environmental conditions. A recent study in the European barn owl (Tyto alba) has shown that natural selection maintains a strong clinal variation in reddish pheomelanin‐based coloration. Studies in the region where phenotypic variation in this owl is the highest in Europe have further demonstrated that dark‐reddish and pale‐reddish owls exploit open and wooded habitats, predate voles and wood mice, and are long‐tailed and short‐tailed, respectively. However, it remains unclear as to whether these traits evolved as a consequence of allopatric evolution of dark colour in northern Europe and white colour in southern Europe, during which owls could have also evolved different morphologies and foraging behaviour. This scenario implies that covariation between coloration and foraging behaviour could be a specificity of the European continent, which is not found in other worldwide‐distributed populations. To investigate this issue, we studied a barn owl population in the Middle East. The results obtained show that, as in Central Europe, dark‐reddish female owls breed more often in the open landscape than their pale‐reddish female conspecifics, their offspring are fed with more voles than Muridae, and they are longer‐winged and longer‐tailed. These findings indicate that, in the barn owl, the association in females between pheomelanin‐based coloration and foraging behaviour and morphology is not restricted to the European continent but may well evolve in sympatry in many barn owl populations worldwide. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 447–454.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

MULTITRAIT,HOST‐ASSOCIATED DIVERGENCE AMONG SETS OF BUTTERFLY POPULATIONS: IMPLICATIONS FOR REPRODUCTIVE ISOLATION AND ECOLOGICAL SPECIATION

When populations use different resources, they tend to diverge in traits that affect performance on those resources. The extent and complexity of divergence that is achieved will depend on gene flow, genetic constraints, and the character of divergent selection. We describe divergent host adaptation among Californian populations of the Melitaeine butterfly, Euphydryas editha. Divergence in seven traits created parallel phenotypic suites, each suite associated with the use of a different host species, either Collinsia torreyi or Pedicularis semibarbata. The suites involved alighting responses of adults (probably to visual stimuli), chemosensory responses to leaf surfaces, vertical positioning of adults and larvae (probably due to geotaxis), partitioning of reproductive effort among clutches, and larval performance. Remarkably, the divergent suites did not occur sympatrically, despite ubiquitous co‐occurrence of the hosts, and we know of only one site where any Collinsia species is used sympatrically with another host. In contrast, E. editha often uses two host genera sympatrically when neither of them is Collinsia. We suggest that adaptation to Collinsia is incompatible with adaptation to other hosts and may generate extrinsic postzygotic reproductive isolation among populations. Despite the apparent rarity of host‐shift‐associated speciation in Melitaeine butterflies, adoption of Collinsia as a host may lead to allopatric ecological speciation.     

Do different rates of gene flow underlie variation in phenotypic and phenological clines in a montane grasshopper community?

Species responses to environmental change are likely to depend on existing genetic and phenotypic variation, as well as evolutionary potential. A key challenge is to determine whether gene flow might facilitate or impede genomic divergence among populations responding to environmental change, and if emergent phenotypic variation is dependent on gene flow rates. A general expectation is that patterns of genetic differentiation in a set of codistributed species reflect differences in dispersal ability. In less dispersive species, we predict greater genetic divergence and reduced gene flow. This could lead to covariation in life‐history traits due to local adaptation, although plasticity or drift could mirror these patterns. We compare genome‐wide patterns of genetic structure in four phenotypically variable grasshopper species along a steep elevation gradient near Boulder, Colorado, and test the hypothesis that genomic differentiation is greater in short‐winged grasshopper species, and statistically associated with variation in growth, reproductive, and physiological traits along this gradient. In addition, we estimate rates of gene flow under competing demographic models, as well as potential gene flow through surveys of phenological overlap among populations within a species. All species exhibit genetic structure along the elevation gradient and limited gene flow. The most pronounced genetic divergence appears in short‐winged (less dispersive) species, which also exhibit less phenological overlap among populations. A high‐elevation population of the most widespread species, Melanoplus sanguinipes, appears to be a sink population derived from low elevation populations. While dispersal ability has a clear connection to the genetic structure in different species, genetic distance does not predict growth, reproductive, or physiological trait variation in any species, requiring further investigation to clearly link phenotypic divergence to local adaptation.     

Local adaptation of a holoparasitic plant,Cuscuta europaea: variation among populations

Locally adapted parasites have higher infectivity and/or fitness on sympatric than on allopatric hosts. We tested local adaptation of a holoparasitic plant, Cuscuta europaea, to its host plant, Urtica dioica. We infected hosts from five sites with holoparasites from the same five sites and measured local adaptation in terms of infectivity and parasite performance (biomass) in a reciprocal cross‐infection experiment. The virulence of the parasite did not differ between sympatric and allopatric hosts. Overall, parasites had higher infectivity on sympatric hosts but infectivity and parasite performance varied among populations. Parasites from one of the populations showed local adaptation in terms of performance, whereas parasites from one of the populations had higher infectivity on allopatric hosts compared with sympatric hosts. This among‐population variation may be explained by random variation in parasite adaptation to host populations or by time‐lagged co‐evolutionary oscillations that lead to fluctuations in the level of local adaptation.     

Divergence and gene flow in the globally distributed blue‐winged ducks

The ability to disperse over long distances can result in a high propensity for colonizing new geographic regions, including uninhabited continents, and lead to lineage diversification via allopatric speciation. However, high vagility can also result in gene flow between otherwise allopatric populations, and in some cases, parapatric or divergence‐with‐gene‐flow models might be more applicable to widely distributed lineages. Here, we use five nuclear introns and the mitochondrial control region along with Bayesian models of isolation with migration to examine divergence, gene flow, and phylogenetic relationships within a cosmopolitan lineage comprising six species, the blue‐winged ducks (genus Anas), which inhabit all continents except Antarctica. We found two primary sub‐lineages, the globally‐distributed shoveler group and the New World blue‐winged/cinnamon teal group. The blue‐winged/cinnamon sub‐lineage is composed of sister taxa from North America and South America, and taxa with parapatric distributions are characterized by low to moderate levels of gene flow. In contrast, our data support strict allopatry for most comparisons within the shovelers. However, we found evidence of gene flow from the migratory, Holarctic northern shoveler A. clypeata and the more sedentary, African Cape shoveler A. smithii into the Australasian shoveler A. rhynchotis, although we could not reject strict allopatry. Given the diverse mechanisms of speciation within this complex, the shovelers and blue‐winged/cinnamon teals can serve as an effective model system for examining how the genome diverges under different evolutionary processes and how genetic variation is partitioned among highly dispersive taxa.     

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.     

Heritability,covariation and natural selection on 24 traits of common evening primrose (Oenothera biennis) from a field experiment

This study explored genetic variation and co‐variation in multiple functional plant traits. Our goal was to characterize selection, heritabilities and genetic correlations among different types of traits to gain insight into the evolutionary ecology of plant populations and their interactions with insect herbivores. In a field experiment, we detected significant heritable variation for each of 24 traits of Oenothera biennis and extensive genetic covariance among traits. Traits with diverse functions formed several distinct groups that exhibited positive genetic covariation with each other. Genetic variation in life‐history traits and secondary chemistry together explained a large proportion of variation in herbivory (r² = 0.73). At the same time, selection acted on lifetime biomass, life‐history traits and two secondary compounds of O. biennis, explaining over 95% of the variation in relative fitness among genotypes. The combination of genetic covariances and directional selection acting on multiple traits suggests that adaptive evolution of particular traits is constrained, and that correlated evolution of groups of traits will occur, which is expected to drive the evolution of increased herbivore susceptibility. As a whole, our study indicates that an examination of genetic variation and covariation among many different types of traits can provide greater insight into the evolutionary ecology of plant populations and plant–herbivore interactions.     

Local divergence of thermal reaction norms among amphibian populations is affected by pond temperature variation

Although temperature variation is known to cause large‐scale adaptive divergence, its potential role as a selective factor over microgeographic scales is less well‐understood. Here, we investigated how variation in breeding pond temperature affects divergence in multiple physiological (thermal performance curve and critical thermal maximum [CT_max]) and life‐history (thermal developmental reaction norms) traits in a network of Rana arvalis populations. The results supported adaptive responses to face two main constraints limiting the evolution of thermal adaptation. First, we found support for the faster–slower model, indicating an adaptive response to compensate for the thermodynamic constraint of low temperatures in colder environments. Second, we found evidence for the generalist–specialist trade‐off with populations from colder and less thermally variable environments exhibiting a specialist phenotype performing at higher rates but over a narrower range of temperatures. By contrast, the local optimal temperature for locomotor performance and CT_max did not match either mean or maximum pond temperatures. These results highlight the complexity of the adaptive multiple‐trait thermal responses in natural populations, and the role of local thermal variation as a selective force driving diversity in life‐history and physiological traits in the presence of gene flow.     

Rapid evolution of thermal plasticity in mountain lake Daphnia populations

Populations at risk of extinction due to climate change may be rescued by adaptive evolution or plasticity. Selective agents, such as introduced predators, may enhance or constrain plastic or adaptive responses to temperature. We tested responses of Daphnia to temperature by collecting populations from lakes across an elevational gradient in the presence and absence of fish predators (long‐term selection). We subsequently grew these populations at two elevations in field mesocosms over two years (short‐term selection), followed by a common‐garden experiment at two temperatures in the lab to measure life‐history traits. Both long‐term and short‐term selection affected traits, suggesting that genetic variation of plasticity within populations enabled individuals to rapidly evolve plasticity in response to high temperature. We found that short‐term selection by high temperature increased plasticity for growth rate in all populations. Fecundity was higher in populations from fishless lakes and body size showed greater plasticity in populations from warm lakes (long‐term selection). Neither body size nor fecundity were affected by short‐term thermal selection. These results demonstrate that plasticity is an important component of the life‐history response of Daphnia, and that genetic variation within populations enabled rapid evolution of plasticity in response to selection by temperature.     

Ecological speciation along an elevational gradient in a tropical passerine bird?

Local adaptation of populations along elevational gradients is well known, but conclusive evidence that such divergence has resulted in the origin of distinct species in parapatry remains lacking. We integrated morphological, vocal, genetic and behavioural data to test predictions pertaining to the hypothesis of parapatric ecological speciation associated with elevation in populations of a tropical montane songbird, the Grey‐breasted Wood‐wren (Henicorhina leucophrys: Troglodytidae), from the Sierra Nevada de Santa Marta, Colombia. We confirmed that two distinct populations exist along the elevational gradient. Phylogenetic analyses tentatively indicate that the two populations are not sister taxa, suggesting they did not differentiate from a single ancestor along the gradient, but rather resulted from separate colonization events. The populations showed marked divergence in morphometrics, vocalizations and genetic variation in mitochondrial and nuclear loci, and little to no evidence of hybridization. Individuals of both populations responded more strongly to their own local songs than to songs from another elevation. Although the two forms do not appear to have differentiated locally in parapatry, morphological and vocal divergence along the elevational gradient is consistent with adaptation, suggesting a possible link between adaptive evolution in morphology and songs and the origin of reproductive isolation via a behavioural barrier to gene flow. The adaptive value of phenotypic differences between populations requires additional study.     

What paths do advantageous alleles take during short‐term evolutionary change?

Combining experimental evolution with whole‐genome resequencing is a promising new strategy for investigating the dynamics of evolutionary change. Published studies that have resequenced laboratory‐selected populations of sexual organisms have typically focused on populations sampled at the end of an evolution experiment. These studies have attempted to associate particular alleles with phenotypic change and attempted to distinguish between different theoretical models of adaptation. However, neither the population used to initiate the experiment nor multiple time points sampled during the evolutionary trajectory are generally available for examination. In this issue of Molecular Ecology, Orozco‐terWengel et al. (2012) take a significant step forward by estimating genome‐wide allele frequencies at the start, 15 generations into and at the end of a 37‐generation Drosophila experimental evolution study. The authors identify regions of the genome that have responded to laboratory selection and describe the temporal dynamics of allele frequency change. They identify two common trajectories for putatively adaptive alleles: alleles either gradually increase in frequency throughout the entire 37 generations or alleles plateau at a new frequency by generation 15. The identification of complex trajectories of alleles under selection contributes to a growing body of literature suggesting that simple models of adaptation, whereby beneficial alleles arise and increase in frequency unimpeded until they become fixed, may not adequately describe short‐term response to selection.     

Evolution of mate choice and the so‐called magic traits in ecological speciation

Non‐random mating provides multiple evolutionary benefits and can result in speciation. Biological organisms are characterised by a myriad of different traits, many of which can serve as mating cues. We consider multiple mechanisms of non‐random mating simultaneously within a unified modelling framework in an attempt to understand better which are more likely to evolve in natural populations going through the process of local adaptation and ecological speciation. We show that certain traits that are under direct natural selection are more likely to be co‐opted as mating cues, leading to the appearance of magic traits (i.e. phenotypic traits involved in both local adaptation and mating decisions). Multiple mechanisms of non‐random mating can interact so that trait co‐evolution enables the evolution of non‐random mating mechanisms that would not evolve alone. The presence of magic traits may suggest that ecological selection was acting during the origin of new species.