Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial 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.

Dynamic Evolution of Repetitive Elements and Chromatin States in Apis mellifera Subspecies

In this study, we elucidate the contribution of repetitive DNA sequences to the establishment of social structures in honeybees (Apis mellifera). Despite recent advancements in understanding the molecular mechanisms underlying the formation of honeybee castes, primarily associated with Notch signaling, the comprehensive identification of specific genomic cis-regulatory sequences remains elusive. Our objective is to characterize the repetitive landscape within the genomes of two honeybee subspecies, namely A. m. mellifera and A. m. ligustica. An observed recent burst of repeats in A. m. mellifera highlights a notable distinction between the two subspecies. After that, we transitioned to identifying differentially expressed DNA elements that may function as cis-regulatory elements. Nevertheless, the expression of these sequences showed minimal disparity in the transcriptome during caste differentiation, a pivotal process in honeybee eusocial organization. Despite this, chromatin segmentation, facilitated by ATAC-seq, ChIP-seq, and RNA-seq data, revealed a distinct chromatin state associated with repeats. Lastly, an analysis of sequence divergence among elements indicates successive changes in repeat states, correlating with their respective time of origin. Collectively, these findings propose a potential role of repeats in acquiring novel regulatory functions.

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.

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.

Genomic evidence for contrasting patterns of host-associated genetic differentiation across shared host-plant species in leaf- and bud-galling sawflies

Resource specialization and ecological speciation arising through host-associated genetic differentiation (HAD) are frequently invoked as an explanation for the high diversity of plant-feeding insects and other organisms with a parasitic lifestyle. While genetic studies have demonstrated numerous examples of HAD in insect herbivores, the rarity of comparative studies means that we still lack an understanding of how deterministic HAD is, and whether patterns of host shifts can be predicted over evolutionary timescales. We applied genome-wide single nucleotide polymorphism and mitochondrial DNA sequence data obtained through genome resequencing to define species limits and to compare host-plant use in population samples of leaf- and bud-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae) collected from seven shared willow (Salicaceae: Salix) host species. To infer the repeatability of long-term cophylogenetic patterns, we also contrasted the phylogenies of the two galler groups with each other as well as with the phylogeny of their Salix hosts estimated based on RADseq data. We found clear evidence for host specialization and HAD in both of the focal galler groups, but also that leaf gallers are more specialized to single host species compared with most bud gallers. In contrast to bud gallers, leaf gallers also exhibited statistically significant cophylogenetic signal with their Salix hosts. The observed discordant patterns of resource specialization and host shifts in two related galler groups that have radiated in parallel across a shared resource base indicate a lack of evolutionary repeatability in the focal system, and suggest that short- and long-term host use and ecological diversification in plant-feeding insects are dominated by stochasticity and/or lineage-specific effects.

Mblk-1/E93, an ecdysone related-transcription factor, targets synaptic plasticity-related genes in the honey bee mushroom bodies

Among hymenopteran insects, aculeate species such as bees, ants, and wasps have enlarged and morphologically elaborate mushroom bodies (MBs), a higher-order brain center in the insect, implying their relationship with the advanced behavioral traits of aculeate species. The molecular bases leading to the acquisition of complicated MB functions, however, remains unclear. We previously reported the constitutive and MB-preferential expression of an ecdysone-signaling related transcription factor, Mblk-1/E93, in the honey bee brain. Here, we searched for target genes of Mblk-1 in the worker honey bee MBs using chromatin immunoprecipitation sequence analyses and found that Mblk-1 targets several genes involved in synaptic plasticity, learning, and memory abilities. We also demonstrated that Mblk-1 expression is self-regulated via Mblk-1-binding sites, which are located upstream of Mblk-1. Furthermore, we showed that the number of the Mblk-1-binding motif located upstream of Mblk-1 homologs increased associated with evolution of hymenopteran insects. Our findings suggest that Mblk-1, which has been focused on as a developmental gene transiently induced by ecdysone, has acquired a novel expression pattern to play a role in synaptic plasticity in honey bee MBs, raising a possibility that molecular evolution of Mblk-1 may have partly contributed to the elaboration of MB function in insects.

Persistent Parental RNAi in the Beetle Tribolium castaneum Involves Maternal Transmission of Long Double-Stranded RNA

Parental RNA interference (pRNAi) is a powerful and widely used method for gene-specific knockdown. Yet in insects its efficacy varies between species, and how the systemic response is transmitted from mother to offspring remains elusive. Using the beetle Tribolium castaneum, an RT-qPCR strategy to distinguish the presence of double-stranded RNA (dsRNA) from endogenous mRNA is reported. It is found that injected dsRNA is directly transmitted into the egg and persists throughout embryogenesis. Despite this depletion of dsRNA from the mother, it is shown that strong pRNAi can persist for months before waning at strain-specific rates. In seeking the receptor proteins for cellular uptake of long dsRNA into the egg, a phylogenomics profiling approach of candidate proteins is also presented. A visualization strategy based on taxonomically hierarchical assessment of orthology clustering data to rapidly assess gene age and copy number changes, refined by sequence-based evidence, is demonstrated. Repeated losses of SID-1-like channel proteins in the arthropods, including wholesale loss in the Heteroptera (true bugs), which are nonetheless highly sensitive to pRNAi, are thereby documented. Overall, practical considerations for insect pRNAi against a backdrop of outstanding questions on the molecular mechanism of dsRNA transmission for long-term, systemic knockdown are elucidated.

Venom profile of the European carpenter bee Xylocopa violacea: Evolutionary and applied considerations on its toxin components

Modern venomics is increasing its focus on hymenopterans such as honeybees, bumblebees, parasitoid wasps, ants and true wasps. However solitary bees remain understudied in comparison and the few available venom studies focus on short melittin-like sequences and antimicrobial peptides. Herein we describe the first comprehensive venom profile of a solitary bee, the violet carpenter bee Xylocopa violacea, by using proteo-transcriptomics. We reveal a diverse and complex venom profile with 43 different protein families identified from dissected venom gland extracts of which 32 are also detected in the defensively injected venom. Melittin and apamin are the most highly secreted components, followed by Phospholipase A2, Icarapin, Secapin and three novel components. Other components, including eight novel protein families, are rather lowly expressed. We further identify multiple forms of apamin-like peptides. The melittin-like sequences of solitary bees separate into two clades, one comprised most sequences from solitary bees including xylopin (the variant in Xylocopa), while sequences from Lasioglossa appear closer related to melittin-like peptides from Bombus (Bombolittins). Our study suggests that more proteo-transcriptomic data from other solitary bees should be complemented with corresponding genome data to fully understand the evolution and complexity of bee venom proteins, and is of a particular need to disentangle the ambiguous phylogenetic relations of short peptides.

Hymenoptera Genome Database: new genomes and annotation datasets for improved go enrichment and orthologue analyses

We report an update of the Hymenoptera Genome Database (HGD; http://HymenopteraGenome.org), a genomic database of hymenopteran insect species. The number of species represented in HGD has nearly tripled, with fifty-eight hymenopteran species, including twenty bees, twenty-three ants, eleven wasps and four sawflies. With a reorganized website, HGD continues to provide the HymenopteraMine genomic data mining warehouse and JBrowse/Apollo genome browsers integrated with BLAST. We have computed Gene Ontology (GO) annotations for all species, greatly enhancing the GO annotation data gathered from UniProt with more than a ten-fold increase in the number of GO-annotated genes. We have also generated orthology datasets that encompass all HGD species and provide orthologue clusters for fourteen taxonomic groups. The new GO annotation and orthology data are available for searching in HymenopteraMine, and as bulk file downloads.

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.

Midgut transcriptome assessment of the cockroach-hunting wasp Ampulex compressa (Apoidea: Ampulicidae)

The emerald jewel wasp Ampulex compressa (Hymenoptera: Ampulicidae) is a solitary wasp that is widely known for its specialized hunting of cockroaches as larvae provision. Adult wasps mainly feed on pollen and nectar, while their larvae feed on the cockroachs’ body, first as ecto- and later as endoparsitoids. Little is known about the expression of digestive, detoxification and stress-response-related genes in the midgut of A. compressa, or about its transcriptional versatility between life stages. To identify gut-biased genes related to digestion, detoxification, and stress response, we explored the midgut transcriptome of lab-reared A. compressa, for both adults and larvae, by focusing on the top 100 significantly up- and down-regulated genes. From the top 100 significantly differentially expressed genes (DEGs), we identified 39 and 36 DEGs putatively related to digestion and detoxification in the adult wasps and larvae, respectively. The two carbohydrases alpha-glucosidase (containing an alpha-amylase domain) and glycosyl hydrolase family 31, as well as the two proteinases chymotrypsin and trypsin, revealed the highest gene diversity. We identified six significant DEGs related to detoxification, which comprise glutathione S-transferase, cytochrome P450s and UDP-glucuronosyltransferase. The gene expression levels that were significantly expressed in both life stages vary strongly between life stages, as found in genes encoding for chymotrypsin and trypsin or glycosyl hydrolases family 31. The number of genes related to alpha-glucosidase, glycosyl hydrolase family 31, and cytochrome P450s was found to be similar across nine reference hymenopteran species, except for the identified glycosyl hydrolase family 31 gene, which was absent in all reference bee species. Phylogenetic analyses of the latter candidate genes revealed that they cluster together with their homologous genes found in the reference hymenopteran species. These identified candidate genes provide a basis for future comparative genomic and proteomic studies on (ontogenetic) dietary transitions in Hymenoptera.

Opsin Evolution in Flower-Visiting Beetles

Flowers have evolved signals that exploit the sensory systems of insect visitors. In the case of visual cues, color signals are thought to have been shaped in large part by the spectral sensitivity of key pollinators, such as hymenopterans. Beetles were some of the first plant pollinators, pre-dating the angiosperm radiation but with the exception of a few well-studied species, the evolution of flower-visiting beetle visual systems is poorly understood. Thus, the ability of beetles to detect and distinguish flower color signals and perhaps their potential role in shaping flower coloration is not well understood. Traditional models of pollinator visual systems often assume a putative tri- or tetrachromatic flower-visitor, as is found in bees, flies and butterflies. Beetles are unique among modern pollinators as ancestrally they did not possess the machinery for trichromatic vision, lacking the blue-sensitive photoreceptor class. Research on the evolution of visual genes responsible for wavelength sensitivity (opsins) has revealed that beetles with putative tri- and tetrachromatic visual systems have evolved independently, along multiple lineages. We explore the evolution of beetle visual genes using newly generated and publicly available RNA-seq data from 25 species with flower associations, including previously unexplored key flower-visitor groups and 20 non-flower visiting relatives. Our findings serve as a resource to inform and guide future studies on beetle-flower interactions, where insight from both signal and receiver is needed to better understand these poorly explored systems.

Genomes of the willow-galling sawflies Euura lappo and Eupontania aestiva (Hymenoptera: Tenthredinidae): a resource for research on ecological speciation, adaptation, and gall induction

Hymenoptera is a hyperdiverse insect order represented by over 153,000 different species. As many hymenopteran species perform various crucial roles for our environments, such as pollination, herbivory, and parasitism, they are of high economic and ecological importance. There are 99 hymenopteran genomes in the NCBI database, yet only five are representative of the paraphyletic suborder Symphyta (sawflies, woodwasps, and horntails), while the rest represent the suborder Apocrita (bees, wasps, and ants). Here, using a combination of 10X Genomics linked-read sequencing, Oxford Nanopore long-read technology, and Illumina short-read data, we assembled the genomes of two willow-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae: Euurina): the bud-galling species Euura lappo and the leaf-galling species Eupontania aestiva. The final assembly for E. lappo is 259.85 Mbp in size, with a contig N50 of 209.0 kbp and a BUSCO score of 93.5%. The E. aestiva genome is 222.23 Mbp in size, with a contig N50 of 49.7 kbp and a 90.2% complete BUSCO score. De novo annotation of repetitive elements showed that 27.45% of the genome was composed of repetitive elements in E. lappo and 16.89% in E. aestiva, which is a marked increase compared to previously published hymenopteran genomes. The genomes presented here provide a resource for inferring phylogenetic relationships among basal hymenopterans, comparative studies on host-related genomic adaptation in plant-feeding insects, and research on the mechanisms of plant manipulation by gall-inducing insects.

Chemoreceptor Diversity in Apoid Wasps and Its Reduction during the Evolution of the Pollen-Collecting Lifestyle of Bees (Hymenoptera: Apoidea)

Chemoreceptors help insects to interact with their environment, to detect and assess food sources and oviposition sites, and to aid in intra- and interspecific communication. In Hymenoptera, species of eusocial lineages possess large chemoreceptor gene repertoires compared with solitary species, possibly because of their additional need to recognize nest-mates and caste. However, a critical piece of information missing so far has been the size of chemoreceptor gene repertoires of solitary apoid wasps. Apoid wasps are a paraphyletic group of almost exclusively solitary Hymenoptera phylogenetically positioned between ant and bee, both of which include eusocial species. We report the chemosensory-related gene repertoire sizes of three apoid wasps: Ampulex compressa, Cerceris arenaria, and Psenulus fuscipennis. We annotated genes encoding odorant (ORs), gustatory, and ionotropic receptors and chemosensory soluble proteins and odorant-binding proteins in transcriptomes of chemosensory tissues of the above three species and in early draft genomes of two species, A. compressa and C. arenaria. Our analyses revealed that apoid wasps possess larger OR repertoires than any bee lineage, that the last common ancestor of Apoidea possessed a considerably larger OR repertoire (∼160) than previously estimated (73), and that the expansion of OR genes in eusocial bees was less extensive than previously assumed. Intriguingly, the evolution of pollen-collecting behavior in the stem lineage of bees was associated with a notable loss of OR gene diversity. Thus, our results support the view that herbivorous Hymenoptera tend to possess smaller OR repertoires than carnivorous, parasitoid, or kleptoparasitic species.