Convergence in sympatric swallowtail butterflies reveals ecological interactions as a key driver of worldwide trait diversification.

Ecological interactions can promote phenotypic diversification in sympatric species. While competition can enhance trait divergence, other ecological interactions may promote convergence in sympatric species. Within butterflies, evolutionary convergences in wing color patterns have been reported between distantly related species, especially in females of palatable species, where mimetic color patterns are promoted by predator communities shared with defended species living in sympatry. Wing color patterns are also often involved in species recognition in butterflies, and divergence in this trait has been reported in closely related species living in sympatry as a result of reproductive character displacement. Here, we investigate the effect of sympatry between species on the convergence vs. divergence of their wing color patterns in relation to phylogenetic distance, focusing on the iconic swallowtail butterflies (family Papilionidae). We developed an unsupervised machine learning-based method to estimate phenotypic distances among wing color patterns of 337 species, enabling us to finely quantify morphological diversity at the global scale among species and allowing us to compute pairwise phenotypic distances between sympatric and allopatric species pairs. We found phenotypic convergence in sympatry, stronger among distantly related species, while divergence was weaker and restricted to closely related males. The convergence was stronger among females than males, suggesting that differential selective pressures acting on the two sexes drove sexual dimorphism. Our results highlight the significant effect of ecological interactions driven by predation pressures on trait diversification in Papilionidae and provide evidence for the interaction between phylogenetic proximity and ecological interactions in sympatry, acting on macroevolutionary patterns of phenotypic diversification.

The Limits of Evolutionary Convergence in Sympatry: Reproductive Interference and Historical Constraints Leading to Local Diversity in Warning Traits.

AbstractMutualistic interactions between defended species represent a striking case of evolutionary convergence in sympatry, driven by the increased protection against predators brought by mimicry in warning traits. However, such convergence is often limited: sympatric defended species frequently display different or imperfectly similar warning traits. The phylogenetic distance between sympatric species may indeed prevent evolution toward the exact same signal. Moreover, warning traits are also involved in mate recognition, so trait convergence might result in heterospecific courtship and mating. Here, we develop a mathematical model to investigate the strength and direction of the evolution of warning traits in defended species with different ancestral traits. Specifically, we determine the effect of phenotypic distances between ancestral trait states of sympatric defended species and of the costs of heterospecific sexual interactions on imperfect mimicry and trait divergence. Our analytical results confirm that reproductive interference and historical constraints limit the convergence of warning traits, leading to either complete divergence or imperfect mimicry. Our model reveals that imperfect mimicry evolves only when ancestral trait values differ between species because of historical constraints and highlights the importance of female and predator discrimination in the evolution of such imperfect mimicry. Our study thus provides new predictions on how reproductive interference interacts with historical constraints and may promote the emergence of novel warning traits, enhancing mimetic diversity.

Dominance mechanisms in supergene alleles controlling butterfly wing pattern variation: insights from gene expression in Heliconius numata.

Loci under balancing selection, where multiple alleles are maintained, offer a relevant opportunity to investigate the role of natural selection in shaping genetic dominance: the high frequency of heterozygotes at these loci has been shown to enable the evolution of dominance among alleles. In the butterfly Heliconius numata, mimetic wing color variations are controlled by an inversion polymorphism of a circa 2 Mb genomic region (supergene P), with strong dominance between sympatric alleles. To test how differences in dominance observed on wing patterns correlate with variations in expression levels throughout the supergene region, we sequenced the complete transcriptome of heterozygotes at the prepupal stage and compared it to corresponding homozygotes. By defining dominance based on non-overlapping ranges of transcript expression between genotypes, we found contrasting patterns of dominance between the supergene and the rest of the genome; the patterns of transcript expression in the heterozygotes were more similar to the expression observed in the dominant homozygotes in the supergene region. Dominance also differed among the three subinversions of the supergene, suggesting possible epistatic interactions among their gene contents underlying dominance evolution. We found the expression pattern of the melanization gene cortex located in the P-region to predict wing pattern phenotype in the heterozygote. We also identify new candidate genes that are potentially involved in mimetic color pattern variations highlighting the relevance of transcriptomic analyses in heterozygotes to pinpoint candidate genes in non-recombining regions.

Evidence of attack deflection suggests adaptive evolution of wing tails in butterflies.

Predation is a powerful selective force shaping many behavioural and morphological traits in prey species. The deflection of predator attacks from vital parts of the prey usually involves the coordinated evolution of prey body shape and colour. Here, we test the deflection effect of hindwing (HW) tails in the swallowtail butterfly Iphiclides podalirius. In this species, HWs display long tails associated with a conspicuous colour pattern. By surveying the wings within a wild population of I. podalirius, we observed that wing damage was much more frequent on the tails. We then used a standardized behavioural assay employing dummy butterflies with real I. podalirius wings to study the location of attacks by great tits Parus major. Wing tails and conspicuous coloration of the HWs were struck more often than the rest of the body by birds. Finally, we characterized the mechanical properties of fresh wings and found that the tail vein was more fragile than the others, suggesting facilitated escape ability of butterflies attacked at this location. Our results clearly support the deflective effect of HW tails and suggest that predation is an important selective driver of the evolution of wing tails and colour pattern in butterflies.

Divergence of climbing escape flight performance in Morpho butterflies living in different microhabitats.

Habitat specialization can influence the evolution of animal movement in promoting divergent locomotor abilities adapted to contrasting environmental conditions, differences in vegetation clutter or predatory communities. While the effect of habitat on the evolution of locomotion and particularly escape performance has been well investigated in terrestrial animals, it remains understudied in flying animals. Here, we investigated whether specialization of Morpho butterfly species into different vertical strata of the Amazonian forest affects the performance of upward escape flight manoeuvres. Using stereoscopic high-speed videography, we compared the climbing flight kinematics of seven Morpho species living either in the forest canopy or in the understory. We show that butterflies from canopy species display strikingly higher climbing speed and steeper ascent angle compared with understory species. Although climbing speed increased with wing speed and angle of attack, the higher climb angle observed in canopy species was best explained by their higher body pitch angle, resulting in more upward-directed aerodynamic thrust forces. Climb angle also scales positively with weight-normalized wing area, and this weight-normalized wing area was higher in canopy species. This shows that a combined divergence in flight behaviour and morphology contributes to the evolution of increased climbing flight abilities in canopy species.

The gene cortex controls mimicry and crypsis in butterflies and moths.

The wing patterns of butterflies and moths (Lepidoptera) are diverse and striking examples of evolutionary diversification by natural selection. Lepidopteran wing colour patterns are a key innovation, consisting of arrays of coloured scales. We still lack a general understanding of how these patterns are controlled and whether this control shows any commonality across the 160,000 moth and 17,000 butterfly species. Here, we use fine-scale mapping with population genomics and gene expression analyses to identify a gene, cortex, that regulates pattern switches in multiple species across the mimetic radiation in Heliconius butterflies. cortex belongs to a fast-evolving subfamily of the otherwise highly conserved fizzy family of cell-cycle regulators, suggesting that it probably regulates pigmentation patterning by regulating scale cell development. In parallel with findings in the peppered moth (Biston betularia), our results suggest that this mechanism is common within Lepidoptera and that cortex has become a major target for natural selection acting on colour and pattern variation in this group of insects.

Origin and persistence of polymorphism in loci targeted by disassortative preference: a general model.

The emergence and persistence of polymorphism within populations generally requires specific regimes of natural or sexual selection. Here, we develop a unified theoretical framework to explore how polymorphism at targeted loci can be generated and maintained by either disassortative mating choice or balancing selection due to, for example, heterozygote advantage. To this aim, we model the dynamics of alleles at a single locus A in a population of haploid individuals, where reproductive success depends on the combination of alleles carried by the parents at locus A. Our theoretical study of the model confirms that the conditions for the persistence of a given level of allelic polymorphism depend on the relative reproductive advantages among pairs of individuals. Interestingly, equilibria with unbalanced allelic frequencies were shown to emerge from successive introduction of mutants. We then investigate the role of the function linking allelic divergence to reproductive advantage on the evolutionary fate of alleles within the population. Our results highlight the significance of the shape of this function for both the number of alleles maintained and their level of genetic divergence. Large number of alleles are maintained with substantial replacement of alleles, when disassortative advantage slowly increases with allelic differentiation . In contrast, few highly differentiated alleles are predicted to be maintained when genetic differentiation has a strong effect on disassortative advantage. These opposite effects predicted by our model explain how disassortative mate choice may lead to various levels of allelic differentiation and polymorphism, and shed light on the effect of mate preferences on the persistence of balanced and unbalanced polymorphism in natural population.

When Do Opposites Attract? A Model Uncovering the Evolution of Disassortative Mating.

AbstractDisassortative mating is a rare form of mate preference that promotes the persistence of polymorphism. While the evolution of assortative mating and its consequences for trait variation and speciation have been extensively studied, the conditions enabling the evolution of disassortative mating are still poorly understood. Mate preferences increase the risk of missing mating opportunities, a cost that can be compensated by a greater fitness of offspring. Heterozygote advantage should therefore promote the evolution of disassortative mating, which maximizes the number of heterozygous offspring. From the analysis of a two-locus diploid model with one locus controlling the mating cue under viability selection and the other locus coding for the level of disassortative preference, we show that heterozygote advantage and negative frequency-dependent viability selection acting at the cue locus promote the evolution of disassortative preferences. We predict conditions of evolution of disassortative mating coherent with selection regimes acting on traits observed in the wild. We also show that disassortative mating generates sexual selection, which disadvantages heterozygotes at the cue locus, limiting the evolution of disassortative preferences. Altogether, our results partially explain why this behavior is rare in natural populations.

Convergence in sympatry: Evolution of blue-banded wing pattern in Morpho butterflies.

Species interactions such as mimicry can promote trait convergence but disentangling this effect from those of shared ecology, evolutionary history, and niche conservatism is often challenging. Here by focusing on wing colour pattern variation within and between three butterfly species living in sympatry in a large proportion of their range, we tested the effect of species interactions on trait diversification. These butterflies display a conspicuous iridescent blue coloration on the dorsal side of their wings and a cryptic brownish colour on the ventral side. Combined with an erratic and fast flight, these colour patterns increase the difficulty of capture by predators and contribute to the high escape abilities of these butterflies. We hypothesize that, beyond their direct contribution to predator escape, these wing patterns can be used as signals of escape abilities by predators, resulting in positive frequency-dependent selection favouring convergence in wing pattern in sympatry. To test this hypothesis, we quantified dorsal wing pattern variations of 723 butterflies from the three species sampled throughout their distribution, including sympatric and allopatric situations and compared the phenotypic distances between species, sex and localities. We detected a significant effect of localities on colour pattern, and higher inter-specific resemblance in sympatry as compared to allopatry, consistent with the hypothesis of local convergence of wing patterns. Our results provide support to the existence of escape mimicry in the wild and stress the importance of estimating trait variation within species to understand trait variation between species, and to a larger extent, trait diversification at the macro-evolutionary scale.

Assessing the Role of Developmental and Environmental Factors in Chemical Defence Variation in Heliconiini Butterflies.

Chemical defences in animals are both incredibly widespread and highly diverse. Yet despite the important role they play in mediating interactions between predators and prey, extensive differences in the amounts and types of chemical compounds can exist between individuals, even within species and populations. Here we investigate the potential role of environment and development on the chemical defences of warningly coloured butterfly species from the tribe Heliconiini, which can both synthesize and sequester cyanogenic glycosides (CGs). We reared 5 Heliconiini species in captivity, each on a single species-specific host plant as larvae, and compared them to individuals collected in the wild to ascertain whether the variation in CG content observed in the field might be the result of differences in host plant availability. Three of these species were reared as larvae on the same host plant, Passiflora riparia, to further test how species, sex, and age affected the type and amount of different defensive CGs, and how they affected the ratio of synthesized to sequestered compounds. Then, focusing on the generalist species Heliconius numata, we specifically explored variation in chemical profiles as a result of the host plant consumed by caterpillars and their brood line, using rearing experiments carried out on two naturally co-occurring host plants with differing CG profiles. Our results show significant differences in both the amount of synthesized and sequestered compounds between butterflies reared in captivity and those collected in the field. We also found a significant effect of species and an effect of sex in some, but not all, species. We show that chemical defences in H. numata continue to increase throughout their life, likely because of continued biosynthesis, and we suggest that variation in the amount of synthesized CGs in this species does not appear to stem from larval host plants, although this warrants further study. Interestingly, we detected a significant effect of brood lines, consistent with heritability influencing CG concentrations in H. numata. Altogether, our results point to multiple factors resulting in chemical defence variation in Heliconiini butterflies and highlight the overlooked effect of synthesis capabilities, which may be genetically determined to some extent.

How do predators generalize warning signals in simple and complex prey communities? Insights from a videogame.

The persistence of distinct warning signals within and between sympatric mimetic communities is a puzzling evolutionary question because selection favours convergence of colour patterns among toxic species. Such convergence is partly shaped by predators' reaction to similar but not identical stimulus (i.e. generalization behaviour), and generalization by predators is likely to be shaped by the diversity of local prey. However, studying generalization behaviour is generally limited to simple variations of prey colour patterns. Here, we used a computer game played by humans as surrogate predators to investigate generalization behaviours in simple (4 morphs) and complex (10 morphs) communities of unprofitable (associated with a penalty) and profitable butterflies. Colour patterns used in the game are observed in the natural populations of unprofitable butterfly species such as Heliconius numata. Analyses of 449 game participants' behaviours show that players avoided unprofitable prey more readily in simple than in complex communities. However, generalization was observed only in players that faced complex communities, enhancing the protection of profitable prey that looked similar to at least one unprofitable morph. Additionally, similarity among unprofitable prey also reduced attack rates only in complex communities. These results are consistent with previous studies using avian predators but artificial colour patterns and suggest that mimicry is more likely to evolve in complex communities where increases in similarity are more likely to be advantageous.

Polymorphism at a mimicry supergene maintained by opposing frequency-dependent selection pressures.

Explaining the maintenance of adaptive diversity within populations is a long-standing goal in evolutionary biology, with important implications for conservation, medicine, and agriculture. Adaptation often leads to the fixation of beneficial alleles, and therefore it erodes local diversity so that understanding the coexistence of multiple adaptive phenotypes requires deciphering the ecological mechanisms that determine their respective benefits. Here, we show how antagonistic frequency-dependent selection (FDS), generated by natural and sexual selection acting on the same trait, maintains mimicry polymorphism in the toxic butterfly Heliconius numata Positive FDS imposed by predators on mimetic signals favors the fixation of the most abundant and best-protected wing-pattern morph, thereby limiting polymorphism. However, by using mate-choice experiments, we reveal disassortative mate preferences of the different wing-pattern morphs. The resulting negative FDS on wing-pattern alleles is consistent with the excess of heterozygote genotypes at the supergene locus controlling wing-pattern variation in natural populations of H. numata The combined effect of positive and negative FDS on visual signals is sufficient to maintain a diversity of morphs displaying accurate mimicry with other local prey, although some of the forms only provide moderate protection against predators. Our findings help understand how alternative adaptive phenotypes can be maintained within populations and emphasize the need to investigate interactions between selective pressures in other cases of puzzling adaptive polymorphism.

Why has transparency evolved in aposematic butterflies? Insights from the largest radiation of aposematic butterflies, the Ithomiini.

Defended species are often conspicuous and this is thought to be an honest signal of defences, i.e. more toxic prey are more conspicuous. Neotropical butterflies of the large Ithomiini tribe numerically dominate communities of chemically defended butterflies and may thus drive the evolution of mimetic warning patterns. Although many species are brightly coloured, most are transparent to some degree. The evolution of transparency from a warning-coloured ancestor is puzzling as it is generally assumed to be involved in concealment. Here, we show that transparent Ithomiini species are indeed less detectable by avian predators (i.e. concealment). Surprisingly, transparent species are not any less unpalatable, and may in fact be more unpalatable than opaque species, the latter spanning a larger range of unpalatability. We put forth various hypotheses to explain the evolution of weak aposematic signals in these butterflies and other cryptic defended prey. Our study is an important step in determining the selective pressures and constraints that regulate the interaction between conspicuousness and unpalatability.

Effects of natural wing damage on flight performance in Morpho butterflies: what can it tell us about wing shape evolution?

Flying insects frequently experience wing damage during their life. Such irreversible alterations of wing shape affect flight performance and ultimately fitness. Insects have been shown to compensate for wing damage through various behavioural adjustments, but the importance of damage location over the wings has scarcely been studied. Using natural variation in wing damage, we tested how the loss of different wing parts affects flight performance. We quantified flight performance in two species of large butterflies, Morpho helenor and Morpho achilles, caught in the wild and displaying large variation in the extent and location of wing damage. We artificially generated more severe wing damage in our sample to contrast natural versus higher magnitude wing loss. Wing shape alteration across our sample was quantified using geometric morphometrics to test the effect of different damage distributions on flight performance. Our results show that impaired flight performance clearly depends on damage location over the wings, pointing to a relative importance of different wing parts for flight. A deteriorated forewing leading edge most critically affected flight performance, specifically decreasing flight speed and the proportion of gliding flight. In contrast, the most frequent natural damage, deteriorated wing margin, had no detectable effect on flight behaviour. Damage located on the hindwings - although having a limited effect on flight - was associated with reduced flight height, suggesting that the forewings and hindwings play different roles in butterfly flight. By contrasting harmless and deleterious consequences of various types of wing damage, our study highlights different selective regimes acting on morphological variations of butterfly wings.

Looking for the right mate in diploid species: How does genetic dominance affect the spatial differentiation of a sexual trait?

Divergence between populations for a given trait can be driven by sexual selection, interacting with migration behaviour. Mating preference for different phenotypes may lead to specific migration behaviour, with departures from populations where the preferred trait is rare. Such preferences can then trigger the emergence and persistence of differentiated populations, even without any local adaptation. However the genetic architecture underlying the trait targeted by mating preference may have a profound impact on population divergence. In particular, dominance between alleles encoding for divergent phenotypes can interfere with the differentiation process. Using a diploid model of a trait determining both mating success and migration rate, we explored differentiation between two connected populations, assuming either co-dominance or strict dominance between alleles. The model assumes that individuals prefer mating with partners displaying the same phenotype and therefore tend to move to the other population when their own phenotype is rare. We show that the emergence of differentiated populations in this diploid moded is limited as compared to results obtained with the same model assuming haploidy. When assuming co-dominance, differentiation arises only when migration is limited compared to the strength of the preference. Such differentiation is less dependent on migration when assuming strict dominance between haplotypes. Dominant alleles frequently invade populations because their phenotype is more frequently expressed, resulting in higher local mating success and a rapid decrease in migration. However, depending on the initial distribution of alleles, this advantage associated with dominance (i.e. Haldane's sieve) may lead to fixation of the dominant allele throughout both populations. Depending on the initial distribution of heterozygotes in the two populations, persistence of polymorphisms within populations can also occur because heterozygotes displaying the predominant phenotype benefit from high mating success. Altogether, our results highlight that heterozygotes' behaviour has a strong impact on population differentiation and highlight the need for diploid models of differentiation and speciation driven by sexual selection.

Genetic architecture and balancing selection: the life and death of differentiated variants.

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

Variation in cyanogenic compounds concentration within a Heliconius butterfly community: does mimicry explain everything?

BACKGROUND: Aposematic species advertise their unpalatability using warning signals such as striking coloration. Given that predators need to sample aposematic prey to learn that they are unprofitable, prey with similar warning signals share the cost of predator learning. This reduction in predation risk drives evolutionary convergence of warning signals among chemically defended prey (Müllerian mimicry). Whether such warning signal convergence is associated to similar defence levels among co-mimics is still an open question that has rarely been tested in wild populations. We quantified variation in cyanide-based (CN) chemical protection in wild caught individuals of eight aposematic Heliconius butterfly species belonging to four sympatric mimicry rings. We then tested for correlations between chemical protection and ecological species-specific traits. RESULTS: We report significant differences in CN concentrations both within and between sympatric species, even when accounting for the phylogeny, and within and between mimicry rings, even after considering inter-specific variation. We found significant correlations between CN concentration and both hostplant specialization and gregarious behaviour in adults and larvae. However, differences in CN concentrations were not significantly linked to mimicry ring abundance, although the two most toxic species did belong to the rarest mimicry ring. CONCLUSIONS: Our results suggest that mimicry can explain the variation in the levels of chemical defence to a certain extent, although other ecological factors are also relevant to the evolution of such variability.

Crossing fitness valleys: empirical estimation of a fitness landscape associated with polymorphic mimicry.

Characterizing fitness landscapes associated with polymorphic adaptive traits enables investigation of mechanisms allowing transitions between fitness peaks. Here, we explore how natural selection can promote genetic mechanisms preventing heterozygous phenotypes from falling into non-adaptive valleys. Polymorphic mimicry is an ideal system to investigate such fitness landscapes, because the direction of selection acting on complex mimetic colour patterns can be predicted by the local mimetic community composition. Using more than 5000 artificial butterflies displaying colour patterns exhibited by the polymorphic Müllerian mimic Heliconius numata, we directly tested the role of wild predators in shaping fitness landscapes. We compared predation rates on mimetic phenotypes (homozygotes at the supergene controlling colour pattern), intermediate phenotypes (heterozygotes), exotic morphs (absent from the local community) and palatable cryptic phenotypes. Exotic morphs were significantly more attacked than local morphs, highlighting predators' discriminatory capacities. Overall, intermediates were attacked twice as much as local homozygotes, suggesting the existence of deep fitness valleys promoting strict dominance and reduced recombination between supergene alleles. By including information on predators' colour perception, we also showed that protection on intermediates strongly depends on their phenotypic similarity to homozygous phenotypes and that ridges exist between similar phenotypes, which may facilitate divergence in colour patterns.

Selection of Valid Reference Genes for Reverse Transcription Quantitative PCR Analysis in Heliconius numata (Lepidoptera: Nymphalidae).

Identifying the genetic basis of adaptive variation is challenging in non-model organisms and quantitative real time PCR. is a useful tool for validating predictions regarding the expression of candidate genes. However, comparing expression levels in different conditions requires rigorous experimental design and statistical analyses. Here, we focused on the neotropical passion-vine butterflies Heliconius, non-model species studied in evolutionary biology for their adaptive variation in wing color patterns involved in mimicry and in the signaling of their toxicity to predators. We aimed at selecting stable reference genes to be used for normalization of gene expression data in RT-qPCR analyses from developing wing discs according to the minimal guidelines described in Minimum Information for publication of Quantitative Real-Time PCR Experiments (MIQE). To design internal RT-qPCR controls, we studied the stability of expression of nine candidate reference genes (actin, annexin, eF1α, FK506BP, PolyABP, PolyUBQ, RpL3, RPS3A, and tubulin) at two developmental stages (prepupal and pupal) using three widely used programs (GeNorm, NormFinder and BestKeeper). Results showed that, despite differences in statistical methods, genes RpL3, eF1α, polyABP, and annexin were stably expressed in wing discs in late larval and pupal stages of Heliconius numata This combination of genes may be used as a reference for a reliable study of differential expression in wings for instance for genes involved in important phenotypic variation, such as wing color pattern variation. Through this example, we provide general useful technical recommendations as well as relevant statistical strategies for evolutionary biologists aiming to identify candidate-genes involved adaptive variation in non-model organisms.

Which cues are sexy? The evolution of mate preference in sympatric species reveals the contrasted effect of adaptation and reproductive interference.

Mate preferences may target traits (a) enhancing offspring adaptation and (b) reducing heterospecific matings. Because similar selective pressures are acting on traits shared by different sympatric species, preference-enhancing offspring adaptation may increase heterospecific mating, in sharp contrast with the classical case of so-called "magic traits." Using a mathematical model, we study which and how many traits will be used during mate choice, when preferences for locally adapted traits increase heterospecific mating. In particular, we study the evolution of preference toward an adaptive versus a neutral trait in sympatric species. We take into account sensory trade-offs, which may limit the emergence of preference for several traits. Our model highlights that the evolution of preference toward adaptive versus neutral traits depends on the selective regimes acting on traits but also on heterospecific interactions. When the costs of heterospecific interactions are high, mate preference is likely to target neutral traits that become a reliable cue limiting heterospecific matings. We show that the evolution of preference toward a neutral trait benefits from a positive feedback loop: The more preference targets the neutral trait, the more it becomes a reliable cue for species recognition. We then reveal the key role of sensory trade-offs and the cost of choosiness favoring the evolution of preferences targeting adaptive traits, rather than traits reducing heterospecific mating. When sensory trade-offs and the cost of choosiness are low, we also show that preferences targeting multiple traits evolve, improving offspring fitness by both transmitting adapted alleles and reducing heterospecific mating. Altogether, our model aims at reconciling "good gene" and reinforcement models to provide general predictions on the evolution of mate preferences within natural communities.