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

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

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

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

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

References

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

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

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

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

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

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

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

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

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

Flowering plants display a staggering diversity of both mating systems and life histories, ranging from almost exclusively selfers to obligate outcrossers, very short-lived annual herbs to super long lived trees. One pervasive pattern that has attracted considerable attention is the tight correlation that is found between mating systems and lifespan [1]. Until recently, theoretical explanations for these patterns have relied on static models exploring the consequences of several non-mutually exclusive important process: levels of inbreeding depression and ability to successfully were center stage. This make sense: successful colonization after long‐distance dispersal is far more likely to happen for self‐compatible than for self‐incompatible individuals in a sexually reproducing species. Furthermore, inbreeding depression (essentially a genetically driven phenomenon) and reproductive insurance are expected to shape the evolution of both mating system and lifespan.
But modelling jointly several processes and how their interplay to shape the evolution of a trait is challenging enough so models for the evolution of mating system tend invariably – for mathematical convenience and tractability – to fix lifespan [2].
However, comparative analysis of between species variations that map traits transitions among sister species in phylogenetic trees reveals a pervasive pattern: frequent transitions from a state outcrossing perennial to selfing annuals. This beg the question: is one transition triggering the other and if so, what comes first or are these transitions happening together? In this work, Lesaffre and Billiard use a very sophisticated machinery developed by Kirkpatrick et al. [3] and consider a general class of so-called modifiers models [4]. They study jointly the evolution of life span and mating system. They do so by using models where different life stages are tracked with life stage having some (fixed for once) amount of inbreeding depression. Their paper is technically demanding, mixing analytics and computer simulations, and along the way generates several important findings that are expected to stimulate further empirical and theoretical studies: (1) pure selfing versus pure outcrossing is the expected stable evolutionary outcomes (despite observation that mixed mating systems can be regularly met in nature), (2) increasing life-span drastically reduces the scope for the evolution of selfing, conversely (3) transition to selfing will also select for shorter life span as a way to mitigate the cumulative effects of inbreeding depression on adult life stages.
As usual there is room for future work, in particular the authors’ model assumes fixed inbreeding depression in the different life stages and this highlights the need for models that explore how inbreeding depression, a pivotal quantity in these models, can itself be molded by both mating system and lifespan. A third-generation of models should be “soon” on the way!

References
[1] Grossenbacher D, Briscoe Runquist R, Goldberg EE, and Brandvain Y. (2015) Geographic range size is predicted by plant mating system. Ecology Letters 18, 706–713. doi: 10.1111/ele.12449
[2] Morgan MT, Schoen DJ, and Bataillon T. (1997) The evolution of self-fertilization in perennials. The American Naturalist 150, 618–638. doi: 10.1086/286085
[3] Kirkpatrick M, Johnson T, and Barton N. (2002) General models of multilocus evolution. Genetics 161, 1727–1750.
[4] Lesaffre, T, and Billiard S. (2019) The joint evolution of lifespan and self-fertilisation. bioRxiv, 420877, ver. 3 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/420877

Heliconius butterflies are famous for their colorful wing patterns acting as a warning of their chemical defenses [1]. Most species are involved in Müllerian mimicry assemblies, as predators learn to associate common wing patterns with unpalatability and preferentially target rare variants. Such positive-frequency dependent selection homogenizes wing patterns at different localities, and in several species, all individuals within a community belong to the same morph [2]. In this respect, H. numata stands out. This species shows stable local polymorphism across multiple localities, with local populations home to up to seven distinct morphs [2]. Although a balance between migration and local positive-frequency dependent selection can allow some degree of local polymorphism, theory suggests that this occurs only when migration is within a narrow window [3].
One factor that potentially enhances local polymorphism in H. numata is disassortative mating. Mate choice assays have in fact revealed that females of this species tend to reject males with the same wing pattern [4]. However the evolution of such mating behavior and its effect on polymorphism remain unclear when selection is locally positive-frequency dependent. Using a mathematical model, Maisonneuve et al. [5] clarify the conditions that favor the evolution of disassortative mating in the complicated system of H. numata. In particular, they investigate whether the genetic basis of wing colour can favor the emergence of disassortative mating. Variation in wing pattern in H. numata is controlled by the supergene P, which is a single genomic region harboring multiple protein coding genes that have ceased to recombine due to chromosomal inversions [6]. If such remarkable genetic configuration allows for the co-adaptation of multiple loci participating to a complex phenotype such as wing color pattern, the absence of recombination can also result in the accumulation of deleterious mutations [7]. In fact, alleles at the P locus have been associated with a recessive genetic load, leading to a fitness advantage for heterozygotes at this locus [8]. Can this fitness advantage to heterozygotes lead to the evolution of disassortative mating? And if so, can such evolution lead to the maintenance of local polymorphism in spite of strong positive frequency-dependent selection?
To investigate these questions, Maisonneuve et al. [5] model evolution at two loci, one is the P locus for wing pattern, and the other influences mating behavior. The population is divided among two connected patches that differ in their butterfly communities, so that different alleles at the P locus are favored by positive frequency-dependent selection in different patches. The different alleles at the P locus are ordered in dominance relationships such that the most dominant over wing color pattern are also those with the highest load. By tracking the dynamics of haplotype frequencies in the population, the authors first show that disassortative mating readily evolves via the invasion of an allele causing females carrying it to reject males that resemble them phenotypically. Such “self-referencing” mechanism of mate choice, however, has never been reported and has been argued to be rare due to its complicated nature [9].
Maisonneuve et al. [5] then compare the evolution of disassortative mating via two alternative mechanisms: attraction and rejection. In these cases, alleles at the mating locus determine attraction to or rejection of specific phenotypes (e.g., under attraction rule, allele “B” encodes attraction to males with phenotype B). With the P and mating loci fully linked, disassortative mating can evolve under all three mechanisms (self-referencing, attraction and rejection), but tends to be less prevalent at equilibrium under attraction rule. This in turn results in the maintenance of less genetic variation under attraction compared to the other mating mechanisms. The loss of variation that occurs under attraction rules is due to a combination of dominance relationships between alleles at the P locus and the searching cost to females in finding rare types of males. When a particular wing pattern, say B, is only expressed in homozygotic form, B males are relatively rare. Females that carry the allele at the mating locus causing them to be attracted to such males then suffer a fitness cost due to lost mating opportunities. This mating allele is therefore purged, and in turn so is the recessive allele for B phenotype at the P locus. Under self-referencing and rejection rules, however, choosy females only reject males of a specific phenotype. They can therefore potentially mate with larger pool of males than females attracted to a single type. As a result, self-referencing and rejection rules are less sensitive to demographic effects and so are more conducive to disassortative mating evolution.
In their final analysis, Maisonneuve et al. [5] investigate the influence of recombination among the P and mating loci. They show that recombination has different effects on disassortative mating evolution depending on the mechanism of mate choice. Under the self-referencing rule, loose linkage leads to higher levels of disassortative mating and polymorphism than when linkage is tight. Under attraction or rejection rule, however, even very limited recombination completely inhibits the evolution of disassortative mating. This is because, with alleles at the mating locus coding for attraction/rejection to specific males, recombination breaks the association between the P and mating loci necessary for disassortative mating. By contrast, disassortative mating via a self-referencing rule does not depend on the linkage among the P and mating loci: females choose males that are different to themselves independently from the alleles they carry at the P locus.
Taken together, Maisonneuve et al.’s analyses [5] show that disassortative mating can readily evolve in a system like H. numata, but that this evolution depends on the genetic architecture of mating behavior. The architectures that are more conducive to the evolution of disassortative mating are: (1) epistatic interactions among the P and mating loci such that females are able to recognize their own phenotype and base their mating decision upon this information (self-referencing rule); and (2) full linkage among the P supergene and a mating locus that triggers rejection of a specific color pattern. While the mechanisms behind disassortative mating remain to be elucidated, assortative mating seems to rely on alleles triggering attraction to specific cues with variation in attraction and cues linked together [10]. These observations support the notion that disassortative mating is due to alleles causing rejection, in tight linkage to the P locus. If so, mating loci would in fact be part of the P supergene, thus controlling not only intricate wing color pattern but also mating behavior.
Beyond the specific system of H. numata, Maisonneuve et al.’s study [5] helps understand the evolution of disassortative mating and its association with the genetic architecture of correlated traits. In particular, Maisonneuve et al. [5] expands the role of supergenes for ecologically relevant traits to mating behavior, further bolstering the relevance of these remarkable genetic elements in the maintenance of variation in complex and elaborate phenotypes.

References

[1] Merrill, R M, K K Dasmahapatra, J W Davey, D D Dell'Aglio, J J Hanly, B Huber, C D Jiggins, et al. (2015). The Diversification of Heliconius butterflies: What Have We Learned in 150 Years? Journal of Evolutionary Biology 28 (8), 1417–38. https://doi.org/10.1111/jeb.12672.
[2] Joron M, IR Wynne, G Lamas, and J Mallet (1999). Variable selection and the coexistence of multiple mimetic forms of the butterfly Heliconius numata. Evolutionary Ecology 13, 721– 754. https://doi.org/10.1023/A:1010875213123
[3] Joron M and Y Iwasa (2005). The evolution of a Müllerian mimic in a spatially distributed community. Journal of Theoretical Biology 237, 87–103. https://doi.org/10.1016/j.jtbi.2005.04.005
[4] Chouteau M, V Llaurens, F Piron-Prunier, and M Joron (2017). Polymorphism at a mimicry su- pergene maintained by opposing frequency-dependent selection pressures. Proceedings of the National Academy of Sciences 114, 8325–8329. https://doi.org/10.1073/pnas.1702482114
[5] Maisonneuve, L, Chouteau, M, Joron, M and Llaurens, V. (2020). Evolution and genetic architecture of disassortative mating at a locus under heterozygote advantage. bioRxiv, 616409, ver. 9 peer-reviewed and recommended by PCI Evolutionary Biology. https://doi.org/10.1101/616409
[6] Joron M, L Frezal, RT Jones, NL Chamberlain, SF Lee, CR Haag, A Whibley, M Becuwe, SW Baxter, L Ferguson, et al. (2011). Chromosomal rearrangements maintain a polymorphic super- gene controlling butterfly mimicry. Nature 477, 203. https://doi.org/10.1038/nature10341
[7] Schwander T, R Libbrecht, and L Keller (2014). Supergenes and Complex Phenotypes.” Current Biology. 24 (7), 288–94. https://doi.org/10.1016/j.cub.2014.01.056.
[8] Jay P, M Chouteau, A Whibley, H Bastide, V Llaurens, H Parrinello, and M Joron (2019). Mutation accumulation in chromosomal inversions maintains wing pattern polymorphism in a butterfly. bioRxiv. https://doi.org/ 10.1101/736504.
[9] Kopp M, MR Servedio, TC Mendelson, RJ Safran, RL Rodrıguez, ME Hauber, EC Scordato, LB Symes, CN Balakrishnan, DM Zonana, et al. (2018). Mechanisms of assortative mating in speciation with gene flow: connecting theory and empirical research. The American Naturalist 191, 1–20. https://doi.org/10.1086/694889
[10] Merrill RM, P Rastas, SH Martin, MC Melo, S Barker, J Davey, WO McMillan, and CD Jiggins (2019). Genetic dissection of assortative mating behavior. PLoS biology 17, e2005902. https://doi.org/10.1371/journal.pbio.2005902

Although most mutations are deleterious, the strongly deleterious ones do not spread in a very large population as their chance of fixation is very small. Another mechanism via which the deleterious mutations can be eliminated is via recombination or sexual reproduction. However, in a finite asexual population, the subpopulation without any deleterious mutation will eventually acquire a deleterious mutation resulting in the reduction of the population size or in other words, an increase in the genetic drift. This, in turn, will lead the population to acquire deleterious mutations at a faster rate eventually leading to a mutational meltdown.

This irreversible (or, at least over some long time scales) accumulation of deleterious mutations is especially relevant to RNA viruses due to their high mutation rate, and while the prior work has dealt with bacteriophages and RNA viruses, the study by Lafforgue et al. [1] makes an interesting contribution to the existing literature by focusing on plants.

In this study, the authors enquire how despite the repeated increase in the strength of genetic drift, how the RNA viruses manage to survive in plants. Following a series of experiments and some numerical simulations, the authors find that as expected, after severe bottlenecks, the fitness of the population decreases significantly. But if the bottlenecks are followed by population expansion, the Muller’s ratchet can be halted due to the genetic diversity generated during population growth. They hypothesize this mechanism as a potential way by which the RNA viruses can survive the mutational meltdown.

As a theoretician, I find this investigation quite interesting and would like to see more studies addressing, e.g., the minimum population growth rate required to counter the potential extinction for a given bottleneck size and deleterious mutation rate. Of course, it would be interesting to see in future work if the hypothesis in this article can be tested in natural populations.

References

[1] Guillaume Lafforgue, Marie Lefebvre, Thierry Michon, Santiago F. Elena (2024) How do plant RNA viruses overcome the negative effect of Muller s ratchet despite strong transmission bottlenecks? bioRxiv, ver. 3 peer-reviewed and recommended by Peer Community In Evolutionary Biology
https://doi.org/10.1101/2023.08.01.550272

Parasitoid wasps have developed different mechanisms to increase their parasitic success, usually at the expense of host survival (Fellowes and Godfray, 2000). Eggs of these insects are deposited inside the juvenile stages of their hosts, which in turn deploy several immune response strategies to eliminate or disable them (Yang et al., 2020). Drosophila melanogaster protects itself against parasitoid attacks through the production of specific elongated haemocytes called lamellocytes which form a capsule around the invading parasite (Lavine and Strand, 2002; Rizki and Rizki, 1992) and the subsequent activation of the phenol-oxidase cascade leading to the release of toxic radicals (Nappi et al., 1995). On the parasitoid side, robust responses have evolved to evade host immune defenses as for example the Drosophila-specific endoparasite Leptopilina boulardi, which releases venom during oviposition that modifies host behaviour (Varaldi et al., 2006) and inhibits encapsulation (Gueguen et al., 2011; Martinez et al., 2012).
Studies have shown that the wasp parasitic capacity is correlated to venom presence and its content (Colinet et al., 2009; Poirié et al., 2014), including that evolution of venom protein composition is driven by different levels of host susceptibility to infection (Cavigliasso et al., 2019). However, it had not been determined to this day, if and how parasitic range can affect venom protein composition and to which extent host specialization requires broad-spectrum factors or a plethora of specialized components.
These outstanding questions are now approached in a study by Cavigliasso and colleagues (Cavigliasso et al., 2021), where they perform experimental evolution of L. boulardi for 9 generations exposing it to different Drosophila host species and genetic backgrounds (two strains of D. melanogaster, D. simulans and D. yakuba). The authors tested whether the parasitic success of each selection regime was host-specific and how they influenced venom composition in parasitoids. For the first part, infection outcomes were assayed for each selection regime when cross-infecting different hosts. To get a finer measurement of the mechanisms under selection, the authors differentiated three phenotypes: overall parasitic success, encapsulation inhibition and escape from capsule. Throughout the course of experimental evolution, only encapsulation inhibition did not show an improved response upon selection on any host. Importantly, the cross-infection scenario revealed a clear specificity to the selected host for each evolved resistance.
As for venom composition, a trend of differential evolution was detected between host species, although a significant part of that was due to a larger differentiation in the D. yakuba regime, which showed a completely different directionality. Importantly, the authors could identify some of the specific proteins targeted by the several selection regimes, whether selected or counter-selected for. Interestingly, the D. yakuba regime is the only case where the key parasitoid protein LbSPNy (Colinet et al., 2009) was not counter-selected and the only regime in which the overall venom composition did not evolve towards the Ism strain, one of the two ancestral strains of L. boulardi used in the study. It is possible that these two results are correlated, since LbSPNy has been described to inhibit activation of the phenoloxidase cascade in D. yakuba and is one of the most abundant proteins in the ISy venom, making it a good target for selection (Colinet et al., 2013). The authors also discuss the possibility that this difference is related to the geographical distribution of the strains of L. boulardi, since each coincide with either D. melanogaster or D. yakuba.
This methodologically broad work by Cavigliasso and colleagues constitutes an important experimental contribution towards the understanding of how parasitoid adaptation to specific hosts is achieved at different phenotypic and mechanistic levels. It provides compelling evidence that venom composition evolves differently in response to specific parasitic ranges, particularly considering the evolutionary difference between the selective hosts. In line with this result, it is also concluded that the majority of venom proteins selected are lineage-specific, although a few broad-spectrum factors could also be detected. 
The question of whether parasitic range can affect venom composition and parasitic success is still open to more contributions. A potentially interesting long-term direction will be to use a similar setup of experimental evolution on the generalist L. heterotoma (Schlenke et al., 2007) . On a more immediate horizon, comparing the venom evolution of both L. heterotoma and L. boulardi under selection with different hosts and under cross-infection scenarios could reveal interesting patterns. The recent sequencing of the L. boulardi genome together with the vast number of studies addressing mechanisms of Drosophila resistance to parasitoid infection, will enable the thorough characterization of the genetic basis of host-parasitoid interactions and the deeper understanding of these ubiquitous and economically-relevant relationships.
 
*This recommendation text has been co-written with Tânia F. Paulo who is not a recommender of PCI Evol Biol

References

Cavigliasso, F., Mathé-Hubert, H., Gatti, J.-L., Colinet, D. and Poirié, M. (2021) Parasitic success and venom composition evolve upon specialization of parasitoid wasps to different host species. bioRxiv, 2020.10.24.353417, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2020.10.24.353417

Cavigliasso, F., Mathé-Hubert, H., Kremmer, L., Rebuf, C., Gatti, J.-L., Malausa, T., Colinet, D., Poiré, M. and  Léne. (2019). Rapid and Differential Evolution of the Venom Composition of a Parasitoid Wasp Depending on the Host Strain. Toxins, 11(629). https://doi.org/10.3390/toxins11110629

Colinet, D., Deleury, E., Anselme, C., Cazes, D., Poulain, J., Azema-Dossat, C., Belghazi, M., Gatti, J. L. and  Poirié, M. (2013). Extensive inter- and intraspecific venom variation in closely related parasites targeting the same host: The case of Leptopilina parasitoids of Drosophila. Insect Biochemistry and Molecular Biology, 43(7), 601–611. https://doi.org/10.1016/j.ibmb.2013.03.010

Colinet, D., Dubuffet, A., Cazes, D., Moreau, S., Drezen, J. M. and  Poirié, M. (2009). A serpin from the parasitoid wasp Leptopilina boulardi targets the Drosophila phenoloxidase cascade. Developmental and Comparative Immunology, 33(5), 681–689. https://doi.org/10.1016/j.dci.2008.11.013

Fellowes, M. D. E. and  Godfray, H. C. J. (2000). The evolutionary ecology of resistance to parasitoids by Drosophila. Heredity, 84(1), 1–8. https://doi.org/10.1046/j.1365-2540.2000.00685.x

Gueguen, G., Rajwani, R., Paddibhatla, I., Morales, J. and  Govind, S. (2011). VLPs of Leptopilina boulardi share biogenesis and overall stellate morphology with VLPs of the heterotoma clade. Virus Research, 160(1–2), 159–165. https://doi.org/10.1016/j.virusres.2011.06.005

Lavine, M. D. and  Strand, M. R. (2002). Insect hemocytes and their role in immunity. Insect Biochemistry and Molecular Biology, 32(10), 1295–1309. https://doi.org/10.1016/S0965-1748(02)00092-9

Martinez, J., Duplouy, A., Woolfit, M., Vavre, F., O’Neill, S. L. and  Varaldi, J. (2012). Influence of the virus LbFV and of Wolbachia in a host-parasitoid interaction. PloS One, 7(4), e35081. https://doi.org/10.1371/journal.pone.0035081

Nappi, A. J., Vass, E., Frey, F. and  Carton, Y. (1995). Superoxide anion generation in Drosophila during melanotic encapsulation of parasites. European Journal of Cell Biology, 68(4), 450–456.

Poirié, M., Colinet, D. and  Gatti, J. L. (2014). Insights into function and evolution of parasitoid wasp venoms. Current Opinion in Insect Science, 6, 52–60. https://doi.org/10.1016/j.cois.2014.10.004

Rizki, T. M. and  Rizki, R. M. (1992). Lamellocyte differentiation in Drosophila larvae parasitized by Leptopilina. Developmental and Comparative Immunology, 16(2–3), 103–110. https://doi.org/10.1016/0145-305X(92)90011-Z

Schlenke, T. A., Morales, J., Govind, S. and  Clark, A. G. (2007). Contrasting infection strategies in generalist and specialist wasp parasitoids of Drosophila melanogaster. PLoS Pathogens, 3(10), 1486–1501. https://doi.org/10.1371/journal.ppat.0030158

Varaldi, J., Petit, S., Boulétreau, M. and  Fleury, F. (2006). The virus infecting the parasitoid Leptopilina boulardi exerts a specific action on superparasitism behaviour. Parasitology, 132(Pt 6), 747–756. https://doi.org/10.1017/S0031182006009930

Yang, L., Qiu, L., Fang, Q., Stanley, D. W. and  Gong‐Yin, Y. (2020). Cellular and humoral immune interactions between Drosophila and its parasitoids. Insect Science. https://doi.org/10.1111/1744-7917.12863