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

Interactions between species are a driving force in evolution. Many organisms host symbiotic partners that live all or part of their life in or on their host. Whether they are mutualistic or parasitic, these symbiotic associations impose strong selective pressures on both partners and affect their evolutionary trajectories. In fine, they can have a significant impact on the diversification patterns of both host and symbiont lineages, with symbiotic lineages sometimes speciating simultaneously with their hosts and/or switching from one host species to another. Long-term associations between species can also result in gene transfers between the involved organisms. Those lateral gene transfers are a source of ecological innovation but can obscure the phylogenetic signals and render the process of phylogenetic reconstructions complex (Lerat et al. 2003).

Methods known as reconciliations explore similarities and differences between phylogenetic trees. They have been widely used to both compare the diversification patterns of hosts and symbionts and identify lateral gene transfers between species. Though the reconciliation approaches used in host/ symbiont and species/ gene phylogenetic studies are identical, they are always applied separately to solve either molecular evolution questions or investigate the evolution of ecological interactions. However, the two questions are often intimately linked and the current interest in multi-level systems (e.g. the holobiont concept) calls for a unique model that will take into account three-level nested organization (gene/symbiont/ host) where both symbiont and genes can transfer among hosts. 

Here Menet and collaborators (2023) provide such a model to produce three-level reconciliations. In order to do so, they extend the two-level reconciliation model implemented in “ALE” software (Szöllősi et al. 2013), one of the most used and proven reconciliation methods. Briefly, given a symbiont gene tree, a symbiont tree and a host tree, as in previous reconciliation models, the symbiont tree is mapped onto the host tree by mixing three types of events: Duplication, Transfer or Loss (DTL), with a possibility of the symbiont evolving temporarily outside the host phylogeny (in a “ghost” host lineage). The gene tree evolves similarly inside the symbiont tree, but horizontal transfers are constrained to symbionts co-occurring within the same host. Joint reconciliation scenarios are reconstructed and DTL event rates and likelihoods are estimated according to the model. As a nice addition, the authors propose a method to infer the symbiont phylogeny through amalgamation from gene trees and a host tree.

The authors then explore the diverse possibilities offered by this method by testing it on both simulated datasets and biological datasets in order to check whether considering three nested levels is worthwhile. They convincingly show that three-level reconciliation has a better capacity to retrieve the symbiont donors and receivers of horizontal gene transfers, probably because transfers are constrained by additional elements relevant to the biological systems. Using, aphids, their obligate endosymbionts, and the symbiont genes involved in their nutritional functions, they identify horizontal gene transfers between aphid symbionts that are missed by two-level reconciliations but detected by expertise (Manzano-Marín et al. 2020). The other dataset presented here is on the human pathogen Helicobacter pylori, which history is supposed to reflect human migration. They use more than 1000 H. pylori gene families, and four populations, and use likelihood computations to compare different hypotheses on the diversification of the host.

In summary, this study is a proof-of-concept of a 3-level reconciliation, where the authors manage to convey the applicability of their framework to many biological systems. Reported complexities, confirmed by reported running times, show that the method is computationally efficient. Without a doubt, the tool presented here will be very useful to evolutionary biologists who want to investigate multi-scale cophylogenies and it will move forward the study of associations between host and symbionts when symbiont genomic data are available.

REFERENCES

Lerat, E., Daubin, V., & Moran, N. A. (2003). From gene trees to organismal phylogeny in prokaryotes: the case of the γ-Proteobacteria. PLoS biology, 1(1), e19.
https://doi.org/10.1371/journal.pbio.0000019
 
Manzano-Marın, A., Coeur d'acier, A., Clamens, A. L., Orvain, C., Cruaud, C., Barbe, V., & Jousselin, E. (2020). Serial horizontal transfer of vitamin-biosynthetic genes enables the establishment of new nutritional symbionts in aphids' di-symbiotic systems. The ISME Journal, 14(1), 259-273.
https://doi.org/10.1038/s41396-019-0533-6

Menet H, Trung AN, Daubin V, Tannier E (2023) Host-symbiont-gene phylogenetic reconciliation. bioRxiv, 2022.07.01.498457, ver. 2 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.07.01.498457

Szöllősi, G. J., Rosikiewicz, W., Boussau, B., Tannier, E., & Daubin, V. (2013). Efficient exploration of the space of reconciled gene trees. Systematic biology, 62(6), 901-912.
https://doi.org/10.1093/sysbio/syt054

The total amount of DNA utilized to store hereditary information varies immensely among eukaryotic organisms. Single copy genome sizes – disregarding differences due to ploidy - differ by more than three orders of magnitude ranging from a few million nucleotides (Mb) to hundreds of billions (Gb). With the ever-increasing availability of fully sequenced genomes we now know that most of the difference is due either to whole genome duplication or to variation in the abundance of repetitive elements. Regarding repetitive elements, the evolutionary forces underlying the large variation 'allowing' more or less elements in a genome remain largely elusive. A tentative correlation between an organism's complexity (however this may be adequately measured) and genome size, the so called C-value paradox [1], has long been dismissed. Studies testing for selection on secondary phenotypic effects associated with genome size (cell size, metabolic rates, nutrient availability) have yielded mixed results. Nonadaptive theories capitalizing on a role of deleterious insertion-deletion mutations and genetic drift as the main drivers have likewise received mixed support [2-3]. Overall, most evidence was derived from analyses across broad taxonomical scales [4-6].

Lefébure and colleagues [7] take a different approach. They confine their considerations to a homogeneous, restricted taxonomical group, isopod crustaceans of the superfamily Aselloidea. This taxonomic focus allows the authors to circumvent many of the confounding factors such as phylogenetic inertia, life history divergence and mutation rate variation that tend to trouble analyses across broad taxonomic timescales. Another important feature of the chosen system is the evolutionary independent transition of habitat use that has occurred at least 11 times. One group of species inhabits subterranean ecosystems (groundwater), another group thrives on surface water. Populations of the former live in low-energy habitats and are expected to be outnumbered by their surface dwelling relatives. Interestingly – and a precondition for the study - the groundwater species have significantly larger genomes (up to 137%). With this unique set-up, the authors are able to investigate the link between genome size and evolutionary forces related to a proxy of long-term population size by removing many of the confounding factors a priori.

Upfront, we learn that the dN/dS ratio is higher in the groundwater species. This may either suggest prevalent positive selection or lower efficacy of purifying selection (relaxed constraint) in the group of species in which population sizes are expected to be low. Using a series of population genetic analyses the authors provide compelling evidence for the latter. Analyses are carefully conducted and include models for estimating the intensity and frequency of purifying and positive selection, the DoS (direction of selection) and α statistic. Next the authors also exclude the possibility that increased dN/dS of the subterranean groundwater species may be due to nonfunctionalization, which may result from the subterranean lifestyle.

Overall, these analyses suggest relaxed constraint in smaller populations as the most plausible alternative to explain increased dN/dS ratios. In addition to the efficacy of selection, the authors estimate the timing of the ecological transition under the rationale that the amount of time a species may have been exposed to the subterranean habitat may reflect long term population sizes. To calibrate the 'colonization clock' they apply a neat trick based on the degree of degeneration of the opsin gene (as vision tends to get lost in these habitats). When finally testing which parameters may explain differences in genome size all factors – ecological status, selection efficiency as measured by dN/dS and colonization time - turned out to be significant predictors. Direct estimates of the short term effective population size Ne from polymorphism data, however, did not correlate with genome size. Ruling out the effect of other co-variates such as body size and growth rate the authors conclude that genome size was overall best predicted by long-term population size change upon habitat shift. In that the authors provide convincing evidence that the increase in genome size is linked to a decrease in long-term reduction of selection efficiency of subterranean species. Assuming a bias for insertion mutations over deletion mutations (which is usually the case in eukaryotes) this result is in agreement with the theory of mutational hazard [4-6]. This theory proposed by Michael Lynch postulates that the accumulation of non-functional DNA has a weak deleterious effect that can only be efficiently opposed by natural selection in species with high Ne.

In conclusion, Lefébure and colleagues provide novel and welcome evidence supporting a 'neutralist' hypothesis of genome size evolution without the need to invoke an adaptive component. Methodologically, the study cautions against the common use of polymorphism-based estimates of Ne which are often obfuscated by transitory demographic change. Instead, alternative measures of selection efficacy linked to long-term population size may serve as better predictors of genome size. We hope that this study will stimulate additional work testing the link between Ne and genome size variation in other taxonomical groups [8-9]. Using genome sequences instead of the transcriptome approach applied here may concomitantly further our understanding of the molecular mechanisms underlying genome size change.

References

[1] Thomas, CA Jr. 1971. The genetic organization of chromosomes. Annual Review of Genetics 5: 237–256. doi: 10.1146/annurev.ge.05.120171.001321

[2] Ågren JA, Greiner S, Johnson MTJ, Wright SI. 2015. No evidence that sex and transposable elements drive genome size variation in evening primroses. Evolution 69: 1053–1062. doi: 10.1111/evo.12627

[3] Bast J, Schaefer I, Schwander T, Maraun M, Scheu S, Kraaijeveld K. 2016. No accumulation of transposable elements in asexual arthropods. Molecular Biology and Evolution 33: 697–706. doi: 10.1093/molbev/msv261

[4] Lynch M. 2007. The Origins of Genome Architecture. Sinauer Associates.

[5] Lynch M, Bobay LM, Catania F, Gout JF, Rho M. 2011. The repatterning of eukaryotic genomes by random genetic drift. Annual Review of Genomics and Human Genetics 12: 347–366. doi: 10.1146/annurev-genom-082410-101412

[6] Lynch M, Conery JS. 2003. The origins of genome complexity. Science 302: 1401–1404. doi: 10.1126/science.1089370

[7] Lefébure T, Morvan C, Malard F, François C, Konecny-Dupré L, Guéguen L, Weiss-Gayet M, Seguin-Orlando A, Ermini L, Der Sarkissian C, Charrier NP, Eme D, Mermillod-Blondin F, Duret L, Vieira C, Orlando L, and Douady CJ. 2017. Less effective selection leads to larger genomes. Genome Research 27: 1016-1028. doi: 10.1101/gr.212589.116

[8] Lower SS, Johnston JS, Stanger-Hall KF, Hjelmen CE, Hanrahan SJ, Korunes K, Hall D. 2017. Genome size in North American fireflies: Substantial variation likely driven by neutral processes. Genome Biolology and Evolution 9: 1499–1512. doi: 10.1093/gbe/evx097

[9] Sessegolo C, Burlet N, Haudry A. 2016. Strong phylogenetic inertia on genome size and transposable element content among 26 species of flies. Biology Letters 12: 20160407. doi: 10.1098/rsbl.2016.0407

Exciting times are afoot for those of us interested in the ecology and evolution of invasive populations. Recent years have seen evolutionary process woven firmly into our understanding of invasions (Miller et al. 2020). This integration has inspired a welter of empirical and theoretical work. We have moved from field observations and verbal models to replicate experiments and sophisticated mathematical models. Progress has been rapid, and we have seen science at its best; an intimate discussion between theory and data.
An area currently under very active development is our understanding of pushed invasions. Here a population spreads through space driven, not by dispersal and growth originating at the leading tip of the invasion, but by dispersal and growth originating deeper in the bulk of the population. These pushed invasions may be quite common – they result when per capita growth and dispersal rates are higher in the bulk of the wave than at the leading tip. They result from a range of well-known phenomena, including Allee effects and density-dependent dispersal (Gandhi et al. 2016; Bîrzu et al. 2019). Pushed invasions travel faster than we would expect given growth and dispersal rates on the leading tip, and they lose genetic diversity more slowly than classical pulled invasions (Roques et al. 2012; Haond et al. 2018; Bîrzu et al. 2019).
Well… in theory, anyway. The theory on pushed waves has momentarily streaked ahead of the empirical work, because empirical systems for studying pushed invasions are rare (though see Gandhi et al. 2016; Gandhi, Korolev, and Gore 2019). In this paper, Dahirel and colleagues (2020) make the argument that we may be able to generate pushed invasions in laboratory systems simply by reducing the connectedness of our experimental landscapes. If true, we might have a simple tool for turning many of our established experimental systems into systems for studying pushed dynamics.
It’s a nice idea, and the paper goes to careful lengths to explore the possibility in their lab system (a parasitoid wasp, Trichogramma). They run experiments on replicate wasp populations comparing strongly- v poorly-connected arrays, and estimate the resulting invasion speeds and rate of diversity loss. They also build a simulation model of the system, allowing them to explore in-silico a range of possible processes underlying their results.
As well as developing these parallel systems, Dahirel and colleagues (2020) go to careful lengths to develop statistical analyses that allow inference on key parameters, and they apply these analyses to both the experimental and simulation data. They have been motivated to apply methods that might be used in both laboratory and field settings to help classify invasions.
Ultimately, they found reasonable evidence that their poorly-connected habitat did induce a pushed dynamic. Their poorly connected invasions travelled faster than they should have if they were pulled, they lost diversity more slowly than the highly connected habitat, and replicates with a higher carrying capacity tended to have higher invasion speeds. All in line with expectations of a pushed dynamic. Interestingly, however, their simulation results suggest that they probably got this perfect result for unexpected reasons. The strong hint is that their poorly-connected habitat induced density dependent dispersal in the wasps. Without this effect, their simulations suggest they should have seen diversity decreasing much more rapidly than it did.
There is a nuanced, thoughtful, and carefully argued discussion about all this in the paper, and it is worth reading. There is much of value in this paper. Theirs is not a perfect empirical system in which all the model assumptions are met and in which huge population sizes make stochastic effects negligible. Here is a system one step closer to the messy reality of biology. The struggle to align this system with new theory has been worth the effort. Not only does it give us hope that we might usefully be able to discriminate between classes of invasions using real-world data, but it hints at a rule that Tolstoy might have expressed this way: all pulled invasions are alike, each pushed invasion is pushed in its own way.

References

Bîrzu, G., Matin, S., Hallatschek, O., and Korolev, K. S. (2019). Genetic drift in range expansions is very sensitive to density dependence in dispersal and growth. Ecology Letters, 22(11), 1817-1827. doi: https://doi.org/10.1111/ele.13364
Dahirel, M., Bertin, A., Haond, M., Blin, A., Lombaert, E., Calcagno, V., Fellous, S., Mailleret, L., Malausa, T., and Vercken, E. (2020). Shifts from pulled to pushed range expansions caused by reduction of landscape connectivity. bioRxiv, 2020.05.13.092775, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. https://doi.org/10.1101/2020.05.13.092775
Gandhi, S. R., Korolev, K. S., and Gore, J. (2019). Cooperation mitigates diversity loss in a spatially expanding microbial population. Proceedings of the National Academy of Sciences, 116(47), 23582-23587. doi: https://doi.org/10.1073/pnas.1910075116
Gandhi, S. R., Yurtsev, E. A., Korolev, K. S., and Gore, J. (2016). Range expansions transition from pulled to pushed waves as growth becomes more cooperative in an experimental microbial population. Proceedings of the National Academy of Sciences, 113(25), 6922-6927. doi: https://doi.org/10.1073/pnas.1521056113
Haond, M., Morel-Journel, T., Lombaert, E., Vercken, E., Mailleret, L. and Roques, L. (2018). When higher carrying capacities lead to faster propagation (2018), bioRxiv, 307322, ver. 4 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/307322
Miller et al. (2020). Eco‐evolutionary dynamics of range expansion. Ecology, 101(10), e03139. doi: https://doi.org/10.1002/ecy.3139
Roques, L., Garnier, J., Hamel, F., and Klein, E. K. (2012). Allee effect promotes diversity in traveling waves of colonization. Proceedings of the National Academy of Sciences, 109(23), 8828-8833. doi: https://doi.org/10.1073/pnas.1201695109

It is a classical result in the virulence evolution literature that treatments decreasing parasite replication within the host should select for higher replication rates, thus driving increased levels of virulence if the two are correlated. There is some evidence for this in vitro but very little in the field. HIV infections in humans offer a unique opportunity to go beyond the simple predictions that treatments should favour more virulent strains because many details of this host-parasite system are known, especially the link between set-point virus load, transmission rate and virulence.

To tackle this question, Herbeck et al. [1] used a detailed individual-based model. This is original because it allows them to integrate existing knowledge from the epidemiology and evolution of HIV (e.g. recent estimates of the ‘heritability’ of set-point virus load from one infection to the next). This detailed model allows them to formulate predictions regarding the effect of different treatment policies; especially regarding the current policy switch away from treatment initiation based on CD4 counts towards universal treatment.

The results show that, perhaps as expected from the theory, treatments based on the level of remaining host target cells (CD4 T cells) do not affect virulence evolution because they do not strongly affect the virulence level that maximizes HIV’s transmission potential. However, early treatments can lead to moderate increase in virulence within several years if coverage is high enough. These results seem quite robust to variation of all the parameters in realistic ranges.

The great step forward in this model is the ability to obtain quantitative prediction regarding how a virus may evolve in response to public health policies. Here the main conclusion is that given our current knowledge in HIV biology, the risk of virulence evolution is perhaps more limited than expected from a direct application of virulence evolution model. Interestingly, the authors also conclude that recently observed increased in HIV virulence [2-3] cannot be explained by the impact of antiretroviral therapy alone; which raises the question about the main mechanism behind this increase. Finally, the authors make the interesting suggestion that “changing virulence is amenable to being monitored alongside transmitted drug resistance in sentinel surveillance”.

References

[1] Herbeck JT, Mittler JE, Gottlieb GS, Goodreau SM, Murphy JT, Cori A, Pickles M, Fraser C. 2016. Evolution of HIV virulence in response to widespread scale up of antiretroviral therapy: a modeling study. Virus Evolution 2:vew028. doi: 10.1093/ve/vew028

[2] Herbeck JT, Müller V, Maust BS, Ledergerber B, Torti C, et al. 2012. Is the virulence of HIV changing? A meta-analysis of trends in prognostic markers of HIV disease progression and transmission. AIDS 26:193-205. doi: 10.1097/QAD.0b013e32834db418

[3] Pantazis N, Porter K, Costagliola D, De Luca A, Ghosn J, et al. 2014. Temporal trends in prognostic markers of HIV-1 virulence and transmissibility: an observational cohort study. Lancet HIV 1:e119-26. doi: 10.1016/s2352-3018(14)00002-2