Not all genetic elements composing genomes are there for the benefit of their carrier. Many have no consequences on fitness, or too mild ones to be eliminated by selection, and thus stem from neutral processes. Many others are indeed the product of selection, but one acting at a different level, increasing the fitness of some elements of the genome only, at the expense of the “organism” as a whole. These can be called selfish genetic elements, and come into a wide variety of flavours [1], illustrating many possible means to cheat with “fair” reproductive processes such as meiosis, and thus get overrepresented in the offspring of their hosts. Producing copies of itself through transposition is one such strategy; a very successful one indeed, explaining a large part of the genomic content of many organisms. Killing non carrier gametes following meiosis in heterozygous carriers is another one. Less know and less common is the ability of some elements to turn heterozygous carriers into homozygous ones, that will thus transmit the selfish elements to all offspring instead of half. This is achieved by nucleic sequences encoding so-called “Homing endonucleases” (HEs). These proteins tend to induce double strand breaks of DNA specifically in regions homologous to their own insertion sites. The recombination machinery is such that the intact homologous region, that is, the one carrying the HE sequence, is then used as a template for the reparation of the break, resulting in the effective conversion of a non-carrier allele into a carrier allele. Such elements can also occur in the mitochondrial genomes of organisms where mitochondria are not strictly transmitted by one parent only, offering mitochondrial HEs some opportunities for “homing” into new non carrier genomes. This is the case in yeasts, where HEs were first reported [2,3].
In this new study, based on genomic experimental data from the fungal maize pathogen Ustilago maydis, Julien Dutheil and colleagues [4] document one possible evolutionary pathway for which little evidence existed before: the passage of a mitochondrial HE into the nuclear genome. The GC content of this region leaves little doubt on its mitochondrial origin, and homologs can indeed be found in the mitochondrial genomes of close relatives. Strangely enough, U. maydis itself does not appear to carry this selfish element in its own mitochondria, suggesting it may have been acquired from a different species, or be subject to a sufficiently rapid turnover to have been recently lost.
Many elements of the story uncovered by this study remain mysterious. How, in the first place, was this HE gene inserted in a nuclear genomic region that shows no apparent homology with its original insertion site, making typical “homing” a not-so-likely explanation? This question may in fact be generalised to many HE systems: is the first insertion into a homing site always the product of a typical homing event, which implies the presence of an homologous template DNA fragment, or can HE genes insert through other means? But then, why specifically in regions that would be targeted by the nuclease they encode? What is the evolutionary fate of this newly inserted element? The new gene may well be on its way to pseudogenisation, as suggested by the truncation of its upper part, precluding its functioning as a HE, and the lack of evidence of selective constraints through dN/dS analysis; but the mutation generated by the insertion event may have phenotypic implications, possibly through the partial truncation of another gene, encoding a helicase. How old is this insertion? The fact that it has accumulated some mutations makes a very recent event rather unlikely, but this insertion has been detected in only one isolate of U. maydis, suggesting it is not so frequent in natural populations.
Whatever the answers to these open questions, that will hopefully be addressed by further work on this system, the present study has revealed that horizontal transmission enlarges the scope of possible evolutionary consequences of HE genes, that may move not only between mitochondrial genomes, but also occasionally into a nucleus.

References

[1] Burt, A., and Trivers, R. (2006). Genes in Conflict: The Biology of Selfish Genetic Elements. Belknap Press.
[2] Coen, D., Deutch, J., Netter, P., Petrochillo, E., and Slonimski, P. (1970). Mitochondrial genetics. I. Methodology and phenomenology. Symposia of the Society for Experimental Biology, 24, 449-496.
[3] Colleaux, L., D’Auriol, L., Betermier, M., Cottarel, G., Jacquier, A., Galibert, F., and Dujon, B. (1986). Universal code equivalent of a yeast mitochondrial intron reading frame is expressed into E. coli as a specific double strand endonuclease. Cell, 44, 521–533. doi: 10.1016/0092-8674(86)90262-X
[4] Dutheil, J. Y., Münch, K., Schotanus, K., Stukenbrock, E. H., and Kahmann, R. (2020). The insertion of a mitochondrial selfish element into the nuclear genome and its consequences. bioRxiv, 787044, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/787044

Awareness has been growing that structural variants in the genome of species play a fundamental role in adaptive evolution and diversification [1]. Here, Le Moan and co-authors [2] report empirical genomic-wide SNP data on the European plaice (Pleuronectes platessa) across a major environmental transmission zone, ranging from the North Sea to the Baltic Sea. Regions of high linkage disequilibrium suggest the presence of two structural variants that appear to have evolved 220 kya. These two putative structural variants show weak signatures of isolation by distance when contrasted against the rest of the genome, but the frequency of the different putative structural variants appears to co-vary in some parts of the studied range with the environment, indicating the involvement of both selective and neutral processes. This study adds to the mounting body of evidence that structural genomic variants harbour significant information that allows species to respond and adapt to the local environmental context.

References

[1] Wellenreuther, M., Mérot, C., Berdan, E., & Bernatchez, L. (2019). Going beyond SNPs: the role of structural genomic variants in adaptive evolution and species diversification. Molecular ecology, 28(6), 1203-1209. doi: 10.1111/mec.15066
[2] Le Moan, A. Bekkevold, D. & Hemmer-Hansen J. (2020). Evolution at two time-frames: ancient and common origin of two structural variants involved in local adaptation of the European plaice (Pleuronectes platessa). bioRxiv, 662577, ver. 5 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/662577

The present-day distribution of extant species is the result of the interplay between their past population demography (e.g., expansion, contraction, isolation, and migration) and adaptation to the environment. Shedding light on the timing and magnitude of key demographic events helps identify potential drivers of such events and interaction of those drivers, such as life history traits and past episodes of environmental shifts.

The understanding of the key factors driving species evolution gives important insights into how the species may respond to changing conditions, which can be particularly relevant for the management of harmful species, such as agricultural pests (e.g. [1]). Meaningful demographic inferences present major challenges. These include formulating evolutionary scenarios fitting species biology and the eco-geographical context and choosing informative molecular markers and accurate quantitative approaches to statistically compare multiple demographic scenarios and estimate the parameters of interest. A further issue comes with result interpretation. Accurately dating the inferred events is far from straightforward since reliable calibration points are necessary to translate the molecular estimates of the evolutionary time into absolute time units (i.e. years). This can be attempted in different ways, such as by using fossil and archaeological records, heterochronous samples (e.g. ancient DNA), and/or mutation rate estimated from independent data (e.g. [2], [3] for review). Nonetheless, most experimental systems rarely meet these conditions, hindering the comprehensive interpretation of results.

The contribution of Chapuis et al. [4] addresses these issues to investigate the recent history of the African insect pest Schistocerca gregaria (desert locust). They apply Approximate Bayesian Computation-Random Forest (ABC-RF) approaches to microsatellite markers. Owing to their fast mutation rate microsatellite markers offer at least two advantages: i) suitability for analyzing recently diverged populations, and ii) direct estimate of the germline mutation rate in pedigree samples. The work of Chapuis et al. [4] benefits of both these advantages, since they have estimates of mutation rate and allele size constraints derived from germline mutations in the species [5].

The main aim of the study is to infer the history of divergence of the two subspecies of the desert locust, which have spatially disjoint distribution corresponding to the dry regions of North and West-South Africa. They first use paleo-vegetation maps to formulate hypotheses about changes in species range since the last glacial maximum. Based on them, they generate 12 divergence models. For the selection of the demographic model and parameter estimation, they apply the recently developed ABC-RF approach, a powerful inferential tool that allows optimizing the use of summary statistics information content, among other advantages [6]. Some methodological novelties are also introduced in this work, such as the computation of the error associated with the posterior parameter estimates under the best scenario. The accuracy of timing estimate is assured in two ways: i) by the use of microsatellite markers with known evolutionary dynamics, as underlined above, and ii) by assessing the divergence time threshold above which posterior estimates are likely to be biased by size homoplasy and limits in allele size range [7]. The best-supported model suggests a recent divergence event of the subspecies of S. gregaria (around 2.6 kya) and a reduction of populations size in one of the subspecies (S. g. flaviventris) that colonized the southern distribution area. As such, results did not support the hypothesis that the southward colonization was driven by the expansion of African dry environments associated with the last glacial maximum, as it has been postulated for other arid-adapted species with similar African disjoint distributions [8]. The estimated time of divergence points at a much more recent origin for the two subspecies, during the late Holocene, in a period corresponding to fairly stable arid conditions similar to current ones [9,10].

Although the authors cannot exclude that their microsatellite data bear limited information on older colonization events than the last one, they bring arguments in favour of alternative explanations. The hypothesis privileged does not involve climatic drivers, but the particularly efficient dispersal behaviour of the species, whose individuals are able to fly over long distances (up to thousands of kilometers) under favourable windy conditions. A single long-distance dispersal event by a few individuals would explain the genetic signature of the bottleneck. There is a growing number of studies in phylogeography in arid regions in the Southern hemisphere, but the impact of past climate changes on the species distribution in this region remains understudied relative to the Northern hemisphere [11,12].

The study presented by Chapuis et al. [4] offers several important insights into demographic changes and the evolutionary history of an agriculturally important pest species in Africa, which could also mirror the history of other organisms in the continent. As the authors point out, there are necessarily some uncertainties associated with the models of past ecosystems and climate, especially for Africa. Interestingly, the authors argue that the information on paleo-vegetation turnover was more informative than climatic niche modeling for the purpose of their study since it made them consider a wider range of bio-geographical changes and in turn a wider range of evolutionary scenarios (see discussion in Supplementary Material). Microsatellite markers have been offering a useful tool in population genetics and phylogeography for decades, but their popularity is perhaps being taken over by single nucleotide polymorphism (SNP) genotyping and whole-genome sequencing (WGS) (the peak year of the number of the publication with “microsatellite” is in 2012 according to PubMed).

This study reaffirms the usefulness of these classic molecular markers to estimate past demographic events, especially when species- and locus-specific microsatellite mutation features are available and a powerful inferential approach is adopted. Nonetheless, there are still hurdles to overcome, such as the limitations in scenario choice associated with the simulation software used (e.g. not allowing for continuous gene flow in this particular case), which calls for further improvement of simulation tools allowing for more flexible modeling of demographic events and mutation patterns. In sum, this work not only contributes to our understanding of the makeup of the African biodiversity but also offers a useful statistical framework, which can be applied to a wide array of species and molecular markers (microsatellites, SNPs, and WGS).

References

[1] Lehmann, P. et al. (2018). Complex responses of global insect pests to climate change. bioRxiv, 425488. doi: https://dx.doi.org/10.1101/425488

[2] Donoghue, P. C., & Benton, M. J. (2007). Rocks and clocks: calibrating the Tree of Life using fossils and molecules. Trends in Ecology & Evolution, 22(8), 424-431. doi: https://dx.doi.org/10.1016/j.tree.2007.05.005

[3] Ho, S. Y., Lanfear, R., Bromham, L., Phillips, M. J., Soubrier, J., Rodrigo, A. G., & Cooper, A. (2011). Time‐dependent rates of molecular evolution. Molecular ecology, 20(15), 3087-3101. doi: https://dx.doi.org/10.1111/j.1365-294X.2011.05178.x

[4] Chapuis, M.-P., Raynal, L., Plantamp, C., Meynard, C. N., Blondin, L., Marin, J.-M. and Estoup, A. (2020). A young age of subspecific divergence in the desert locust Schistocerca gregaria, inferred by ABC Random Forest. bioRxiv, 671867, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: https://dx.doi.org/10.1101/671867

5] Chapuis, M.-P., Plantamp, C., Streiff, R., Blondin, L., & Piou, C. (2015). Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations. Molecular ecology, 24(24), 6107-6119. doi: https://dx.doi.org/10.1111/mec.13465

[6] Raynal, L., Marin, J. M., Pudlo, P., Ribatet, M., Robert, C. P., & Estoup, A. (2018). ABC random forests for Bayesian parameter inference. Bioinformatics, 35(10), 1720-1728. doi: https://dx.doi.org/10.1093/bioinformatics/bty867

[7] Estoup, A., Jarne, P., & Cornuet, J. M. (2002). Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Molecular ecology, 11(9), 1591-1604. doi: https://dx.doi.org/10.1046/j.1365-294X.2002.01576.x

[8] Moodley, Y. et al. (2018). Contrasting evolutionary history, anthropogenic declines and genetic contact in the northern and southern white rhinoceros (Ceratotherium simum). Proceedings of the Royal Society B, 285(1890), 20181567. doi: https://dx.doi.org/10.1098/rspb.2018.1567

[9] Kröpelin, S. et al. (2008). Climate-driven ecosystem succession in the Sahara: the past 6000 years. science, 320(5877), 765-768. doi: https://dx.doi.org/10.1126/science.1154913

[10] Maley, J. et al. (2018). Late Holocene forest contraction and fragmentation in central Africa. Quaternary Research, 89(1), 43-59. doi: https://dx.doi.org/10.1017/qua.2017.97

[11] Beheregaray, L. B. (2008). Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere. Molecular Ecology, 17(17), 3754-3774. doi: https://dx.doi.org/10.1111/j.1365-294X.2008.03857.x

[12] Dubey, S., & Shine, R. (2012). Are reptile and amphibian species younger in the Northern Hemisphere than in the Southern Hemisphere?. Journal of evolutionary biology, 25(1), 220-226. doi: https://dx.doi.org/10.1111/j.1420-9101.2011.02417.x

*****

A video about this preprint is available here:

The science of naming species (taxonomy) has been renewed with the developments of molecular sequencing, digitization of museum specimens, and novel analytical tools. However, naming species can be highly subjective, sometimes considered as an art [1], because it is based on human-based criteria that vary among taxonomists. Nonetheless, taxonomists often argue that species names are hypotheses, which are therefore testable and refutable as new evidence is provided. This challenge comes with a more and more recognized and critical need for rigorously delineated species not only for producing accurate species inventories, but more importantly many questions in evolutionary biology (e.g. speciation), ecology (e.g. ecosystem structure and functioning), conservation biology (e.g. targeting priorities) or biogeography (e.g. diversification processes) depend in part on those species inventories and our knowledge of species [2-3]. Inaccurate species boundaries or diversity estimates may lead us to deliver biased answers to those questions, exactly as phylogenetic trees must be reconstructed rigorously and analyzed critically because they are a first step toward discussing broader questions [2-3]. In this context, biological diversity needs to be studied from multiple and complementary perspectives requiring the collaboration of morphologists, molecular biologists, biogeographers, and modelers [4-5]. Integrative taxonomy has been proposed as a solution to tackle the challenge of delimiting species [2], especially in highly diverse and undocumented groups of organisms.
In an elegant study that harbors all the characteristics of an integrative approach, Alda et al. [6] tackle the delimitation of species within the snail genus Galba (Lymnaeidae). Snails of this genus represent a peculiar case study for species delineation with a long and convoluted taxonomic history in which previous works recognized a number of species ranging from 4 to 30. The confusion is likely due to a loose morphology (labile shell features and high plasticity), which makes the identification and naming of species very unstable and likely subjective. An integrative taxonomic approach was needed. After two decades of taxon sampling and visits of type localities, the authors present an impressively dense taxon sampling at a global scale for the genus, which includes all described species. When it comes to delineate species, taxon sampling is often the key if we want to embrace the genetic and morphological diversity. Molecular data was obtained for several types of markers (microsatellites and DNA sequences for four genes), which were combined to morphology of shell and of internal organs, and to geographic distribution. All the data are thoroughly analyzed with cutting-edge methods starting from Bayesian phylogenetic reconstructions using multispecies coalescent models, followed by models of species delimitation based on the molecular specimen-level phylogeny, and then Bayesian divergence time estimates. They also used probabilistic models of ancestral state estimation to infer the ancestral phenotypic state of the Galba ancestors.
Their numerous phylogenetic and delimitation analyses allow to redefine the species boundaries that indicate that the genus Galba comprises six species. Interestingly, four of these species are morphologically cryptic and likely constitute species with extensive genetic diversity and widespread geographic distribution. The other two species have more geographically restricted distributions and exhibit an alternative morphology that is more phylogenetically derived than the cryptic one. Although further genomic studies would be required to strengthen some species status, this novel delimitation of Galba species has important implications for our understanding of convergence and morphological stasis, or the role for stabilizing selection in amphibious habitats; topics that are rarely addressed with invertebrate groups. For instance, in terms of macroevolutionary history, it is striking that an invertebrate clade of that age (22 million years ago) has only given birth to six species today. Including 30 (ancient taxonomy) or 6 (integrative taxonomy) species in a similar amount of evolutionary time does not tell us the same story when studying the diversification processes [7]. Here, Alda et al. [6] present a convincing case study that should foster similar studies following their approach, which will provide stimulating perspectives for testing the concepts of species and their effects on evolutionary biology.

References

[1] Ohl, M. (2018). The art of naming. MIT Press.
[2] Dayrat, B. (2005). Towards integrative taxonomy. Biological Journal of the Linnean Society, 85(3), 407–415. doi: 10.1111/j.1095-8312.2005.00503.x
[3] De Queiroz, K. (2007). Species concepts and species delimitation. Systematic Biology, 56(6), 879–886. doi: 10.1080/10635150701701083
[4] Padial, J. M., Miralles, A., De la Riva, I., and Vences, M. (2010). The integrative future of taxonomy. Frontiers in Zoology, 7(1), 16. doi: 10.1186/1742-9994-7-16
[5] Schlick-Steiner, B. C., Steiner, F. M., Seifert, B., Stauffer, C., Christian, E., and Crozier, R. H. (2010). Integrative taxonomy: A multisource approach to exploring biodiversity. Annual Review of Entomology, 55(1), 421–438. doi: 10.1146/annurev-ento-112408-085432
[6] Alda, P. et al. (2019). Systematics and geographical distribution of Galba species, a group of cryptic and worldwide freshwater snails. BioRxiv, 647867, v3 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/647867
[7] Ruane, S., Bryson, R. W., Pyron, R. A., and Burbrink, F. T. (2014). Coalescent species delimitation in milksnakes (Genus Lampropeltis) and impacts on phylogenetic comparative analyses. Systematic Biology, 63(2), 231–250. doi: 10.1093/sysbio/syt099

Is evolution adaptive? Not if there is no variation for natural selection to work with. Theory predicts that how fast a population can adapt to a new environment can be limited by the supply of new mutations coming into it. This supply, in turn, depends on two things: how often mutations occur and in how many individuals. If there are few mutations, or few individuals in whom they can originate, individuals will be mostly identical in their DNA, and natural selection will be impotent.
This theoretical prediction has been hard to test. The rate at which new mutations arise in a population can be manipulated experimentally, and some work has shown that the fitness of a population increases more rapidly if more new mutations appear per generation, lending support to the mutation-limitation hypothesis [1]. However, the question remains whether this limitation has played a role in the history of life over the evolutionary timescale. Maybe all natural populations are so large, the mutation rate so high, and/or the environment changes so slowly, that any novel variant required for adaptation is already there when selection starts to act? Some recent work does suggest that when strong selection begins to favor a certain phenotype, multiple distinct genetic variants producing this phenotype spread; this is what has happened, for instance, at the origin of insecticide resistance in wild populations of Drosophila melanogaster [2] or lactose persistence in humans [3]. In many other cases, though, adaptations seem to originate through a single mutation event, suggesting that the time needed for this mutation to arise may be important.
To complicate things, adaptation is hard to quantify. It leaves a trace in differences between individuals of the same species as well as of different species. However, this trace is often masked or confounded by other processes, including natural selection disfavoring newly arising deleterious variants, interference from selection acting at linked sites, and changes in population size. In 1991, McDonald and Kreitman [4] have come up with a method to infer the rate of adaptation in the presence of strong negative selection, and later work has developed upon it to control for some of the other confounders. Still, the method is data-intensive, and previous attempts to employ it to compare the rates of adaptation between species have yielded somewhat contradictory results.
The new paper by Rousselle et al. recommended by PCI Evol Biol [5] fills this gap. The authors use published data as well as their own newly generated dataset to analyze, in a McDonald and Kreitman-like framework, both closely and distantly related species. Importantly, these comparisons cover species with very different polymorphism levels, spanning two orders of magnitude of difference levels.
So is adaptation in fact limited by supply of new mutations? The answer is, it depends. It does indeed seem that the species with a lower level of polymorphism adapt at a lower rate, consistent with the mutation-limitation hypothesis. However, this only is true for those groups of species in which the variability is low. Therefore, if a population is very small or the mutation rate very low, there may be in fact not enough mutations to secure its need to adapt.
In more polymorphic species, and in comparisons of distant species, the data hint instead at the opposite relationship: the rate of adaptations declines with variability. This is consistent with a different explanation: when a population is small, it needs to adapt more frequently, repairing the weakly deleterious mutations that can’t be prevented by selection under small population sizes.
There are quite a few problems small populations have to deal with. Some of them are ecological: e.g., small numbers make populations more vulnerable to stochastic fluctuations in size or sex ratio. Others, however, are genetic. Small populations are prone to inbreeding depression and have an increased rate of genetic drift, leading to spread of deleterious alleles. Indeed, selection against deleterious mutations is less efficient when populations are small, and less numerable species accumulate more of such mutations over the course of evolution [6]. The work by Rouselle et al. [5] suggests that small populations also face an additional burden: a reduced ability to adapt.
Has the rate of adaptation in our own species also been limited by our deficit of diversity? The data hints at this. Homo sapiens, as well as the two other studied extinct representatives of the genus Homo, Neanderthals and Denisovans, belong to the domain of relatively low polymorphism levels, where an increase in polymorphism matters for the rate at which adaptive substitutions accumulate. Perhaps, if our ancestors were more numerous or more mutable, they would have been able to get themselves out of trouble, and there would be multiple human species still alive rather than just one.

References

[1] G, J. A., Visser, M. de, Zeyl, C. W., Gerrish, P. J., Blanchard, J. L., and Lenski, R. E. (1999). Diminishing Returns from Mutation Supply Rate in Asexual Populations. Science, 283(5400), 404–406. doi: 10.1126/science.283.5400.404
[2] Karasov, T., Messer, P. W., and Petrov, D. A. (2010). Evidence that Adaptation in Drosophila Is Not Limited by Mutation at Single Sites. PLOS Genetics, 6(6), e1000924. doi: 10.1371/journal.pgen.1000924
[3] Jones, B. L., Raga, T. O., Liebert, A., Zmarz, P., Bekele, E., Danielsen, E. T., Olsen, A. K., Bradman, N., Troelsen, J. T., and Swallow, D. M. (2013). Diversity of Lactase Persistence Alleles in Ethiopia: Signature of a Soft Selective Sweep. The American Journal of Human Genetics, 93(3), 538–544. doi: 10.1016/j.ajhg.2013.07.008
[4] McDonald, J. H., and Kreitman, M. (1991). Adaptive protein evolution at the Adh locus in Drosophila. Nature, 351(6328), 652–654. doi: 10.1038/351652a0
[5] Rousselle, M., Simion, P., Tilak, M. K., Figuet, E., Nabholz, B., and Galtier, N. (2019). Is adaptation limited by mutation? A timescale-dependent effect of genetic diversity on the adaptive substitution rate in animals. BioRxiv, 643619, ver 4 peer-reviewed and recommended by Peer Community In Evolutionary Biology. doi: 10.1101/643619
[6] Popadin, K., Polishchuk, L. V., Mamirova, L., Knorre, D., and Gunbin, K. (2007). Accumulation of slightly deleterious mutations in mitochondrial protein-coding genes of large versus small mammals. Proceedings of the National Academy of Sciences, 104(33), 13390–13395. doi: 10.1073/pnas.0701256104