Monnin et al. [1] here studied how Drosophila populations are affected when exposed to a high virulent endosymbiotic wMelPop Wolbachia strain and why virulent vertically transmitting endosymbionts persist in nature. This virulent wMelPop strain has been described to be a blocker of dengue and other arboviral infections in arthropod vector species, such as Aedes aegypti. Whereas it can thus function as a mutualistic symbiont, it here acts as an antagonist along the mutualism-antagonism continuum symbionts operate. The wMelPop strain is not a natural occurring strain in Drosophila melanogaster and thus the start of this experiment can be seen as a novel host-pathogen association. Through experimental evolution of 17 generations, the authors studied how high temperature affects wMelPop Wolbachia virulence and Drosophila melanogaster survival. The authors used Drosophila strains that were selected for late reproduction, given that this should favor evolution to a lower virulence. Assumptions for this hypothesis are not given in the manuscript here, but it can indeed be assumed that energy that is assimilated to symbiont tolerance instead of reproduction may lead to reduced virulence evolution. This has equally been suggested by Reyserhove et al. [2] in a dynamics energy budget model tailored to Daphnia magna virulence evolution upon a viral infection causing White fat Cell disease, reconstructing changing environments through time.
Contrary to their expectations for vertically transmitting symbionts, the authors did not find a reduction in wMelPop Wolbachia virulence during the course of the experimental evolution experiment under high temperature. Important is what this learns for virulence evolution, also for currently horizontal transmitting disease epidemics (such as COVID-19). It mainly reflects that evolution of virulence for new host-pathogen associations is difficult to predict and that it may take multiple generations before optimal levels of virulence are reached [3,4]. These optimal levels of virulence will depend on trade-offs with other life history traits of the symbiont, but also on host demography, host heterogeneity, amongst others [5,6]. Multiple microbial interactions may affect the outcome of virulence evolution [7]. Given that no germ-free individuals were used, it can be expected that other components of the Drosophila microbiome may have played a role in the virulence evolution. In most cases, microbiota have been described as defensive or protective for virulent symbionts [8], but they may also have stimulated the high levels of virulence. Especially, given that upon higher temperatures, Wolbachia growth may have been increased, host metabolic demands increased [9], host immune responses affected and microbial communities changed [10]. This may have resulted in increased competitive interactions to retrieve host resources, sustaining high virulence levels of the symbiont.
A nice asset of this study is that the phenotypic results obtained in the experimental evolution set-up were linked with wMelPop density measurement and octomom copy number quantifications. Octomom is a specific 8-n genes region of the Wolbachia genome responsible for wMelPop virulence, so there is a link between the phenotypic and molecular functions of the involved symbiont. The authors found that density, octomom copy number and virulence were correlated to each other. An important note the authors address in their discussion is that, to exclude the possibility that octomom copy number has an effect on density, and density on virulence, the effect of these variables should be assessed independently of temperature and age. The obtained results are a valuable contribution to the ongoing debate on the relationship between wMelPop octomom copy number, density and virulence.

References

[1] Monnin, D., Kremer, N., Michaud, C., Villa, M., Henri, H., Desouhant, E. and Vavre, F. (2020) Experimental evolution of virulence and associated traits in a Drosophila melanogaster – Wolbachia symbiosis. bioRxiv, 2020.04.26.062265, ver. 4 peer-reviewed and recommended by PCI Evol Biol. doi: https://doi.org/10.1101/2020.04.26.062265
[2] Reyserhove, L., Samaey, G., Muylaert, K., Coppé, V., Van Colen, W., and Decaestecker, E. (2017). A historical perspective of nutrient change impact on an infectious disease in Daphnia. Ecology, 98(11), 2784-2798. doi: https://doi.org/10.1002/ecy.1994
[3] Ebert, D., and Bull, J. J. (2003). Challenging the trade-off model for the evolution of virulence: is virulence management feasible?. Trends in microbiology, 11(1), 15-20. doi: https://doi.org/10.1016/S0966-842X(02)00003-3
[4] Houwenhuyse, S., Macke, E., Reyserhove, L., Bulteel, L., and Decaestecker, E. (2018). Back to the future in a petri dish: Origin and impact of resurrected microbes in natural populations. Evolutionary Applications, 11(1), 29-41. doi: https://doi.org/10.1111/eva.12538
[5] Day, T., and Gandon, S. (2007). Applying population‐genetic models in theoretical evolutionary epidemiology. Ecology Letters, 10(10), 876-888. doi: https://doi.org/10.1111/j.1461-0248.2007.01091.x
[6] Alizon, S., Hurford, A., Mideo, N., and Van Baalen, M. (2009). Virulence evolution and the trade‐off hypothesis: history, current state of affairs and the future. Journal of evolutionary biology, 22(2), 245-259. doi: https://doi.org/10.1111/j.1420-9101.2008.01658.x
[7] Alizon, S., de Roode, J. C., and Michalakis, Y. (2013). Multiple infections and the evolution of virulence. Ecology letters, 16(4), 556-567. doi: https://doi.org/10.1111/ele.12076
[8] Decaestecker, E., and King, K. (2019). Red queen dynamics. Reference module in earth systems and environmental sciences, 3, 185-192. doi: https://doi.org/10.1016/B978-0-12-409548-9.10550-0
[9] Kirk, D., Jones, N., Peacock, S., Phillips, J., Molnár, P. K., Krkošek, M., and Luijckx, P. (2018). Empirical evidence that metabolic theory describes the temperature dependency of within-host parasite dynamics. PLoS biology, 16(2), e2004608. doi: https://doi.org/10.1371/journal.pbio.2004608
[10] Frankel-Bricker, J., Song, M. J., Benner, M. J., and Schaack, S. (2019). Variation in the microbiota associated with Daphnia magna across genotypes, populations, and temperature. Microbial ecology, 1-12. doi: https://doi.org/10.1007/s00248-019-01412-9

150 years after Darwin published ‘Descent of man and selection in relation to sex’ (Darwin, 1871), the evolutionary mechanism that he laid out in his treatise continues to fascinate us. Sexual selection is responsible for some of the most spectacular traits among animals, and plants, and it appeals to our interest in all things reproductive and sexual (Bell, 1982). In addition, sexual selection poses some of the more intractable problems in evolutionary biology: Its realm encompasses traits that are subject to markedly different selection pressures, particularly when distinct, yet associated, traits tend to be associated with males, e.g. courtship signals, and with females, e.g. preferences (cf. Ah-King & Ahnesjo, 2013). While separate, such traits cannot evolve independently of each other (Arnqvist & Rowe, 2005), and complex feedback loops and correlations between them are predicted (Greenfield et al., 2014). Traditionally, sexual selection has been modelled under simplifying assumptions, and quantitative genetic approaches that avoided evolutionary dynamics have prevailed. New computing methods may be able to free the field from these constraints, and a trio of theoreticians (Chevalier, De Coligny & Labonne 2020) describe here a novel application of a ‘demo-genetic agent (or individual) based model’, a mouthful hereafter termed DG-ABM, for arriving at a holistic picture of the sexual selection trajectory. The application is built on the premise that traits, e.g. courtship, preference, gamete investment, competitiveness for mates, can influence the genetic architecture, e.g. correlations, of those traits. In turn, the genetic architecture can influence the expression and evolvability of the traits. Much of this influence occurs via demographic features, i.e. social environment, generated by behavioral interactions during sexual advertisement, courtship, mate guarding, parental care, post-mating dispersal, etc.
The authors provide a lengthy verbal description of their model, specifying the genomic and behavioral parameters that can be set and how a ‘run’ may be initialized. There is a link to an internet site where users can then enter their own parameter values and begin exploring hypotheses. Back in the article several simulations illustrate simple tests; e.g. how gamete investment and preference jointly evolve given certain survival costs. One obvious test would have been the preference – courtship genetic correlation that represents the core of Fisherian runaway selection, and it is regrettable that it was not examined under a range of demographic parameters. As presented the author’s DG-ABM appears particularly geared toward mating systems in ‘higher’ vertebrates, where couples form during a discrete mating season and are responsible for most reproduction. It is not clear how applicable the model could be to a full range of mating systems and nuances, including those in arthropods and other invertebrates as well as plants.
What is the likely value of the DG-ABM for sexual selection researchers? We will not be able to evaluate its potential impact until readers with specialized understanding of a question and taxon begin exploring and comparing their results with prior expectations. Of course, lack of congruence with earlier predictions would not invalidate the model. Hopefully, some of these specialists will have opportunities for comparing results with pertinent empirical data.

References

Ah-King, M. and Ahnesjo, I. 2013. The ‘sex role’ concept: An overview and evaluation Evolutionary Biology, 40, 461-470. doi: https://doi.org/10.1007/s11692-013-9226-7
Arnqvist, G. and Rowe, L. 2005. Sexual Conflict. Princeton University Press, Princeton. doi: https://doi.org/10.1515/9781400850600
Bell, G. 1982. The Masterpiece of Nature: The Evolution and Genetics of Sexuality. University of California Press, Berkeley.
Chevalier, L., De Coligny, F. and Labonne, J. (2020) A demogenetic individual based model for the evolution of traits and genome architecture under sexual selection. bioRxiv, 2020.04.01.014514, ver. 4 peer-reviewed and recommended by PCI Evol Biol. doi: https://doi.org/10.1101/2020.04.01.014514
Darwin, C. 1871. The Descent of Man and Selection in Relation to Sex. J. Murray, London.
Greenfield, M.D., Alem, S., Limousin, D. and Bailey, N.W. 2014. The dilemma of Fisherian sexual selection: Mate choice for indirect benefits despite rarity and overall weakness of trait-preference genetic correlation. Evolution, 68, 3524-3536. doi: https://doi.org/10.1111/evo.12542

One can only hope that one day, we will be able to evaluate how the ecological complexity surrounding natural populations affects their ability to adapt. This is more like a long term quest than a simple scientific aim. Many steps are heading in the right direction. This paper by Monroe and colleagues (2021) is one of them.
Many ecological and genetic mechanisms shape the evolutionary potential of phenotypic trait variation and many of them involve environmental heterogeneity (Pujol et al 2018). To date, we cannot look into these ecological and genetic mechanisms without oversimplifying their effects. We often look into trait variation one trait at a time albeit the variation of multiple phenotypic traits is often linked at the genetic or environmental level. As a consequence, we put our conclusions at risk by not accounting for the reciprocal impacts of trait changes upon each other (Teplitsky et al 2014). We also usually restrict the study of a continuous gradient of environmental conditions to a few conditions because it would otherwise be impossible to model its environmental effect. As a consequence, we miss the full picture of the continuous often nonlinear phenotypic plastic response. Whether the trait undergo threshold effect changes thereby remains obscured to us. Collectively, these issues impede our ability to understand how selection shapes the ecological strategy of organisms under variable environments.
In this paper, Monroe and colleagues (2021) propose an original approach that raised to these two challenges. They analysed phenotypic plastic changes in response to a continuous environment in a multidimensional trait space, namely the response of Brachypodium plant developmental and physiological traits to a continuous gradient of soil moisture. They used dry down experimental treatments to produce the continuous soil moisture gradient and compared the plant capacity to use water between annual B. distachyon and perennial B. sylvaticum. Their results revealed the best mathematical functions that model the nonlinear curvature of the continuous plastic response of Brachypodium plants. This work reinforces our view that nonlinear plastic responses can result in greater or lesser trait values at any stage of the environmental gradient that were unexpected on the basis of linear predictors (Gienapp and Brommer 2014). Their findings also imply that different threshold responses characterize different genotypes. These could otherwise have been missed by a classical approach. By shedding light on unforeseen interactions between traits that make their correlation vary along the nonlinear response, they were able to describe more accurately Brachypodium ecological strategies and the changes in evolutionary constraints along the soil moisture gradient.
Their empirical approach allows to test what environmental conditions maximises the opportunity for selection to shape trait variation. For example, it revealed unforeseen divergence in potentially adaptive mechanisms or life history strategies – and not just trait values – between annual and perennial species of Brachypodium. Behind every environmental variation of the constraints to the future evolutionary change of multiple traits, we can expect that the evolutionary history of the populations shaped their trait genetic correlations. Investigating the nonlinear signature of adaptive evolution across continuous environments will get us into uncharted territory.
Our ability to predict the adaptive potential of species is limited. With their approach of continuous environmental gradients beyond linearity, Monroe and collaborators (2021) improve our understanding of plant phenotypic responses and open a brand new range of exciting developments. As they mention: "the opportunity for scaling up" their approach is big. To illustrate this prospect, I can easily think of an example: the quantitative genetic random regression model. This model allows to use any degree of genetic relatedness in a wild population to estimate the genetic variation of phenotypic plastic reaction norms (Nussey et al 2007, Pujol and Galaud 2013). However, in this approach, only a few modalities of the environmental gradient are used to model nonlinear phenotypic plastic responses. From there, it is rather intuitive. Combining the best of these two approaches (continuity of genetic relatedness in the wild & continuity of environmental gradient in experiments) could open ground breaking new perspectives in research.

References

Gienapp P. & J.E. Brommer. 2014. Evolutionary dynamics in response to climate change. In: Charmentier A, Garant D, Kruuk LEB, editors. Quantitative genetics in the wild. Oxford: Oxford University Press, Oxford. pp. 254–273. doi: https://doi.org/10.1093/acprof:oso/9780199674237.003.0015
Monroe, J. G., Cai, H., and Des Marais, D. L. (2020). Trait plasticity and covariance along a continuous soil moisture gradient. bioRxiv, 2020.02.17.952853, ver. 5 peer-reviewed and recommended by PCI Evol Biol. doi: https://doi.org/10.1101/2020.02.17.952853
Pujol et al. (2018). The missing response to selection in the wild. Trends in ecology & evolution, 33(5), 337-346. doi: https://doi.org/10.1016/j.tree.2018.02.007
Pujol, B., and Galaud, J. P. (2013). A practical guide to quantifying the effect of genes underlying adaptation in a mixed genomics and evolutionary ecology approach. Botany Letters, 160(3-4), 197-204. doi: https://doi.org/10.1080/12538078.2013.799045
Nussey, D. H., Wilson, A. J., and Brommer, J. E. (2007). The evolutionary ecology of individual phenotypic plasticity in wild populations. Journal of evolutionary biology, 20(3), 831-844. doi: https://doi.org/10.1111/j.1420-9101.2007.01300.x
Teplitsky et al. (2014). Assessing multivariate constraints to evolution across ten long-term avian studies. PLoS One, 9(3), e90444. doi: https://doi.org/10.1371/journal.pone.0090444

Biological invasions are natural experiments, and often show that evolution can affect dynamics in important ways [1-3]. While we often think of invasions as a conservation problem stemming from anthropogenic introductions [4,5], biological invasions are much more commonplace than this, including phenomena as diverse as natural range shifts, the spread of novel pathogens, and the growth of tumors. A major question across all these settings is which set of traits determine the ability of a population to invade new space [6,7]. Traits such as: increased growth or reproductive rate, dispersal ability and ability to defend from predation often show large evolutionary shifts across invasion history [1,6,8]. Are such multi-trait shifts driven by selection on multiple traits, or a correlated response by multiple traits to selection on one? Resolving this question is important for both theoretical and practical reasons [9,10]. But despite the importance of this issue, it is not easy to perform the necessary manipulative experiments [9].
Foucaud et al. [11] tackled this issue by performing experimental evolution on source populations of the invasive ladybug Harmonia axyridis. The authors tested if selection on a single trait could generate correlated responses in other life history traits. Specifically, they used experimental evolution to impose divergent selection on female mass, and reproductive timing. After ten generations, they found that selection for weight did not affect almost any other life history trait. However, nine generations of selection for faster reproduction led to correlated phenotypic changes in developmental, reproduction and survival rate of populations, although not always in the direction we might have expected. Despite this correlated response, none of their selected lines were able to fully recapitulate the trait shifts seen in natural invasions of this species. This implies that selection during natural invasions is operating on multiple traits; a finding in agreement with our growing understanding of how selection acts during introduction and invasion [12,13].
Populations undergoing a colonization process may also be subject to a multitude of different selective pressures [14,15]. The authors expanded their work in this direction by testing whether food availability alters the observed correlations between life history traits. The pervasiveness of genotype by environment interactions observed also points to a role for multiple selective pressures in shaping the suite of life-history shifts observed in wild ladybug populations. The work from Foucaud and colleagues [11] adds to a small but growing list of important studies that use experimental evolution to investigate how life-history traits evolve, and how they evolve during invasions in particular.

References

[1] Sakai, A.K., Allendorf, F.W., Holt, J.S. et al. (2001). The population biology of invasive species. Annual review of ecology and systematics, 32(1), 305-332. doi: 10.1146/annurev.ecolsys.32.081501.114037
[2] Hairston Jr, N. G., Ellner, S. P., Geber, M. A., Yoshida, T. and Fox, J. A. (2005). Rapid evolution and the convergence of ecological and evolutionary time. Ecology letters, 8(10), 1114-1127. doi: 10.1111/j.1461-0248.2005.00812.x
[3] Chuang, A. and Peterson, C. R. (2016). Expanding population edges: theories, traits, and trade‐offs. Global change biology, 22(2), 494-512. doi: 10.1111/gcb.13107
[4] Whitney, K. D. and Gabler, C. A. (2008). Rapid evolution in introduced species,‘invasive traits’ and recipient communities: challenges for predicting invasive potential. Diversity and Distributions, 14(4), 569-580. doi: 10.1111/j.1472-4642.2008.00473.x
[5] Catullo, R. A., Llewelyn, J., Phillips, B. L. and Moritz, C. C. (2019). The Potential for Rapid Evolution under Anthropogenic Climate Change. Current Biology, 29(19), R996-R1007. doi: 10.1016/j.cub.2019.08.028
[6] Suarez, A. V. and Tsutsui, N. D. (2008). The evolutionary consequences of biological invasions. Molecular Ecology, 17(1), 351-360. doi: 10.1111/j.1365-294X.2007.03456.x
[7] Deforet, M., Carmona-Fontaine, C., Korolev, K. S. and Xavier, J. B. (2019). Evolution at the edge of expanding populations. The American Naturalist, 194(3), 291-305. doi: 10.1086/704594
[8] Phillips, B. L., Brown, G. P., and Shine, R. (2010). Life‐history evolution in range‐shifting populations. Ecology, 91(6), 1617-1627. doi: 10.1890/09-0910.1
[9] Colautti, R. I. and Lau, J. A. (2015). Contemporary evolution during invasion: evidence for differentiation, natural selection, and local adaptation. Molecular ecology, 24(9), 1999-2017. doi: 10.1111/mec.13162
[10] Szűcs, M., Melbourne, B. A., Tuff, T., Weiss‐Lehman, C. and Hufbauer, R. A. (2017). Genetic and demographic founder effects have long‐term fitness consequences for colonising populations. Ecology Letters, 20(4), 436-444. doi: 10.1111/ele.12743
[11] Foucaud, J., Hufbauer, R. A., Ravigné, V., Olazcuaga, L., Loiseau, A., Ausset, A., Wang, S., Zang, L.-S., Lemenager, N., Tayeh, A., Weyna, A., Gneux, P., Bonnet, E., Dreuilhe, V., Poutout, B., Estoup, A. and Facon, B. (2020). How do invasion syndromes evolve? An experimental evolution approach using the ladybird Harmonia axyridis. bioRxiv, 849968 ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/849968
[12] Simons, A. M. (2003). Invasive aliens and sampling bias. Ecology Letters, 6(4), 278-280. doi: 10.1046/j.1461-0248.2003.00430.x
[13] Phillips, B. L. and Perkins, T. A. (2019). Spatial sorting as the spatial analogue of natural selection. Theoretical Ecology, 12(2), 155-163. doi: 10.1007/s12080-019-0412-9
[14] Lavergne, S. and Molofsky, J. (2007). Increased genetic variation and evolutionary potential drive the success of an invasive grass. Proceedings of the National Academy of Sciences, 104(10), 3883-3888. doi: 10.1073/pnas.0607324104
[15] Moran, E. V. and Alexander, J. M. (2014). Evolutionary responses to global change: lessons from invasive species. Ecology Letters, 17(5), 637-649. doi: 10.1111/ele.12262

Heliconius butterflies are famous for their colorful wing patterns acting as a warning of their chemical defenses [1]. Most species are involved in Müllerian mimicry assemblies, as predators learn to associate common wing patterns with unpalatability and preferentially target rare variants. Such positive-frequency dependent selection homogenizes wing patterns at different localities, and in several species, all individuals within a community belong to the same morph [2]. In this respect, H. numata stands out. This species shows stable local polymorphism across multiple localities, with local populations home to up to seven distinct morphs [2]. Although a balance between migration and local positive-frequency dependent selection can allow some degree of local polymorphism, theory suggests that this occurs only when migration is within a narrow window [3].
One factor that potentially enhances local polymorphism in H. numata is disassortative mating. Mate choice assays have in fact revealed that females of this species tend to reject males with the same wing pattern [4]. However the evolution of such mating behavior and its effect on polymorphism remain unclear when selection is locally positive-frequency dependent. Using a mathematical model, Maisonneuve et al. [5] clarify the conditions that favor the evolution of disassortative mating in the complicated system of H. numata. In particular, they investigate whether the genetic basis of wing colour can favor the emergence of disassortative mating. Variation in wing pattern in H. numata is controlled by the supergene P, which is a single genomic region harboring multiple protein coding genes that have ceased to recombine due to chromosomal inversions [6]. If such remarkable genetic configuration allows for the co-adaptation of multiple loci participating to a complex phenotype such as wing color pattern, the absence of recombination can also result in the accumulation of deleterious mutations [7]. In fact, alleles at the P locus have been associated with a recessive genetic load, leading to a fitness advantage for heterozygotes at this locus [8]. Can this fitness advantage to heterozygotes lead to the evolution of disassortative mating? And if so, can such evolution lead to the maintenance of local polymorphism in spite of strong positive frequency-dependent selection?
To investigate these questions, Maisonneuve et al. [5] model evolution at two loci, one is the P locus for wing pattern, and the other influences mating behavior. The population is divided among two connected patches that differ in their butterfly communities, so that different alleles at the P locus are favored by positive frequency-dependent selection in different patches. The different alleles at the P locus are ordered in dominance relationships such that the most dominant over wing color pattern are also those with the highest load. By tracking the dynamics of haplotype frequencies in the population, the authors first show that disassortative mating readily evolves via the invasion of an allele causing females carrying it to reject males that resemble them phenotypically. Such “self-referencing” mechanism of mate choice, however, has never been reported and has been argued to be rare due to its complicated nature [9].
Maisonneuve et al. [5] then compare the evolution of disassortative mating via two alternative mechanisms: attraction and rejection. In these cases, alleles at the mating locus determine attraction to or rejection of specific phenotypes (e.g., under attraction rule, allele “B” encodes attraction to males with phenotype B). With the P and mating loci fully linked, disassortative mating can evolve under all three mechanisms (self-referencing, attraction and rejection), but tends to be less prevalent at equilibrium under attraction rule. This in turn results in the maintenance of less genetic variation under attraction compared to the other mating mechanisms. The loss of variation that occurs under attraction rules is due to a combination of dominance relationships between alleles at the P locus and the searching cost to females in finding rare types of males. When a particular wing pattern, say B, is only expressed in homozygotic form, B males are relatively rare. Females that carry the allele at the mating locus causing them to be attracted to such males then suffer a fitness cost due to lost mating opportunities. This mating allele is therefore purged, and in turn so is the recessive allele for B phenotype at the P locus. Under self-referencing and rejection rules, however, choosy females only reject males of a specific phenotype. They can therefore potentially mate with larger pool of males than females attracted to a single type. As a result, self-referencing and rejection rules are less sensitive to demographic effects and so are more conducive to disassortative mating evolution.
In their final analysis, Maisonneuve et al. [5] investigate the influence of recombination among the P and mating loci. They show that recombination has different effects on disassortative mating evolution depending on the mechanism of mate choice. Under the self-referencing rule, loose linkage leads to higher levels of disassortative mating and polymorphism than when linkage is tight. Under attraction or rejection rule, however, even very limited recombination completely inhibits the evolution of disassortative mating. This is because, with alleles at the mating locus coding for attraction/rejection to specific males, recombination breaks the association between the P and mating loci necessary for disassortative mating. By contrast, disassortative mating via a self-referencing rule does not depend on the linkage among the P and mating loci: females choose males that are different to themselves independently from the alleles they carry at the P locus.
Taken together, Maisonneuve et al.’s analyses [5] show that disassortative mating can readily evolve in a system like H. numata, but that this evolution depends on the genetic architecture of mating behavior. The architectures that are more conducive to the evolution of disassortative mating are: (1) epistatic interactions among the P and mating loci such that females are able to recognize their own phenotype and base their mating decision upon this information (self-referencing rule); and (2) full linkage among the P supergene and a mating locus that triggers rejection of a specific color pattern. While the mechanisms behind disassortative mating remain to be elucidated, assortative mating seems to rely on alleles triggering attraction to specific cues with variation in attraction and cues linked together [10]. These observations support the notion that disassortative mating is due to alleles causing rejection, in tight linkage to the P locus. If so, mating loci would in fact be part of the P supergene, thus controlling not only intricate wing color pattern but also mating behavior.
Beyond the specific system of H. numata, Maisonneuve et al.’s study [5] helps understand the evolution of disassortative mating and its association with the genetic architecture of correlated traits. In particular, Maisonneuve et al. [5] expands the role of supergenes for ecologically relevant traits to mating behavior, further bolstering the relevance of these remarkable genetic elements in the maintenance of variation in complex and elaborate phenotypes.

References

[1] Merrill, R M, K K Dasmahapatra, J W Davey, D D Dell'Aglio, J J Hanly, B Huber, C D Jiggins, et al. (2015). The Diversification of Heliconius butterflies: What Have We Learned in 150 Years? Journal of Evolutionary Biology 28 (8), 1417–38. https://doi.org/10.1111/jeb.12672.
[2] Joron M, IR Wynne, G Lamas, and J Mallet (1999). Variable selection and the coexistence of multiple mimetic forms of the butterfly Heliconius numata. Evolutionary Ecology 13, 721– 754. https://doi.org/10.1023/A:1010875213123
[3] Joron M and Y Iwasa (2005). The evolution of a Müllerian mimic in a spatially distributed community. Journal of Theoretical Biology 237, 87–103. https://doi.org/10.1016/j.jtbi.2005.04.005
[4] Chouteau M, V Llaurens, F Piron-Prunier, and M Joron (2017). Polymorphism at a mimicry su- pergene maintained by opposing frequency-dependent selection pressures. Proceedings of the National Academy of Sciences 114, 8325–8329. https://doi.org/10.1073/pnas.1702482114
[5] Maisonneuve, L, Chouteau, M, Joron, M and Llaurens, V. (2020). Evolution and genetic architecture of disassortative mating at a locus under heterozygote advantage. bioRxiv, 616409, ver. 9 peer-reviewed and recommended by PCI Evolutionary Biology. https://doi.org/10.1101/616409
[6] Joron M, L Frezal, RT Jones, NL Chamberlain, SF Lee, CR Haag, A Whibley, M Becuwe, SW Baxter, L Ferguson, et al. (2011). Chromosomal rearrangements maintain a polymorphic super- gene controlling butterfly mimicry. Nature 477, 203. https://doi.org/10.1038/nature10341
[7] Schwander T, R Libbrecht, and L Keller (2014). Supergenes and Complex Phenotypes.” Current Biology. 24 (7), 288–94. https://doi.org/10.1016/j.cub.2014.01.056.
[8] Jay P, M Chouteau, A Whibley, H Bastide, V Llaurens, H Parrinello, and M Joron (2019). Mutation accumulation in chromosomal inversions maintains wing pattern polymorphism in a butterfly. bioRxiv. https://doi.org/ 10.1101/736504.
[9] Kopp M, MR Servedio, TC Mendelson, RJ Safran, RL Rodrıguez, ME Hauber, EC Scordato, LB Symes, CN Balakrishnan, DM Zonana, et al. (2018). Mechanisms of assortative mating in speciation with gene flow: connecting theory and empirical research. The American Naturalist 191, 1–20. https://doi.org/10.1086/694889
[10] Merrill RM, P Rastas, SH Martin, MC Melo, S Barker, J Davey, WO McMillan, and CD Jiggins (2019). Genetic dissection of assortative mating behavior. PLoS biology 17, e2005902. https://doi.org/10.1371/journal.pbio.2005902