Koutsovoulos et al. [1] have generated and analysed the first population genomic dataset in root-knot nematode Meloidogyne incognita. Why is this interesting? For two major reasons. First, M. incognita has been documented to be apomictic, i.e., to lack any form of sex. This is a trait of major evolutionary importance, with implications on species adaptive potential. The study of genome evolution in asexuals is fascinating and has the potential to inform on the forces governing the evolution of sex and recombination. Even small amounts of sex, however, are sufficient to restore most of the population genetic properties of true sexuals [2]. Because rare events of sex can remain undetected in the field, to confirm asexuality in M. incognita using genomic data is an important step. The second reason why M. incognita is of interest is that this nematode is one of the most harmful pests currently living on earth. M. incognita feeds on the roots of many cultivated plants, including tomato, bean, and cotton, and has been of major agricultural importance for decades. A number of races were defined based on host specificity. These have played a key role in attempts to control the dynamic of M. incognita populations via crop rotations. Races and management strategies so far lack any genetic basis, hence the second major interest of this study.
The authors newly sequenced the full genome of eleven strains from Brazil and added nine already available samples from Africa and North-America. They report that, in all likelihood, M. incognita is indeed a purely asexual species. This is supported by (i) the confirmation that the genome is in its major part haploid, and (ii) a spectacularly high level of linkage disequilibrium, which does not decline with genetic distance between loci at a 100kb scale. The absence of sex and recombination is associated in M. incognita with a remarkably low amount of genetic diversity - one order of magnitude less than in typical sexual nematodes - and an heavy load of deleterious mutations, as measured by the ratio of non-synonymous (=amino-acid changing) to synonymous (=amino-acid conservative) diversity in coding sequences. The other important result of this study is that the population substructure in M. incognita is in no way related to host races or geography. The tree genetic clusters that are identified include strains from several continents and feeding on a diversity of host plants.
The implications of this work are numerous. First, the results suggest that M. incognita is an ancient asexual. Asexuality, which was here demonstrated via linkage disequilibrium analysis, must be ancient enough for diploidy (or, in this case, maybe triploidy) to have been lost - i.e., formerly homologous chromosomes have accumulated enough mutations to be assembled as distinct entities. So we are not talking about a highly successful clone having recently spread the world - rather a long-term obligate parthenogen. Asexual organisms are deprived of the source of genetic variation offered by recombination, which is why asexuality is thought to be an evolutionary dead-end. Long-term asexuals are uncommon and even the most famous ones, bdelloid rotifers, are suspected to experience between-individual genetic transfers [3]. M. incognita is apparently a true 'evolutionary scandal', and as such deserves particular attention from molecular evolutionary geneticists.
The lack of any host race effect on the genetic diversity of M. incognita is another important finding. So-called 'races' have largely contributed to shape researchers' view of the structure of the species so far. This study demonstrates that a mental effort is now needed to forget about races, and consider host-specificity for what it is - a phenotypic trait. This result implies that many host shifts must have independently occurred in the three M. incognita genetic lineages, suggesting an arms race between plants and nematodes, which in the absence of sex and recombination must be entirely mutation-driven on the nematode side. Genes functionally involved in the arms race might therefore be expected to have experienced convergent evolution, if distinct M. incognita lineages have adopted the same solutions to overcome plant defenses. The present study paves the way for such a genome scan. The authors rightly discuss that the strong adaptive potential of M. incognita, at least in terms of host shift, despite no sex and tiny amounts of genetic diversity, is a paradox that would deserve to be further investigated.

References

[1] Koutsovoulos, G. D., Marques, E., Arguel, M. J., Duret, L., Machado, A. C. Z., Carneiro, R. M. D. G., Kozlowski, D. K., Bailly-Bechet, M., Castagnone-Sereno, P., Albuquerque, E. V., & Danchin, E. G. J. (2019). Population genomics supports clonal reproduction and multiple gains and losses of parasitic abilities in the most devastating nematode plant pest. bioRxiv, 362129, ver. 5, peer-reviewed and recommended by Peer Community in Evolutionary Biology. doi: 10.1101/362129
[2] Hartfield, M. (2016). Evolutionary genetic consequences of facultative sex and outcrossing. Journal of evolutionary biology, 29(1), 5-22. doi: 10.1111/jeb.12770
[3] Debortoli, N., Li, X., Eyres, I., Fontaneto, D., Hespeels, B., Tang, C. Q., Flot, J. F. & Van Doninck, K. (2016). Genetic exchange among bdelloid rotifers is more likely due to horizontal gene transfer than to meiotic sex. Current Biology, 26(6), 723-732. doi: 10.1016/j.cub.2016.01.031

Flowering plants display a staggering diversity of both mating systems and life histories, ranging from almost exclusively selfers to obligate outcrossers, very short-lived annual herbs to super long lived trees. One pervasive pattern that has attracted considerable attention is the tight correlation that is found between mating systems and lifespan [1]. Until recently, theoretical explanations for these patterns have relied on static models exploring the consequences of several non-mutually exclusive important process: levels of inbreeding depression and ability to successfully were center stage. This make sense: successful colonization after long‐distance dispersal is far more likely to happen for self‐compatible than for self‐incompatible individuals in a sexually reproducing species. Furthermore, inbreeding depression (essentially a genetically driven phenomenon) and reproductive insurance are expected to shape the evolution of both mating system and lifespan.
But modelling jointly several processes and how their interplay to shape the evolution of a trait is challenging enough so models for the evolution of mating system tend invariably – for mathematical convenience and tractability – to fix lifespan [2].
However, comparative analysis of between species variations that map traits transitions among sister species in phylogenetic trees reveals a pervasive pattern: frequent transitions from a state outcrossing perennial to selfing annuals. This beg the question: is one transition triggering the other and if so, what comes first or are these transitions happening together? In this work, Lesaffre and Billiard use a very sophisticated machinery developed by Kirkpatrick et al. [3] and consider a general class of so-called modifiers models [4]. They study jointly the evolution of life span and mating system. They do so by using models where different life stages are tracked with life stage having some (fixed for once) amount of inbreeding depression. Their paper is technically demanding, mixing analytics and computer simulations, and along the way generates several important findings that are expected to stimulate further empirical and theoretical studies: (1) pure selfing versus pure outcrossing is the expected stable evolutionary outcomes (despite observation that mixed mating systems can be regularly met in nature), (2) increasing life-span drastically reduces the scope for the evolution of selfing, conversely (3) transition to selfing will also select for shorter life span as a way to mitigate the cumulative effects of inbreeding depression on adult life stages.
As usual there is room for future work, in particular the authors’ model assumes fixed inbreeding depression in the different life stages and this highlights the need for models that explore how inbreeding depression, a pivotal quantity in these models, can itself be molded by both mating system and lifespan. A third-generation of models should be “soon” on the way!

References
[1] Grossenbacher D, Briscoe Runquist R, Goldberg EE, and Brandvain Y. (2015) Geographic range size is predicted by plant mating system. Ecology Letters 18, 706–713. doi: 10.1111/ele.12449
[2] Morgan MT, Schoen DJ, and Bataillon T. (1997) The evolution of self-fertilization in perennials. The American Naturalist 150, 618–638. doi: 10.1086/286085
[3] Kirkpatrick M, Johnson T, and Barton N. (2002) General models of multilocus evolution. Genetics 161, 1727–1750.
[4] Lesaffre, T, and Billiard S. (2019) The joint evolution of lifespan and self-fertilisation. bioRxiv, 420877, ver. 3 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/420877

The ability of a population to adapt to a new niche is an important phenomenon in evolutionary biology. The colonisation of a new volcanic island by plant species; the colonisation of a host treated by antibiotics by a-resistant strain; the Ebola virus transmitting from bats to humans and spreading epidemically in Western Africa, are all examples of a population invading a new niche, adapting and eventually establishing in this new environment.

Adaptation to a new niche can be studied using source-sink models. In the original environment —the “source”—, the population enjoys a positive growth-rate and is self-sustaining, while in the new environment —the “sink”— the population has a negative growth rate and is able to sustain only by the continuous influx of migrants from the source. Understanding the dynamics of adaptation to the sink environment is challenging from a theoretical standpoint, because it requires modelling the demography of the sink as well as the transient dynamics of adaptation. Moreover, local selection in the sink and immigration from the source create distributions of genotypes that complicate the use of many common mathematical approaches.

In their paper, Lavigne et al. [1], develop a new deterministic model of adaptation to a harsh sink environment in an asexual species. The fitness of an individual is maximal when a number of phenotypes are tuned to an optimal value, and declines monotonously as phenotypes are further away from this optimum. This model —called Fisher’s Geometric Model— generates a GxE interaction for fitness because the phenotypic optimum in the sink environment is distinct from that in the source environment [2]. The authors circumvent mathematical difficulties by developing an original approach based on tracking the deterministic dynamics of the cumulant generating function of the fitness distribution in the sink. They derive a number of important results on the dynamics of adaptation to the sink:

• From the point where immigration from the source to the sink starts, four phases of adaptation are observed. After a short transient phase (phase 1), a migration-selection balance is reached in the sink (phase 2). After a while, thanks to the immigration of rare adapted migrants and mutation in the sink, a small fraction of the sink population exhibits a close-to-optimal phenotype. This small adapted fraction grows in frequency and mean fitness rapidly increases in the sink (phase 3). Finally, the population settles around the sink optimum (phase 4) and, hurray, the sink is now a source!

• Interestingly, in this model the evolutionary dynamics do not depend on the immigration rate. In other words, adaptation will proceed at the same rate regardless of how many immigrants invade the sink. This is because the impact of immigration on adaptation depends on the rate of immigration relative to the sink density. This ratio is actually independent of immigration in a model where the sink is initially empty, migration from the sink back to the source is negligible and without density-dependence in the sink.

• In this model, mutation is a double-edged sword. Adapted phenotypes emerge from new mutations, and under this effect alone a higher mutation rate would translate into a shorter time to establishment in the sink. However, mutations may also have deleterious effects by displacing the phenotype away from the optimum. This mutation load will be greater when individuals need to simultaneously tune a large number of phenotypes. As a consequence of these two effects of mutations, time to establishment is minimal for an intermediate mutation rate. This result emerges from Fisher’s Geometric Model, but may hold more generally for biologically plausible fitness landscapes where mutations generates both beneficial (allowing adaptation to the sink) and deleterious genotypes.

• Lastly, in Fisher’s Geometric Model, the time to establishment increases superlinearly with harshness of the sink when the sink is too harsh, and establishment may occur only after a very long time. In these harsh sinks, the adapted genotypes are very few and increase very slowly in frequency, making the second phase of adaptation much longer. Thus, and as a direct consequence of Fisher’s Geometric Model, adding a “stepping stone” intermediate environment would allow faster adaptation to the extreme environment.

In conclusion, this theoretical work presents a method based on Fisher’s Geometric Model and the use of cumulant generating functions to resolve some aspects of adaptation to a sink environment. It generates a number of theoretical predictions for the adaptive colonisation of a sink by an asexual species with some standing genetic variation. It will be a fascinating task to examine whether these predictions hold in experimental evolution systems: will we observe the four phases of the dynamics of mean fitness in the sink environment? Will the rate of adaptation indeed be independent of the immigration rate? Is there an optimal rate of mutation for adaptation to the sink? Such critical tests of the theory will greatly improve our understanding of adaptation to novel environments.

References

[1] Lavigne, F., Martin, G., Anciaux, Y., Papaïx, J., and Roques, L. (2019). When sinks become sources: adaptive colonization in asexuals. bioRxiv, 433235, ver. 5 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/433235
[2] Martin, G., and Lenormand, T. (2006). A general multivariate extension of Fisher's geometrical model and the distribution of mutation fitness effects across species. Evolution, 60, 893-907. doi: 10.1111/j.0014-3820.2006.tb01169.x

Recent advances in next-generation sequencing are allowing us to uncover the evolution of sex chromosomes in non-model organisms. This study [1] represents an example of this application to birds of two Sylvioidea species from the genus Zosterops (commonly known as white-eyes). The study is exemplary in the amount and types of data generated and in the thoroughness of the analysis applied. Both male and female genomes were sequenced to allow the authors to identify sex-chromosome specific scaffolds. These data were augmented by generating the transcriptome (RNA-seq) data set. The findings after the analysis of these extensive data are intriguing: neoZ and neoW chromosome scaffolds and their breakpoints were identified. Novel sex chromosome formation appears to be accompanied by translocation events. The timing of formation of novel sex chromosomes was identified using molecular dating and appears to be relatively recent. Yet first signatures of distinct evolutionary patterns of sex chromosomes vs. autosomes could be already identified. These include the accumulation of transposable elements and changes in GC content. The changes in GC content could be explained by biased gene conversion and altered recombination landscape of the neo sex chromosomes. The authors also study divergence and diversity of genes located on the neo sex chromosomes. Here their findings appear to be surprising and need further exploration. The neoW chromosome already shows unique patterns of divergence and diversity at protein-coding genes as compared with genes on either neoZ or autosomes. In contrast, the genes on the neoZ chromosome do not display divergence or diversity patterns different from those for autosomes. This last observation is puzzling and I believe should be explored in further studies. Overall, this study significantly advances our knowledge of the early stages of sex chromosome evolution in vertebrates, provides an example of how such a study could be conducted in other non-model organisms, and provides several avenues for future work.

References

[1] Leroy T., Anselmetti A., Tilak M.K., Bérard S., Csukonyi L., Gabrielli M., Scornavacca C., Milá B., Thébaud C. and Nabholz B. (2019). A bird’s white-eye view on neo-sex chromosome evolution. bioRxiv, 505610, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/505610

Sex chromosomes are special in the genome because they are often highly differentiated over much of their lengths and marked by degenerative evolution of their gene content. Understanding why sex chromosomes differentiate requires deciphering the forces driving their recombination patterns. Suppression of recombination may be subject to selection, notably because of functional effects of locking together variation at different traits, as well as longer-term consequences of the inefficient purge of deleterious mutations, both of which may contribute to patterns of differentiation [1]. As an example, male and female functions may reveal intrinsic antagonisms over the optimal genotypes at certain genes or certain combinations of interacting genes. As a result, selection may favour the recruitment of rearrangements blocking recombination and maintaining the association of sex-antagonistic allele combinations with the sex-determining locus.
The hypothesis that sexually antagonistic selection might drive recombination suppression along the sex chromosomes is not new, but there are surprisingly few studies examining this empirically [1]. Support mainly comes from the study of guppy populations Poecilia reticulata in which the level of sexual dimorphism (notably due to male ornaments, subject to sexual selection) varies among populations, and was found to correlate with the length of the non-recombining region on the sex chromosome [2]. But the link is not always that clear. For instance in the fungus Microbotryum violaceum, the mating type loci is characterized by adjacent segments with recombination suppression, despite the near absence of functional differentiation between mating types [3].
In this study, Shearn and colleagues [4] explore the patterns of recombination suppression on the sex chromosomes of primates. X and Y chromosomes are strongly differentiated, except in a small region where they recombine with each other, the pseudoautosomal region (PAR). In the clade of apes and monkeys, including humans, large rearrangements have extended the non recombining region stepwise, eroding the PAR. Could this be driven by sexually antagonistic selection in a clade showing strong sexual differentiation?
To evaluate this idea, Shearn et al. have compared the structure of recombination in apes and monkeys to their sister clade with lower levels of sexual dimorphism, the lemurs and the lorises. If sexual antagonism was important in shaping recombination suppression, and assuming lower measures of sexual dimorphism reflect lower sexual antagonism [5], then lemurs and lorises would be predicted to show a shorter non-recombining region than apes and monkeys.
Lemurs and lorises were terra incognita in terms of genomic research on the sex chromosomes, so Shearn et al. have sequenced the genomes of males and females of different species. To assess whether sequences came from a recombining or non-recombining segment, they used coverage information in males vs females to identify sequences on the X whose copy on the Y is absent or too divergent to map, indicating long-term differentiation (absence of recombination). This approach reveals that the two lineages have undergone different recombination dynamics since they split from their common ancestor: regions which have undergone further structural rearrangements extending the non-recombining region in apes and monkeys, have continued to recombine normally in lemurs and lorises. Consistent with the prediction, macroevolutionary variation in the differentiation of males and females is indeed accompanied by variation in the size of the non-recombining region on the sex chromosome.
Sex chromosomes are excellent examples of how genomes are shaped by selection. By directly exploring recombination patterns on the sex chromosome across all extant primate groups, this study comes as a nice addition to the short series of empirical studies evaluating whether sexual antagonism may drive certain aspects of genome structure. The sexual selection causing sometimes spectacular morphological or behavioural differences between sexes in many animals may be the visible tip of the iceberg of all the antagonisms that characterise male vs. female functions generally [5]. Further research should bring insight into how different flavours or intensities of antagonistic selection can contribute to shape genome variation.

References

[1] Charlesworth D (2017) Evolution of recombination rates between sex chromosomes. Philosophical Transactions of the Royal Society B: Biological Sciences, 372, 20160456. https://doi.org/10.1098/rstb.2016.0456
[2] Wright AE, Darolti I, Bloch NI, Oostra V, Sandkam B, Buechel SD, Kolm N, Breden F, Vicoso B, Mank JE (2017) Convergent recombination suppression suggests role of sexual selection in guppy sex chromosome formation. Nature Communications, 8, 14251. https://doi.org/10.1038/ncomms14251
[3] Branco S, Badouin H, Vega RCR de la, Gouzy J, Carpentier F, Aguileta G, Siguenza S, Brandenburg J-T, Coelho MA, Hood ME, Giraud T (2017) Evolutionary strata on young mating-type chromosomes despite the lack of sexual antagonism. Proceedings of the National Academy of Sciences, 114, 7067–7072. https://doi.org/10.1073/pnas.1701658114
[4] Shearn R, Wright AE, Mousset S, Régis C, Penel S, Lemaitre J-F, Douay G, Crouau-Roy B, Lecompte E, Marais GAB (2020) Evolutionary stasis of the pseudoautosomal boundary in strepsirrhine primates. bioRxiv, 445072. https://doi.org/10.1101/445072
[5] Connallon T, Clark AG (2014) Evolutionary inevitability of sexual antagonism. Proceedings of the Royal Society B: Biological Sciences, 281, 20132123. https://doi.org/10.1098/rspb.2013.2123