Quantitative trait loci involved in the reproductive success of a parasitoid wasp

Hot and cold waves decrease sperm production and bias sex ratio in the parasitoid wasp Cotesia typhae (Hymenoptera, Braconidae)

Parasitoid wasps are haplodiploid, meaning that sperm stored by egg laying females are only used to produce daughters. Thus, the sex ratio of the offspring depends on the availability of sperm after mating. In these insects, males are sensitive to temperature at the pupal stage. This stress leads to subfertility due to a drastic reduction in the number of sperm produced and transferred to females. Experiments were conducted under controlled conditions on the parasitoid wasp Cotesia typhae (Hymenoptera, Braconidae), a natural enemy of the invading pest Sesamia nonagrioides (Lepidoptera, Noctuidae). At 25-27 °C, sperm production was measured for 7 days, and found to reach a plateau at the third day of adult life. It leads to a final amount around 25,000 sperm per male. A male can successfully inseminate at least 10 females, producing predominantly female offspring. Sperm production decreased significantly after 1 day of pupal exposure to heat at 34 or 36 °C and 7 days of cold at 0, 5 or 10 °C. This highlights that both cold and heat are stressful. After mating with one male treated at 10 or 34 °C, females store fewer sperm than the control, and produce fewer daughters. The sex ratio of the offspring is male biased when males experienced temperature stresses during development, like other parasitoid wasps. In the field, C. typhae populations would be affected by heat and cold, at least at the pupal stage. This lowers overwintering risk in case this biological agent was introduced in Europe. This risk is both economical, as companies seek to establish costly continuous production to sell beneficial insects, and ecological as the introduced population would not settle in the ecosystem. Lastly, the transport and storage of this insect of agronomic interest would need to consider temperature variations to ensure successful application.

Evolution of parasitoid host preference and performance in response to an invasive host acting as evolutionary trap

The invasion of a novel host species can create a mismatch in host choice and offspring survival (performance) when native parasitoids attempt to exploit the invasive host without being able to circumvent its resistance mechanisms. Invasive hosts can therefore act as evolutionary trap reducing parasitoids' fitness and this may eventually lead to their extinction. Yet, escape from the trap can occur when parasitoids evolve behavioral avoidance or a physiological strategy compatible with the trap host, resulting in either host‐range expansion or a complete host‐shift. We developed an individual based model to investigate which conditions promote parasitoids to evolve behavioral preference that matches their performance, including host‐trap avoidance, and which conditions lead to adaptations to the unsuitable hosts. The model was inspired by solitary endo‐parasitoids attacking larval host stages. One important aspect of these conditions was reduced host survival during incompatible interaction, where a failed parasitization attempt by a parasitoid resulted not only in death of her offspring but also in host killing. This non‐reproductive host mortality had a strong influence on the likelihood of establishment of novel host–parasitoid relationship, in some cases constraining adaptation to the trap host species. Moreover, our model revealed that host‐search efficiency and genetic variation in host‐preference play a key role in the likelihood that parasitoids will include the suboptimal host in their host range, or will evolve behavioral avoidance resulting in specialization and host‐range conservation, respectively. Hence, invasive species might change the evolutionary trajectory of native parasitoid species, which is important for predicting biocontrol ability of native parasitoids towards novel hosts.

Assemblies of the genomes of parasitic wasps using meta-assembly and scaffolding with genetic linkage

Safe, effective biological-control introductions against invasive pests depend on narrowly host-specific natural enemies with the ability to adapt to a changing environment. As part of a project on the genetic architectures of these traits, we assembled and annotated the genomes of two aphid parasitoids, Aphelinus atriplicis and Aphelinus certus. We report here several assemblies of A. atriplicis made with Illumina and PacBio data, which we combined into a meta-assembly. We scaffolded the meta-assembly with markers from a genetic map of hybrids between A. atriplicis and A. certus. We used this genetic-linkage scaffolded (GLS) assembly of A. atriplicis to scaffold a de novo assembly of A. certus. The de novo assemblies of A. atriplicis differed in contiguity, and the meta-assembly of these assemblies was more contiguous than the best de novo assembly. Scaffolding with genetic-linkage data allowed chromosomal-level assembly of the A. atriplicis genome and scaffolding a de novo assembly of A. certus with this GLS assembly, greatly increased the contiguity of the A. certus assembly to the point where it was also at the chromosomal-level. However, completeness of the A. atriplicis assembly, as measured by percent complete, single-copy BUSCO hymenopteran genes, varied little among de novo assemblies and was not increased by meta-assembly or genetic scaffolding. Furthermore, the greater contiguity of the meta-assembly and GLS assembly had little or no effect on the numbers of genes identified, the proportions with homologs or functional annotations. Increased contiguity of the A. certus assembly provided modest improvement in assembly completeness, as measured by percent complete, single-copy BUSCO hymenopteran genes. The total genic sequence increased, and while the number of genes declined, gene length increased, which together suggest greater accuracy of gene models. More contiguous assemblies provide uses other than gene annotation, for example, identifying the genes associated with quantitative trait loci and understanding of chromosomal rearrangements associated with speciation.

Genome-Wide Patterns of Bracovirus Chromosomal Integration into Multiple Host Tissues during Parasitism

Bracoviruses are domesticated viruses found in parasitic wasp genomes. They are composed of genes of nudiviral origin that are involved in particle production and proviral segments containing virulence genes that are necessary for parasitism success. During particle production, proviral segments are amplified and individually packaged as DNA circles in nucleocapsids. These particles are injected by parasitic wasps into host larvae together with their eggs. Bracovirus circles of two wasp species were reported to undergo chromosomal integration in parasitized host hemocytes, through a conserved sequence named the host integration motif (HIM). Here, we used bulk Illumina sequencing to survey integrations of Cotesia typhae bracovirus circles in the DNA of its host, the maize corn borer (Sesamia nonagrioides), 7 days after parasitism. First, assembly and annotation of a high-quality genome for C. typhae enabled us to characterize 27 proviral segments clustered in proviral loci. Using these data, we characterized large numbers of chromosomal integrations (from 12 to 85 events per host haploid genome) for all 16 bracovirus circles containing a HIM. Integrations were found in four S. nonagrioides tissues and in the body of a caterpillar in which parasitism had failed. The 12 remaining circles do not integrate but are maintained at high levels in host tissues. Surprisingly, we found that HIM-mediated chromosomal integration in the wasp germ line has occurred accidentally at least six times during evolution. Overall, our study furthers our understanding of wasp-host genome interactions and supports HIM-mediated chromosomal integration as a possible mechanism of horizontal transfer from wasps to their hosts. IMPORTANCE Bracoviruses are endogenous domesticated viruses of parasitoid wasps that are injected together with wasp eggs into wasp host larvae during parasitism. Several studies have shown that some DNA circles packaged into bracovirus particles become integrated into host somatic genomes during parasitism, but the phenomenon has never been studied using nontargeted approaches. Here, we use bulk Illumina sequencing to systematically characterize and quantify bracovirus circle integrations that occur in four tissues of the Mediterranean corn borer (Sesamia nonagrioides) during parasitism by the Cotesia typhae wasp. Our analysis reveals that all circles containing a HIM integrate at substantial levels (from 12 to 85 integrations per host cell, in total) in all tissues, while other circles do not integrate. In addition to shedding new light on wasp-bracovirus-host interactions, our study supports HIM-mediated chromosomal integration of bracovirus as a possible source of wasp-to-host horizontal transfer, with long-term evolutionary consequences.

Quantitative trait locus analysis of parasitoid counteradaptation to symbiont-conferred resistance

Insect hosts and parasitoids are engaged in an intense struggle of antagonistic coevolution. Infection with heritable bacterial endosymbionts can substantially increase the resistance of aphids to parasitoid wasps, which exerts selection on parasitoids to overcome this symbiont-conferred protection (counteradaptation). Experimental evolution in the laboratory has produced counteradapted populations of the parasitoid wasp Lysiphlebus fabarum. These populations can parasitize black bean aphids (Aphis fabae) protected by the bacterial endosymbiont Hamiltonella defensa, which confers high resistance against L. fabarum. We used two experimentally evolved parasitoid populations to study the genetic architecture of the counteradaptation to symbiont-conferred resistance by QTL analysis. With simple crossing experiments, we showed that the counteradaptation is a recessive trait depending on the maternal genotype. Based on these results, we designed a customized crossing scheme to genotype a mapping population phenotyped for the ability to parasitize Hamiltonella-protected aphids. Using 1835 SNP markers obtained by ddRAD sequencing, we constructed a high-density linkage map consisting of six linkage groups (LGs) with an overall length of 828.3 cM and an average marker spacing of 0.45 cM. We identified a single QTL associated with the counteradaptation to Hamiltonella in L. fabarum on linkage group 2. Out of 120 genes located in this QTL, several genes encoding putative venoms may represent candidates for counteradaptation, as parasitoid wasps inject venoms into their hosts during oviposition.