Genome of the pincer wasp Gonatopus flavifemur reveals unique venom evolution and a dual adaptation to parasitism and predation

An expanded odorant-binding protein mediates host cue detection in the parasitic wasp Baryscapus dioryctriae basis of the chromosome-level genome assembly analysis

BACKGROUND

Baryscapus dioryctriae (Chalcidodea: Eulophidae) is a parasitic wasp that parasitizes the pupae of many Pyralidae members and has been used as a biological control agent against Dioryctria pests of pinecones.

RESULTS

This B. dioryctriae assembly has a genome size of 485.5 Mb with a contig N50 of 2.17 Mb, and scaffolds were assembled onto six chromosomes using Hi-C analysis, significantly increasing the scaffold N50 to 91.17 Mb, with more than 96.13% of the assembled bases located on chromosomes, and an analysis revealed that 94.73% of the BUSCO gene set. A total of 54.82% (279.27 Mb) of the assembly was composed of repetitive sequences and 24,778 protein-coding genes were identified. Comparative genomic analysis demonstrated that the chemosensory perception, genetic material synthesis, and immune response pathways were primarily enriched in the expanded genes. Moreover, the functional characteristics of an odorant-binding protein (BdioOBP45) with ovipositor-biased expression identified from the expanded olfactory gene families were investigated by the fluorescence competitive binding and RNAi assays, revealing that BdioOBP45 primarily binds to the D. abietella-induced volatile compounds, suggesting that this expanded OBP is likely involved in locating female wasp hosts and highlighting a direction for future research.

CONCLUSIONS

Taken together, this work not only provides new genomic sequences for the Hymenoptera systematics, but also the high-quality chromosome-level genome of B. dioryctriae offers a valuable foundation for studying the molecular, evolutionary, and parasitic processes of parasitic wasps.

The Chromosome-level Genome Provides Insights into the Evolution and Adaptation of Extreme Aggression

Extremely aggressive behavior, as the special pattern, is rare in most species and characteristic as contestants severely injured or killed ending the combat. Current studies of extreme aggression are mainly from the perspectives of behavioral ecology and evolution, while lacked the aspects of molecular evolutionary biology. Here, a high-quality chromosome-level genome of the parasitoid Anastatus disparis was provided, in which the males exhibit extreme mate-competition aggression. The integrated multiomics analysis highlighted that neurotransmitter dopamine overexpression, energy metabolism (especially from lipid), and antibacterial activity are likely major aspects of evolutionary formation and adaptation for extreme aggression in A. disparis. Conclusively, our study provided new perspectives for molecular evolutionary studies of extreme aggression as well as a valuable genomic resource in Hymenoptera.

The state of parasitoid wasp genomics

Parasitoid wasps represent a group of parasitic insects with high species diversity that have played a pivotal role in biological control and evolutionary studies. Over the past 20 years, developments in genomics have greatly enhanced our understanding of the biology of these species. Technological leaps in sequencing have facilitated the improvement of genome quality and quantity, leading to the availability of hundreds of parasitoid wasp genomes. Here, we summarize recent progress in parasitoid wasp genomics, focusing on the evolution of genome size (GS) and the genomic basis of several key traits. We also discuss the contributions of genomics in studying venom evolution and endogenization of viruses. Finally, we advocate for increased sequencing and functional research to better understand parasitoid biology and enhance biological control.

A taxon-rich and genome-scale phylogeny of Opisthokonta

Ancient divergences within Opisthokonta-a major lineage that includes organisms in the kingdoms Animalia, Fungi, and their unicellular relatives-remain contentious. To assess progress toward a genome-scale Opisthokonta phylogeny, we conducted the most taxon rich phylogenomic analysis using sets of genes inferred with different orthology inference methods and established the geological timeline of Opisthokonta diversification. We also conducted sensitivity analysis by subsampling genes or taxa from the full data matrix based on filtering criteria previously shown to improve phylogenomic inference. We found that approximately 85% of internal branches were congruent across data matrices and the approaches used. Notably, the use of different orthology inference methods was a substantial contributor to the observed incongruence: analyses using the same set of orthologs showed high congruence of 97% to 98%, whereas different sets of orthologs resulted in somewhat lower congruence (87% to 91%). Examination of unicellular Holozoa relationships suggests that the instability observed across varying gene sets may stem from weak phylogenetic signals. Our results provide a comprehensive Opisthokonta phylogenomic framework that will be useful for illuminating ancient evolutionary episodes concerning the origin and diversification of the 2 major eukaryotic kingdoms and emphasize the importance of investigating effects of orthology inference on phylogenetic analyses to resolve ancient divergences.

Host and venom evolution in parasitoid wasps: does independently adapting to the same host shape the evolution of the venom gland transcriptome?

BACKGROUND

Venoms have repeatedly evolved over 100 occasions throughout the animal tree of life, making them excellent systems for exploring convergent evolutionary novelty. Growing evidence supports that venom evolution is predominantly driven by prey or host-related selection pressures, and the expression patterns of venom glands reflect adaptive evolution. However, it remains elusive whether the evolution of expression patterns in venom glands is likewise a convergent evolution driven by their prey/host species.

RESULTS

We utilized parasitoid wasps that had independently adapted to Drosophila hosts as models to investigate the convergent evolution of venom gland transcriptomes in 19 hymenopteran species spanning ~ 200 million years of evolution. Comparative transcriptome analysis reveals that the global expression patterns among the venom glands of Drosophila parasitoid wasps do not achieve higher similarity compared to non-Drosophila parasitoid wasps. Further evolutionary analyses of expression patterns at the single gene, orthogroup, and Gene Ontology (GO) term levels indicate that some orthogroups/GO terms show correlation with the Drosophila parasitoid wasps. However, these groups rarely include genes highly expressed in venom glands or putative venom genes in the Drosophila parasitoid wasps.

CONCLUSIONS

Our study suggests that convergent evolution may not play a predominant force shaping gene expression levels in the venom gland of the Drosophila parasitoid wasps, offering novel insights into the co-evolution between venom and prey/host.

Identification and Functional Analysis of Adipokinetic Hormone Receptor in Ostrinia furnacalis Guenée Larvae Parasitized by Macrocentrus cingulum

As a typical G protein-coupled receptor, the adipokinetic hormone receptor (AKHR) has seven transmembrane domains (TMDs), and its structure and function are similar to the gonadotropin-releasing hormone receptor (GnRHR) in vertebrates. However, there is a dearth of information on other components of the AKHR signaling pathway and how it functions in the interaction between insect hosts and parasitoids. In this study, we cloned and analyzed the multifunctional Ostrinia furnacalis AKHR (OfAKHR) cDNA (GenBank accession number MF797868). OfAKHR has a 2206 bp full-length cDNA, which includes an open reading frame containing 1194 bp. OfAKHR contains the typical seven TMDs, and a "DRY" motif. OfAKHR has the highest relative expression in the fat body and the fifth instar larvae. The results revealed that ApoLpⅢ, PPO2, GS, TPS, Cecropin, and Moricin decreased the transcription levels from 48 to 72 h after the knockdown of OfAKHR expression by dsOfAKHR injection in the fourth instar O. furnacalis larvae. The parasitization of Macrocentrus cingulum selectively upregulated the expression levels of nutrition metabolism and immune-related genes in parasitized O. furnacalis larvae, stimulated lysozyme activity, and obviously raised the concentrations of triglyceride and trehalose in the hemolymph of O. furnacalis larvae. However, they inhibited the activities of PO and trehalase. This study is conducive to a deeper cognition of the roles of OfAKHR in nutrition and immune homeostasis, coevolution, and coexistence between parasitic wasps and hosts. It also sheds light on the potential as the target of pest control reagents.

Comparative transcriptome and proteome reveal the unique genes and proteins of female parasitic wasps, Lysiphlebia japonica Ashmead

BACKGROUND

Lysiphlebia japonica Ashmead (Hymenoptera, Braconidae) is an endophagous parasitoid wasp and its host, Aphis gossypii Glover (Hemiptera, Aphididae) is a major cotton pest. L. japonica affects the growth and fatty acid metabolism of cotton aphids after parasitization and has been widely used as a biocontrol agent. However, there are currently few reports about the molecular characteristics of L. japonica, especially the differences between male and female.

RESULTS

In this study, using transcriptome and proteome analysis of the abdomen of female and male parasitic wasps, respectively, we obtained a total of 27,169 DEGs and 1,194 DEPs, then a total of 909 positively correlated high-expression proteins and genes were obtained by combined omics analysis. Subsequently, 20 differentially expressed abdomen specific proteins were selected for validation by RT-qPCR and Multiple Reaction Monitoring (MRM) protein verification. The result of RT-qPCR demonstrated that all 20 genes were highly expressed in the abdomen of females, and five target proteins with unique peptide fragments and identification profiles were identified by MRM, which were venom protease, tropomyosin, lipase member I, venom serine carboxypeptidase and calreticulin, respectively.

CONCLUSION

Overall, these results provided molecular resources for the differences between males and females in L. japonica and the screened 20 abdomen specific proteins were verified to demonstrate the validity of the data, which offered important molecular data resources for further studies on the related functional genes of parasitic wasps and the mechanism of parasitoids regulating the growth of aphids. © 2023 Society of Chemical Industry.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

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.

Comprehensive isoform-level analysis reveals the contribution of alternative isoforms to venom evolution and repertoire diversity

Animal venom systems have emerged as valuable models for investigating how novel polygenic phenotypes may arise from gene evolution by varying molecular mechanisms. However, a significant portion of venom genes produce alternative mRNA isoforms that have not been extensively characterized, hindering a comprehensive understanding of venom biology. In this study, we present a full-length isoform-level profiling workflow integrating multiple RNA sequencing technologies, allowing us to reconstruct a high-resolution transcriptome landscape of venom genes in the parasitoid wasp Pteromalus puparum Our findings demonstrate that more than half of the venom genes generate multiple isoforms within the venom gland. Through mass spectrometry analysis, we confirm that alternative splicing contributes to the diversity of venom proteins, acting as a mechanism for expanding the venom repertoire. Notably, we identified seven venom genes that exhibit distinct isoform usages between the venom gland and other tissues. Furthermore, evolutionary analyses of venom serpin3 and orcokinin further reveal that the co-option of an ancient isoform and a newly evolved isoform, respectively, contributes to venom recruitment, providing valuable insights into the genetic mechanisms driving venom evolution in parasitoid wasps. Overall, our study presents a comprehensive investigation of venom genes at the isoform level, significantly advancing our understanding of alternative isoforms in venom diversity and evolution and setting the stage for further in-depth research on venoms.

Molecular Characterization and Functional Analysis of the Dipeptidyl Peptidase IV from Venom of the Ectoparasitoid Scleroderma guani

Dipeptidyl peptidase IV (DPPIV) is a proline-specific serine peptidase that remains poorly investigated in terms of venom composition. Here, we describe the molecular characteristics and possible functions of DPPIV as a major venom component of the ant-like bethylid ectoparasitoid, Scleroderma guani, named SgVnDPPIV. The SgVnDPPIV gene was cloned, which encodes a protein with the conserved catalytic triads and substrate binding sites of mammalian DPPIV. This venom gene is highly expressed in the venom apparatus. Recombinant SgVnDPPIV, produced in Sf9 cells using the baculovirus expression system, has high enzymatic activity, which can be efficiently inhibited by vildagliptin and sitagliptin. Functional analysis revealed that SgVnDPPIV affects genes related to detoxification, lipid synthesis and metabolism, response to stimuli, and ion exchange in pupae of Tenebrio molitor, an envenomated host of S. guani. The present work contributes towards understanding the role of venom DPPIV involved in the interaction between parasitoid wasp and its host.

Functional and Proteomic Insights into Aculeata Venoms

Aculeate hymenopterans use their venom for a variety of different purposes. The venom of solitary aculeates paralyze and preserve prey without killing it, whereas social aculeates utilize their venom in defence of their colony. These distinct applications of venom suggest that its components and their functions are also likely to differ. This study investigates a range of solitary and social species across Aculeata. We combined electrophoretic, mass spectrometric, and transcriptomic techniques to characterize the compositions of venoms from an incredibly diverse taxon. In addition, in vitro assays shed light on their biological activities. Although there were many common components identified in the venoms of species with different social behavior, there were also significant variations in the presence and activity of enzymes such as phospholipase A2s and serine proteases and the cytotoxicity of the venoms. Social aculeate venom showed higher presence of peptides that cause damage and pain in victims. The venom-gland transcriptome from the European honeybee (Apis mellifera) contained highly conserved toxins which match those identified by previous investigations. In contrast, venoms from less-studied taxa returned limited results from our proteomic databases, suggesting that they contain unique toxins.

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.

Genomic analyses provide insights into the genome evolution and environmental adaptation of the tobacco moth Ephestia elutella

Ephestia elutella is a major pest responsible for significant damage to stored tobacco over many years. Here, we conduct a comparative genomic analysis on this pest, aiming to explore the genetic bases of environmental adaptation of this species. We find gene families associated with nutrient metabolism, detoxification, antioxidant defense and gustatory receptors are expanded in the E. elutella genome. Detailed phylogenetic analysis of P450 genes further reveals obvious duplications in the CYP3 clan in E. elutella compared to the closely related species, the Indianmeal moth Plodia interpunctella. We also identify 229 rapidly evolving genes and 207 positively selected genes in E. elutella, respectively, and highlight two positively selected heat shock protein 40 (Hsp40) genes. In addition, we find a number of species-specific genes related to diverse biological processes, such as mitochondria biology and development. These findings advance our understanding of the mechanisms underlying processes of environmental adaptation on E. elutella and will enable the development of novel pest management strategies.

Genomic signatures associated with maintenance of genome stability and venom turnover in two parasitoid wasps

Parasitoid wasps are rapidly developing as a model for evolutionary biology. Here we present chromosomal genomes of two Anastatus wasps, A. japonicus and A. fulloi, and leverage these genomes to study two fundamental questions-genome size evolution and venom evolution. Anastatus shows a much larger genome than is known among other wasps, with unexpectedly recent bursts of LTR retrotransposons. Importantly, several genomic innovations, including Piwi gene family expansion, ubiquitous Piwi expression profiles, as well as transposable element-piRNA coevolution, have likely emerged for transposable element silencing to maintain genomic stability. Additionally, we show that the co-option evolution arose by expression shifts in the venom gland plays a dominant role in venom turnover. We also highlight the potential importance of non-venom genes that are coexpressed with venom genes during venom evolution. Our findings greatly advance the current understanding of genome size evolution and venom evolution, and these genomic resources will facilitate comparative genomics studies of insects in the future.

Transposons and non-coding regions drive the intrafamily differences of genome size in insects

Genome size (GS) can vary considerably between phylogenetically close species, but the landscape of GS changes in insects remain largely unclear. To better understand the specific evolutionary factors that determine GS in insects, we examined flow cytometry-based published GS data from 1,326 insect species, spanning 700 genera, 155 families, and 21 orders. Model fitting showed that GS generally followed an Ornstein-Uhlenbeck adaptive evolutionary model in Insecta overall. Ancestral reconstruction indicated a likely GS of 1,069 Mb, suggesting that most insect clades appeared to undergo massive genome expansions or contractions. Quantification of genomic components in 56 species from nine families in four insect orders revealed that the proliferation of transposable elements contributed to high variation in GS between close species, such as within Coleoptera. This study sheds lights on the pattern of GS variation in insects and provides a better understanding of insect GS evolution.

Genome of the parasitoid wasp Cotesia chilonis sheds light on amino acid resource exploitation

Background

A fundamental feature of parasitism is the nutritional exploitation of host organisms by their parasites. Parasitoid wasps lay eggs on arthropod hosts, exploiting them for nutrition to support larval development by using diverse effectors aimed at regulating host metabolism. However, the genetic components and molecular mechanisms at the basis of such exploitation, especially the utilization of host amino acid resources, remain largely unknown. To address this question, here, we present a chromosome-level genome assembly of the parasitoid wasp Cotesia chilonis and reconstruct its amino acid biosynthetic pathway.

Results

Analyses of the amino acid synthetic pathway indicate that C. chilonis lost the ability to synthesize ten amino acids, which was confirmed by feeding experiments with amino acid-depleted media. Of the ten pathways, nine are known to have been lost in the common ancestor of animals. We find that the ability to synthesize arginine was also lost in C. chilonis because of the absence of two key genes in the arginine synthesis pathway. Further analyses of the genomes of 72 arthropods species show that the loss of arginine synthesis is common in arthropods. Metabolomic analyses by UPLC-MS/MS reveal that the temporal concentrations of arginine, serine, tyrosine, and alanine are significantly higher in host (Chilo suppressalis) hemolymph at 3 days after parasitism, whereas the temporal levels of 5-hydroxylysine, glutamic acid, methionine, and lysine are significantly lower. We sequence the transcriptomes of a parasitized host and non-parasitized control. Differential gene expression analyses using these transcriptomes indicate that parasitoid wasps inhibit amino acid utilization and activate protein degradation in the host, likely resulting in the increase of amino acid content in host hemolymph.

Conclusions

We sequenced the genome of a parasitoid wasp, C. chilonis, and revealed the features of trait loss in amino acid biosynthesis. Our work provides new insights into amino acid exploitation by parasitoid wasps, and this knowledge can specifically be used to design parasitoid artificial diets that potentially benefit mass rearing of parasitoids for pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01313-3.

Comparative Genomics Sheds Light on the Convergent Evolution of Miniaturized Wasps

Miniaturization has occurred in many animal lineages, including insects and vertebrates, as a widespread trend during animal evolution. Among Hymenoptera, miniaturization has taken place in some parasitoid wasp lineages independently, and may have contributed to the diversity of species. However, the genomic basis of miniaturization is little understood. Diverged approximately 200 Ma, Telenomus wasps (Platygastroidea) and Trichogramma wasps (Chalcidoidea) have both evolved to a highly reduced body size independently, representing a paradigmatic example of convergent evolution. Here, we report a high-quality chromosomal genome of Telenomus remus, a promising candidate for controlling Spodoptera frugiperda, a notorious pest that has recently caused severe crop damage. The T. remus genome (129 Mb) is characterized by a low density of repetitive sequence and a reduction of intron length, resulting in the shrinkage of genome size. We show that hundreds of genes evolved faster in two miniaturized parasitoids Trichogramma pretiosum and T. remus. Among them, 38 genes exhibit extremely accelerated evolutionary rates in these miniaturized wasps, possessing diverse functions in eye and wing development as well as cell size control. These genes also highlight potential roles in body size regulation. In sum, our analyses uncover a set of genes with accelerated evolutionary rates in Tri. pretiosum and T. remus, which might be responsible for their convergent adaptations to miniaturization, and thus expand our understanding on the evolutionary basis of miniaturization. Additionally, the genome of T. remus represents the first genome resource of superfamily Platygastroidea, and will facilitate future studies of Hymenoptera evolution and pest control.