Comparative genomics of the miniature wasp and pest control agent Trichogramma pretiosum

Knockdown of vitellogenin receptor based on minute insect RNA interference methods affects the initial mature egg load in the pest natural enemy Trichogramma dendrolimi

Vitellogenin receptor (VgR) plays a crucial role in oogenesis by mediating endocytosis of vitellogenin and a portion of the yolk proteins in many insect species. However, the function of VgR in minute parasitoid wasps is largely unknown. Here, we applied Trichogramma dendrolimi, a minute egg parasitoid, as a study model to investigate the function of VgR in parasitoids. We developed RNA interference (RNAi) methods based on microinjection of prepupae in T. dendrolimi. RNAi employs nanomaterial branched amphipathic peptide capsules (BAPC) as a carrier for double-stranded RNA (dsRNA), significantly enhancing delivery efficiency. Also, artificial hosts without medium were used to culture the injected prepupae in vitro. Utilizing these methods, we found that ovarian growth was disrupted after knockdown of TdVgR, as manifested by the suppressed development of the ovariole and the inhibition of nurse cell internalization by oocytes. Also, the initial mature egg load in the ovary was significantly reduced. Notably, the parasitic capacity of the female adult with ovarian dysplasia was significantly decreased, possibly resulting from the low availability of mature eggs. Moreover, ovarian dysplasia in T. dendrolimi caused by VgR deficiency are conserved despite feeding on different hosts. The results confirmed a critical role of TdVgR in the reproductive ability of T. dendrolimi and provided a reference for gene functional studies in minute insects.

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.

Posterior concentration of Wolbachia during the early embryogenesis of the host dynamically shapes the tissue tropism of Wolbachia in host Trichogramma wasps

Introduction

The bacterial endosymbiont, Wolbachia spp. induce thelytokous parthenogenesis in certain parasitoid wasps, such as the egg parasitoid wasps Trichogramma spp. To complete the cycle of vertical transmission, Wolbachia displays efficient transovarial transmission by targeting the reproductive tissues and often exhibits strong tissue-specific tropism in their host.

Method

The present study aimed to describe the basic Wolbachia distribution patterns that occur during the development of Wolbachia-infected, thelytokous Trichogramma dendrolimi, and T. pretiosum. We used fluorescence in situ hybridization (FISH) to investigate Wolbachia signal dynamics during early embryogenesis (from 30 to 120 min). Wolbachia titers and distributions from the embryo to adult stages of Trichogramma after early embryogenesis were detected by absolute quantitative polymerase chain reaction (AQ-PCR) and FISH. The symmetry ratios (SR) of the Wolbachia signals were calculated using the SR odds ratios in the anterior and posterior parts of the host. The SR was determined to describe Wolbachia tropism during early embryogenesis and various developmental stages of Trichogramma.

Results

Wolbachia was concentrated in the posterior part of the embryo during early embryogenesis and the various developmental stages of both T. dendrolimi and T. pretiosum. Wolbachia density increased with the number of nuclei and the initial mitotic division frequency during early embryogenesis. The total Wolbachia titer increased with postembryogenesis development in both T. dendrolimi and T. pretiosum. However, the Wolbachia densities relative to body size were significantly lower at the adult and pupal stages than they were at the embryonic stage.

Discussion

The present work revealed that posterior Wolbachia concentration during early host embryogenesis determined Wolbachia localization in adult wasps. By this mechanism, Wolbachia exhibits efficient vertical transmission across generations by depositing only female Wolbachia-infected offspring. The results of this study describe the dynamics of Wolbachia during the development of their Trichogramma host. The findings of this investigation helped clarify Wolbachia tropism in Trichogramma wasps.

Parthenogenesis in dipterans: a genetic perspective

Parthenogenesis has been documented in almost every phylum of animals, and yet this phenomenon is largely understudied. It has particular importance in dipterans since some parthenogenetic species are also disease vectors and agricultural pests. Here, we present a catalogue of parthenogenetic dipterans, although it is likely that many more remain to be identified, and we discuss how their developmental biology and interactions with diverse environments may be linked to different types of parthenogenetic reproduction. We discuss how the advances in genetics and genomics have identified chromosomal loci associated with parthenogenesis. In particular, a polygenic cause of facultative parthenogenesis has been uncovered in Drosophila mercatorum, allowing the corresponding genetic variants to be tested for their ability to promote parthenogenesis in another species, Drosophila melanogaster. This study probably identifies just one of many routes that could be followed in the evolution of parthenogenesis. We attempt to account for why the phenomenon has evolved so many times in the dipteran order and why facultative parthenogenesis appears particularly prevalent. We also discuss the significance of coarse genomic changes, including non-disjunction, aneuploidy, and polyploidy and how, together with changes to specific genes, these might relate to both facultative and obligate parthenogenesis in dipterans and other parthenogenetic animals.

Finding an egg in a haystack: variation in chemical cue use by egg parasitoids of herbivorous insects

Egg parasitoids of herbivorous insects use an interplay of short- and long-range chemical cues emitted by hosts and host plants to find eggs to parasitize. Volatile compounds that attract egg parasitoids can be identified via behavioral assays and used to manipulate parasitoid behavior in the field for biological control of herbivorous pests. However, how and when a particular cue will be used varies over the life of an individual, as well as at and below species level. Future research should expand taxonomic coverage to explore variation in chemical cue use in more natural, dynamic settings. More nuanced understanding of the variability of egg parasitoid host-finding strategies will aid in disentangling the underlying genetics and further enhancing biological control.

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.

Comparative studies of genomic and epigenetic factors influencing transcriptional variation in two insect species

Different genes show different levels of expression variability. For example, highly expressed genes tend to exhibit less expression variability. Genes whose promoters have TATA box and initiator motifs tend to have increased expression variability. On the other hand, DNA methylation of transcriptional units, or gene body DNA methylation, is associated with reduced gene expression variability in many species. Interestingly, some insect lineages, most notably Diptera including the canonical model insect Drosophila melanogaster, have lost DNA methylation. Therefore, it is of interest to determine whether genomic features similarly influence gene expression variability in lineages with and without DNA methylation. We analyzed recently generated large-scale data sets in D. melanogaster and honey bee (Apis mellifera) to investigate these questions. Our analysis shows that increased gene expression levels are consistently associated with reduced expression variability in both species, while the presence of TATA box is consistently associated with increased gene expression variability. In contrast, initiator motifs and gene lengths have weak effects limited to some data sets. Importantly, we show that a sequence characteristics indicative of gene body DNA methylation is strongly and negatively associate with gene expression variability in honey bees, while it shows no such association in D. melanogaster. These results suggest the evolutionary loss of DNA methylation in some insect lineages has reshaped the molecular mechanisms concerning the regulation of gene expression variability.

Laboratory evaluation of the compatibility of Beauveria bassiana with the egg parasitoid Trichogramma dendrolimi (Hymenoptera: Trichogrammatidae) for joint application against the oriental fruit moth Grapholita molesta (Lepidoptera: Tortricidae)

BACKGROUND

The entomopathogenic fungus Beauveria bassiana and the egg parasitoid Trichogramma dendrolimi can both contribute to biological control of the global fruit borer Grapholita molesta. To derive insights for optimizing their combined application in augmentation programs, we assayed fungal pathogenicity to both G. molesta (eggs, larvae and adults) and adult parasitoids, and assessed wasp acquisition and transmission of fungal spores following their emergence from B. bassiana-treated host eggs.

RESULTS

Concentrations of 10⁸ conidia mL^-1 of B. bassiana and higher caused virtually complete mortality of G. molesta larvae, but less than 30% mortality of eggs, and female moths surviving treatments had reduced fecundity. More than 80% of adult T. dendrolimi survived B. bassiana treatments of 10⁷ conidia mL^-1 , and more than 60% survived 10⁹ conidia mL^-1 , with surviving females achieving reproductive success equivalent to controls. Parasitoid females preferred healthy eggs over B. bassiana-infected ones in choice tests; wasps emerging from treated eggs carried about 10⁴ conidia each and contributed an additional 11% host mortality in trials combining both agents.

CONCLUSION

Our results indicate that combined applications of B. bassiana and T. dendrolimi can have complimentary impacts on G. molesta, the wasps compensating for low fungal pathogenicity to eggs and their progeny potentially aiding in subsequent disease transmission. Although foraging wasps tended to avoid infected eggs, negative interactions between the two agents might be further mitigated by timing B. bassiana applications from late instar larva to early moth stage, and T. dendrolimi releases several days later to coincide with peak oviposition periods. © 2022 Society of Chemical Industry.

Whole-genome sequencing analysis and protocol for RNA interference of the endoparasitoid wasp Asobara japonica

Asobara japonica is an endoparasitic wasp that parasitizes Drosophila flies. It synthesizes various toxic components in the venom gland and injects them into host larvae during oviposition. To identify and characterize these toxic components for enabling parasitism, we performed the whole-genome sequencing (WGS) and devised a protocol for RNA interference (RNAi) with A. japonica. Because it has a parthenogenetic lineage due to Wolbachia infection, we generated a clonal strain from a single wasp to obtain highly homogenous genomic DNA. The WGS analysis revealed that the estimated genome size was 322 Mb with a heterozygosity of 0.132%. We also performed RNA-seq analyses for gene annotation. Based on the qualified WGS platform, we cloned ebony-Aj, which encodes the enzyme N-β-alanyl dopamine synthetase, which is involved in melanin production. The microinjection of double-stranded RNA (dsRNA) targeting ebony-Aj led to body colour changes in adult wasps, phenocopying ebony-Dm mutants. Furthermore, we identified putative venom genes as a target of RNAi, confirming that dsRNA injection-based RNAi specifically suppressed the expression of the target gene in wasp adults. Taken together, our results provide a powerful genetic toolkit for studying the molecular mechanisms of parasitism.

Mitochondrial Genomes of Two Asexual Trichogramma (Hymenoptera: Trichogrammatidae) Strains and Comparison with Their Sexual Relatives

Simple Summary

Sexual reproduction is dominant in animals, while asexual lineages are rare and evolutionarily short-lived. However, sexual reproduction has substantial costs, such as male production, inputs to courtship and mating, increased risk of predator exposure, and sexually transmitted diseases. A large body of theories has been proposed to explain the paradox of sex. One favored explanation is that asexuals are more likely to accumulate a greater number of deleterious mutations, known as Muller’s ratchet. Trichogramma is a genus of egg parasitoid wasps and is widely used as a biological control agent for agricultural and forest pests. With asexual lineages in at least 16 species, Trichogramma provides an excellent model to investigate the causes and consequences of asexual reproduction. In this study, we sequenced and assembled two asexual Trichogramma mitogenomes, representing two divergent origins of asexual reproduction. The asexual T. pretiosum is induced by the endosymbiont Wolbachia, while T. cacoeciae presumably originates from interspecific hybridization. To test Muller’s ratchet hypothesis, we compared these two asexual mitogenomes with their sexual relatives and found no association between asexual reproduction and mutation accumulation. This study provides a basis for further investigation into mitochondrial evolution and asexual reproduction in Trichogramma.

Abstract

Despite its substantial costs, sexual reproduction dominates in animals. One popular explanation for the paradox of sex is that asexual reproduction is more likely to accumulate deleterious mutations than sexual reproduction. To test this hypothesis, we compared the mitogenomes of two asexual wasp strains, Trichogramma cacoeciae and T. pretiosum, to their sexual relatives. These two asexual strains represent two different transition mechanisms in Trichogramma from sexual to asexual reproduction. Asexual T. pretiosum is induced by Wolbachia, while T. cacoeciae presumably originated from interspecific hybridization. We sequenced and assembled complete mitochondrial genomes of asexual T. cacoeciae and T. pretiosum. Compared to four sexual relatives, we found no evidence of higher mutation accumulation in asexual Trichogramma mitogenomes than in their sexual relatives. We also did not detect any relaxed selection in asexual Trichogramma mitogenomes. In contrast, the intensified selection was detected in Nad1 and Nad4 of the asexual T. pretiosum mitogenome, suggesting more purifying selection. In summary, no higher mitochondrial mutation accumulation was detected in these two asexual Trichogramma strains. This study provides a basis for further investigating mitochondrial evolution and asexual reproduction in Trichogramma.

Crystal Structures of Wolbachia CidA and CidB Reveal Determinants of Bacteria-induced Cytoplasmic Incompatibility and Rescue

Cytoplasmic incompatibility (CI) results when Wolbachia bacteria-infected male insects mate with uninfected females, leading to embryonic lethality. “Rescue” of viability occurs if the female harbors the same Wolbachia strain. CI is caused by linked pairs of Wolbachia genes called CI factors (CifA and CifB). The co-evolution of CifA-CifB pairs may account in part for the incompatibility patterns documented in insects infected with different Wolbachia strains, but the molecular mechanisms remain elusive. Here, we use X-ray crystallography and AlphaFold to analyze the CI factors from Wolbachia strain wMel called CidA^wMel and CidB^wMel. Substituting CidA^wMel interface residues with those from CidA^wPip (from strain wPip) enables the mutant protein to bind CidB^wPip and rescue CidB^wPip-induced yeast growth defects, supporting the importance of CifA-CifB interaction in CI rescue. Sequence divergence in CidA^wPip and CidB^wPip proteins affects their pairwise interactions, which may help explain the complex incompatibility patterns of mosquitoes infected with different wPip strains.

Wolbachia induced cytoplasmic incompatibility (CI) is caused by linked pairs of genes named cifA and cifB. Here, authors show that the residues at interfaces of the CidA-CidB complex is crucial for their binding and contribute to the diversity of CI.

Comparative Genomics Sheds Light on the Convergent Evolution of Miniaturized Wasps

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

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.

Olfactory Sensilla and Olfactory Genes in the Parasitoid Wasp Trichogramma pretiosum Riley (Hymenoptera: Trichogrammatidae)

Simple Summary

Parasitic wasps are the major natural enemies of many organisms, and therefore they are broadly used in the biological control of numerous agricultural and horticultural pests. For example, Trichogramma pretiosum Riley (Hymenoptera: Trichogrammatidae) is a tiny natural egg parasitoid of various agricultural pest insects, including Plutella xylostella, Helicoverpa armigera, Spodoptera frugiperda and Ectomyelois ceratoniae. However, how T. pretiosum seek and localise host insect eggs is still not clear. The olfactory system is critical in guiding insect behaviours, including mating, feeding and oviposition, which play pivotal roles in the interactions between parasitoid wasps and their hosts. This project aimed to investigate T. pretiosum major olfactory tissue (antennae) and the olfactory genes, including odorant binding proteins (OBPs) and odorant receptors (ORs). T. pretiosum adult antennae were examined under scanning electron microscopy, and four types of olfactory sensilla were observed. Using T. pretiosum genome, 22 OBPs and 105 ORs were identified, which were further compared with olfactory genes of other Hymenoptera insect species. The expression patterns of OBPs between T. pretiosum male and female adults were examined to identify female- or male-specific OBPs. This study enriches our knowledge of T. pretiosum olfactory system and will help better use it in the integrated pest management (IPM) for many insect pest species.

Abstract

Trichogramma pretiosum Riley (Hymenoptera: Trichogrammatidae) is a tiny natural egg parasitoid of several agricultural pest insects, which has been widely used in the biological control for Plutella xylostella, Helicoverpa armigera, Spodoptera frugiperda and Ectomyelois ceratoniae. However, limited studies have been conducted on T. pretiosum olfactory system, which is critical in regulating insect behaviours. In this study, T. pretiosum adult antennae were investigated under ascanning electron microscopy (SEM). Four types of olfactory sensilla were observed, including chaetica sensilla (CS), trichoid sensilla (TS), faleate sensilla (FS) and placoid sensilla (PS). Using T. pretiosum genome, 22 putative odorant binding proteins (OBPs) and 105 odorant receptors (ORs) were identified, which were further compared with olfactory genes of Apis mellifera, Nasonia vitripennis and Diachasma alloeum. The expression patterns of OBPs between T. pretiosum male and female adults were examined by quantitative real time PCR (qRT-PCR) approaches. Three female-specific OBPs (TpreOBP19, TpreOBP15 and TpreOBP3) were identified, which may play crucial roles in T. pretiosum host-seeking and oviposition behaviours. This study enriches our knowledge of T. pretiosum olfactory genes and improves our understanding of its olfactory system.

Transitions to asexuality and evolution of gene expression in Artemia brine shrimp

While sexual reproduction is widespread among many taxa, asexual lineages have repeatedly evolved from sexual ancestors. Despite extensive research on the evolution of sex, it is still unclear whether this switch represents a major transition requiring major molecular reorganization, and how convergent the changes involved are. In this study, we investigated the phylogenetic relationship and patterns of gene expression of sexual and asexual lineages of Eurasian Artemia brine shrimp, to assess how gene expression patterns are affected by the transition to asexuality. We find only a few genes that are consistently associated with the evolution of asexuality, suggesting that this shift may not require an extensive overhauling of the meiotic machinery. While genes with sex-biased expression have high rates of expression divergence within Eurasian Artemia, neither female- nor male-biased genes appear to show unusual evolutionary patterns after sexuality is lost, contrary to theoretical expectations.

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.

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.

Lineage and Parent-of-Origin Effects in DNA Methylation of Honey Bees (Apis mellifera) Revealed by Reciprocal Crosses and Whole-Genome Bisulfite Sequencing

Parent-of-origin methylation arises when the methylation patterns of a particular allele are dependent on the parent it was inherited from. Previous work in honey bees has shown evidence of parent-of-origin-specific expression, yet the mechanisms regulating such pattern remain unknown in honey bees. In mammals and plants, DNA methylation is known to regulate parent-of-origin effects such as genomic imprinting. Here, we utilize genotyping of reciprocal European and Africanized honey bee crosses to study genome-wide allele-specific methylation patterns in sterile and reproductive individuals. Our data confirm the presence of allele-specific methylation in honey bees in lineage-specific contexts but also importantly, though to a lesser degree, parent-of-origin contexts. We show that the majority of allele-specific methylation occurs due to lineage rather than parent-of-origin factors, regardless of the reproductive state. Interestingly, genes affected by allele-specific DNA methylation often exhibit both lineage and parent-of-origin effects, indicating that they are particularly labile in terms of DNA methylation patterns. Additionally, we re-analyzed our previous study on parent-of-origin-specific expression in honey bees and found little association with parent-of-origin-specific methylation. These results indicate strong genetic background effects on allelic DNA methylation and suggest that although parent-of-origin effects are manifested in both DNA methylation and gene expression, they are not directly associated with each other.

A chromosome-level genome assembly of the parasitoid wasp, Cotesia glomerata (Hymenoptera: Braconidae)

Hymenopterans make up about 20% of all animal species, but most are poorly known and lack high-quality genomic resources. One group of important, yet under-studied hymenopterans, are parasitoid wasps in the family Braconidae. Among this under-studied group are braconid wasps in the genus Cotesia; a clade of ~1,000 species routinely used in studies of physiology, ecology, biological control, and genetics. However, our ability to understand these organisms has been hindered by a lack of genomic resources. We helped bridge this gap by generating a high-quality genome assembly for the parasitoid wasp, Cotesia glomerata (Braconidae; Microgastrinae). We generated this assembly using multiple sequencing technologies, including Oxford Nanopore, whole-genome shotgun sequencing, and 3-D chromatin contact information (Hi-C). Our assembly is one of the most contiguous, complete, and publicly available hymenopteran genomes, represented by 3,355 scaffolds with a scaffold N50 of ~28Mb and a BUSCO score of ~99%. Given the genome sizes found in closely related species, our genome assembly was ~50% larger than expected, which was apparently induced by runaway amplification of three types of repetitive elements: simple repeats, Long Terminal Repeats (LTRs), and Long Interspersed Nuclear Elements (LINEs). This assembly is another step forward for genomics across this hyper-diverse, yet understudied, order of insects. The assembled genomic data and metadata files are publicly available via Figshare (https://doi.org/10.6084/m9.figshare.13010549).

Jekyll or Hyde? The genome (and more) of Nesidiocoris tenuis, a zoophytophagous predatory bug that is both a biological control agent and a pest

Nesidiocoris tenuis (Reuter) is an efficient predatory biological control agent used throughout the Mediterranean Basin in tomato crops but regarded as a pest in northern European countries. From the family Miridae, it is an economically important insect yet very little is known in terms of genetic information and no genomic or transcriptomic studies have been published. Here, we use a linked-read sequencing strategy on a single female N. tenuis. From this, we assembled the 355 Mbp genome and delivered an ab initio, homology-based and evidence-based annotation. Along the way, the bacterial "contamination" was removed from the assembly. In addition, bacterial lateral gene transfer (LGT) candidates were detected in the N. tenuis genome. The complete gene set is composed of 24 688 genes; the associated proteins were compared to other hemipterans (Cimex lectularis, Halyomorpha halys and Acyrthosiphon pisum). We visualized the genome using various cytogenetic techniques, such as karyotyping, CGH and GISH, indicating a karyotype of 2n = 32. Additional analyses include the localization of 18S rDNA and unique satellite probes as well as pooled sequencing to assess nucleotide diversity and neutrality of the commercial population. This is one of the first mirid genomes to be released and the first of a mirid biological control agent.

Genomic Features of Parthenogenetic Animals

Evolution without sex is predicted to impact genomes in numerous ways. Case studies of individual parthenogenetic animals have reported peculiar genomic features that were suggested to be caused by their mode of reproduction, including high heterozygosity, a high abundance of horizontally acquired genes, a low transposable element load, or the presence of palindromes. We systematically characterized these genomic features in published genomes of 26 parthenogenetic animals representing at least 18 independent transitions to asexuality. Surprisingly, not a single feature was systematically replicated across a majority of these transitions, suggesting that previously reported patterns were lineage-specific rather than illustrating the general consequences of parthenogenesis. We found that only parthenogens of hybrid origin were characterized by high heterozygosity levels. Parthenogens that were not of hybrid origin appeared to be largely homozygous, independent of the cellular mechanism underlying parthenogenesis. Overall, despite the importance of recombination rate variation for the evolution of sexual animal genomes, the genome-wide absence of recombination does not appear to have had the dramatic effects which are expected from classical theoretical models. The reasons for this are probably a combination of lineage-specific patterns, the impact of the origin of parthenogenesis, and a survivorship bias of parthenogenetic lineages.

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.

Amino acid synthesis loss in parasitoid wasps and other hymenopterans

Insects utilize diverse food resources which can affect the evolution of their genomic repertoire, including leading to gene losses in different nutrient pathways. Here, we investigate gene loss in amino acid synthesis pathways, with special attention to hymenopterans and parasitoid wasps. Using comparative genomics, we find that synthesis capability for tryptophan, phenylalanine, tyrosine, and histidine was lost in holometabolous insects prior to hymenopteran divergence, while valine, leucine, and isoleucine were lost in the common ancestor of Hymenoptera. Subsequently, multiple loss events of lysine synthesis occurred independently in the Parasitoida and Aculeata. Experiments in the parasitoid Cotesia chilonis confirm that it has lost the ability to synthesize eight amino acids. Our findings provide insights into amino acid synthesis evolution, and specifically can be used to inform the design of parasitoid artificial diets for pest control.

Hybrid Genome Assembly and Evidence-Based Annotation of the Egg Parasitoid and Biological Control Agent Trichogramma brassicae

Trichogramma brassicae (Bezdenko) are egg parasitoids that are used throughout the world as biological control agents and in laboratories as model species. Despite this ubiquity, few genetic resources exist beyond COI, ITS2, and RAPD markers. Aided by a Wolbachia infection, a wild-caught strain from Germany was reared for low heterozygosity and sequenced in a hybrid de novo strategy, after which several assembling strategies were evaluated. The best assembly, derived from a DBG2OLC-based pipeline, yielded a genome of 235 Mbp made up of 1,572 contigs with an N50 of 556,663 bp. Following a rigorous ab initio-, homology-, and evidence-based annotation, 16,905 genes were annotated and functionally described. As an example of the utility of the genome, a simple ortholog cluster analysis was performed with sister species T. pretiosum, revealing over 6000 shared clusters and under 400 clusters unique to each species. The genome and transcriptome presented here provides an essential resource for comparative genomics of the commercially relevant genus Trichogramma, but also for research into molecular evolution, ecology, and breeding of T. brassicae.

Distinct epigenomic and transcriptomic modifications associated with Wolbachia-mediated asexuality

Wolbachia are maternally transmitted intracellular bacteria that induce a range of pathogenic and fitness-altering effects on insect and nematode hosts. In parasitoid wasps of the genus Trichogramma, Wolbachia infection induces asexual production of females, thus increasing transmission of Wolbachia. It has been hypothesized that Wolbachia infection accompanies a modification of the host epigenome. However, to date, data on genome-wide epigenomic changes associated with Wolbachia are limited, and are often confounded by background genetic differences. Here, we took sexually reproducing Trichogramma free of Wolbachia and introgressed their genome into a Wolbachia-infected cytoplasm, converting them to Wolbachia-mediated asexuality. Wolbachia was then cured from replicates of these introgressed lines, allowing us to examine the genome-wide effects of wasps newly converted to asexual reproduction while controlling for genetic background. We thus identified gene expression and DNA methylation changes associated with Wolbachia-infection. We found no overlaps between differentially expressed genes and differentially methylated genes, indicating that Wolbachia-infection associated DNA methylation change does not directly modulate levels of gene expression. Furthermore, genes affected by these mechanisms exhibit distinct evolutionary histories. Genes differentially methylated due to the infection tended to be evolutionarily conserved. In contrast, differentially expressed genes were significantly more likely to be unique to the Trichogramma lineage, suggesting host-specific transcriptomic responses to infection. Nevertheless, we identified several novel aspects of Wolbachia-associated DNA methylation changes. Differentially methylated genes included those involved in oocyte development and chromosome segregation. Interestingly, Wolbachia-infection was associated with higher levels of DNA methylation. Additionally, Wolbachia infection reduced overall variability in gene expression, even after accounting for the effect of DNA methylation. We also identified specific cases where alternative exon usage was associated with DNA methylation changes due to Wolbachia infection. These results begin to reveal distinct genes and molecular pathways subject to Wolbachia induced epigenetic modification and/or host responses to Wolbachia-infection.

Wolbachia is an extremely common endosymbiotic infection of arthropods and nematodes. One of the reasons why Wolbachia can so successfully infect diverse species is the bacterium’s ability to profoundly alter the reproductive behavior of its host. It has been proposed that Wolbachia may modify host’s epigenetic programs to alter its reproductive behavior. However, it has been difficult to study how epigenetic programs change with Wolbachia infection, due to the confounding effects of genetic backgrounds. Here, we studied host transcriptome and epigenome changes associated with Wolbachia infection in a homogenous genetic background, by carrying out an innovative introgression scheme. By doing so, we show, for the first time, high-resolution molecular consequences of intracellular infection and offer insights into epigenetic and transcriptomic regulation of invertebrates.

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

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

Gene content evolution in the arthropods

BACKGROUND

Arthropods comprise the largest and most diverse phylum on Earth and play vital roles in nearly every ecosystem. Their diversity stems in part from variations on a conserved body plan, resulting from and recorded in adaptive changes in the genome. Dissection of the genomic record of sequence change enables broad questions regarding genome evolution to be addressed, even across hyper-diverse taxa within arthropods.

RESULTS

Using 76 whole genome sequences representing 21 orders spanning more than 500 million years of arthropod evolution, we document changes in gene and protein domain content and provide temporal and phylogenetic context for interpreting these innovations. We identify many novel gene families that arose early in the evolution of arthropods and during the diversification of insects into modern orders. We reveal unexpected variation in patterns of DNA methylation across arthropods and examples of gene family and protein domain evolution coincident with the appearance of notable phenotypic and physiological adaptations such as flight, metamorphosis, sociality, and chemoperception.

CONCLUSIONS

These analyses demonstrate how large-scale comparative genomics can provide broad new insights into the genotype to phenotype map and generate testable hypotheses about the evolution of animal diversity.

Sex biased expression and co-expression networks in development, using the hymenopteran Nasonia vitripennis

Sexual dimorphism requires regulation of gene expression in developing organisms. These developmental differences are caused by differential expression of genes and isoforms. The effect of expressing a gene is also influenced by which other genes are simultaneously expressed (functional interactions). However, few studies have described how these processes change across development. We compare the dynamics of differential expression, isoform switching and functional interactions in the sexual development of the model parasitoid wasp Nasonia vitripennis, a system that permits genome wide analysis of sex bias from early embryos to adults. We find relatively little sex-bias in embryos and larvae at the gene level, but several sub-networks show sex-biased functional interactions in early developmental stages. These networks provide new candidates for hymenopteran sex determination, including histone modification. In contrast, sex-bias in pupae and adults is driven by the differential expression of genes. We observe sex-biased isoform switching consistently across development, but mostly in genes that are already differentially expressed. Finally, we discover that sex-biased networks are enriched by genes specific to the Nasonia clade, and that those genes possess the topological properties of key regulators. These findings suggest that regulators in sex-biased networks evolve more rapidly than regulators of other developmental networks.

A partial genome assembly of the miniature parasitoid wasp, Megaphragma amalphitanum

Body size reduction, also known as miniaturization, is an important evolutionary process that affects a number of physiological and phenotypic traits and helps animals conquer new ecological niches. However, this process is poorly understood at the molecular level. Here, we report genomic and transcriptomic features of arguably the smallest known insect-the parasitoid wasp, Megaphragma amalphitanum (Hymenoptera: Trichogrammatidae). In contrast to expectations, we find that the genome and transcriptome sizes of this parasitoid wasp are comparable to other members of the Chalcidoidea superfamily. Moreover, compared to other chalcid wasps the gene content of M. amalphitanum is remarkably conserved. Intriguingly, we observed significant changes in M. amalphitanum transposable element dynamics over time, in which an initial burst was followed by suppression of activity, possibly due to a recent reinforcement of the genome defense machinery. Overall, while the M. amalphitanum genomic data reveal certain features that may be linked to the unusual biological properties of this organism, miniaturization is not associated with a large decrease in genome complexity.

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

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

Chemical cues and genetic divergence in insects on plants: conceptual cross pollination between mutualistic and antagonistic systems

Cascading or reciprocal genetic diversification of herbivores, parasitoids, and pollinators can track chemotypic variation in host resources, and can lead to non-overlapping communities. Because plants simultaneously interact with both pollinators and herbivores, models investigating the genetic divergence of antagonistic herbivores and mutualistic pollinators should be merged in order to study how both processes interact using a common conceptual and methodological approach. We expect insects to mediate divergence in many systems, with outcomes depending on the level of pollinator or herbivore specialisation, and the relative selective pressures they impose. Applying approaches widely used to study insect pollinators, for example genomic tools and integration of behavioural, genetic and chemical data, to both pollinators and herbivores in the same system will facilitate our understanding of patterns of genetic divergence across multiple interacting species.

Sex loss in insects: causes of asexuality and consequences for genomes

Boasting a staggering diversity of reproductive strategies, insects provide attractive models for the comparative study of the causes and consequences of transitions to asexuality. We provide an overview of some contemporary studies of reproductive systems in insects and compile an initial database of asexual insect genome resources. Insect systems have already yielded some important insights into various mechanisms by which sex is lost, including genetic, endosymbiont-mediated, and hybridization. Studies of mutation and substitution after loss of sex provide the strongest empirical support for hypothesized effects of asexuality, whereas there is mixed evidence for ecological hypotheses such as increased parasite load and altered niche breadth in asexuals. Most hypotheses have been explored in a select few taxa (e.g. stick insects, aphids), such that much of the great taxonomic breadth of insects remain understudied. Given the variation in the proximate causes of asexuality in insects, we argue for expanding the taxonomic breadth of study systems. Despite some challenges for investigating sex in insects, the increasing cost-effectiveness of genomic sequencing makes data generation for closely-related asexual and sexual lineages increasingly feasible.

No limits: Breaking constraints in insect miniaturization

Small arthropods are not simply scaled-down versions of their larger closest relatives, as changes in morphology and functional characters are largely governed by scaling laws. These same scaling laws set strict limits to size change toward smaller sizes. The evolution of extreme miniaturized forms involves the breaking of these constraints, by means of design innovations that allow evolutionary change to evade the limits posed by scaling laws. Here we review several cases studies in insects and other arthropods that illustrate this evolutionary path. We examine morphologies commonly recurring in miniaturized forms but not exclusive to them, morphologies exclusive to miniaturized forms and novel functional solutions supported by unconventional morphologies. We also discuss miniaturization and its evolvability taking into consideration arthropod postembryonic development and modular body organization. The modification of features commonly supposed not to change appears as a recurring pattern in arthropod miniaturization.

Genomic Landscape of Methylation Islands in Hymenopteran Insects

Recent genome-wide DNA methylation analyses of insect genomes accentuate an intriguing contrast compared with those in mammals. In mammals, most CpGs are heavily methylated, with the exceptions of clusters of hypomethylated sites referred to as CpG islands. In contrast, DNA methylation in insects is localized to a small number of CpG sites. Here, we refer to clusters of methylated CpGs as “methylation islands (MIs),” and investigate their characteristics in seven hymenopteran insects with high-quality bisulfite sequencing data. Methylation islands were primarily located within gene bodies. They were significantly overrepresented in exon–intron boundaries, indicating their potential roles in splicing. Methylated CpGs within MIs exhibited stronger evolutionary conservation compared with those outside of MIs. Additionally, genes harboring MIs exhibited higher and more stable levels of gene expression compared with those that do not harbor MIs. The effects of MIs on evolutionary conservation and gene expression are independent and stronger than the effect of DNA methylation alone. These results indicate that MIs may be useful to gain additional insights into understanding the role of DNA methylation in gene expression and evolutionary conservation in invertebrate genomes.