Genetic markers are used for in modern population genetics/genomics to uncover the past neutral and selective history of population and species. Besides Single Nucleotide Polymorphisms (SNPs) obtained from whole genome data, microsatellites (or Short Tandem Repeats, SSR) have been common markers of choice in numerous population genetics studies of non-model species with large sample sizes [1]. Microsatellites can be used to uncover and draw inference of the past population demography (e.g. expansion, decline, bottlenecks…), population split, population structure and gene flow, but also life history traits and modes of reproduction (e.g. [2,3]). These markers are widely used in conservation genetics [4] or to study parasites or disease vectors [5]. Microsatellites do show higher mutation rate than SNPs increasing, on the one hand, the statistical power to infer recent events (for example crop domestication, [2,3]), while, on the other hand, decreasing their statistical power over longer time scales due to homoplasy [6].
To perform such analyses, however, an excellent and reliable quality of data is required. As emphasized in the article by De Meeûs et al. [7] three main issues do bias the observed heterozygosity at microsatellites: null alleles, short allele dominance (SAD) and stuttering. These originates from poor PCR amplification. As a result, an excess of homozygosity is observed at the microsatellite loci leading to overestimation of the variation statistics FIS and FST as well as increased linage disequilibrium (LD). For null alleles, several methods and software do help to reduce the bias, and in the present study, De Meeûs et al. [7] propose a way to tackle issues with SAD and stuttering.
The authors study a dataset consisting of 387 samples from 61 subsamples genotyped at nine loci of the species Ixodes scapularis, i.e. ticks transmitting the Lyme disease. Based on correlation methods and FST, FIS they can uncover null alleles and SAD. Stuttering is detected by evaluating the heterozygote deficit between alleles displaying a single repeat difference. Without correction, six loci are affected by one of these amplification problems generating a large deficit of heterozygotes (measured by significant FIS and FST) remaining so after correction for the false discovery rate (FDR). These results would be classically interpreted as a strong Wahlund effect and/or selection at several loci.
After correcting for null alleles, the authors apply two novel corrections: 1) a re-examination of the chromatograms reveals previously disregarded larger alleles thus decreasing SAD, and 2) pooling alleles close in size decreasing stuttering. The corrected dataset shows then a significant excess of heterozygotes as could be expected in a dioecious species with strong population structure. The FDR correction removes then the significant excess of homozygotes and LD between pairs of loci. FST on the cured dataset is used to demonstrate the strong population structure and small effective subpopulation sizes. This is confirmed by a clustering analysis using discriminant analysis of principal components (DAPC).
While based on a specific dataset of ticks from different populations sampled across the USA, the generality of the authors’ approach is presented in Figure 6 in which they provide a step by step flowchart to cure microsatellite datasets from null alleles, SAD and stuttering. Several criteria based on FIS, FST and LD between loci are used as decision keys in the flowchart. An excel file is also provided as help for the curation steps. This study and the proposed methodology are thus extremely useful for all population geneticists working on non-model species with large number of samples genotyped at microsatellite markers. The method not only allows more accurate estimates of heterozygosity but also prevents the thinning of datasets due to the removal of problematic loci. As a follow-up and extension of this work, an exhaustive simulation study could investigate the influence of these data quality issues on past demographic and population structure inference under a wide range of scenarios. This would allow to quantify the current biases in the literature and the robustness of the methodology devised by De Meeûs et al. [7].

References

[1] Jarne, P., and Lagoda, P. J. (1996). Microsatellites, from molecules to populations and back. Trends in ecology & evolution, 11(10), 424-429. doi: 10.1016/0169-5347(96)10049-5
[2] Cornille, A., Giraud, T., Bellard, C., Tellier, A., Le Cam, B., Smulders, M. J. M., Kleinschmit, J., Roldan-Ruiz, I. and Gladieux, P. (2013). Postglacial recolonization history of the E uropean crabapple (Malus sylvestris M ill.), a wild contributor to the domesticated apple. Molecular Ecology, 22(8), 2249-2263. doi: 10.1111/mec.12231
[3] Parat, F., Schwertfirm, G., Rudolph, U., Miedaner, T., Korzun, V., Bauer, E., Schön C.-C. and Tellier, A. (2016). Geography and end use drive the diversification of worldwide winter rye populations. Molecular ecology, 25(2), 500-514. doi: 10.1111/mec.13495
[4] Broquet, T., Ménard, N., & Petit, E. (2007). Noninvasive population genetics: a review of sample source, diet, fragment length and microsatellite motif effects on amplification success and genotyping error rates. Conservation Genetics, 8(1), 249-260. doi: 10.1007/s10592-006-9146-5
[5] Koffi, M., De Meeûs, T., Séré, M., Bucheton, B., Simo, G., Njiokou, F., Salim, B., Kaboré, J., MacLeod, A., Camara, M., Solano, P., Belem, A. M. G. and Jamonneau, V. (2015). Population genetics and reproductive strategies of African trypanosomes: revisiting available published data. PLoS neglected tropical diseases, 9(10), e0003985. doi: 10.1371/journal.pntd.0003985
[6] 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: 10.1046/j.1365-294X.2002.01576.x
[7] De Meeûs, T., Chan, C. T., Ludwig, J. M., Tsao, J. I., Patel, J., Bhagatwala, J., and Beati, L. (2019). Deceptive combined effects of short allele dominance and stuttering: an example with Ixodes scapularis, the main vector of Lyme disease in the USA. bioRxiv, 622373, ver. 4 peer-reviewed and recommended by Peer Community In Evolutionary Biology. doi: 10.1101/622373

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

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

The study of sexual selection goes as far as Darwin himself. Since then, elaborate theories concerning both intra- and inter-sexual sexual have been developed, and elegant experiments have been designed to test this body of theory. It may thus come as a surprise that the community is still debating on the correct way to measure simple components of sexual selection, such as the Bateman gradient (i.e., the covariance between the number of matings and the number of offspring)[1,2], or to quantify complex behaviours such as mate choice (the non-random choice of individuals with particular characters as mates)[3,4] and their consequences.
One difficulty in the study of sexual selection is evaluating the consequences of non-random mating. Indeed, when non-random mating is observed in a population, it is often difficult to establish whether such mating pattern leads to i) sexual selection per se (selection pressures favouring certain phenotypes), and/or ii) the non-random association of parental genes in their offspring or not. These two processes differ. In particular, assortative (and disassortative) mating can shape genetic covariances without leading to changes in gene frequencies in the population. Their distinction matters because these two processes lead to different evolutionary outcomes, which can have large ripple effects in the evolution of sexual behaviours, sexual ornamentation, and speciation.
In his paper, entitled “Multi-model inference of non-random mating from an information theoretic approach” [5], Carvajal-Rodríguez tackled this issue. The author generated a simple model in which the consequences of non-random mating can be inferred from information on the population frequencies before and after mating. The procedure is as follows: from the initial population frequencies of phenotypes (or genotypes) of both sexes, the model generates predictions on the frequencies after mating, assuming that particular mating patterns have occurred. This leads to different predictions for the phenotypic (or genotypic) frequencies after mating. The particular mating pattern leading to the best fit with the real frequencies is then identified via a model selection procedure (performing model averaging to combine different mating patterns is also possible).
This study builds on a framework introduced by Carvajal-Rodríguez’s colleagues [6] and encompasses later methodological developments involving the author himself [7]. Compared to early work, the new method proposed by the author builds on the relationship between mating pattern and information [8] to distinguish among scenarios that would lead to non-random mating due to different underlying processes, using simple model selection criterion such as the AICc.
The great asset of the proposed method is that it can be applied to the study of natural populations in which the study of mate choice and sexual selection is notoriously difficult. In the manuscript, the procedure is tested on a population of marine gastropods (Littorina saxatilis). This allows the reader to grasp how the method can be applied to a real system. In fact, anyone can try out the method thanks to the freely available software InfoMating programmed by the author. One important assumption underlying the current method is that the frequencies of unmated individuals do not change during the mating season. If this is not the case, the reader may refer to another publication of the same author which relaxes this assumption [9]. These papers are both instrumental for empiricists interested in testing sexual selection theory.

References

[1] Bateman, A. J. (1948). Intra-sexual selection in Drosophila. Heredity, 2(3), 349-368. doi: 10.1038/hdy.1948.21
[2] Jones, A. G. (2009). On the opportunity for sexual selection, the Bateman gradient and the maximum intensity of sexual selection. Evolution: International Journal of Organic Evolution, 63(7), 1673-1684. doi: 10.1111/j.1558-5646.2009.00664.x
[3] Andersson, M., & Simmons, L. W. (2006). Sexual selection and mate choice. Trends in ecology & evolution, 21(6), 296-302. doi: 10.1016/j.tree.2006.03.015
[4] Kuijper, B., Pen, I., & Weissing, F. J. (2012). A guide to sexual selection theory. Annual Review of Ecology, Evolution, and Systematics, 43, 287-311. doi: 10.1146/annurev-ecolsys-110411-160245
[5] Carvajal-Rodríguez, A. (2019). Multi-model inference of non-random mating from an information theoretic approach. bioRxiv, 305730, ver. 5 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/305730
[6] Rolán‐Alvarez, E., & Caballero, A. (2000). Estimating sexual selection and sexual isolation effects from mating frequencies. Evolution, 54(1), 30-36. doi: 10.1111/j.0014-3820.2000.tb00004.x
[7] Carvajal-Rodríguez, A., & Rolan-Alvarez, E. (2006). JMATING: a software for the analysis of sexual selection and sexual isolation effects from mating frequency data. BMC Evolutionary Biology, 6(1), 40. doi: 10.1186/1471-2148-6-40
[8] Carvajal-Rodríguez, A. (2018). Non-random mating and information theory. Theoretical population biology, 120, 103-113. doi: 10.1016/j.tpb.2018.01.003
[9] Carvajal-Rodríguez, A. (2019). A generalization of the informational view of non-random mating: Models with variable population frequencies. Theoretical population biology, 125, 67-74. doi: 10.1016/j.tpb.2018.12.004

Artificial selection experiments are key experiments in evolutionary biology. The demonstration that application of selective pressure across multiple generations results in heritable phenotypic changes is a tangible and reproducible proof of the evolution by natural selection.
Artificial selection experiments are used to evaluate the joint effects of selection on multiple traits, their genetic covariances and differences in responses in different environments. Most studies on artificial selection experiments report and base their analyses on phenotypic changes [1]. More recently, changes in allele frequency and other patterns of molecular genetic diversity have been used to identify genomic locations where selection has had an effect. However, so far the changes in gene expression have not been in the focus of artificial selection experiment studies (see [2] for an example though).
In plants, one of the most famous artificial selection experiments is the Illinois Corn Experiment where maize (Zea mays) is selected for oil and protein content [3], but in addition, similar experiments have been conducted also for other traits in maize. In Saclay divergent selection experiment [4] two maize inbred lines (F252 and MBS847) have been selected for early and late flowering for 13 generations, resulting in two week difference in flowering time.
In ”Transcriptomic response to divergent selection for flowering time in maize reveals convergence and key players of the underlying gene regulatory network ” [5] Maud Tenaillon and her coworkers study the gene expression differences among these two independently selected maize populations. Their experiments cover two years in field conditions and they use samples of shoot apical meristem at three different developmental stages: vegetative, transitioning and reproductive. They use RNA-seq transcriptome level differences and qRT-PCR for gene expression pattern investigation. The work is continuation to earlier genetic and phenotypic studies on the same material [4, 6].
The reviewers and I agree that dataset is unique and its major benefit is that it has been obtained from field conditions similar to those that species may face under natural setting during selection. Their tissue sampling is supported by flowering time phenotypic observations and covers the developmental transition stage, making a good effort to identify key transcriptional and phenotypic changes and their timing affected by selection.
Tenaillon et al. [5] identify more than 2000 genes that are differentially expressed among early and late flowering populations. Expectedly, they are enriched for known flowering time genes. As they point out, differential expression of thousands of genes does not mean that they all were independently affected by selection, but rather that the whole transcriptional network has shifted, possibly due to just few upstream or hub-genes. Also, the year-to-year variation had smaller effect in gene expression compared to developmental stage or genetic background, possibly indicating selection for stability across environmental fluctuation for such an important phenotype as flowering time.
Another noteworthy observation is that they find convergent patterns of transcriptional changes among the two selected lines. 115 genes expression patterns are shifted due to selection in both genetic backgrounds. This convergent pattern can be a result of either selection on standing variation or de novo mutations. The data does not allow testing which process is underlying the observed convergence. However, their results show that this is an interesting future question that can be addressed using genotype and gene expression data from the same ancestral and derived material and possibly their hybrids.

References

[1] Hill, W. G., & Caballero, A. (1992). Artificial selection experiments. Annual Review of Ecology and Systematics, 23(1), 287-310. doi: 10.1146/annurev.es.23.110192.001443
[2] Konczal, M., Babik, W., Radwan, J., Sadowska, E. T., & Koteja, P. (2015). Initial molecular-level response to artificial selection for increased aerobic metabolism occurs primarily through changes in gene expression. Molecular biology and evolution, 32(6), 1461-1473. doi: 10.1093/molbev/msv038
[3] Moose, S. P., Dudley, J. W., & Rocheford, T. R. (2004). Maize selection passes the century mark: a unique resource for 21st century genomics. Trends in plant science, 9(7), 358-364. doi: 10.1016/j.tplants.2004.05.005
[4] Durand, E., Tenaillon, M. I., Ridel, C., Coubriche, D., Jamin, P., Jouanne, S., Ressayre, A., Charcosset, A. and Dillmann, C. (2010). Standing variation and new mutations both contribute to a fast response to selection for flowering time in maize inbreds. BMC evolutionary biology, 10(1), 2. doi: 10.1186/1471-2148-10-2
[5] Tenaillon, M. I., Seddiki, K., Mollion, M., Le Guilloux, M., Marchadier, E., Ressayre, A. and Dillmann C. (2019). Transcriptomic response to divergent selection for flowering time in maize reveals convergence and key players of the underlying gene regulatory network. BioRxiv, 461947 ver. 5 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/461947
[6] Durand, E., Tenaillon, M. I., Raffoux, X., Thépot, S., Falque, M., Jamin, P., Bourgais A., Ressayre, A. and Dillmann, C. (2015). Dearth of polymorphism associated with a sustained response to selection for flowering time in maize. BMC evolutionary biology, 15(1), 103. doi: 10.1186/s12862-015-0382-5