A new approach for investigating venom function applied to venom calreticulin in a parasitoid wasp

Transcriptome Analysis Reveals the Venom Genes of the Ectoparasitoid Habrobracon hebetor (Hymenoptera: Braconidae)

The ectoparasitoid Habrobracon hebetor (Hymenoptera: Braconidae) exhibits a broad parasitic capability towards various lepidopteran pests, with venom serving as a crucial virulent factor ensuring successful parasitization and subsequent host mortality. Analyzing the constituents of its venom is essential for elucidating the mechanisms underlying efficient host killing by this parasitoid and for exploring potentially functional venom proteins. Through a transcriptomic analysis, a total of 34 venom proteins were identified within the venom of H. hebetor, encompassing known components such as serine protease, metalloproteinase, esterase, and serine protease inhibitors commonly present in parasitoid venoms. Unique components like paralytic protein and ion transport peptide-like were identified, possibly specific to certain parasitoids, along with novel proteins with uncharacterized functions. Spatial gene expression profiling of the identified venom proteins using transcriptomic data, corroborated by quantitative PCR validation for 13 randomly selected proteins, revealed abundant expression levels in the venom apparatus, affirming them as genuine venom components. Notably, the paralytic protein exhibited prominent expression, with the highest FPKM (fragments per kilobase of transcript per million fragments mapped) value of 24,704.87 in the venom apparatus, indicative of its significant role in successful parasitism by H. hebetor. The identification of these venom proteins establishes a foundation for the further exploration of bioactive agents for pest management strategies.

Immune system modulation & virus transmission during parasitism identified by multi-species transcriptomics of a declining insect biocontrol system

BACKGROUND

The Argentine stem weevil (ASW, Listronotus bonariensis) is a significant pasture pest in Aotearoa New Zealand, primarily controlled by the parasitoid biocontrol agent Microctonus hyperodae. Despite providing effective control of ASW soon after release, M. hyperodae parasitism rates have since declined significantly, with ASW hypothesised to have evolved resistance to its biocontrol agent. While the parasitism arsenal of M. hyperodae has previously been investigated, revealing many venom components and an exogenous novel DNA virus Microctonus hyperodae filamentous virus (MhFV), the effects of said arsenal on gene expression in ASW during parasitism have not been examined. In this study, we performed a multi-species transcriptomic analysis to investigate the biology of ASW parasitism by M. hyperodae, as well as the decline in efficacy of this biocontrol system.

RESULTS

The transcriptomic response of ASW to parasitism by M. hyperodae involves modulation of the weevil's innate immune system, flight muscle components, and lipid and glucose metabolism. The multispecies approach also revealed continued expression of venom components in parasitised ASW, as well as the transmission of MhFV to weevils during parasitism and some interrupted parasitism attempts. Transcriptomics did not detect a clear indication of parasitoid avoidance or other mechanisms to explain biocontrol decline.

CONCLUSIONS

This study has expanded our understanding of interactions between M. hyperodae and ASW in a biocontrol system of critical importance to Aotearoa-New Zealand's agricultural economy. Transmission of MhFV to ASW during successful and interrupted parasitism attempts may link to a premature mortality phenomenon in ASW, hypothesised to be a result of a toxin-antitoxin system. Further research into MhFV and its potential role in ASW premature mortality is required to explore whether manipulation of this viral infection has the potential to increase biocontrol efficacy in future.

A serpin gene from a parasitoid wasp disrupts host immunity and exhibits adaptive alternative splicing

Alternative splicing (AS) is a major source of protein diversity in eukaryotes, but less is known about its evolution compared to gene duplication (GD). How AS and GD interact is also largely understudied. By constructing the evolutionary trajectory of the serpin gene PpSerpin-1 (Pteromalus puparum serpin 1) in parasitoids and other insects, we found that both AS and GD jointly contribute to serpin protein diversity. These two processes are negatively correlated and show divergent features in both protein and regulatory sequences. Parasitoid wasps exhibit higher numbers of serpin protein/domains than nonparasitoids, resulting from more GD but less AS in parasitoids. The potential roles of AS and GD in the evolution of parasitoid host-effector genes are discussed. Furthermore, we find that PpSerpin-1 shows an exon expansion of AS compared to other parasitoids, and that several isoforms are involved in the wasp immune response, have been recruited to both wasp venom and larval saliva, and suppress host immunity. Overall, our study provides an example of how a parasitoid serpin gene adapts to parasitism through AS, and sheds light on the differential features of AS and GD in the evolution of insect serpins and their associations with the parasitic life strategy.

Aphidius ervi venom regulates Buchnera contribution to host nutritional suitability

The association between the pea aphid, Acyrthosiphon pisum (Harris) (Homoptera: Aphididae), and the endophagous parasitoid wasp Aphidius ervi Haliday (Hymenoptera: Braconidae) offers a unique model system for studying the molecular mechanisms underlying the complex interactions between the parasitoid, its host and the associated primary symbiont. Here, we investigate in vivo the functional role of the most abundant component of A. ervi venom, Ae-γ-glutamyl transpeptidase (Ae-γ-GT), which is known to induce host castration. Microinjections of double-stranded RNA into A. ervi pupae stably knocked down Ae-γ-GT1 and Ae-γ-GT2 paralogue genes in newly emerged females. These females were used to score the phenotypic changes both in parasitized hosts and in the parasitoid's progeny, as affected by a venom blend lacking Ae-γ-GT. Ae-γ-GT gene silencing enhanced growth both of host and parasitoid, supported by a higher load of the primary bacterial symbiont Buchnera aphidicola. Emerging adults showed a reduced survival and fecundity, suggesting a trade-off with body size. This demonstrates in vivo the primary role of Ae-γ-GT in host ovary degeneration and suggests that this protein counterbalances the proliferation of Buchnera likely triggered by other venom components. Our study provides a new approach to unravelling the complexity of aphid parasitoid venom in vivo, and sheds light on a novel role for Ae-γ-GT in host regulation.

The venom composition and parthenogenesis mechanism of the parasitoid wasp Microctonus hyperodae, a declining biocontrol agent

A biocontrol system in New Zealand using the endoparasitoid Microctonus hyperodae is failing, despite once being one of the most successful examples of classical biocontrol worldwide. Though it is of significant economic importance as a control agent, little is known about the genetics of M. hyperodae. In this study, RNA-seq was used to characterise two key traits of M. hyperodae in this system, the venom, critical for the initial success of biocontrol, and the asexual reproduction mode, which influenced biocontrol decline. Expanded characterisation of M. hyperodae venom revealed candidates involved in manipulating the host environment to source nutrition for the parasitoid egg, preventing a host immune response against the egg, as well as two components that may stimulate the host's innate immune system. Notably lacking from the venom-specific expression list was calreticulin, as it also had high expression in the ovaries. In-situ hybridisation revealed this ovarian expression was localised to the follicle cells, which may result in the deposition of calreticulin into the egg exochorion. Investigating the asexual reproduction of M. hyperodae revealed core meiosis-specific genes had conserved expression patterns with the highest expression in the ovaries, suggesting M. hyperodae parthenogenesis involves meiosis and that the potential for sexual reproduction may have been retained. Upregulation of genes involved in endoreduplication provides a potential mechanism for the restoration of diploidy in eggs after meiosis.

Effect of Larval Diets on House Fly (Diptera: Muscidae) Production and Host Suitability for a Pupal House Fly Parasitoid (Spalangia gemina) (Hymenoptera: Pteromalidae)

A larval medium for house fly (Musca domestica L.) suitable for house fly parasitoid production was developed from locally available and inexpensive food ingredients. Biological parameters and life table parameters were estimated for house flies treated with five different diets. It was found that percentage survival of 1st-3rd instar larvae of house fly, percentage of pupation, percentage of hatching adults, growth index, and life table parameters (net reproductive rate [R0] the cohort generation time [T], intrinsic rate of increase [r], and finite rate of increase [λ]) were significantly highest on larvae produced on diet 5 (composed of rice barn and chicken feed). Diet composition and performance of house fly larvae in larval medium are discussed, based on the results. The results of Spalangia gemina Bouček (Hymenoptera, Pteromalidae) parasitization on house fly pupae revealed that the number of total parasitized pupa and number of parasitoids hatched from pupa reared with diet 5 (rice bran and chicken feed) were greater than with diet 1 (rice bran, powdered milk, dry yeast, fish meal, soybean meal, and chicken feed). Additionally, S. gemina offspring from diet 5 treatment had a higher proportion of females. Our results indicated that pupal size resulting from larval diet was an important factor for parasitization.

A venom protein, Kazal-type serine protease inhibitor, of ectoparasitoid Pachycrepoideus vindemiae inhibits the hemolymph melanization of host Drosophila melanogaster

Parasitic wasps inject various virulence factors into the host insects while laying eggs, among which the venom proteins, one of the key players in host insect/parasitoid relationships, act in host cellular and humoral immune regulation to ensure successful development of wasp progeny. Although the investigations into actions of venom proteins are relatively ample in larval parasitoids, their regulatory mechanisms have not been thoroughly understood in pupal parasitoids. Here, we identified a venom protein, Kazal-type serine protease inhibitor, in the pupal ectoparasitoid Pachycrepoideus vindemiae (PvKazal). Sequence analysis revealed that PvKazal is packed by a signal peptide and a highly conserved "Kazal" domain. Quantitative polymerase chain reaction analysis recorded a higher transcript level of PvKazal in the venom apparatus relative to that in the carcass, and the PvKazal messenger RNA level appeared to reach a peak on day 5 posteclosion. Recombinant PvKazal strongly inhibited the hemolymph melanization of host Drosophila melanogaster. Additionally, the heterologous expression of PvKazal in transgenic Drosophila reduced the crystal cell numbers and blocked the melanization of host pupal hemolymph. Our present work underlying the roles of PvKazal undoubtedly increases the understanding of venom-mediated host-parasitoid crosstalk.

Identification and Comparative Analysis of Venom Proteins in a Pupal Ectoparasitoid, Pachycrepoideus vindemmiae

Parasitoid wasps inject venom containing complex bioactive compounds to regulate the immune response and development of host arthropods and sometime paralyze host arthropods. Although extensive studies have been conducted on the identification of venom proteins in larval parasitoids, relatively few studies have examined the pupal parasitoids. In our current study, a combination of transcriptomic and proteomic methods was used to identify 64 putative venom proteins from Pachycrepoideus vindemmiae, an ectoparasitoid of Drosophila. Expression analysis revealed that 20 tested venom proteins have 419-fold higher mean expression in the venom apparatus than in other wasp tissues, indicating their specialization to venom. Comparisons of venom proteins from P. vindemmiae and other five species spanning three parasitoid families detected a core set of "ancient" orthologs in Pteromalidae. Thirty-five venom proteins of P. vindemmiae were assigned to the orthologous groups by reciprocal best matches with venoms of other pteromalids, while the remaining 29 were not. Of the 35 categories, twenty-seven have orthologous relationships with Nasonia vitripennis venom proteins and 25 with venoms of Pteromalus puparum. More distant relationships detected that five and two venom proteins of P. vindemmiae are orthologous with venoms of two Figitidae parasitoids and a Braconidae representative, respectively. Moreover, twenty-two venoms unique to P. vindemmiae were also detected, indicating considerable interspecific variation of venom proteins in parasitoids. Phylogenetic reconstruction based on a set of single-copy genes clustered P. vindemmiae with P. puparum, N. vitripennis, and other members of the family Pteromalidae. These findings provide strong evidence that P. vindemmiae venom proteins are well positioned for future functional and evolutionary studies.

Immune Suppressive Extracellular Vesicle Proteins of Leptopilina heterotoma Are Encoded in the Wasp Genome

Leptopilina heterotoma are obligate parasitoid wasps that develop in the body of their Drosophila hosts. During oviposition, female wasps introduce venom into the larval hosts' body cavity. The venom contains discrete, 300 nm-wide, mixed-strategy extracellular vesicles (MSEVs), until recently referred to as virus-like particles. While the crucial immune suppressive functions of L. heterotoma MSEVs have remained undisputed, their biotic nature and origin still remain controversial. In recent proteomics analyses of L. heterotoma MSEVs, we identified 161 proteins in three classes: conserved eukaryotic proteins, infection and immunity related proteins, and proteins without clear annotation. Here we report 246 additional proteins from the L. heterotoma MSEV proteome. An enrichment analysis of the entire proteome supports vesicular nature of these structures. Sequences for more than 90% of these proteins are present in the whole-body transcriptome. Sequencing and de novo assembly of the 460 Mb-sized L. heterotoma genome revealed 90% of MSEV proteins have coding regions within the genomic scaffolds. Altogether, these results explain the stable association of MSEVs with their wasps, and like other wasp structures, their vertical inheritance. While our results do not rule out a viral origin of MSEVs, they suggest that a similar strategy for co-opting cellular machinery for immune suppression may be shared by other wasps to gain advantage over their hosts. These results are relevant to our understanding of the evolution of figitid and related wasp species.

Functional Characterization of a Venom Protein Calreticulin in the Ectoparasitoid Pachycrepoideus vindemiae

Venom proteins act in the immunological interactions between parasitoids and their host insects. The effect of venom proteins on host immunity is not fully understood in pupal parasitoids. We identified the functions of a venom protein, calreticulin (PvCRT), in the pupal ectoparasitoid Pachycrepoideus vindemiae. Here, we report that PvCRT features a signal peptide and two conserved "calreticulin" domains. Multiple sequence alignments show that PvCRT shares 83.54% amino acid identity with CRT from both Pteromalus puparum and Nasonia vitripennis, which infers a close relationship among these three species. Using qPCR analysis, we found a lower expression level of PvCRT (0.27-fold) in the venom apparatus compared to the corresponding carcass. Immunohistochemical localization revealed that PvCRT was ubiquitously expressed in venom gland. The expression of the PvCRT gene in Drosophila transgenic lines via the UAS/Gal4 binary expression system reduced the self-encapsulation phenotype of tu(1)Sz1 mutants. Additionally, studies on humoral immunity indicate that PvCRT does not affect the antimicrobial immune responses of the host. This work on an ectoparasitoid will increase our understanding of venom-mediated host-parasitoid interactions.

The Pupal Ectoparasitoid Pachycrepoideus vindemmiae Regulates Cellular and Humoral Immunity of Host Drosophila melanogaster

The immunological interaction between Drosophila melanogaster and its larval parasitoids has been thoroughly investigated, however, little is known about the interaction between the host and its pupal parasitoids. Pachycrepoideus vindemmiae, a pupal ectoparasitoid of D. melanogaster, injects venom into its host while laying eggs on the puparium, which regulates host immunity and interrupts host development. To resist the invasion of parasitic wasps, various immune defense strategies have been developed in their hosts as a consequence of co-evolution. In this study, we mainly focused on the host immunomodulation by P. vindemmiae and thoroughly investigated cellular and humoral immune response, including cell adherence, cell viability, hemolymph melanization and the Toll, Imd, and JAK/STAT immune pathways. Our results indicated that venom had a significant inhibitory effect on lamellocyte adherence and induced plasmatocyte cell death. Venom injection and in vitro incubation strongly inhibited hemolymph melanization. More in-depth investigation revealed that the Toll and Imd immune pathways were immediately activated upon parasitization, followed by the JAK/STAT pathway, which was activated within the first 24 h post-parasitism. These regulatory effects were further validated by qPCR. Our present study manifested that P. vindemmiae regulated the cellular and humoral immune system of host D. melanogaster in many aspects. These findings lay the groundwork for studying the immunological interaction between D. melanogaster and its pupal parasitoid.

Evaluating the evolution and function of the dynamic Venom Y protein in ectoparasitoid wasps

Venom of the parasitoid wasp Nasonia vitripennis changes the metabolism and gene expression in its fly host Sarcophaga bullata to induce developmental arrest, suppression of the immune response and various other venom effects. Yet, the venom of ectoparasitoid wasps has not been fully characterized. A major component of N. vitripennis venom is an uncharacterized, high-expressing protein referred to as Venom Y. Here we describe the evolutionary history and possible functions of this venom protein. We found that Venom Y is a relatively young gene that has duplicated to form two distinct paralogue groups. A copy of Venom Y has been recruited as a venom protein in at least five wasp species. Functional analysis found that Venom Y affects detoxification and immunity genes in envenomated fly hosts. Many of these genes are fat-body specific, suggesting that Venom Y may have a targeted effect on fat body tissue. We also show that Venom Y may mitigate negative effects of other venom proteins. Finally, protein sequencing indicates that Venom Y is post-translationally modified. This study contributes to elucidating parasitoid venom by using RNA interference knockdown to investigate venom protein function in the context of the whole venom cocktail.

Genome and Ontogenetic-Based Transcriptomic Analyses of the Flesh Fly, Sarcophaga bullata

The flesh fly, Sarcophaga bullata, is a widely-used model for examining the physiology of insect diapause, development, stress tolerance, neurobiology, and host-parasitoid interactions. Flies in this taxon are implicated in myiasis (larval infection of vertebrates) and feed on carrion, aspects that are important in forensic studies. Here we present the genome of S. bullata, along with developmental- and reproduction-based RNA-Seq analyses. We predict 15,768 protein coding genes, identify orthology in relation to closely related flies, and establish sex and developmental-specific gene sets based on our RNA-Seq analyses. Genomic sequences, predicted genes, and sequencing data sets have been deposited at the National Center for Biotechnology Information. Our results provide groundwork for genomic studies that will expand the flesh fly’s utility as a model system.

Combination of label-free quantitative proteomics and transcriptomics reveals intraspecific venom variation between the two strains of Tetrastichus brontispae, a parasitoid of two invasive beetles

The venom apparatus is a conserved organ in parasitoids that shows adaptations correlated with life-style diversification. Combining transcriptomics and label-free quantitative proteomics, here we explored the venom apparatus components of the endoparasitoid Tetrastichus brontispae (Eulophidae), and provide a comparison of the venom apparatus proteomes between its two closely related strains, T. brontispae-Octodonta nipae (Tb-On) and T. brontispae-Brontispa longissima (Tb-Bl). Tb-Bl targets the B. longissima pupa as its habitual host. However, Tb-On is an experimental derivative of Tb-Bl, which has been exposed to the O. nipae pupa as host consecutively for over 40 generation. Results showed that approximately 1505 venom proteins were identified in the T. brontispae venom apparatus. The extracts contained novel venom proteins, such as 4-coumarate-CoA ligase 4. A comparative venom proteome analysis revealed that significant quantitative and qualitative differences in venom composition exist between the two strains; although the most abundant venom proteins were shared between them. The differentially produced proteins were mainly enriched in fatty acid biosynthesis and melanotic encapsulation response. Six of these enriched proteins presented increased levels in Tb-On, and this result was validated by parallel reaction monitoring (PRM) analysis. Overall, our data reveal that venom composition can evolve quickly and respond to host selection.

Venom is beneficial but not essential for development and survival of Nasonia

Parasitoid wasps sting and inject venom in arthropod hosts, which alters host metabolism and development while keeping the host alive for several days, presumably to induce benefits for the parasitoid young.Here we investigate the consequences of host envenomation on development and fitness of wasp larvae in the ectoparasitoid Nasonia vitripennis, by comparing wasps reared on live unstung, previously stung, and cold-killed hosts. Developmental arrest and suppression of host response to larvae are major venom effects that occur in both stung and cold-killed hosts, but not unstung hosts; while cold-killed hosts lack venom effects that require a living host. Thus, cold-killed hosts mimic some of the effects of venom, but not others.Eggs placed on live unstung hosts have significantly higher mortality during development, however successfully developing wasps from these hosts have similar lifetime fecundity to wasps from cold-killed or stung hosts. Therefore, although venom is beneficial, it is not required for wasp survival.While wasps developing on cold-killed versus stung hosts have similar fitness, multiple generations of rearing on cold-killed hosts results in significant fitness reductions of wasps.We conclude that the largest benefits of venom are induction of host developmental arrest and suppression of host response to larva (e.g. immune responses), although more subtle benefits may accrue across generations, or under stressful conditions.

Functional characterization of the transcriptional regulatory elements of three highly expressed constitutive genes in the jewel wasp, Nasonia vitripennis

The jewel wasp, Nasonia vitripennis Ashmead (Hymenoptera: Pteromalidae), is an easily reared parasitoid that is providing an ever increasingly malleable model for examining the biology and genetics of Hymenoptera. Utilizing genomic and transcriptome resources, 5' upstream transcriptional regulatory sequences (TREs) from three highly expressed genes were identified and cloned. Criteria for TRE selection included the presence of an adjacent gene 5' of the translation initiation site. One gene was methylated whereas the other two were nonmethylated. Each TRE, heat-shock protein 70 (hsp70), activator of 90 kDa hsp ATPase protein 1 (hsp90A), and lipid storage droplet surface-binding protein 1 (lsdp) was linked with enhanced green fluorescent protein (EGFP) coding sequence and cloned into both pDP9e somatic and piggyBac germline transformation vectors. EGFP expression patterns under control of each TRE were compared with patterns of DsRed fluorescence produced from the transformation vector cassette. Functional activity of each TRE was observed in cultured Spodoptera frugiperda 9 (Sf9) cells and Drosophila melanogaster as well as in N. vitripennis embryos demonstrating that all three sequences had functional transcriptional regulatory activity in three different insect orders. Identification and functional characterization of these three TREs will provide critical and necessary resources for further genetic analyses of N. vitripennis, Hymenoptera and other insects.

The Evolution of Venom by Co-option of Single-Copy Genes

The classic model for the evolution of novel gene function is through gene duplication followed by evolution of a new function by one of the copies (neofunctionalization) [1, 2]. However, other modes have also been found, such as novel genes arising from non-coding DNA, chimeric fusions, and lateral gene transfers from other organisms [3-7]. Here we use the rapid turnover of venom genes in parasitoid wasps to study how new gene functions evolve. In contrast to the classic gene duplication model, we find that a common mode of acquisition of new venom genes in parasitoid wasps is co-option of single-copy genes from non-venom progenitors. Transcriptome and proteome sequencing reveal that recruitment and loss of venom genes occur primarily by rapid cis-regulatory expression evolution in the venom gland. Loss of venom genes is primarily due to downregulation of expression in the gland rather than gene death through coding sequence degradation. While the majority of venom genes have specialized expression in the venom gland, recent losses of venom function occur primarily among genes that show broader expression in development, suggesting that they can more readily switch functional roles. We propose that co-option of single-copy genes may be a common but relatively understudied mechanism of evolution for new gene functions, particularly under conditions of rapid evolutionary change.

Protein Discovery: Combined Transcriptomic and Proteomic Analyses of Venom from the Endoparasitoid Cotesia chilonis (Hymenoptera: Braconidae)

Many species of endoparasitoid wasps provide biological control services in agroecosystems. Although there is a great deal of information on the ecology and physiology of host/parasitoid interactions, relatively little is known about the protein composition of venom and how specific venom proteins influence physiological systems within host insects. This is a crucial gap in our knowledge because venom proteins act in modulating host physiology in ways that favor parasitoid development. Here, we identified 37 possible venom proteins from the polydnavirus-carrying endoparasitoid Cotesia chilonis by combining transcriptomic and proteomic analyses. The most abundant proteins were hydrolases, such as proteases, peptidases, esterases, glycosyl hydrolase, and endonucleases. Some components are classical parasitoid venom proteins with known functions, including extracellular superoxide dismutase 3, serine protease inhibitor and calreticulin. The venom contains novel proteins, not recorded from any other parasitoid species, including tolloid-like proteins, chitooligosaccharidolytic β-N-acetylglucosaminidase, FK506-binding protein 14, corticotropin-releasing factor-binding protein and vascular endothelial growth factor receptor 2. These new data generate hypotheses and provide a platform for functional analysis of venom components.

Laterally Transferred Gene Recruited as a Venom in Parasitoid Wasps

Parasitoid wasps use venom to manipulate the immunity and metabolism of their host insects in a variety of ways to provide resources for their offspring. Yet, how genes are recruited and evolve to perform venom functions remain open questions. A recently recognized source of eukaryotic genome innovation is lateral gene transfer (LGT). Glycoside hydrolase family 19 (GH19) chitinases are widespread in bacteria, microsporidia, and plants where they are used in nutrient acquisition or defense, but have previously not been known in metazoans. In this study, a GH19 chitinase LGT is described from the unicellular microsporidia/Rozella clade into parasitoid wasps of the superfamily Chalcidoidea, where it has become recruited as a venom protein. The GH19 chitinase is present in 15 species of chalcidoid wasps representing four families, and phylogenetic analysis indicates that it was laterally transferred near or before the origin of Chalcidoidea (∼95 Ma). The GH19 chitinase gene is highly expressed in the venom gland of at least seven species, indicating a role in the complex host manipulations performed by parasitoid wasp venom. RNAi knockdown in the model parasitoid Nasonia vitripennis reveals that-following envenomation-the GH19 chitinase induces fly hosts to upregulate genes involved in an immune response to fungi. A second, independent LGT of GH19 chitinase from microsporidia into mosquitoes was also found, also supported by phylogenetic reconstructions. Besides these two LGT events, GH19 chitinase is not found in any other sequenced animal genome, or in any fungi outside the microsporidia/Rozella clade.