Optimal sex ratios have been topic of extensive studies so far. Fisherian 1:1 proportions of males and females are known to be optimal in most (diploid) organisms, but many deviations from this golden rule are observed. These deviations not only attract a lot of attention from evolutionary biologists but also from population ecologists as they eventually determine long-term population growth. Because sex ratios are tightly linked to fitness, they can be under strong selection or plastic in response to changing demographic conditions. Hamilton [1] pointed out that an equality of the sex ratio breaks down when there is local competition for mates. Competition for mates can be considered as a special case of local resource competition. In short, this theory predicts females to adjust their offspring sex ratio conditional on cues indicating the level of local mate competition that their sons will experience. When cues indicate high levels of LMC mothers should invest more resources in the production of daughters to maximise their fitness, while offspring sex ratios should be closer to 50:50 when cues indicate low levels of LMC.
In isolated populations, Macke et al. [2] found sex ratio to evolve fast in response to changes in population sex-structure in the spider mite Tetranychus urticae. Spider mites are becoming top-models in evolutionary biology because of their easy housekeeping, fast generation times and well-studied genome [3]. The species is known to respond fast to changes in relatedness and kin-structure by changing its mating strategy [4], but also dispersal [5]. Sex ratio adjustments are likely mediated by differential investments in egg size, with small eggs possibly experiencing lower chances of fertilization, and thus to develop in haploid males [4].
Alison Duncan and colleagues [6] asked the question whether sex ratios change plastically in response to changes in the local population structure. They additionally questioned whether maternal effects could drive changes in sex-allocation of spider mite mothers. Indeed, theory predicts that if environmental changes are predictable across generations, intergenerational plasticity might be more adaptive than intragenerational plasticity [7]. Especially in spatially structured and highly dynamics populations, female spider mites may experience highly variable demographic conditions from one generation to another. During range expansions, spatial variation in local relatedness and inbreeding are documented to change and to impact eco-evolutionary trajectories as well (e.g. [8]).
Duncan et al. [6] specifically investigate whether the offspring sex ratio of T. urticae females changes in response to 1) the current number of females in the same patch, 2) the number of females in the patches of their mothers and 3) their relatedness to their mate. They surprisingly find the maternal environment to be more important than the actual experienced sex-ratio conditions. These insights thus show the maternal environment to be a reliable predictor of LMC experienced by grand-children. Maternal effects have been found to impact many traits, but this study is the first to convincingly demonstrate maternal effects in sex allocation. It therefore provides an alternative explanation of the apparent fast evolved responses under constant demographic conditions [2], and adds evidence to the importance of non-genetic trait changes for adaptation towards changing demographic and environmental conditions.

References

[1] Hamilton, W. D. (1967). Extraordinary Sex Ratios. Science, 156(3774), 477–488. doi: 10.1126/science.156.3774.477
[2] Macke, E., Magalhães, S., Bach, F., & Olivieri, I. (2011). Experimental evolution of reduced sex ratio adjustment under local mate competition. Science, 334(6059), 1127–1129. doi: 10.1126/science.1212177
[3] Grbić, M., Van Leeuwen, T., Clark, R. M., et al. (2011). The genome of Tetranychus urticae reveals herbivorous pest adaptations. Nature, 479(7374), 487–492. doi: 10.1038/nature10640
[4] Macke, E., Magalhães, S., Bach, F., & Olivieri, I. (2012). Sex-ratio adjustment in response to local mate competition is achieved through an alteration of egg size in a haplodiploid spider mite. Proceedings of the Royal Society B: Biological Sciences, 279(1747), 4634–4642. doi: 10.1098/rspb.2012.1598
[5] Bitume, E. V., Bonte, D., Ronce, O., Bach, F., Flaven, E., Olivieri, I., & Nieberding, C. M. (2013). Density and genetic relatedness increase dispersal distance in a subsocial organism. Ecology Letters, 16(4), 430–437. doi: 10.1111/ele.12057
[6] Duncan, A., Marinosci, C., Devaux, C., Lefèvre, S., Magalhães, S., Griffin, J., Valente, A., Ronce, O., Olivieri, I. (2018). Transgenerational cues about local mate competition affect offspring sex ratios in the spider mite Tetranychus urticae. BioRxiv, 240127, ver. 3. doi: 10.1101/240127
[7] Petegem, K. V., Moerman, F., Dahirel, M., Fronhofer, E. A., Vandegehuchte, M. L., Leeuwen, T. V., Wybouw, N., Stoks, R., Bonte, D. (2018). Kin competition accelerates experimental range expansion in an arthropod herbivore. Ecology Letters, 21(2), 225–234. doi: 10.1111/ele.12887
[8] Marshall, D. J., & Uller, T. (2007). When is a maternal effect adaptive? Oikos, 116(12), 1957–1963. doi: 10.1111/j.2007.0030-1299.16203.x

Duplications, mutations and recombination may be considered the main sources of genomic variation and evolution. In addition, sexual recombination is essential in purging deleterious mutations and allowing advantageous allelic combinations to occur (Glémin et al. 2019). However, in parthenogenetic asexual organisms, variation cannot be explained by sexual recombination, and other mechanisms must account for it. Although it is known that transposable elements (TE) may influence on genome structure and gene expression patterns, their role as a primary source of genomic variation and rapid adaptability has received less attention. An important role of TE on adaptive genome evolution has been documented for fungal phytopathogens (Faino et al. 2016), suggesting that TE activity might explain the evolutionary dynamics of this type of organisms.
The phytopathogen nematode Meloidogyne incognita is one of the worst agricultural pests in warm climates (Savary et al. 2019). This species, as well as other root-knot nematodes (RKN), shows a wide geographical distribution range infecting diverse groups of plants. Although allopolyploidy may have played an important role on the wide adaptation of this phytopathogen, it may not explain by itself the rapid changes required to overcome plant resistance in a few generations. Paradoxically, M. incognita reproduces asexually via mitotic parthenogenesis (Trudgill and Blok 2001; Castagnone-Sereno and Danchin 2014) and only few single nucleotide variations were identified between different host races isolates (Koutsovoulos et al. 2020). Therefore, this is an interesting model to explore other sources of genomic variation such the TE activity and its role on the success and adaptability of this phytopathogen.
To address these questions, Kozlowski et al. (2020) estimated the TE mobility across 12 geographical isolates that presented phenotypic variations in Meloidogyne incognita, concluding that recent activity of TE in both genic and regulatory regions might have given rise to relevant functional differences between genomes. This was the first estimation of TE activity as a mechanism probably involved in genome plasticity of this root-knot nematode. This study also shed light on evolutionary mechanisms of asexual organisms with an allopolyploid origin. These authors re-annotated the 185 Mb triploid genome of M. incognita for TE content analysis using stringent filters (Kozlowski 2020a), and estimated activity by their distribution using a population genomics approach including isolates from different crops and locations. Canonical TE represented around 4.7% of the M. incognita genome of which mostly correspond to TIR (Terminal Inverted Repeats) and MITEs (Miniature Inverted repeat Transposable Elements) followed by Maverick DNA transposons and LTR (Long Terminal Repeats) retrotransposons. The result that most TE found were represented by DNA transposons is similar to the previous studies with the nematode species model Caenorhabditis elegans (Bessereau 2006; Kozlowski 2020b) and other nematodes as well. Canonical TE annotations were highly similar to their consensus sequences containing transposition machinery when TE are autonomous, whereas no genes involved in transposition were found in non-autonomous ones. These findings suggest recent activity of TE in the M. incognita genome. Other relevant result was the significant variation in TE presence frequencies found in more than 3,500 loci across isolates, following a bimodal distribution within isolates. However, variation in TE frequencies was low to moderate between isolates recapitulating the phylogenetic signal of isolates DNA sequences polymorphisms. A detailed analysis of TE frequencies across isolates allowed identifying polymorphic TE loci, some of which might be neo-insertions mostly of TIRs and MITEs (Kozlowski 2020c). Interestingly, the two thirds of the fixed neo-insertions were located in coding regions or in regulatory regions impacting expression of specific genes in M. incognita. Future research on proteomics is needed to evaluate the functional impact that these insertions have on adaptive evolution in M. incognita. In this line, this pioneer research of Kozlowski et al. (2020) is a first step that is also relevant to remark the role that allopolyploidy and reproduction have had on shaping nematode genomes.

References

[1] Bessereau J-L. 2006. Transposons in C. elegans. WormBook. 10.1895/wormbook.1.70.1
[2] Castagnone-Sereno P, Danchin EGJ. 2014. Parasitic success without sex - the nematode experience. J. Evol. Biol. 27:1323-1333. 10.1111/jeb.12337
[3] Faino L, Seidl MF, Shi-Kunne X, Pauper M, Berg GCM van den, Wittenberg AHJ, Thomma BPHJ. 2016. Transposons passively and actively contribute to evolution of the two-speed genome of a fungal pathogen. Genome Res. 26:1091-1100. 10.1101/gr.204974.116
[4] Glémin S, François CM, Galtier N. 2019. Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species. In: Anisimova M, editor. Evolutionary Genomics: Statistical and Computational Methods. Methods in Molecular Biology. New York, NY: Springer. p. 331-369. 10.1007/978-1-4939-9074-0_11
[5] Koutsovoulos GD, Marques E, Arguel M-J, Duret L, Machado ACZ, Carneiro RMDG, Kozlowski DK, Bailly-Bechet M, Castagnone-Sereno P, Albuquerque EVS, et al. 2020. Population genomics supports clonal reproduction and multiple independent gains and losses of parasitic abilities in the most devastating nematode pest. Evol. Appl. 13:442-457. 10.1111/eva.12881
[6] Kozlowski D. 2020a. Transposable Elements prediction and annotation in the M. incognita genome. Portail Data INRAE. 10.15454/EPTDOS
[7] Kozlowski D. 2020b. Transposable Elements prediction and annotation in the C. elegans genome. Portail Data INRAE. 10.15454/LQCIW0
[8] Kozlowski D. 2020c. TE polymorphisms detection and analysis with PopoolationTE2. Portail Data INRAE. 10.15454/EWJCT8
[9] Kozlowski DK, Hassanaly-Goulamhoussen R, Da Rocha M, Koutsovoulos GD, Bailly-Bechet M, Danchin EG (2020) Transposable Elements are an evolutionary force shaping genomic plasticity in the parthenogenetic root-knot nematode Meloidogyne incognita. bioRxiv, 2020.04.30.069948, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. 10.1101/2020.04.30.069948
[10] Savary S, Willocquet L, Pethybridge SJ, Esker P, McRoberts N, Nelson A. 2019. The global burden of pathogens and pests on major food crops. Nat. Ecol. Evol. 3:430-439. 10.1038/s41559-018-0793-y
[11] Trudgill DL, Blok VC. 2001. Apomictic, polyphagous root-knot nematodes: exceptionally successful and damaging biotrophic root pathogens. Annu Rev Phytopathol 39:53-77. 10.1146/annurev.phyto.39.1.53

Local adaptation, the higher fitness a population achieves in its local “home” environment relative to other environments is a crucial phase in the divergence of populations, and as such both generates and maintains diversity. Local adaptation is enhanced by selection and genetic variation in the relevant traits, and decreased by gene flow and genetic drift.

Demonstrating local adaptation is laborious, and is typically done with a reciprocal transplant design [1], documenting repeated geographic clines [e.g. 2, 3] also provides strong evidence of local adaptation. Even when well documented, it is often unknown which aspects of the environment impose selection. Indeed, differences in environment between different sites that are measured during studies of local adaptation explain little of the variance in the degree of local adaptation [4]. This poses a problem to population management. Given climate change and habitat destruction, understanding the environmental drivers of local adaptation can be crucially important to conducting successful assisted migration or targeted gene flow.

In this manuscript, Macdonald et al. [5] propose a means of identifying which aspects of the environment select for local adaptation without conducting a reciprocal transplant experiment. The idea is that the strength of relationships between traits and environmental variables that are due to plastic responses to the environment will not be influenced by gene flow, but the strength of trait-environment relationships that are due to local adaptation should decrease with gene flow. This then can be used to reduce the somewhat arbitrary list of environmental variables on which data are available down to a targeted list more likely to drive local adaptation in specific traits. To perform such an analysis requires three things: 1) measurements of traits of interest in a species across locations, 2) an estimate of gene flow between locations, which can be replaced with a biologically meaningful estimate of how well connected those locations are from the point of view of the study species, and 3) data on climate and other environmental variables from across a species’ range, many of which are available on line.

Macdonald et al. [5] demonstrate their approach using a skink (Lampropholis coggeri). They collected morphological and physiological data on individuals from multiple populations. They estimated connectivity among those locations using information on habitat suitability and dispersal potential [6], and gleaned climatic data from available databases and the literature. They find that two physiological traits, the critical minimum and maximum temperatures, show the strongest signs of local adaptation, specifically local adaptation to annual mean precipitation, precipitation of the driest quarter, and minimum annual temperature. These are then aspects of skink phenotype and skink habitats that could be explored further, or could be used to provide background information if migration efforts, for example for genetic rescue [7] were initiated. The approach laid out has the potential to spark a novel genre of research on local adaptation. It its simplest form, knowing that local adaptation is eroded by gene flow, it is intuitive to consider that if connectivity reduces the strength of the relationship between an environmental variable and a trait, that the trait might be involved in local adaptation. The approach is less intuitive than that, however – it relies not connectivity per-se, but the interaction between connectivity and different environmental variables and how that interaction alters trait-environment relationships. The authors lay out a number of useful caveats and potential areas that could use further development. It will be interesting to see how the community of evolutionary biologists responds.

References

[1] Blanquart F, Kaltz O, Nuismer SL and Gandon S. 2013. A practical guide to measuring local adaptation. Ecology Letters, 16: 1195-1205. doi: 10.1111/ele.12150

[2] Huey RB, Gilchrist GW, Carlson ML, Berrigan D and Serra L. 2000. Rapid evolution of a geographic cline in size in an introduced fly. Science, 287: 308-309. doi: 10.1126/science.287.5451.308

[3] Milesi P, Lenormand T, Lagneau C, Weill M and Labbé P. 2016. Relating fitness to long-term environmental variations in natura. Molecular Ecology, 25: 5483-5499. doi: 10.1111/mec.13855

[4] Hereford, J. 2009. A quantitative survey of local adaptation and fitness trade-offs. The American Naturalist 173: 579-588. doi: 10.1086/597611

[5] Macdonald SL, Llewelyn J and Phillips BL. 2017. Using connectivity to identify climatic drivers of local adaptation. bioRxiv, ver. 4 of October 4, 2017. doi: 10.1101/145169

[6] Macdonald SL, Llewelyn J, Moritz C and Phillips BL. 2017. Peripheral isolates as sources of adaptive diversity under climate change. Frontiers in Ecology and Evolution, 5:88. doi: 10.3389/fevo.2017.00088

[7] Whiteley AR, Fitzpatrick SW, Funk WC and Tallmon DA. 2015. Genetic rescue to the rescue. Trends in Ecology & Evolution, 30: 42-49. doi: 10.1016/j.tree.2014.10.009

Many organisms across the Tree of life have the ability to produce seeds, eggs, cysts, or spores, that can remain dormant for several generations before hatching. This widespread adaptive trait in bacteria, fungi, plants and animals, has a significant impact on the ecology, population dynamics and population genetics of species that express it (Evans and Dennehy 2005).

In population genetics, and despite the recognition of its evolutionary importance in many empirical studies, few theoretical models have been developed to characterize the evolutionary consequences of this trait on the level and distribution of neutral genetic diversity (see, e.g., Kaj et al. 2001; Vitalis et al. 2004), and also on the dynamics of selected alleles (see, e.g., Živković and Tellier 2018). However, due to the complexity of the interactions between evolutionary forces in the presence of dormancy, the fate of selected mutations in their genomic environment is not yet fully understood, even from the most recently developed models.

In a comprehensive article, Korfmann et al. (2023) aim to fill this gap by investigating the effect of germ banking on the probability of (and time to) fixation of beneficial mutations, as well as on the shape of the selective sweep in their vicinity. To this end, Korfmann et al. (2023) developed and released their own forward-in-time simulator of genome-wide data, including neutral and selected polymorphisms, that makes use of Kelleher et al.’s (2018) tree sequence toolkit to keep track of gene genealogies.

The originality of Korfmann et al.’s (2023) study is to provide new quantitative results for the effect of dormancy on the time to fixation of positively selected mutations, the shape of the genomic landscape in the vicinity of these mutations, and the temporal dynamics of selective sweeps. Their major finding is the prediction that germ banking creates narrower signatures of sweeps around positively selected sites, which are detectable for increased periods of time (as compared to a standard Wright-Fisher population).

The availability of Korfmann et al.’s (2023) code will allow a wider range of parameter values to be explored, to extend their results to the particularities of the biology of many species. However, as they chose to extend the haploid coalescent model of Kaj et al. (2001), further development is needed to confirm the robustness of their results with a more realistic diploid model of seed dormancy.

REFERENCES

Evans, M. E. K., and J. J. Dennehy (2005) Germ banking: bet-hedging and variable release from egg and seed dormancy. The Quarterly Review of Biology, 80(4): 431-451. https://doi.org/10.1086/498282

Kaj, I., S. Krone, and M. Lascoux (2001) Coalescent theory for seed bank models. Journal of Applied Probability, 38(2): 285-300. https://doi.org/10.1239/jap/996986745

Kelleher, J., K. R. Thornton, J. Ashander, and P. L. Ralph (2018) Efficient pedigree recording for fast population genetics simulation. PLoS Computational Biology, 14(11): e1006581. https://doi.org/10.1371/journal.pcbi.1006581

Korfmann, K., D. Abu Awad, and A. Tellier (2023) Weak seed banks influence the signature and detectability of selective sweeps. bioRxiv, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.04.26.489499

Vitalis, R., S. Glémin, and I. Olivieri (2004) When genes go to sleep: the population genetic consequences of seed dormancy and monocarpic perenniality. American Naturalist, 163(2): 295-311. https://doi.org/10.1086/381041

Živković, D., and A. Tellier (2018). All but sleeping? Consequences of soil seed banks on neutral and selective diversity in plant species. Mathematical Modelling in Plant Biology, 195-212. https://doi.org/10.1007/978-3-319-99070-5_10

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