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

The gene regulatory mechanisms shaping the heterogeneity of venom production in the Cape coral snake

BACKGROUND

Venoms and their associated glands and delivery structures have evolved numerous times among animals. Within these venom systems, the molecular, cellular, and morphological components interact and co-evolve to generate distinct, venom phenotypes that are increasingly recognized as models for studying adaptive evolution. However, toxins are often unevenly distributed across venom-producing tissues in patterns that are not necessarily adaptive but instead likely result from constraints associated with protein secretion.

RESULTS

We generate a high-quality draft genome of the Cape coral snake (Aspidelaps lubricus) and combine analyses of venom gland single-cell RNA-seq data with spatial venom gland in situ toxin distributions. Our results reveal that while different toxin families are produced by distinct populations of cells, toxin expression is fine-tuned by regulatory modules that result in further specialization of toxin production within each cell population. We also find that the evolution of regulatory elements closely mirrors the evolution of their associated toxin genes, resulting in spatial association of closely related and functionally similar toxins in the venom gland. While this compartmentalization is non-adaptive, the modularity of the underlying regulatory network likely facilitated the repeated evolution of defensive venom in spitting cobras.

CONCLUSIONS

Our results provide new insight into the variability of toxin regulation across snakes, reveal the molecular mechanisms underlying the heterogeneous toxin production in snake venom glands, and provide an example of how constraints can result in non-adaptive character states that appear to be adaptive, which may nevertheless facilitate evolutionary innovation and novelty.

Genomics of insect natural enemies in agroecosystems

Currently, a wealth of genomic data are now accessible for numerous insect natural enemies, serving as valuable resources that deepen our understanding of the genetic basis of biocontrol traits in these organisms. We summarize the current state of genome sequencing and highlight candidate genes related to biocontrol traits that hold promise for genetic improvement. We also review the recent population genomic studies in biological control and the discovery of potential insecticidal genes in parasitoid wasps. Collectively, current genomic works have shown the powerful ability to identify candidate genes responsible for desirable traits or promising effectors. However, further functional study is necessary to gain a mechanistic understanding of these genes, and future efforts are also needed to develop suitable approaches to translate genomic insights into field applications.

Integrative multi-omics analysis reveals the contribution of neoVTX genes to venom diversity of Synanceia verrucosa

BACKGROUND

Animal venom systems are considered as valuable model for investigating the molecular mechanisms underlying phenotypic evolution. Stonefish are the most venomous and dangerous fish because of severe human envenomation and occasionally fatalities, whereas the genomic background of their venom has not been fully explored compared with that in other venomous animals.

RESULTS

In this study, we followed modern venomic pipelines to decode the Synanceia verrucosa venom components. A catalog of 478 toxin genes was annotated based on our assembled chromosome-level genome. Integrative analysis of the high-quality genome, the transcriptome of the venom gland, and the proteome of crude venom revealed mechanisms underlying the venom complexity in S. verrucosa. Six tandem-duplicated neoVTX subunit genes were identified as the major source for the neoVTX protein production. Further isoform sequencing revealed massive alternative splicing events with a total of 411 isoforms demonstrated by the six genes, which further contributed to the venom diversity. We then characterized 12 dominantly expressed toxin genes in the venom gland, and 11 of which were evidenced to produce the venom protein components, with the neoVTX proteins as the most abundant. Other major venom proteins included a presumed CRVP, Kuntiz-type serine protease inhibitor, calglandulin protein, and hyaluronidase. Besides, a few of highly abundant non-toxin proteins were also characterized and they were hypothesized to function in housekeeping or hemostasis maintaining roles in the venom gland. Notably, gastrotropin like non-toxin proteins were the second highest abundant proteins in the venom, which have not been reported in other venomous animals and contribute to the unique venom properties of S. verrucosa.

CONCLUSIONS

The results identified the major venom composition of S. verrucosa, and highlighted the contribution of neoVTX genes to the diversity of venom composition through tandem-duplication and alternative splicing. The diverse neoVTX proteins in the venom as lethal particles are important for understanding the adaptive evolution of S. verrucosa. Further functional studies are encouraged to exploit the venom components of S. verrucosa for pharmaceutical innovation.

Parasitoid speciation and diversification

Parasitoid wasps may well be the most species-rich animal group on Earth, and host-parasitoid interactions may thereby be one of the most common types of species interactions. Understanding the major mechanisms underlying diversification in parasitoids should be a high priority, not the least in order to predict consequences from high extinction rates currently observed. The two major hypotheses explaining host-associated diversification are the escape-and-radiate hypothesis and the oscillation hypothesis, where the former assumes that key innovations are major drivers of radiation bursts, whereas the latter rather assumes that diversification depends on processes acting on the standing genetic variation that influences host use. This paper reviews the recent literature on parasitoid speciation in light of these major hypotheses to identify potential key innovations and host use variability underlying diversification. The paper also calls upon recent theoretical advances from a similar system, plant-butterfly interactions, to provide shortcuts in the development of theories explaining the high diversity of parasitoid wasps.

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.

Drosophila are hosts to the first described parasitoid wasp of adult flies

Parasitoid wasps are exceptionally diverse and use specialized adaptations capable of manipulating the physiology and behaviour of host organisms1. In more than two centuries since the first records of Drosophila-parasitizing wasps, nearly 200 described and provisional parasitoid species of drosophilids have been identified2. These include endoparasitoids and ectoparasitoids, as well as species attacking larval and pupal hosts3. Despite a deep history of research attention and remarkable biodiversity, a wasp species that attacks and develops inside the adult stage of a fly host has not been described previously. Here we report the discovery of a wasp species that infects the adult stage of fruit flies in the genus Drosophila, including one of the most deeply studied model organisms in biology, Drosophila melanogaster. Notably, this wasp can be easily collected from backyard fly baits and has a broad geographic distribution throughout the eastern USA. We document its life history and unique host interactions, including egg-laying into and larval emergence from adult flies, and provide protocols to raise wasps from wild-caught host flies. Our results emphasize the need for ongoing research investment in insect biodiversity and systematics. As parasitoid research continues to uncover unusual biology and supports fundamental mechanistic insights into immunity4, metabolism5, ecology6, evolution7-9 and behaviour10-12, we anticipate that this wasp's association with the laboratory model organism, D. melanogaster, will provide new research opportunities across the life sciences.

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.

Coordinated molecular and ecological adaptations underlie a highly successful parasitoid

The success of an organism depends on the molecular and ecological adaptations that promote its beneficial fitness. Parasitoids are valuable biocontrol agents for successfully managing agricultural pests, and they have evolved diversified strategies to adapt to both the physiological condition of hosts and the competition of other parasitoids. Here, we deconstructed the parasitic strategies in a highly successful parasitoid, Trichopria drosophilae, which parasitizes a broad range of Drosophila hosts, including the globally invasive species D. suzukii. We found that T. drosophilae had developed specialized venom proteins that arrest host development to obtain more nutrients via secreting tissue inhibitors of metalloproteinases (TIMPs), as well as a unique type of cell-teratocytes-that digest host tissues for feeding by releasing trypsin proteins. In addition to the molecular adaptations that optimize nutritional uptake, this pupal parasitoid has evolved ecologically adaptive strategies including the conditional tolerance of intraspecific competition to enhance parasitic success in older hosts and the obligate avoidance of interspecific competition with larval parasitoids. Our study not only demystifies how parasitoids weaponize themselves to colonize formidable hosts but also provided empirical evidence of the intricate coordination between the molecular and ecological adaptations that drive evolutionary success.

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.

The genome of the ant Tetramorium bicarinatum reveals a tandem organization of venom peptides genes allowing the prediction of their regulatory and evolutionary profiles

BACKGROUND

Venoms have evolved independently over a hundred times in the animal kingdom to deter predators and/or subdue prey. Venoms are cocktails of various secreted toxins, whose origin and diversification provide an appealing system for evolutionary researchers. Previous studies of the ant venom of Tetramorium bicarinatum revealed several Myrmicitoxin (MYRTX) peptides that gathered into seven precursor families suggesting different evolutionary origins. Analysis of the T. bicarinatum genome enabling further genomic approaches was necessary to understand the processes underlying the evolution of these myrmicitoxins.

RESULTS

Here, we sequenced the genome of Tetramorium bicarinatum and reported the organisation of 44 venom peptide genes (vpg). Of the eleven chromosomes that make up the genome of T. bicarinatum, four carry the vpg which are organized in tandem repeats. This organisation together with the ML evolutionary analysis of vpg sequences, is consistent with evolution by local duplication of ancestral genes for each precursor family. The structure of the vpg into two or three exons is conserved after duplication events while the promoter regions are the least conserved parts of the vpg even for genes with highly identical sequences. This suggests that enhancer sequences were not involved in duplication events, but were recruited from surrounding regions. Expression level analysis revealed that most vpg are highly expressed in venom glands, although one gene or group of genes is much more highly expressed in each family. Finally, the examination of the genomic data revealed that several genes encoding transcription factors (TFs) are highly expressed in the venom glands. The search for binding sites (BS) of these TFs in the vpg promoters revealed hot spots of GATA sites in several vpg families.

CONCLUSION

In this pioneering investigation on ant venom genes, we provide a high-quality assembly genome and the annotation of venom peptide genes that we think can fosters further genomic research to understand the evolutionary history of ant venom biochemistry.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality

BACKGROUND

Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected.

RESULTS

Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification.

CONCLUSIONS

Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

The genome sequence of a chalcid wasp, Gastracanthus pulcherrimus (Westwood, 1833)

We present a genome assembly from an individual female Gastracanthus pulcherrimus (a chalcid wasp; Arthropoda; Insecta; Hymenoptera; Pteromalidae). The genome sequence is 1,010.0 megabases in span. Most of the assembly is scaffolded into 5 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 24.4 kilobases in length.

Morphological, functional, compositional and transcriptional constraints shape the distinct venom profiles of the assassin bug Sycanus croceovittatus

Predatory bugs employ a salivary venom apparatus to generate complex venoms for capturing and digesting prey. The venom apparatus consists of different glands for the production of distinct venom sets, but the underlying mechanisms behind this process remain poorly understood. Here we present a comprehensive analysis of the morphological, functional, compositional and transcriptional characteristics of venoms derived from posterior main gland (PMG), anterior main gland (AMG), and accessory gland (AG) of the assassin bug Sycanus croceovittatus. Structural observations revealed the intricate constructions of the venom apparatus, enabling the production and storage of three distinct venom sets in anatomically varied glands and allowing them to be modulated in a context-dependent manner upon utilization. There were remarkable differences in the biological activities exhibited by PMG, AMG, and AG venoms. Proteotranscriptomic analysis demonstrated that these venoms displayed compositional heterogeneity at both the quantity and variety levels of proteins. Transcriptional profiles of the identified venom proteins revealed gland-specific or biased expression patterns. These findings indicate that the divergence in venom profiles among different glands arises from morphological, functional, compositional and transcriptional constraints on the venom apparatus, reflecting remarkable morphogenesis and regulatory gene networks responsible for the compartmentalized production of venom proteins in different glands.

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.

Whole Genome Duplication and Gene Evolution in the Hyperdiverse Venomous Gastropods

The diversity of venomous organisms and the toxins they produce have been increasingly investigated, but taxonomic bias remains important. Neogastropods, a group of marine predators representing almost 22% of the known gastropod diversity, evolved a wide range of feeding strategies, including the production of toxins to subdue their preys. However, whether the diversity of these compounds is at the origin of the hyperdiversification of the group and how genome evolution may correlate with both the compounds and species diversities remain understudied. Among the available gastropods genomes, only eight, with uneven quality assemblies, belong to neogastropods. Here, we generated chromosome-level assemblies of two species belonging to the Tonnoidea and Muricoidea superfamilies (Monoplex corrugatus and Stramonita haemastoma). The two obtained high-quality genomes had 3 and 2.2 Gb, respectively, and 92-89% of the total assembly conformed 35 pseudochromosomes in each species. Through the analysis of syntenic blocks, Hox gene cluster duplication, and synonymous substitutions distribution pattern, we inferred the occurrence of a whole genome duplication event in both genomes. As these species are known to release venom, toxins were annotated in both genomes, but few of them were found in homologous chromosomes. A comparison of the expression of ohnolog genes (using transcriptomes from osphradium and salivary glands in S. haemastoma), where both copies were differentially expressed, showed that most of them had similar expression profiles. The high quality of these genomes makes them valuable reference in their respective taxa, facilitating the identification of genome-level processes at the origin of their evolutionary success.

Horizontal gene transfer underlies the painful stings of asp caterpillars (Lepidoptera: Megalopygidae)

Larvae of the genus Megalopyge (Lepidoptera: Zygaenoidea: Megalopygidae), known as asp or puss caterpillars, produce defensive venoms that cause severe pain. Here, we present the anatomy, chemistry, and mode of action of the venom systems of caterpillars of two megalopygid species, the Southern flannel moth Megalopyge opercularis and the black-waved flannel moth Megalopyge crispata. We show that megalopygid venom is produced in secretory cells that lie beneath the cuticle and are connected to the venom spines by canals. Megalopygid venoms consist of large aerolysin-like pore-forming toxins, which we have named megalysins, and a small number of peptides. The venom system differs markedly from those of previously studied venomous zygaenoids of the family Limacodidae, suggestive of an independent origin. Megalopygid venom potently activates mammalian sensory neurons via membrane permeabilization and induces sustained spontaneous pain behavior and paw swelling in mice. These bioactivities are ablated by treatment with heat, organic solvents, or proteases, indicating that they are mediated by larger proteins such as the megalysins. We show that the megalysins were recruited as venom toxins in the Megalopygidae following horizontal transfer of genes from bacteria to the ancestors of ditrysian Lepidoptera. Megalopygids have recruited aerolysin-like proteins as venom toxins convergently with centipedes, cnidarians, and fish. This study highlights the role of horizontal gene transfer in venom evolution.

Multi-Omic Identification of Venom Proteins Collected from Artificial Hosts of a Parasitoid Wasp

Habrobracon hebetor is a parasitoid wasp capable of infesting many lepidopteran larvae. It uses venom proteins to immobilize host larvae and prevent host larval development, thus playing an important role in the biocontrol of lepidopteran pests. To identify and characterize its venom proteins, we developed a novel venom collection method using an artificial host (ACV), i.e., encapsulated amino acid solution in paraffin membrane, allowing parasitoid wasps to inject venom. We performed protein full mass spectrometry analysis of putative venom proteins collected from ACV and venom reservoirs (VRs) (control). To verify the accuracy of proteomic data, we also collected venom glands (VGs), Dufour's glands (DGs) and ovaries (OVs), and performed transcriptome analysis. In this paper, we identified 204 proteins in ACV via proteomic analysis; compared ACV putative venom proteins with those identified in VG, VR, and DG via proteome and transcriptome approaches; and verified a set of them using quantitative real-time polymerase chain reaction. Finally, 201 ACV proteins were identified as potential venom proteins. In addition, we screened 152 and 148 putative venom proteins identified in the VG transcriptome and the VR proteome against those in ACV, and found only 26 and 25 putative venom proteins, respectively, were overlapped with those in ACV. Altogether, our data suggest proteome analysis of ACV in combination with proteome-transcriptome analysis of other organs/tissues will provide the most comprehensive identification of true venom proteins in parasitoid wasps.

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.

Genome assembly of the ectoparasitoid wasp Theocolax elegans

The ectoparasitoid wasp Theocolax elegans is a cosmopolitan and generalist pteromalid parasitoid of several major storage insect pests, and can effectively suppress a host population in warehouses. However, little molecular information about this wasp is currently available. In this study, we assembled the genome of T. elegans using PacBio long-read sequencing, Illumina sequencing, and Hi-C methods. The genome assembly is 662.73 Mb in length with contig and scaffold N50 values of 1.15 Mb and 88.8 Mb, respectively. The genome contains 56.4% repeat sequences and 23,212 protein-coding genes were annotated. Phylogenomic analyses revealed that T. elegans diverged from the lineage leading to subfamily Pteromalinae (Nasonia vitripennis and Pteromalus puparum) approximately 110.5 million years ago. We identified 130 significantly expanded gene families, 34 contracted families, 248 fast-evolving genes, and 365 positively selected genes in T. elegans. Additionally, 260 olfactory receptors and 285 venom proteins were identified. This genome assembly provides valuable genetic bases for future investigations on evolution, molecular biology and application of T. elegans.

Genome Assembly and Comparative Analysis of the Egg Parasitoid Wasp Trichogramma dendrolimi Shed Light on the Composition and Evolution of Olfactory Receptors and Venoms

Trichogramma dendrolimi is one of the most successfully industrialized Trichogramma species used to control agricultural and forestry pests in China. However, the molecular mechanisms underlying its host recognition and parasitism remain largely unknown, partially due to the limited genome information of this parasitoid wasp. Here, we present a high-quality de novo assembly of T. dendrolimi through a combination of Illumina and PacBio sequencing technologies. The final assembly had a length of 215.2 Mb and contains 316 scaffolds with a scaffold N50 size of 1.41 Mb. Repetitive sequences with a length of 63.4 Mb and 12,785 protein-coding genes were identified. Significantly expanded gene families were identified to be involved in the development and regulatory processes, while remarkably contracted gene families were involved in the transport processes in T. dendrolimi. The olfactory and venom-associated genes were identified in T. dendrolimi and 24 other hymenopteran species, using uniform methods combining BLAST and HMM profiling. The identified venom genes of T. dendrolimi were enriched in antioxidant activity, tricarboxylic acid cycle, response to oxidative stress and cell redox homeostasis. Our study provides an important resource for comparative genomics and functional studies to interpret the molecular mechanisms underlying host recognition and parasitism of Trichogramma species.

Proteinaceous Venom Expression of the Yellow Meadow Ant, Lasius flavus (Hymenoptera: Formicidae)

Ants are one of the important groups of venomous animals with about 14,000 described species. Studies so far focused on the discovery of venom proteins are only available for limited stinging ants, and the proteinaceous compositions of the stingless ants are completely unknown. Here, we used the transcriptomic approach to identify venom components from the yellow meadow ant, Lasius flavus, a stingless ant. The transcriptomic analysis yielded an extraordinary simplicity of the venom expression profile, with 17 venom proteins, such as phospholipase B, odorant binding protein, and apolipoprotein D. Ten of them were discovered as novel toxins for future functional investigations. Quantitative real time PCR analysis revealed that genes encoding the identified venom proteins display exclusively or highly expression profiles in venom glands, validating them as venom compositions. Our findings contribute to the understanding of the evolutional diversity of toxins between stinging and stingless ants.

Spatial Venomics─Cobra Venom System Reveals Spatial Differentiation of Snake Toxins by Mass Spectrometry Imaging

Among venomous animals, toxic secretions have evolved as biochemical weapons associated with various highly specialized delivery systems on many occasions. Despite extensive research, there is still limited knowledge of the functional biology of most animal toxins, including their venom production and storage, as well as the morphological structures within sophisticated venom producing tissues that might underpin venom modulation. Here, we report on the spatial exploration of a snake venom gland system by matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI), in combination with standard proteotranscriptomic approaches, to enable in situ toxin mapping in spatial intensity maps across a venom gland sourced from the Egyptian cobra (Naja haje). MALDI-MSI toxin visualization on the elapid venom gland reveals a high spatial heterogeneity of different toxin classes at the proteoform level, which may be the result of physiological constraints on venom production and/or storage that reflects the potential for venom modulation under diverse stimuli.

Co-option of a carotenoid cleavage dioxygenase gene (CCD4a) into the floral symmetry gene regulatory network contributes to the polymorphic floral shape-color combinations in Chrysanthemum sensu lato

Morphological novelties, including formation of trait combinations, may result from de novo gene origination and/or co-option of existing genes into other developmental contexts. A variety of shape-color combinations of capitular florets occur in Chrysanthemum and its allies. We hypothesized that co-option of a carotenoid cleavage dioxygenase gene into the floral symmetry gene network would generate a white zygomorphic ray floret. We tested this hypothesis in an evolutionary context using species in Chrysanthemum sensu lato, a monophyletic group with diverse floral shape-color combinations, based on morphological investigation, interspecific crossing, molecular interaction and transgenic experiments. Our results showed that white color was significantly associated with floret zygomorphy. Specific expression of the carotenoid cleavage dioxygenase gene CCD4a in marginal florets resulted in white color. Crossing experiments between Chrysanthemum lavandulifolium and Ajania pacifica indicated that expression of CCD4a is trans-regulated. The floral symmetry regulator CYC2g can activate expression of CCD4a with a dependence on TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING (TCP) binding element 8 on the CCD4a promoter. Based on all experimental findings, we propose that gene co-option of carotenoid degradation into floral symmetry regulation, and the subsequent dysfunction or loss of either CYC2g or CCD4a, may have led to evolution of capitular shape-color patterning in Chrysanthemum sensu lato.

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.

Mining Amphibian and Insect Transcriptomes for Antimicrobial Peptide Sequences with rAMPage

Antibiotic resistance is a global health crisis increasing in prevalence every day. To combat this crisis, alternative antimicrobial therapeutics are urgently needed. Antimicrobial peptides (AMPs), a family of short defense proteins, are produced naturally by all organisms and hold great potential as effective alternatives to small molecule antibiotics. Here, we present rAMPage, a scalable bioinformatics discovery platform for identifying AMP sequences from RNA sequencing (RNA-seq) datasets. In our study, we demonstrate the utility and scalability of rAMPage, running it on 84 publicly available RNA-seq datasets from 75 amphibian and insect species—species known to have rich AMP repertoires. Across these datasets, we identified 1137 putative AMPs, 1024 of which were deemed novel by a homology search in cataloged AMPs in public databases. We selected 21 peptide sequences from this set for antimicrobial susceptibility testing against Escherichia coli and Staphylococcus aureus and observed that seven of them have high antimicrobial activity. Our study illustrates how in silico methods such as rAMPage can enable the fast and efficient discovery of novel antimicrobial peptides as an effective first step in the strenuous process of antimicrobial drug development.

A chromosome-level genome assembly of the pollinating fig wasp Valisia javana

Fig wasp has always been thought the species-specific pollinator for their host fig (Moraceae, Ficus) and constitute a model system with its host to study co-evolution and co-speciation. The availability of a high-quality genome will help to further reveal the mechanisms underlying these characteristics. Here, we present a high-quality chromosome-level genome for Valisa javana developed by a combination of PacBio long-read and Illumina short-read. The assembled genome size is 296.34 Mb from 13 contigs with a contig N50 length of 26.76 kb. Comparative genomic analysis revealed expanded and positively selected genes related to biological features that aid fig wasps living in syconium of its highly specific host. Protein-coding genes associated with chemosensory, detoxification and venom genes were identified. Several differentially expressed genes in transcriptome data of V. javana between odor-stimulated samples and the controls have been identified in some olfactory signal transduction pathways, e.g. olfactory transduction, cAMP, cGMP-PKG, Calcim, Ras and Rap1. This study provides a valuable genomic resource for a fig wasp, and sheds insight into further revealing the mechanisms underlying their adaptive traits to their hosts in different places and co-speciation with their host.

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.

Research gaps and new insights in the evolution of Drosophila seminal fluid proteins

While the striking effects of seminal fluid proteins (SFPs) on females are fairly conserved among Diptera, most SFPs lack detectable homologues among the SFP repertoires of phylogenetically distant species. How such a rapidly changing proteome conserves functions across taxa is a fascinating question. However, this and other pivotal aspects of SFPs' evolution remain elusive because discoveries on these proteins have been mainly restricted to the model Drosophila melanogaster. Here, we provide an overview of the current knowledge on the inter-specific divergence of the SFP repertoire in Drosophila and compile the increasing amount of relevant genomic information from multiple species. Capitalizing on the accumulated knowledge in D. melanogaster, we present novel sets of high-confidence SFP candidates and transcription factors presumptively involved in regulating the expression of SFPs. We also address open questions by performing comparative genomic analyses that failed to support the existence of many conserved SFPs shared by most dipterans and indicated that gene co-option is the most frequent mechanism accounting for the origin of Drosophila SFP-coding genes. We hope our update establishes a starting point to integrate further data and thus widen the understanding of the intricate evolution of these proteins.

Genetic, morphometric, and molecular analyses of interspecies differences in head shape and hybrid developmental defects in the wasp genus Nasonia

Males in the parasitoid wasp genus Nasonia have distinct, species-specific, head shapes. The availability of fertile hybrids among the species, along with obligate haploidy of males, facilitates analysis of complex gene interactions in development and evolution. Previous analyses showed that both the divergence in head shape between Nasonia vitripennis and Nasonia giraulti, and the head-specific developmental defects of F2 haploid hybrid males, are governed by multiple changes in networks of interacting genes. Here, we extend our understanding of the gene interactions that affect morphogenesis in male heads. Use of artificial diploid male hybrids shows that alleles mediating developmental defects are recessive, while there are diverse dominance relationships among other head shape traits. At the molecular level, the sex determination locus doublesex plays a major role in male head shape differences, but it is not the only important factor. Introgression of a giraulti region on chromsome 2 reveals a recessive locus that causes completely penetrant head clefting in both males and females in a vitripennis background. Finally, a third species (N. longicornis) was used to investigate the timing of genetic changes related to head morphology, revealing that most changes causing defects arose after the divergence of N. vitripennis from the other species, but prior to the divergence of N. giraulti and N. longicornis from each other. Our results demonstrate that developmental gene networks can be dissected using interspecies crosses in Nasonia, and set the stage for future fine-scale genetic dissection of both head shape and hybrid developmental defects.

Long-Read Assembly and Annotation of the Parasitoid Wasp Muscidifurax raptorellus, a Biological Control Agent for Filth Flies

The parasitoid wasp Muscidifurax raptorellus (Hymenoptera: Pteromalidae) is a gregarious species that has received extensive attention for its potential in biological pest control against house fly, stable fly, and other filth flies. It has a high reproductive capacity and can be reared easily. However, genome assembly is not available for M. raptorellus or any other species in this genus. Previously, we assembled a complete circular mitochondrial genome with a length of 24,717 bp. Here, we assembled and annotated a high-quality nuclear genome of M. raptorellus, using a combination of long-read (104× genome coverage) and short-read (326× genome coverage) sequencing technologies. The assembled genome size is 314 Mbp in 226 contigs, with a 97.9% BUSCO completeness score and a contig N50 of 4.67 Mb, suggesting excellent continuity of this assembly. Our assembly builds the foundation for comparative and evolutionary genomic analysis in the genus of Muscidifurax and possible future biocontrol applications.

The emergence of ecotypes in a parasitoid wasp: a case of incipient sympatric speciation in Hymenoptera?

Background

To understand which reproductive barriers initiate speciation is a major question in evolutionary research. Despite their high species numbers and specific biology, there are only few studies on speciation in Hymenoptera. This study aims to identify very early reproductive barriers in a local, sympatric population of Nasonia vitripennis (Walker 1836), a hymenopterous parasitoid of fly pupae. We studied ecological barriers, sexual barriers, and the reduction in F1-female offspring as a postmating barrier, as well as the population structure using microsatellites.

Results

We found considerable inbreeding within female strains and a population structure with either three or five subpopulation clusters defined by microsatellites. In addition, there are two ecotypes, one parasitizing fly pupae in bird nests and the other on carrion. The nest ecotype is mainly formed from one of the microsatellite clusters, the two or four remaining microsatellite clusters form the carrion ecotype. There was slight sexual isolation and a reduction in F1-female offspring between inbreeding strains from the same microsatellite clusters and the same ecotypes. Strains from different microsatellite clusters are separated by a reduction in F1-female offspring. Ecotypes are separated only by ecological barriers.

Conclusions

This is the first demonstration of very early reproductive barriers within a sympatric population of Hymenoptera. It demonstrates that sexual and premating barriers can precede ecological separation. This indicates the complexity of ecotype formation and highlights the general need for more studies within homogenous populations for the identification of the earliest barriers in the speciation process.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01938-y.

A Short Review of the Venoms and Toxins of Spider Wasps (Hymenoptera: Pompilidae)

Parasitoid wasps represent the plurality of venomous animals, but have received extremely little research in proportion to this taxonomic diversity. The lion’s share of investigation into insect venoms has focused on eusocial hymenopterans, but even this small sampling shows great promise for the development of new active substances. The family Pompilidae is known as the spider wasps because of their reproductive habits which include hunting for spiders, delivering a paralyzing sting, and entombing them in burrows with one of the wasp’s eggs to serve as food for the developing larva. The largest members of this family, especially the tarantula hawks of the genus Pepsis, have attained notoriety for their large size, dramatic coloration, long-term paralysis of their prey, and incredibly painful defensive stings. In this paper we review the existing research regarding the composition and function of pompilid venoms, discuss parallels from other venom literatures, identify possible avenues for the adaptation of pompilid toxins towards human purposes, and future directions of inquiry for the field.

Neofunctionalization of an ancient domain allows parasites to avoid intraspecific competition by manipulating host behaviour

Intraspecific competition is a major force in mediating population dynamics, fuelling adaptation, and potentially leading to evolutionary diversification. Among the evolutionary arms races between parasites, one of the most fundamental and intriguing behavioural adaptations and counter-adaptations are superparasitism and superparasitism avoidance. However, the underlying mechanisms and ecological contexts of these phenomena remain underexplored. Here, we apply the Drosophila parasite Leptopilina boulardi as a study system and find that this solitary endoparasitic wasp provokes a host escape response for superparasitism avoidance. We combine multi-omics and in vivo functional studies to characterize a small set of RhoGAP domain-containing genes that mediate the parasite's manipulation of host escape behaviour by inducing reactive oxygen species in the host central nervous system. We further uncover an evolutionary scenario in which neofunctionalization and specialization gave rise to the novel role of RhoGAP domain in avoiding superparasitism, with an ancestral origin prior to the divergence between Leptopilina specialist and generalist species. Our study suggests that superparasitism avoidance is adaptive for a parasite and adds to our understanding of how the molecular manipulation of host behaviour has evolved in this system.

Tissue-specific gene expression shows a cynipid wasp repurposes oak host gene networks to create a complex and novel parasite-specific organ

Every organism on Earth depends on interactions with other organisms to survive. In each of these interactions, an organism must utilize the limited toolbox of genes and proteins it possesses to successfully manipulate or cooperate with another species, but it can also co-opt the genome machinery of its partner to expand its available tools. Insect-induced plant galls are an extreme example of this, wherein an insect hijacks the plant's genome to direct the initiation and development of galls consisting of plant tissue. However, previous transcriptomic studies have not evaluated individual tissues within a gall to determine the full extent to which a galling insect manipulates its host plant. Here we demonstrate that the cynipid wasp Dryocosmus quercuspalustris creates a complex parasite-specific organ from red oak tissue via massive changes in host gene expression. Our results show that the gall wasp is not merely modifying oak leaf tissue but creating extensive changes in gene expression between galled and ungalled tissue (differential expression in 28% of genes) and distinct gall tissue types (20% of genes). The outer gall tissue shows increases in various plant defence systems, which is consistent with its predicted functional role of protecting the wasp larva. The inner larval capsule shows suppression of large parts of the plant innate immune system and evidence for the wasp utilizing the plant's RNA interference mechanisms, which may be a potential mechanism for the wasp's control on gall growth.

Rapid Genomic Evolution Drives the Diversification of Male Reproductive Genes in Dung Beetles

The molecular basis for the evolution of novel phenotypes is a central question in evolutionary biology. In recent years, dung beetles have emerged as models for novel trait evolution as they possess distinct precopulatory traits such as sexually dimorphic horns on their head and thorax. Here, we use functional and evolutionary genomics to investigate the origins and the evolution of postcopulatory reproductive traits in male dung beetles. Male ejaculates that underlie postcopulatory sexual selection are excellent candidates to study novel trait evolution as they are complex, fast evolving, and often highly divergent in insects. We assemble de novo transcriptomes of male accessory glands and testes of a widespread dung beetle, Catharsius molossus, and we perform an evolutionary analysis of closely and distantly related insect genomes. Our results show there is rapid innovation at the genomic level even among closely related dung beetles. Genomic expansion and contraction drive the divergence of male reproductive traits and their functions. The birth of scores of completely novel reproductive genes is reinforced by the recruitment of these genes for high expression in male reproductive tissues, especially in the accessory glands. We find that male accessory glands of C. molossus are specialized for secretory function and express female, egg, and embryo-related genes as well as serine protease inhibitors, whilst the testes are specialized for spermatogenesis and sperm function. Finally, we touch upon putative functions of these evolutionary novelties using structure-function analysis as these proteins bear no homology to any other known proteins.

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

Background

Hymenoptera comprise extremely diverse insect species with extensive variation in their life histories. The Dryinidae, a family of solitary wasps of Hymenoptera, have evolved innovations that allow them to hunt using venom and a pair of chelae developed from the fore legs that can grasp prey. Dryinidae larvae are also parasitoids of Auchenorrhyncha, a group including common pests such as planthoppers and leafhoppers. Both of these traits make them effective and valuable for pest control, but little is yet known about the genetic basis of its dual adaptation to parasitism and predation.

Results

We sequenced and assembled a high-quality genome of the dryinid wasp Gonatopus flavifemur, which at 636.5 Mb is larger than most hymenopterans. The expansion of transposable elements, especially DNA transposons, is a major contributor to the genome size enlargement. Our genome-wide screens reveal a number of positively selected genes and rapidly evolving proteins involved in energy production and motor activity, which may contribute to the predatory adaptation of dryinid wasp. We further show that three female-biased, reproductive-associated yellow genes, in response to the prey feeding behavior, are significantly elevated in adult females, which may facilitate the egg production. Venom is a powerful weapon for dryinid wasp during parasitism and predation. We therefore analyze the transcriptomes of venom glands and describe specific expansions in venom Idgf-like genes and neprilysin-like genes. Furthermore, we find the LWS2-opsin gene is exclusively expressed in male G. flavifemur, which may contribute to partner searching and mating.

Conclusions

Our results provide new insights into the genome evolution, predatory adaptation, venom evolution, and sex-biased genes in G. flavifemur, and present genomic resources for future in-depth comparative analyses of hymenopterans that may benefit pest control.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01081-6.

Utilisation of compounds from venoms in drug discovery

Difficult drug targets are becoming the normal course of business in drug discovery, sometimes due to large interacting surfaces or only small differences in selectivity regions. For these, a different approach is merited: compounds lying somewhere between the small molecule and the large antibody in terms of many properties including stability, biodistribution and pharmacokinetics. Venoms have evolved over millions of years to be complex mixtures of stable molecules derived from other somatic molecules, the stability comes from the pressure to be ready for delivery at a moment's notice. Snakes, spiders, scorpions, jellyfish, wasps, fish and even mammals have evolved independent venom systems with complex mixtures in their chemical arsenal. These venom-derived molecules have been proven to be useful tools, such as for the development of antihypotensive angiotensin converting enzyme (ACE) inhibitors and have also made successful drugs such as Byetta® (Exenatide), Integrilin® (Eptifibatide) and Echistatin. Only a small percentage of the available chemical space from venoms has been investigated so far and this is growing. In a new era of biological therapeutics, venom peptides present opportunities for larger target engagement surface with greater stability than antibodies or human peptides. There are challenges for oral absorption and target engagement, but there are venom structures that overcome these and thus provide substrate for engineering novel molecules that combine all desired properties. Venom researchers are characterising new venoms, species, and functions all the time, these provide great substrate for solving the challenges presented by today's difficult targets.

The evolutionary dynamics of venom toxins made by insects and other animals

Animal venoms are recognised as unique biological systems in which to study molecular evolution. Venom use has evolved numerous times among the insects, and insects today use venom to capture prey, defend themselves from predators, or to subdue and modulate host responses during parasitism. However, little is known about most insect venom toxins or the mode and tempo by which they evolve. Here, I review the evolutionary dynamics of insect venom toxins, and argue that insects offer many opportunities to examine novel aspects of toxin evolution. The key questions addressed are: How do venomous animals evolve from non-venomous animals, and how does this path effect the composition and pharmacology of the venom? What genetic processes (gene duplication, co-option, neofunctionalisation) are most important in toxin evolution? What kinds of selection pressures are acting on toxin-encoding genes and their cognate targets in envenomated animals? The emerging evidence highlights that venom composition and pharmacology adapts quickly in response to changing selection pressures resulting from new ecological interactions, and that such evolution occurs through a stunning variety of genetic mechanisms. Insects offer many opportunities to investigate the evolutionary dynamics of venom toxins due to their evolutionary history rich in venom-related adaptations, and their quick generation time and suitability for culture in the laboratory.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Coadaptation between host genome and microbiome under long-term xenobiotic-induced selection

One of the most difficult experimental challenges today is testing the evolutionary dynamics shaping complex host-microbiome interactions. We investigated host-microbiome codiversification in response to xenobiotic-induced selection using an experimental evolution approach. To this end, we exposed the parasitoid wasp Nasonia vitripennis to sublethal concentrations of the widely used herbicide atrazine for 85 generations. Our results reveal that atrazine exposure not only mediated adaptive changes within the microbiome, which conferred host resistance to atrazine toxicity, but also exerted selective pressure on the host genome and altered host gene expression and immune response. Furthermore, microbiome transplant experiments reveal a decreased survival of adults from the control population after exposure to the evolved microbiome of the atrazine-exposed population, while no such decrease occurred in the reciprocal transplant. These results indicate that xenobiotic-induced selection mediated host-microbiome coadaptation, ultimately leading to a new host genome-microbiome equilibrium.

Production, composition, and mode of action of the painful defensive venom produced by a limacodid caterpillar, Doratifera vulnerans

Venoms have evolved independently several times in Lepidoptera. Limacodidae is a family with worldwide distribution, many of which are venomous in the larval stage, but the composition and mode of action of their venom is unknown. Here, we use imaging technologies, transcriptomics, proteomics, and functional assays to provide a holistic picture of the venom system of a limacodid caterpillar, Doratifera vulnerans Contrary to dogma that defensive venoms are simple in composition, D. vulnerans produces a complex venom containing 151 proteinaceous toxins spanning 59 families, most of which are peptides <10 kDa. Three of the most abundant families of venom peptides (vulnericins) are 1) analogs of the adipokinetic hormone/corazonin-related neuropeptide, some of which are picomolar agonists of the endogenous insect receptor; 2) linear cationic peptides derived from cecropin, an insect innate immune peptide that kills bacteria and parasites by disrupting cell membranes; and 3) disulfide-rich knottins similar to those that dominate spider venoms. Using venom fractionation and a suite of synthetic venom peptides, we demonstrate that the cecropin-like peptides are responsible for the dominant pain effect observed in mammalian in vitro and in vivo nociception assays and therefore are likely to cause pain after natural envenomations by D. vulnerans Our data reveal convergent molecular evolution between limacodids, hymenopterans, and arachnids and demonstrate that lepidopteran venoms are an untapped source of novel bioactive peptides.

Variation in Parasitoid Virulence of Tetrastichus brontispae during the Targeting of Two Host Beetles

In host-parasitoid interactions, antagonistic relationship drives parasitoids to vary in virulence in facing different hosts, which makes these systems excellent models for stress-induced evolutionary studies. Venom compositions varied between two strains of Tetrastichus brontispae, Tb-Bl and Tb-On. Tb-Bl targets Brontispa longissima pupae as hosts, and Tb-On is a sub-population of Tb-Bl, which has been experimentally adapted to a new host, Octodonta nipae. Aiming to examine variation in parasitoid virulence of the two strains toward two hosts, we used reciprocal injection experiments to compare effect of venom/ovarian fluids from the two strains on cytotoxicity, inhibition of immunity and fat body lysis of the two hosts. We found that Tb-Onvenom was more virulent towards plasmatocyte spreading, granulocyte function and phenoloxidase activity than Tb-Blvenom. Tb-Blovary was able to suppress encapsulation and phagocytosis in both hosts; however, Tb-Onovary inhibition targeted only B. longissima. Our data suggest that the venom undergoes rapid evolution when facing different hosts, and that the wasp has good evolutionary plasticity.

A chimeric gene paternally instructs female sex determination in the haplodiploid wasp Nasonia

Various primary signals direct insect sex determination. In hymenopteran insects, the presence of a paternal genome is needed to initiate female development. When absent, uniparental haploid males develop. We molecularly and functionally identified the instructor sex-determination gene, wasp overruler of masculinization (wom), of the haplodiploid wasp Nasonia vitripennis This gene contains a P53-like domain coding region and arose by gene duplication and genomic rearrangements. Maternal silencing of wom results in male development of haploid embryos. Upon fertilization, early zygotic transcription from the paternal wom allele is initiated, followed by a timely zygotic expression of transformer (tra), leading to female development. Wom is an instructor gene with a parent-of-origin effect in sex determination.

Two novel venom proteins underlie divergent parasitic strategies between a generalist and a specialist parasite

Parasitoids are ubiquitous in natural ecosystems. Parasitic strategies are highly diverse among parasitoid species, yet their underlying genetic bases are poorly understood. Here, we focus on the divergent adaptation of a specialist and a generalist drosophilid parasitoids. We find that a novel protein (Lar) enables active immune suppression by lysing the host lymph glands, eventually leading to successful parasitism by the generalist. Meanwhile, another novel protein (Warm) contributes to a passive strategy by attaching the laid eggs to the gut and other organs of the host, leading to incomplete encapsulation and helping the specialist escape the host immune response. We find that these diverse parasitic strategies both originated from lateral gene transfer, followed with duplication and specialization, and that they might contribute to the shift in host ranges between parasitoids. Our results increase our understanding of how novel gene functions originate and how they contribute to host adaptation.

How the Toxin got its Toxicity

Venom systems are functional and ecological traits, typically used by one organism to subdue or deter another. A predominant subset of their constituent molecules—“toxins”—share this ecological function and are therefore molecules that mediate interactions between organisms. Such molecules have been referred to as “exochemicals.” There has been debate within the field of toxinology concerning the evolutionary pathways leading to the “recruitment” of a gene product for a toxic role within venom. We review these discussions and the evidence interpreted in support of alternate pathways, along with many of the most popular models describing the origin of novel molecular functions in general. We note that such functions may arise with or without gene duplication occurring and are often the consequence of a gene product encountering a novel “environment,” i.e., a range of novel partners for molecular interaction. After stressing the distinction between “activity” and “function,” we describe in detail the results of a recent study which reconstructed the evolutionary history of a multigene family that has been recruited as a toxin and argue that these results indicate that a pluralistic approach to understanding the origin of novel functions is advantageous. This leads us to recommend that an expansive approach be taken to the definition of “neofunctionalization”—simply the origins of a novel molecular function by any process—and “recruitment”—the “weaponization” of a molecule via the acquisition of a toxic function in venom, by any process. Recruitment does not occur at the molecular level or even at the level of gene expression, but only when a confluence of factors results in the ecological deployment of a physiologically active molecule as a toxin. Subsequent to recruitment, the evolutionary regime of a gene family may shift into a more dynamic form of “birth-and-death.” Thus, recruitment leads to a form of “downwards causation,” in which a change at the ecological level at which whole organisms interact leads to a change in patterns of evolution at the genomic level.

Comparative analysis reveals the expansion of mitochondrial DNA control region containing unusually high G-C tandem repeat arrays in Nasonia vitripennis

Insect mitochondrial DNA (mtDNA) ranges from 14 to 19 kbp, and the size difference is attributed to the AT-rich control region. Jewel wasps have a parasitoid lifestyle, which may affect mitochondria function and evolution. We sequenced, assembled, and annotated mitochondrial genomes in Nasonia and outgroup species. Gene composition and order are conserved within Nasonia, but they differ from other parasitoids by two large inversion events that were not reported before. We observed a much higher substitution rate relative to the nuclear genome and mitochondrial introgression between N. giraulti and N. oneida, which is consistent with previous studies. Most strikingly, N. vitripennis mtDNA has an extremely long control region (7665 bp), containing twenty-nine 217 bp tandem repeats and can fold into a super-cruciform structure. In contrast to tandem repeats commonly found in other mitochondria, these high-copy repeats are highly conserved (98.7% sequence identity), much longer in length (approximately 8 Kb), extremely GC-rich (50.7%), and CpG-rich (percent CpG 19.4% vs. 1.1% in coding region), resulting in a 23 kbp mtDNA beyond the typical size range in insects. These N. vitripennis-specific mitochondrial repeats are not related to any known sequences in insect mitochondria. Their evolutionary origin and functional consequences warrant further investigations.

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.

Polymorphism and Divergence of Novel Gene Expression Patterns in Drosophila melanogaster

Transcriptomes may evolve by multiple mechanisms, including the evolution of novel genes, the evolution of transcript abundance, and the evolution of cell, tissue, or organ expression patterns. Here, we focus on the last of these mechanisms in an investigation of tissue and organ shifts in gene expression in Drosophila melanogaster. In contrast to most investigations of expression evolution, we seek to provide a framework for understanding the mechanisms of novel expression patterns on a short population genetic timescale. To do so, we generated population samples of D. melanogaster transcriptomes from five tissues: accessory gland, testis, larval salivary gland, female head, and first-instar larva. We combined these data with comparable data from two outgroups to characterize gains and losses of expression, both polymorphic and fixed, in D. melanogaster We observed a large number of gain- or loss-of-expression phenotypes, most of which were polymorphic within D. melanogaster Several polymorphic, novel expression phenotypes were strongly influenced by segregating cis-acting variants. In support of previous literature on the evolution of novelties functioning in male reproduction, we observed many more novel expression phenotypes in the testis and accessory gland than in other tissues. Additionally, genes showing novel expression phenotypes tend to exhibit greater tissue-specific expression. Finally, in addition to qualitatively novel expression phenotypes, we identified genes exhibiting major quantitative expression divergence in the D. melanogaster lineage.

Genome Report: Whole Genome Sequence and Annotation of the Parasitoid Jewel Wasp Nasonia giraulti Laboratory Strain RV2X[u]

Jewel wasps in the genus of Nasonia are parasitoids with haplodiploidy sex determination, rapid development and are easy to culture in the laboratory. They are excellent models for insect genetics, genomics, epigenetics, development, and evolution. Nasonia vitripennis (Nv) and N. giraulti (Ng) are closely-related species that can be intercrossed, particularly after removal of the intracellular bacterium Wolbachia, which serve as a powerful tool to map and positionally clone morphological, behavioral, expression and methylation phenotypes. The Nv reference genome was assembled using Sanger, PacBio and Nanopore approaches and annotated with extensive RNA-seq data. In contrast, Ng genome is only available through low coverage resequencing. Therefore, de novo Ng assembly is in urgent need to advance this system. In this study, we report a high-quality Ng assembly using 10X Genomics linked-reads with 670X sequencing depth. The current assembly has a genome size of 259,040,977 bp in 3,160 scaffolds with 38.05% G-C and a 98.6% BUSCO completeness score. 97% of the RNA reads are perfectly aligned to the genome, indicating high quality in contiguity and completeness. A total of 14,777 genes are annotated in the Ng genome, and 72% of the annotated genes have a one-to-one ortholog in the Nv genome. We reported 5 million Ng-Nv SNPs which will facility mapping and population genomic studies in Nasonia In addition, 42 Ng-specific genes were identified by comparing with Nv genome and annotation. This is the first de novo assembly for this important species in the Nasonia model system, providing a useful new genomic toolkit.