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

References

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

A classic problem in evolutionary biology is to understand the genetic variance in fitness. The simplest hypothesis is that variation exists, even in well-adapted populations, as a result of the balance between mutational input and selective elimination. This variation causes a reduction in mean fitness, known as the mutation load. Though mutation load is difficult to quantify empirically, indirect evidence of segregating genetic variation in fitness is often readily obtained by comparing the fitness of inbred and outbred offspring, i.e., by measuring inbreeding depression. Mutation-selection balance models have been studied as a means of understanding the genetic variance in fitness, mutation load, and inbreeding depression. Since their inception, such models have increased in sophistication, allowing us to ask these questions under more realistic and varied scenarios. The new theoretical work by Abu Awad and Roze [1] is a substantial step forward in understanding how arbitrary levels of self-fertilization affect variation, load and inbreeding depression under mutation-selection balance.
It has never been entirely clear how selfing should affect these population genetic properties in a multi-locus model. From the single-locus perspective, selfing increases homozygosity, which allows for more efficient purging leading to a prediction of less variance and lower load. On the other hand, selfing directly and indirectly affects several types of multilocus associations, which tend to make selection less efficient. Though this is certainly not the first study to consider mutation-selection balance in species with selfing (e.g., [2-5]), it is perhaps the most biologically realistic. The authors consider a model where n traits are under stabilizing selection and where each locus affects an arbitrary subset of these traits. As others have argued [6-7], this type of fitness landscape model “naturally” gives rise to dominance and epistatic effects. Abu Awad and Roze [1] thoroughly investigate this model both with analytical approximations and stochastic simulations (incorporating the effects of drift).
Their analysis reveals three major parameter regimes. The first regime occurs under low mutation rates, when segregating deleterious alleles are sufficiently rare across the genome that multi-locus genetic associations (disequilibria) can be ignored. As expected, in this regime, increased selfing facilitates purging, thereby leading to less standing genetic variation, lower load and less inbreeding depression.
In the second regime, mutation rates are higher and segregating deleterious alleles are more common. Though the effects of multilocus genetic associations cannot be ignored, Abu Awad and Roze [1] show that a good approximation can be obtained by considering only two-locus associations (ignoring the multitude of higher order associations). This is where the sophistication of their analysis yields the greatest insights. Their analysis shows that two different types of interlocus associations are important. First, selfing directly generates identity disequilibrium (correlation in homozygosity between two loci) that occurs because individuals produced through outbreeding tend to be heterozygous across multiple loci whereas individuals produced by selfing tend to be homozygous across multiple loci. These correlations reduce the efficiency of selection when deleterious effects are partially recessive [5]. Second, selfing indirectly affects traditional linkage disequilibrium. Epistatic selection resulting from the fitness landscape generates negative linkage disequilibrium between alleles at different loci that cause the same direction of deviation in a trait from its optimum. Because selfing reduces the effective rate of recombination, linkage disequilibrium reaches higher levels. Because selection tends to generate compensatory combinations of alleles, partially masking their deleterious effects, these associations also make purging less efficient. Their analysis shows the strength of the effect from identity disequilibrium scales with U, the genome-wide rate of deleterious mutations, but the effect of linkage disequilibrium scales with U/n because with more traits (higher n) two randomly chosen alleles are less likely to affect the same trait and so be subject to epistatic selection. Together, the effects of multilocus associations increase the load and can, in some cases, cause the load to increase as selfing increase from moderate to high levels.
However, their analytical approximations become inaccurate under conditions when the number of epistatically interacting segregating mutations (proportional to U/n) becomes large relative to the effective recombination rate (dependent on outcrossing and recombination rates). In this third regime, higher order genetic associations become important. In the limit of no recombination, model behaves as if the whole genome is a single locus with a very large number of alleles, becoming equivalent to previous studies [2-3].
The study by Abu Awad and Roze [1] helps us better understand the “simplest” explanation for genetic variance in fitness—mutation-selection balance—in a model of considerable complexity involving multiple traits under stabilizing selection, which ‘naturally’ allows for pleiotropy and epistasis. Their model tends to confirm the classic prediction of lower variation in fitness, less load, and inbreeding depression in species with higher levels of selfing. However, their careful analysis provides a clearer picture of how (and by how much) epistasis and selfing affect key population genetic properties.

References

[1] Abu Awad D and Roze D. 2017. Effects of partial selfing on the equilibrium genetic variance, mutation load and inbreeding depression under stabilizing selection. bioRxiv, 180000, ver. 4 of 17th November 2017. doi: 10.1101/180000

[2] Lande R. 1977. The influence of the mating system on the maintenance of genetic variability in polygenic characters. Genetics 86: 485–498.

[3] Charlesworth D and Charlesworth B. 1987. Inbreeding depression and its evolutionary consequences. Annual Review of Ecology and Systematics. 18: 237–268. doi: 10.1111/10.1146/annurev.es.18.110187.001321

[4] Lande R and Porcher E. 2015. Maintenance of quantitative genetic variance under partial self-fertilization, with implications for the evolution of selfing. Genetics 200: 891–906. doi: 10.1534/genetics.115.176693

[5] Roze D. 2015. Effects of interference between selected loci on the mutation load, inbreeding depression, and heterosis. Genetics 201: 745–757. doi: 10.1534/genetics.115.178533

[6] Martin G and Lenormand T. 2006. A general multivariate extension of Fisher's geometrical model and the distribution of mutation fitness effects across species. Evolution 60: 893–907. doi: 10.1111/j.0014-3820.2006.tb01169.x

[7] Martin G, Elena SF and Lenormand T. 2007. Distributions of epistasis in microbes fit predictions from a fitness landscape model. Nature Genetics 39: 555–560. doi: 10.1038/ng1998

Viruses establish intimate relationships with the cells they infect. The virocell is a novel entity, different from the original host cell and beyond the mere combination of viral and cellular genetic material. In these close encounters, viral and cellular genomes often hybridise, combine, recombine, merge and excise. Such chemical promiscuity leaves genomics scars that can be passed on to descent, in the form of deletions or duplications and, importantly, insertions and back and forth exchange of genetic material between viruses and their hosts.
In this preprint [1], Di Giovanni and coworkers report the identification of 13 genes present in the extant genomes of members of the Leptopilina wasp genus, bearing sound signatures of having been horizontally acquired from an ancestral virus. Importantly the authors identify Leptopilina boulardi filamentous virus (LbFV) as an extant relative of the ancestral virus that served as donor for the thirteen horizontally transferred genes. While pinpointing genes with a likely possible viral origin in eukaryotic genomes is only relatively rare, identifying an extant viral lineage related to the ancestral virus that continues to infect an extant relative of the ancestral host is remarkable. But the amazing evolutionary history of the Leptopilina hosts and these filamentous viruses goes beyond this shared genes. These wasps are endoparasitoids of Drosophila larvae, the female wasp laying the eggs inside the larvae and simultaneously injecting venom that hinders the immune response. The composition of the venoms is complex, varies between wasp species and also between individuals within a species, but a central component of all these venoms are spiked structures that vary in morphology, symmetry and size, often referred to as virus-like particles (VLPs).
In this preprint, the authors convincingly show that the expression pattern in the Leptopilina wasps of the thirteen genes identified to have been horizontally acquired from the LbFV ancestor coincides with that of the production of VLPs in the female wasp venom gland. Based on this spatio-temporal match, the authors propose that these VLPs have a viral origin. The data presented in this preprint will undoubtedly stimulate further research on the composition, function, origin, evolution and diversity of these VLP structures, which are highly debated (see for instance [2] and [3]).

References

[1] Di Giovanni, D., Lepetit, D., Boulesteix, M., Ravallec, M., & Varaldi, J. (2018). A behavior-manipulating virus relative as a source of adaptive genes for parasitoid wasps. bioRxiv, 342758, ver. 5 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/342758
[2] Poirié, M., Colinet, D., & Gatti, J. L. (2014). Insights into function and evolution of parasitoid wasp venoms. Current Opinion in Insect Science, 6, 52-60. doi: 10.1016/j.cois.2014.10.004
[3] Heavner, M. E., Ramroop, J., Gueguen, G., Ramrattan, G., Dolios, G., Scarpati, M., ... & Govind, S. (2017). Novel organelles with elements of bacterial and eukaryotic secretion systems weaponize parasites of Drosophila. Current Biology, 27(18), 2869-2877. doi: 10.1016/j.cub.2017.08.019

Monnin et al. [1] here studied how Drosophila populations are affected when exposed to a high virulent endosymbiotic wMelPop Wolbachia strain and why virulent vertically transmitting endosymbionts persist in nature. This virulent wMelPop strain has been described to be a blocker of dengue and other arboviral infections in arthropod vector species, such as Aedes aegypti. Whereas it can thus function as a mutualistic symbiont, it here acts as an antagonist along the mutualism-antagonism continuum symbionts operate. The wMelPop strain is not a natural occurring strain in Drosophila melanogaster and thus the start of this experiment can be seen as a novel host-pathogen association. Through experimental evolution of 17 generations, the authors studied how high temperature affects wMelPop Wolbachia virulence and Drosophila melanogaster survival. The authors used Drosophila strains that were selected for late reproduction, given that this should favor evolution to a lower virulence. Assumptions for this hypothesis are not given in the manuscript here, but it can indeed be assumed that energy that is assimilated to symbiont tolerance instead of reproduction may lead to reduced virulence evolution. This has equally been suggested by Reyserhove et al. [2] in a dynamics energy budget model tailored to Daphnia magna virulence evolution upon a viral infection causing White fat Cell disease, reconstructing changing environments through time.
Contrary to their expectations for vertically transmitting symbionts, the authors did not find a reduction in wMelPop Wolbachia virulence during the course of the experimental evolution experiment under high temperature. Important is what this learns for virulence evolution, also for currently horizontal transmitting disease epidemics (such as COVID-19). It mainly reflects that evolution of virulence for new host-pathogen associations is difficult to predict and that it may take multiple generations before optimal levels of virulence are reached [3,4]. These optimal levels of virulence will depend on trade-offs with other life history traits of the symbiont, but also on host demography, host heterogeneity, amongst others [5,6]. Multiple microbial interactions may affect the outcome of virulence evolution [7]. Given that no germ-free individuals were used, it can be expected that other components of the Drosophila microbiome may have played a role in the virulence evolution. In most cases, microbiota have been described as defensive or protective for virulent symbionts [8], but they may also have stimulated the high levels of virulence. Especially, given that upon higher temperatures, Wolbachia growth may have been increased, host metabolic demands increased [9], host immune responses affected and microbial communities changed [10]. This may have resulted in increased competitive interactions to retrieve host resources, sustaining high virulence levels of the symbiont.
A nice asset of this study is that the phenotypic results obtained in the experimental evolution set-up were linked with wMelPop density measurement and octomom copy number quantifications. Octomom is a specific 8-n genes region of the Wolbachia genome responsible for wMelPop virulence, so there is a link between the phenotypic and molecular functions of the involved symbiont. The authors found that density, octomom copy number and virulence were correlated to each other. An important note the authors address in their discussion is that, to exclude the possibility that octomom copy number has an effect on density, and density on virulence, the effect of these variables should be assessed independently of temperature and age. The obtained results are a valuable contribution to the ongoing debate on the relationship between wMelPop octomom copy number, density and virulence.

References

[1] Monnin, D., Kremer, N., Michaud, C., Villa, M., Henri, H., Desouhant, E. and Vavre, F. (2020) Experimental evolution of virulence and associated traits in a Drosophila melanogaster – Wolbachia symbiosis. bioRxiv, 2020.04.26.062265, ver. 4 peer-reviewed and recommended by PCI Evol Biol. doi: https://doi.org/10.1101/2020.04.26.062265
[2] Reyserhove, L., Samaey, G., Muylaert, K., Coppé, V., Van Colen, W., and Decaestecker, E. (2017). A historical perspective of nutrient change impact on an infectious disease in Daphnia. Ecology, 98(11), 2784-2798. doi: https://doi.org/10.1002/ecy.1994
[3] Ebert, D., and Bull, J. J. (2003). Challenging the trade-off model for the evolution of virulence: is virulence management feasible?. Trends in microbiology, 11(1), 15-20. doi: https://doi.org/10.1016/S0966-842X(02)00003-3
[4] Houwenhuyse, S., Macke, E., Reyserhove, L., Bulteel, L., and Decaestecker, E. (2018). Back to the future in a petri dish: Origin and impact of resurrected microbes in natural populations. Evolutionary Applications, 11(1), 29-41. doi: https://doi.org/10.1111/eva.12538
[5] Day, T., and Gandon, S. (2007). Applying population‐genetic models in theoretical evolutionary epidemiology. Ecology Letters, 10(10), 876-888. doi: https://doi.org/10.1111/j.1461-0248.2007.01091.x
[6] Alizon, S., Hurford, A., Mideo, N., and Van Baalen, M. (2009). Virulence evolution and the trade‐off hypothesis: history, current state of affairs and the future. Journal of evolutionary biology, 22(2), 245-259. doi: https://doi.org/10.1111/j.1420-9101.2008.01658.x
[7] Alizon, S., de Roode, J. C., and Michalakis, Y. (2013). Multiple infections and the evolution of virulence. Ecology letters, 16(4), 556-567. doi: https://doi.org/10.1111/ele.12076
[8] Decaestecker, E., and King, K. (2019). Red queen dynamics. Reference module in earth systems and environmental sciences, 3, 185-192. doi: https://doi.org/10.1016/B978-0-12-409548-9.10550-0
[9] Kirk, D., Jones, N., Peacock, S., Phillips, J., Molnár, P. K., Krkošek, M., and Luijckx, P. (2018). Empirical evidence that metabolic theory describes the temperature dependency of within-host parasite dynamics. PLoS biology, 16(2), e2004608. doi: https://doi.org/10.1371/journal.pbio.2004608
[10] Frankel-Bricker, J., Song, M. J., Benner, M. J., and Schaack, S. (2019). Variation in the microbiota associated with Daphnia magna across genotypes, populations, and temperature. Microbial ecology, 1-12. doi: https://doi.org/10.1007/s00248-019-01412-9

Occasional events of asexual reproduction in otherwise sexual taxa have been documented since a long time. Accounts range from observations of offspring development from unfertilized eggs in Drosophila to rare offspring production by isolated females in lizards and birds (e.g., Stalker 1954, Watts et al 2006, Ryder et al. 2021). Many more such cases likely await documentation, as rare events are inherently difficult to observe. These rare events of asexual reproduction are often associated with low offspring fitness (“tychoparthenogenesis”), and have mostly been discarded in the evolutionary literature as reproductive accidents without evolutionary significance. Recently, however, there has been an increased interest in the details of evolutionary transitions from sexual to asexual reproduction (e.g., Archetti 2010, Neiman et al.2014, Lenormand et al. 2016), because these details may be key to understanding why successful transitions are rare, why they occur more frequently in some groups than in others, and why certain genetic mechanisms of ploidy maintenance or ploidy restoration are more often observed than others. In this context, the hypothesis has been formulated that regular or even obligate asexual reproduction may evolve from these rare events of asexual reproduction (e.g., Schwander et al. 2010).

A new study by Capdevielle Dulac et al. (2022) now investigates this question in a parasitoid wasp, highlighting also the fact that what is considered rare or occasional may differ from one system to the next. The results show “rare” parthenogenetic production of diploid daughters occurring at variable frequencies (from zero to 2 %) in different laboratory strains, as well as in a natural population. They also demonstrate parthenogenetic production of female offspring in both virgin females and mated ones, as well as no reduced fecundity of parthenogenetically produced offspring. These findings suggest that parthenogenetic production of daughters, while still being rare, may be a more regular and less deleterious reproductive feature in this species than in other cases of occasional asexuality. Indeed, haplodiploid organisms, such as this parasitoid wasp have been hypothesized to facilitate evolutionary transitions to asexuality (Neimann et al. 2014, Van Der Kooi et al. 2017). First, in haploidiploid organisms, females are diploid and develop from normal, fertilized eggs, but males are haploid as they develop parthenogenetically from unfertilized eggs. This means that, in these species, fertilization is not necessarily needed to trigger development, thus removing one of the constraints for transitions to obligate asexuality (Engelstädter 2008, Vorburger 2014). Second, spermatogenesis in males occurs by a modified meiosis that skips the first meiotic division (e.g., Ferree et al. 2019). Haploidiploid organisms may thus have a potential route for an evolutionary transition to obligate parthenogenesis that is not available to organisms: The pathways for the modified meiosis may be re-used for oogenesis, which might result in unreduced, diploid eggs. Third, the particular species studied here regularly undergoes inbreeding by brother-sister mating within their hosts. Homozygosity, including at the sex determination locus (Engelstädter 2008), is therefore expected to have less negative effects in this species compared to many other, non-inbreeding haplodipoids (see also Little et al. 2017). This particular species may therefore be less affected by loss of heterozygosity, which occurs in a fashion similar to self-fertilization under many forms of non-clonal parthenogenesis. 

Indeed, the study also addresses the mechanisms underlying parthenogenesis in the species. Surprisingly, the authors find that parthenogenetically produced females are likely produced by two distinct genetic mechanisms. The first results in clonality (maintenance of the maternal genotype), whereas the second one results in a loss of heterozygosity towards the telomeres, likely due to crossovers occurring between the centromeres and the telomeres. Moreover, bacterial infections appear to affect the propensity of parthenogenesis but are unlikely the primary cause. Together, the finding suggests that parthenogenesis is a variable trait in the species, both in terms of frequency and mechanisms. It is not entirely clear to what degree this variation is heritable, but if it is, then these results constitute evidence for low-frequency existence of variable and heritable parthenogenesis phenotypes, that is, the raw material from which evolutionary transitions to more regular forms of parthenogenesis may occur.

References

Archetti M (2010) Complementation, Genetic Conflict, and the Evolution of Sex and Recombination. Journal of Heredity, 101, S21–S33. https://doi.org/10.1093/jhered/esq009

Capdevielle Dulac C, Benoist R, Paquet S, Calatayud P-A, Obonyo J, Kaiser L, Mougel F (2022) Spontaneous parthenogenesis in the parasitoid wasp Cotesia typhae: low frequency anomaly or evolving process? bioRxiv, 2021.12.13.472356, ver. 6 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.12.13.472356

Engelstädter J (2008) Constraints on the evolution of asexual reproduction. BioEssays, 30, 1138–1150. https://doi.org/10.1002/bies.20833

Ferree PM, Aldrich JC, Jing XA, Norwood CT, Van Schaick MR, Cheema MS, Ausió J, Gowen BE (2019) Spermatogenesis in haploid males of the jewel wasp Nasonia vitripennis. Scientific Reports, 9, 12194. https://doi.org/10.1038/s41598-019-48332-9

van der Kooi CJ, Matthey-Doret C, Schwander T (2017) Evolution and comparative ecology of parthenogenesis in haplodiploid arthropods. Evolution Letters, 1, 304–316. https://doi.org/10.1002/evl3.30

Lenormand T, Engelstädter J, Johnston SE, Wijnker E, Haag CR (2016) Evolutionary mysteries in meiosis. Philosophical Transactions of the Royal Society B: Biological Sciences, 371, 20160001. https://doi.org/10.1098/rstb.2016.0001

Little CJ, Chapuis M-P, Blondin L, Chapuis E, Jourdan-Pineau H (2017) Exploring the relationship between tychoparthenogenesis and inbreeding depression in the Desert Locust, Schistocerca gregaria. Ecology and Evolution, 7, 6003–6011. https://doi.org/10.1002/ece3.3103

Neiman M, Sharbel TF, Schwander T (2014) Genetic causes of transitions from sexual reproduction to asexuality in plants and animals. Journal of Evolutionary Biology, 27, 1346–1359. https://doi.org/10.1111/jeb.12357

Ryder OA, Thomas S, Judson JM, Romanov MN, Dandekar S, Papp JC, Sidak-Loftis LC, Walker K, Stalis IH, Mace M, Steiner CC, Chemnick LG (2021) Facultative Parthenogenesis in California Condors. Journal of Heredity, 112, 569–574. https://doi.org/10.1093/jhered/esab052

Schwander T, Vuilleumier S, Dubman J, Crespi BJ (2010) Positive feedback in the transition from sexual reproduction to parthenogenesis. Proceedings of the Royal Society B: Biological Sciences, 277, 1435–1442. https://doi.org/10.1098/rspb.2009.2113

Stalker HD (1954) Parthenogenesis in Drosophila. Genetics, 39, 4–34. https://doi.org/10.1093/genetics/39.1.4

Vorburger C (2014) Thelytoky and Sex Determination in the Hymenoptera: Mutual Constraints. Sexual Development, 8, 50–58. https://doi.org/10.1159/000356508

Watts PC, Buley KR, Sanderson S, Boardman W, Ciofi C, Gibson R (2006) Parthenogenesis in Komodo dragons. Nature, 444, 1021–1022. https://doi.org/10.1038/4441021a