Nearly 8 million people are suffering from Chagas disease in the Americas. The etiological agent, Trypanosoma cruzi, is mainly transmitted by triatomine bugs, also known as kissing or vampire bugs, which suck blood and transmit the parasite through their feces. Among these triatomine species, Rhodnius prolixus is considered the main vector, and many studies have focused on characterizing its biology, physiology, ecology and evolution. 

Interestingly, given that Rhodnius species feed almost exclusively on blood, their diet is unbalanced, and the insects can lack nutrients and vitamins that they cannot synthetize themself, such as B-vitamins. In all insects feeding exclusively on blood, symbioses with microbes producing B-vitamins (mainly biotin, riboflavin and folate) have been widely described (see review in Duron and Gottlieb 2020) and are critical for insect development and reproduction. These co-evolved relationships between blood feeders and nutritional symbionts could now be considered to develop new control methods, by targeting the ‘Achille’s heel’ of the symbiotic association (i.e., transfer of nutrient and / or control of nutritional symbiont density). But for this, it is necessary to better characterize the relationships between triatomines and their symbionts. 

R. prolixus is known to be associated with several symbionts. The extracellular gut symbiont Rhodococcus rhodnii, which reaches high bacterial densities and is almost fixed in R. prolixus populations, appears to be a nutritional symbiont under many blood sources. This symbiont can provide B-vitamins such as biotin (B7), niacin (B3), thiamin (B1), pyridoxin (B6) or riboflavin (B2) and can play an important role in the development and the reproduction of R. prolixus (Pachebat et al. (2013) and see review in Salcedo-Porras et al. (2020)). This symbiont is orally acquired through egg smearing, ensuring the fidelity of transmission of the symbiont from mother to offspring. However, as recently highlighted by Tobias et al. (2020) and Gilliland et al. (2022), other gut microbes could also participate to the provision of B-vitamins, and R. rhodnii could additionally provide metabolites (other than B-vitamins) increasing bug fitness. In the study from Filée et al., the authors focused on Wolbachia, an intracellular, maternally inherited bacterium, known to be a nutritional symbiont in other blood-sucking insects such as bedbugs (Nikoh et al. 2014), and its potential role in vitamin provision in triatomine bugs. 

After screening 17 different triatomine species from the 3 phylogenetic groups prolixus, pallescens and pictipes, they first show that Wolbachia symbionts are widely distributed in the different Rhodnius species. Contrary to R. rhodnii that were detected in all samples, Wolbachia prevalence was patchy and rarely fixed. The authors then sequenced, assembled, and compared 13 Wolbachia genomes from the infected Rhodnius species. They showed that all Wolbachia are phylogenetically positioned in the supergroup F that contains wCle (the Wolbachia from bedbugs). In addition, 8 Wolbachia strains (out of 12) encode a biotin operon under strong purifying selection, suggesting the preservation of the biological function and the metabolic potential of Wolbachia to supplement biotin in their Rhodnius host. From the study of insect genomes, the authors also evidenced several horizontal transfers of genes from Wolbachia to Rhodnius genomes, which suggests a complex evolutionary interplay between vampire bugs and their intracellular symbiont. 

This nice piece of work thus provides valuable information to the fields of multiple partners / nutritional symbioses and Wolbachia research. Dual symbioses described in insects feeding on unbalanced diets generally highlight a certain complementarity between symbionts that ensure the whole nutritional complementation. The study presented by Filée et al. leads rather to consider the impact of multiple symbionts with different lifestyles and transmission modes in the provision of a specific nutritional benefit (here, biotin). Because of the low prevalence of Wolbachia in certain species, a “ménage à trois” scenario would rather be replaced by an “open couple”, where the host relationship with new symbiotic partners (more or less stable at the evolutionary timescale) could provide benefits in certain ecological situations. The results also support the potential for Wolbachia to evolve rapidly along a continuum between parasitism and mutualism, by acquiring operons encoding critical pathways of vitamin biosynthesis.

References

Duron O. and Gottlieb Y. (2020) Convergence of Nutritional Symbioses in Obligate Blood Feeders. Trends in Parasitology 36(10):816-825. https://doi.org/10.1016/j.pt.2020.07.007

Filée J., Agésilas-Lequeux K., Lacquehay L., Bérenger J.-M., Dupont L., Mendonça V., Aristeu da Rosa J. and Harry M. (2023) Wolbachia genomics reveals a potential for a nutrition-based symbiosis in blood-sucking Triatomine bugs. bioRxiv, 2022.09.06.506778, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.09.06.506778

Gilliland C.A. et al. (2022) Using axenic and gnotobiotic insects to examine the role of different microbes on the development and reproduction of the kissing bug Rhodnius prolixus (Hemiptera: Reduviidae). Molecular Ecology. https://doi.org/10.1111/mec.16800

Nikoh et al. (2014) Evolutionary origin of insect–Wolbachia nutritional mutualism. PNAS. 111(28):10257-10262. https://doi.org/10.1073/pnas.1409284111

Pachebat, J.A. et al. (2013). Draft genome sequence of Rhodococcus rhodnii strain LMG5362, a symbiont of Rhodnius prolixus (Hemiptera, Reduviidae, Triatominae), the principle vector of Trypanosoma cruzi. Genome Announc. 1(3):e00329-13. https://doi.org/10.1128/genomea.00329-13

Salcedo-Porras N., et al. (2020). The role of bacterial symbionts in Triatomines: an evolutionary perspective. Microorganisms. 8:1438. https://doi.org/10.3390%2Fmicroorganisms8091438

Tobias N.J., Eberhard F.E., Guarneri A.A. (2020) Enzymatic biosynthesis of B-complex vitamins is supplied by diverse microbiota in the Rhodnius prolixus anterior midgut following Trypanosoma cruzi infection. Computational and Structural Biotechnology Journal. 3395-3401. https://doi.org/10.1016/j.csbj.2020.10.031 

The article “Drivers of genomic landscapes of differentiation across Populus divergence gradient” by Shang et al. describes an amazing dataset where genomic variations among 21 pairs of diverging poplar species are compared. Such comparisons are still quite rare and are needed to shed light on the processes shaping genomic divergence along the speciation gradient. Relying on two hundred whole-genome resequenced samples from 8 species that diverged from 1.3 to 4.8 million years ago, the authors aim at identifying the key factors involved in the genomic differentiation between species. They carried out a wide range of robust statistical tests aiming at characterizing the genomic differentiation along the genome of these species pairs. They highlight in particular the role of linked selection and gene flow in shaping the divergence along the genomes of species pairs. They also confirm the significance of introgression among species with a net divergence larger than the upper boundaries of the grey zone of speciation previously documented in animals (da from 0.005 to 0.02, Roux et al. 2016). Because these findings pave the way to research about the genomic mechanisms associated with speciation in species with allopatric and parapatric distributions, I warmingly recommend this article.

References

Roux C, Fraïsse C, Romiguier J, Anciaux Y, Galtier N, Bierne N (2016) Shedding Light on the Grey Zone of Speciation along a Continuum of Genomic Divergence. PLOS Biology, 14, e2000234. https://doi.org/10.1371/journal.pbio.2000234

Shang H, Rendón-Anaya M, Paun O, Field DL, Hess J, Vogl C, Liu J, Ingvarsson PK, Lexer C, Leroy T (2023) Drivers of genomic landscapes of differentiation across Populus divergence gradient. bioRxiv, 2021.08.26.457771, ver. 5 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.08.26.457771

The causes and consequences of the evolution of sexual reproduction are a major topic in evolutionary biology. With advances in sequencing technology, it becomes possible to compare sexual chromosomes across species and infer the neutral and selective processes shaping polymorphism at these chromosomes. Most sex and mating-type chromosomes exhibit an absence of recombination in large genomic regions around the animal, plant or fungal sex-determining genes. This suppression of recombination likely occurred in several time steps generating stepwise increasing genomic regions starting around the sex-determining genes. This mechanism generates so-called evolutionary strata of differentiation between sex chromosomes (Nicolas et al., 2004, Bergero and Charlesworth, 2009, Hartmann et al. 2021). The evolution of extended regions of recombination suppression is also documented on mating-type chromosomes in fungi (Hartmann et al., 2021) and around supergenes (Yan et al., 2020, Jay et al., 2021). The exact reason and evolutionary mechanisms for this phenomenon are still, however, debated.

Two hypotheses are proposed: 1) sexual antagonism (Charlesworth et al., 2005), which, nevertheless, does explain the observed occurrence of the evolutionary strata, and 2) the sheltering of deleterious alleles by inversions carrying a lower load than average in the population (Charlesworth and Wall, 1999, Antonovics and Abrams, 2004). In the latter, the mechanism is as follows. A genetic inversion or a suppressor of recombination in cis may exhibit some overdominance behaviour. The inversion exhibiting less recessive deleterious mutations (compared to others at the same locus) may increase in frequency, before at higher frequency occurring at the homozygous state, expressing its genetic load. However, the inversion may never be at the homozygous state if it is genetically linked to a gene in a permanently heterozygous state. The inversion can then be advantageous and may reach fixation at the sex chromosome (Charlesworth and Wall, 1999, Antonovics and Abrams, 2004, Jay et al., 2022). These selective mechanisms promote thus the suppression of recombination around the sex-determining gene, and recessive deleterious mutations are permanently sheltered. This hypothesis is corroborated by the sheltering of deleterious mutations observed around loci under balancing selection (Llaurens et al. 2009, Lenz et al. 2016) and around mating-type genes in fungi and supergenes (Jay et al. 2021, Jay et al., 2022).

In this present theoretical study, Tezenas et al. (2022) analyse how linkage to a necessarily heterozygous fungal mating type locus influences the persistence/extinction time of a new mutation at a second selected locus. This mutation can either be deleterious and recessive, or overdominant. There is arbitrary linkage between the two loci, and sexual reproduction occurs either between 1) gametes of different individuals (outcrossing), or 2) by selfing with gametes originating from the same (intra-tetrad) or different (inter-tetrad) tetrads produced by that individual. Note, here, that the mating-type gene does not prevent selfing. The authors study the initial stochastic dynamics of the mutation using a multi-type branching process (and simulations when analytical results cannot be obtained) to compute the extinction time of the deleterious mutation. The main result is that the presence of a mating-type locus always decreases the purging probability and increases the purging time of the mutations under selfing. Ultimately, deleterious mutations can indeed accumulate near the mating-type locus over evolutionary time scales. In a nutshell, high selfing or high intra-tetrad mating do increase the sheltering effect of the mating-type locus. In effect, the outcome of sheltering of deleterious mutations depends on two opposing mechanisms: 1) a higher selfing rate induces a greater production of homozygotes and an increased effect of the purging of deleterious mutations, while 2) a higher intra-tetrad selfing rate (or linkage with the mating-type locus) generates heterozygotes which have a small genetic load (and are favoured). The authors also show that rare events of extremely long maintenance of deleterious mutations can occur.

The authors conclude by highlighting the manifold effect of selfing which reduces the effective population size and thus impairs the efficiency of selection and increases the mutational load, while also favouring the purge of deleterious homozygous mutations. Furthermore, this study emphasizes the importance of studying the maintenance and accumulation of deleterious mutations in the vicinity of heterozygous loci (e.g. under balancing selection) in selfing species.

References

Antonovics J, Abrams JY (2004) Intratetrad Mating and the Evolution of Linkage Relationships. Evolution, 58, 702–709. https://doi.org/10.1111/j.0014-3820.2004.tb00403.x

Bergero R, Charlesworth D (2009) The evolution of restricted recombination in sex chromosomes. Trends in Ecology & Evolution, 24, 94–102. https://doi.org/10.1016/j.tree.2008.09.010

Charlesworth D, Morgan MT, Charlesworth B (1990) Inbreeding Depression, Genetic Load, and the Evolution of Outcrossing Rates in a Multilocus System with No Linkage. Evolution, 44, 1469–1489. https://doi.org/10.1111/j.1558-5646.1990.tb03839.x

Charlesworth D, Charlesworth B, Marais G (2005) Steps in the evolution of heteromorphic sex chromosomes. Heredity, 95, 118–128. https://doi.org/10.1038/sj.hdy.6800697

Charlesworth B, Wall JD (1999) Inbreeding, heterozygote advantage and the evolution of neo–X and neo–Y sex chromosomes. Proceedings of the Royal Society of London. Series B: Biological Sciences, 266, 51–56. https://doi.org/10.1098/rspb.1999.0603

Hartmann FE, Duhamel M, Carpentier F, Hood ME, Foulongne-Oriol M, Silar P, Malagnac F, Grognet P, Giraud T (2021) Recombination suppression and evolutionary strata around mating-type loci in fungi: documenting patterns and understanding evolutionary and mechanistic causes. New Phytologist, 229, 2470–2491. https://doi.org/10.1111/nph.17039

Jay P, Chouteau M, Whibley A, Bastide H, Parrinello H, Llaurens V, Joron M (2021) Mutation load at a mimicry supergene sheds new light on the evolution of inversion polymorphisms. Nature Genetics, 53, 288–293. https://doi.org/10.1038/s41588-020-00771-1

Jay P, Tezenas E, Véber A, Giraud T (2022) Sheltering of deleterious mutations explains the stepwise extension of recombination suppression on sex chromosomes and other supergenes. PLOS Biology, 20, e3001698. https://doi.org/10.1371/journal.pbio.3001698

Lenz TL, Spirin V, Jordan DM, Sunyaev SR (2016) Excess of Deleterious Mutations around HLA Genes Reveals Evolutionary Cost of Balancing Selection. Molecular Biology and Evolution, 33, 2555–2564. https://doi.org/10.1093/molbev/msw127

Llaurens V, Gonthier L, Billiard S (2009) The Sheltered Genetic Load Linked to the S Locus in Plants: New Insights From Theoretical and Empirical Approaches in Sporophytic Self-Incompatibility. Genetics, 183, 1105–1118. https://doi.org/10.1534/genetics.109.102707

Nicolas M, Marais G, Hykelova V, Janousek B, Laporte V, Vyskot B, Mouchiroud D, Negrutiu I, Charlesworth D, Monéger F (2004) A Gradual Process of Recombination Restriction in the Evolutionary History of the Sex Chromosomes in Dioecious Plants. PLOS Biology, 3, e4. https://doi.org/10.1371/journal.pbio.0030004

Tezenas E, Giraud T, Véber A, Billiard S (2022) The fate of recessive deleterious or overdominant mutations near mating-type loci under partial selfing. bioRxiv, 2022.10.07.511119, ver. 2 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.10.07.511119

Yan Z, Martin SH, Gotzek D, Arsenault SV, Duchen P, Helleu Q, Riba-Grognuz O, Hunt BG, Salamin N, Shoemaker D, Ross KG, Keller L (2020) Evolution of a supergene that regulates a trans-species social polymorphism. Nature Ecology & Evolution, 4, 240–249. https://doi.org/10.1038/s41559-019-1081-1

Studying the effects of speciation, hybridisation, and evolutionary outcomes following reproduction from divergent populations is a major research area in evolutionary genetics [1]. There are two phenomena that have been the focus of contemporary research. First, a classic concept is the formation of ‘Bateson-Dobzhansky-Muller’ incompatibilities (BDMi) [2–4] that negatively affect hybrid fitness. Here, two diverging populations accumulate mutations over time that are unique to that subpopulation. If they subsequently meet, then these mutations might negatively interact, leading to a loss in fitness or even a complete lack of reproduction. BDMi formation can be complex, involving multiple genes and the fitness changes can depend on the direction of introgression [5]. Second, such secondary contact can instead lead to heterosis, where offspring are fitter than their parental progenitors [6].

Understanding which outcomes are likely to arise require one to know the potential fitness effects of mutations underlying reproductive isolation, to determine whether they are likely to reduce or enhance fitness when hybrids are formed. This is far from an easy task, as it requires one to track mutations at several loci, along with their effects, across a fitness landscape.

The work of De Sanctis et al. [7] neatly fills in this knowledge gap, by creating a general mathematical framework for describing the consequences of a cross from two divergent populations. The derivations are based on Fisher’s Geometric Model, which is widely used to quantify selection acting on a general fitness landscape that is affected by several biological traits [8,9], and has previously been used in theoretical studies of hybridisation [10–12]. By doing so, they are able to decompose how divergence at multiple loci affects offspring fitness through both additive and dominance effects.

A key result arising from their analyses is demonstrating how offspring fitness can be captured by two main functions. The first one is the ‘net effect of evolutionary change’ that, broadly defined, measures how phenotypically divergent two populations are. The second is the ‘total amount of evolutionary change’, which reflects how many mutations contribute to divergence and the effect sizes captured by each of them. The authors illustrate these measurements using simulations covering different scenarios, demonstrating how different parental states can lead to similar fitness outcomes. They also propose experimental methods to measure the underlying mutational effects.

This study neatly demonstrates how complex genetic phenomena underlying hybridisation can be captured using fairly simple mathematical formulae. This powerful approach will thus open the door for future research to investigate hybridisation in more detail, whether it is by expanding on these theoretical models or using the elegant outcomes to quantify fitness effects in experiments.

References

1. Coyne JA, Orr HA. Speciation. Sunderland, Mass: Sinauer Associates; 2004.
2. Bateson W, Seward A. Darwin and modern science. Heredity and variation in modern lights. 1909;85: 101. https://doi.org/10.1017/CBO9780511693953.007
3. Dobzhansky T. Genetics and the Origin of Species. Columbia university press; 1937.
4. Muller HJ. Isolating mechanisms, evolution and temperature. Biol Symp. 1942;6: 71-125.
5. Fraïsse C, Elderfield JAD, Welch JJ. The genetics of speciation: are complex incompatibilities easier to evolve? J Evol Biol. 2014;27: 688-699. https://doi.org/10.1111/jeb.12339
6. Birchler JA, Yao H, Chudalayandi S, Vaiman D, Veitia RA. Heterosis. The Plant Cell. 2010;22: 2105-2112. https://doi.org/10.1105/tpc.110.076133
7. De Sanctis B, Schneemann H, Welch JJ. How does the mode of evolutionary divergence affect reproductive isolation? bioRxiv. 2022. 2022.03.08.483443 version 4. https://doi.org/10.1101/2022.03.08.483443 
8. Fisher RA. The genetical theory of natural selection. Oxford: The Clarendon Press; 1930. https://doi.org/10.5962/bhl.title.27468 
9. Tenaillon O. The Utility of Fisher's Geometric Model in Evolutionary Genetics. Annu Rev Ecol Evol Syst. 2014;45: 179-201. https://doi.org/10.1146/annurev-ecolsys-120213-091846
10. Barton NH. The role of hybridization in evolution. Molecular Ecology. 2001;10: 551-568. https://doi.org/10.1046/j.1365-294x.2001.01216.x 
11. Chevin L-M, Decorzent G, Lenormand T. Niche Dimensionality and The Genetics of Ecological Speciation. Evolution. 2014;68: 1244-1256. https://doi.org/10.1111/evo.12346 
12. Fraïsse C, Gunnarsson PA, Roze D, Bierne N, Welch JJ. The genetics of speciation: Insights from Fisher's geometric model. Evolution. 2016;70: 1450-1464. https://doi.org/10.1111/evo.12968

Evolutionary genomics, and population genetics in particular, aim to decipher the respective influence of neutral and selective forces shaping genetic polymorphism in a species/population. This is a much-needed requirement before scanning genome data for footprints of species adaptation to their biotic and abiotic environment (Johri et al. 2022). In general, we would like to quantify the proportion of the genome evolving neutrally and under selective (positive, balancing and negative) pressures (Kern and Hahn 2018, Johri et al. 2021). We thus need to understand patterns of linked selection along the genome, that is how the distribution of genetic polymorphisms is shaped by selected sites and the recombination landscape. The present contribution by Pavinato et al. (2022) provides an additional method in the population genomics toolbox to quantify the extent of linked positive and negative selection using temporal data.

The availability of genomics data for model and non-model species has led to improvement of the modeling framework for demography and selection (Johri et al. 2022), but also new inference methods making use of the full genome data based on the Sequential Markovian Coalescent (SMC, Li and Durbin 2011), Approximate Bayesian Computation (ABC, Jay et al. 2019), ABC and machine learning (Pudlo et al. 2016, Raynal et al. 2019) or Deep Learning (Sanchez et al. 2021). These methods are based on one sample in time and the use of the coalescent theory to reconstruct the past (demographic) history. However, it is also possible to obtain for many species temporal data sampled over several time points. For species with short generation time (in experimental evolution or monitored populations), one can sample a population every couple of generations as exemplified with Drosophila melanogaster (Bergland et al. 2010). For species with longer generation times that cannot be easily regularly sampled in time, it becomes possible to sequence available specimens from museums (e.g. Cridland et al. 2018) or ancient DNA samples. Methods using temporal data are based on the classical population genomics assumption that demography (migration, population subdivision, population size changes) leaves a genome-wide signal, while selection leaves a localized signal in the close vicinity of the causal mutation. Several methods do assess the demography of a population (change in effective population size, Ne, in time) using temporal data (e.g. Jorde and Ryman 2007) which can be used to calibrate the detection of loci under strong positive selection (Foll et al. 2014). Recently Buffalo and Coop (2020) used genome-wide covariance between allele frequency changes across time samples (and across replicates) to quantify the effects of linked selection over short timescales. 

In the present contribution, Pavinato et al. (2022) make use of temporal data to draw the joint estimation of demographic and selective parameters using a simulation-based method (ABC-Random Forests). This study by Pavinato et al. (2022) builds a framework allowing to infer the census size of the population in time (N) separately from the effect of genetic drift, which is determined by change in effective population size (Ne) in time, as well estimates of genome-wide parameters of selection. In a nutshell, the authors use a forward simulator and summarize genome data by genomic windows using classic statistics (nucleotide diversity, Tajima’s D, FST, heterozygosity) between time samples and for each sample. They specifically use the distributions (higher moments) of these statistics among all windows. The authors combine as input for the ABC-RF, vectors of summary statistics, model parameters and five latent variables: Ne, the ratio Ne/N, the number of beneficial mutations under strong selection, the average selection coefficient of strongly selected mutations, and the average substitution load. Indeed, the authors are interested in three different types of selection components: 1) the adaptive potential of a population which is estimated as the population mutation rate of beneficial mutations (θb), 2) the number of mutations under strong selection (irrespective of whether they reached fixation or not), and 3) the overall population fitness which is a function of the genetic load. In other words, the novelty of this method is not to focus on the detection of loci under selection, but to infer key parameters/distributions summarizing the genome-wide signal of demography and (positive and negative) selection. As a proof of principle, the authors then apply their method to a dataset of feral populations of honey bees (Apis mellifera) collected in California across many years and recovered from Museum samples (Cridland et al. 2018). The approach yields estimates of Ne which are on the same order of magnitude of previous estimates in hymenopterans, and the authors discuss why the different populations show various values of Ne and N which can be explained by different history of admixture with wild but also domesticated lineages of bees.

This study focuses on quantifying the genome-wide joint footprints of demography, and strong positive and negative selection to determine which proportion of the genome evolves neutrally or not. Further application of this method can be anticipated, for example, to study species with ecological and life-history traits which generate discrepancies between census size and Ne, for example for plants with selfing or seed banking (Sellinger et al. 2020), and for which the genome-wide effect of linked selection is not fully understood.

References

Johri P, Aquadro CF, Beaumont M, Charlesworth B, Excoffier L, Eyre-Walker A, Keightley PD, Lynch M, McVean G, Payseur BA, Pfeifer SP, Stephan W, Jensen JD (2022) Recommendations for improving statistical inference in population genomics. PLOS Biology, 20, e3001669. https://doi.org/10.1371/journal.pbio.3001669

Kern AD, Hahn MW (2018) The Neutral Theory in Light of Natural Selection. Molecular Biology and Evolution, 35, 1366–1371. https://doi.org/10.1093/molbev/msy092

Johri P, Riall K, Becher H, Excoffier L, Charlesworth B, Jensen JD (2021) The Impact of Purifying and Background Selection on the Inference of Population History: Problems and Prospects. Molecular Biology and Evolution, 38, 2986–3003. https://doi.org/10.1093/molbev/msab050

Pavinato VAC, Mita SD, Marin J-M, Navascués M de (2022) Joint inference of adaptive and demographic history from temporal population genomic data. bioRxiv, 2021.03.12.435133, ver. 6 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.03.12.435133

Li H, Durbin R (2011) Inference of human population history from individual whole-genome sequences. Nature, 475, 493–496. https://doi.org/10.1038/nature10231

Jay F, Boitard S, Austerlitz F (2019) An ABC Method for Whole-Genome Sequence Data: Inferring Paleolithic and Neolithic Human Expansions. Molecular Biology and Evolution, 36, 1565–1579. https://doi.org/10.1093/molbev/msz038

Pudlo P, Marin J-M, Estoup A, Cornuet J-M, Gautier M, Robert CP (2016) Reliable ABC model choice via random forests. Bioinformatics, 32, 859–866. https://doi.org/10.1093/bioinformatics/btv684

Raynal L, Marin J-M, Pudlo P, Ribatet M, Robert CP, Estoup A (2019) ABC random forests for Bayesian parameter inference. Bioinformatics, 35, 1720–1728. https://doi.org/10.1093/bioinformatics/bty867

Sanchez T, Cury J, Charpiat G, Jay F (2021) Deep learning for population size history inference: Design, comparison and combination with approximate Bayesian computation. Molecular Ecology Resources, 21, 2645–2660. https://doi.org/10.1111/1755-0998.13224

Bergland AO, Behrman EL, O’Brien KR, Schmidt PS, Petrov DA (2014) Genomic Evidence of Rapid and Stable Adaptive Oscillations over Seasonal Time Scales in Drosophila. PLOS Genetics, 10, e1004775. https://doi.org/10.1371/journal.pgen.1004775

Cridland JM, Ramirez SR, Dean CA, Sciligo A, Tsutsui ND (2018) Genome Sequencing of Museum Specimens Reveals Rapid Changes in the Genetic Composition of Honey Bees in California. Genome Biology and Evolution, 10, 458–472. https://doi.org/10.1093/gbe/evy007

Jorde PE, Ryman N (2007) Unbiased Estimator for Genetic Drift and Effective Population Size. Genetics, 177, 927–935. https://doi.org/10.1534/genetics.107.075481

Foll M, Shim H, Jensen JD (2015) WFABC: a Wright–Fisher ABC-based approach for inferring effective population sizes and selection coefficients from time-sampled data. Molecular Ecology Resources, 15, 87–98. https://doi.org/10.1111/1755-0998.12280

Buffalo V, Coop G (2020) Estimating the genome-wide contribution of selection to temporal allele frequency change. Proceedings of the National Academy of Sciences, 117, 20672–20680. https://doi.org/10.1073/pnas.1919039117

Sellinger TPP, Awad DA, Moest M, Tellier A (2020) Inference of past demography, dormancy and self-fertilization rates from whole genome sequence data. PLOS Genetics, 16, e1008698. https://doi.org/10.1371/journal.pgen.1008698