In their perspective article, F Pouyet and KJ Gilbert (2020), propose an interesting overview of all the processes that sculpt patterns of molecular evolution. This well documented article covers most (if not all) important facets of the recurrent debate that has marked the history of molecular evolution: the relative importance of natural selection and neutral processes (i.e. genetic drift). I particularly enjoyed reading this review, that instead of taking a clear position on the debate, catalogs patiently every pieces of information that can help understand how patterns we observed at the genome level, can be understood from a selectionnist point of view, from a neutralist one, and, to quote their title, from "everything in between". The review covers the classical objects of interest in population genetics (genetic drift, selection, demography and structure) but also describes several genomic processes (meiotic drive, linked selection, gene conversion and mutation processes) that obscure the interpretation of these population processes. The interplay between all these processes is very complex (to say the least) and have resulted in many cases in profound confusions while analyzing data. It is always very hard to fully acknowledge our ignorance and we have many times payed the price of model misspecifications. This review has the grand merit to improve our awareness in many directions. Being able to cover so many aspects of a wide topic, while expressing them simply and clearly, connecting concepts and observations from distant fields, is an amazing "tour de force". I believe this article constitutes an excellent up-to-date introduction to the questions and problems at stake in the field of molecular evolution and will certainly also help established researchers by providing them a stimulating overview supported with many relevant references.

References

[1] Pouyet F, Gilbert KJ (2020) Towards an improved understanding of molecular evolution: the relative roles of selection, drift, and everything in between. arXiv:1909.11490 [q-bio]. ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. url:https://arxiv.org/abs/1909.11490

This manuscript [1] investigates the evolutionary history of the DAN gene family—a group of genes important for embryonic development of limbs, kidneys, and left-right axis speciation. This gene family has also been implicated in a number of diseases, including cancer and nephropathies. DAN genes have been associated with the inhibition of the bone morphogenetic protein (BMP) signaling pathway. Despite this detailed biochemical and functional knowledge and clear importance for development and disease, evolution of this gene family has remained understudied. The diversification of this gene family was investigated in all major groups of vertebrates. The monophyly of the gene members belonging to this gene family was confirmed. A total of five clades were delineated, and two novel lineages were discovered. The first lineage was only retained in cephalochordates (amphioxus), whereas the second one (GREM3) was retained by cartilaginous fish, holostean fish, and coelanth. Moreover, the patterns of chromosomal synteny in the chromosomal regions harboring DAN genes were investigated. Additionally, the authors reconstructed the ancestral gene repertoires and studied the differential retention/loss of individual gene members across the phylogeny. They concluded that the ancestor of gnathostome vertebrates possessed eight DAN genes that underwent differential retention during the evolutionary history of this group. During radiation of vertebrates, GREM1, GREM2, SOST, SOSTDC1, and NBL1 were retained in all major vertebrate groups. At the same time, GREM3, CER1, and DAND5 were differentially lost in some vertebrate lineages. At least two DAN genes were present in the common ancestor of vertebrates, and at least three DAN genes were present in the common ancestor of chordates. Therefore the patterns of retention and diversification in this gene family appear to be complex. Evolutionary slowdown for the DAN gene family was observed in mammals, suggesting selective constraints. Overall, this article puts the biomedical importance of the DAN family in the evolutionary perspective.

References

[1] Opazo JC, Hoffmann FG, Zavala K, Edwards SV (2020) Evolution of the DAN gene family in vertebrates. bioRxiv, 794404, ver. 3 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/794404

*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.](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](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](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](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](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](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.](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. ](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](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](https://doi.org/10.1371/journal.pbio.2005902)

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

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