WaspAtlas: a Nasonia vitripennis gene database and analysis platform

Decoding the genetic and chemical basis of sexual attractiveness in parasitic wasps

Attracting and securing potential mating partners is of fundamental importance for reproduction. Therefore, signaling sexual attractiveness is expected to be tightly coordinated in communication systems synchronizing senders and receivers. Chemical signaling has permeated through all taxa of life as the earliest and most widespread form of communication and is particularly prevalent in insects. However, it has been notoriously difficult to decipher how exactly information related to sexual signaling is encoded in complex chemical profiles. Similarly, our knowledge of the genetic basis of sexual signaling is very limited and usually restricted to a few case studies with comparably simple pheromonal communication mechanisms. The present study jointly addresses these two knowledge gaps by characterizing two fatty acid synthase genes that most likely evolved by tandem gene duplication and that simultaneously impact sexual attractiveness and complex chemical surface profiles in parasitic wasps. Gene knockdown in female wasps dramatically reduces their sexual attractiveness coinciding with a drastic decrease in male courtship and copulation behavior. Concordantly, we found a striking shift of methyl-branching patterns in the female surface pheromonal compounds, which we subsequently demonstrate to be the main cause for the greatly reduced male mating response. Intriguingly, this suggests a potential coding mechanism for sexual attractiveness mediated by specific methyl-branching patterns in complex cuticular hydrocarbon (CHC) profiles. So far, the genetic underpinnings of methyl-branched CHCs are not well understood despite their high potential for encoding information. Our study sheds light on how biologically relevant information can be encoded in complex chemical profiles and on the genetic basis of sexual attractiveness.

Effect of usinggreen fluorescent proteindouble-stranded RNA as non-target negative control inNasonia vitripennisRNA interference assays

RNA interference (RNAi) is a technique used in many insects to study gene function. However, prior research suggests possible off-target effects when usingGreen Fluorescent Protein(GFP) sequence as a non-target control. We used a transcriptomic approach to study the effect ofGFPRNAi (GFP-i) inNasonia vitripennis, a widely used parasitoid wasp model system. Our study identified 3.4% of total genes being differentially expressed in response toGFP-i. A subset of these genes appears involved in microtubule and sperm functions.In silicoanalysis identified 17 potential off-targets, of which only one was differentially expressed afterGFP-i. We suggest the primary cause for differential expression afterGFP-i is the non-specific activation of the RNAi machinery at the injection site, and a potentially disturbed spermatogenesis. Still, we advise that any RNAi study involving the genes deregulated in this study, exercises caution in drawing conclusions and uses a different non-target control.

Genomics of sex allocation in the parasitoid wasp Nasonia vitripennis

BACKGROUND

Whilst adaptive facultative sex allocation has been widely studied at the phenotypic level across a broad range of organisms, we still know remarkably little about its genetic architecture. Here, we explore the genome-wide basis of sex ratio variation in the parasitoid wasp Nasonia vitripennis, perhaps the best studied organism in terms of sex allocation, and well known for its response to local mate competition.

RESULTS

We performed a genome-wide association study (GWAS) for single foundress sex ratios using iso-female lines derived from the recently developed outbred N. vitripennis laboratory strain HVRx. The iso-female lines capture a sample of the genetic variation in HVRx and we present them as the first iteration of the Nasonia vitripennis Genome Reference Panel (NVGRP 1.0). This panel provides an assessment of the standing genetic variation for sex ratio in the study population. Using the NVGRP, we discovered a cluster of 18 linked SNPs, encompassing 9 annotated loci associated with sex ratio variation. Furthermore, we found evidence that sex ratio has a shared genetic basis with clutch size on three different chromosomes.

CONCLUSIONS

Our approach provides a thorough description of the quantitative genetic basis of sex ratio variation in Nasonia at the genome level and reveals a number of inter-related candidate loci underlying sex allocation regulation.

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.

Insect genomes: progress and challenges

In the wake of constant improvements in sequencing technologies, numerous insect genomes have been sequenced. Currently, 1219 insect genome-sequencing projects have been registered with the National Center for Biotechnology Information, including 401 that have genome assemblies and 155 with an official gene set of annotated protein-coding genes. Comparative genomics analysis showed that the expansion or contraction of gene families was associated with well-studied physiological traits such as immune system, metabolic detoxification, parasitism and polyphagy in insects. Here, we summarize the progress of insect genome sequencing, with an emphasis on how this impacts research on pest control. We begin with a brief introduction to the basic concepts of genome assembly, annotation and metrics for evaluating the quality of draft assemblies. We then provide an overview of genome information for numerous insect species, highlighting examples from prominent model organisms, agricultural pests and disease vectors. We also introduce the major insect genome databases. The increasing availability of insect genomic resources is beneficial for developing alternative pest control methods. However, many opportunities remain for developing data-mining tools that make maximal use of the available insect genome resources. Although rapid progress has been achieved, many challenges remain in the field of insect genomics.

Regulatory and sequence evolution in response to selection for improved associative learning ability in Nasonia vitripennis

Background

Selection acts on the phenotype, yet only the genotype is inherited. While both the phenotypic and genotypic response to short-term selection can be measured, the link between these is a major unsolved problem in evolutionary biology, in particular for complex behavioural phenotypes.

Results

Here we characterize the genomic and the transcriptomic basis of associative learning ability in the parasitic wasp Nasonia vitripennis and use gene network analysis to link the two. We artificially selected for improved associative learning ability in four independent pairs of lines and identified signatures of selection across the genome. Allele frequency diverged consistently between the selected and control lines in 118 single nucleotide polymorphisms (SNPs), clustering in 51 distinct genomic regions containing 128 genes. The majority of SNPs were found in regulatory regions, suggesting a potential role for gene expression evolution. We therefore sequenced the transcriptomes of selected and control lines and identified 36 consistently differentially expressed transcripts with large changes in expression. None of the differentially expressed genes also showed sequence divergence as a result of selection. Instead, gene network analysis showed many of the genes with consistent allele frequency differences and all of the differentially expressed genes to cluster in a single co-expression network. At a functional level, both genomic and transcriptomic analyses implicated members of gene networks known to be involved in neural plasticity and cognitive processes.

Conclusions

Taken together, our results reveal how specific cognitive abilities can readily respond to selection via a complex interplay between regulatory and sequence evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5310-9) contains supplementary material, which is available to authorized users.

Decay of Sexual Trait Genes in an Asexual Parasitoid Wasp

Trait loss is a widespread phenomenon with pervasive consequences for a species’ evolutionary potential. The genetic changes underlying trait loss have only been clarified in a small number of cases. None of these studies can identify whether the loss of the trait under study was a result of neutral mutation accumulation or negative selection. This distinction is relatively clear-cut in the loss of sexual traits in asexual organisms. Male-specific sexual traits are not expressed and can only decay through neutral mutations, whereas female-specific traits are expressed and subject to negative selection. We present the genome of an asexual parasitoid wasp and compare it to that of a sexual lineage of the same species. We identify a short-list of 16 genes for which the asexual lineage carries deleterious SNP or indel variants, whereas the sexual lineage does not. Using tissue-specific expression data from other insects, we show that fifteen of these are expressed in male-specific reproductive tissues. Only one deleterious variant was found that is expressed in the female-specific spermathecae, a trait that is heavily degraded and thought to be under negative selection in L. clavipes. Although the phenotypic decay of male-specific sexual traits in asexuals is generally slow compared with the decay of female-specific sexual traits, we show that male-specific traits do indeed accumulate deleterious mutations as expected by theory. Our results provide an excellent starting point for detailed study of the genomics of neutral and selected trait decay.

Maternal provision of transformer-2 is required for female development and embryo viability in the wasp Nasonia vitripennis

In insect sex determination a primary signal starts the genetic sex determination cascade that, in most insect orders, is subsequently transduced down the cascade by a transformer (tra) ortholog. Only a female-specifically spliced tra mRNA yields a functional TRA-protein that forms a complex with TRA2, encoded by a transformer-2 (tra2) ortholog, to act as a sex specific splicing regulator of the downstream transcription factors doublesex (dsx) and fruitless (fru). Here, we identify the tra2 ortholog of the haplodiploid parasitoid wasp N. vitripennis (Nv-tra2) and confirm its function in N. vitripennis sex determination. Knock down of Nv-tra2 by parental RNA interference (pRNAi) results in complete sex reversal of diploid offspring from female to male, indicating the requirement of Nv-tra2 for female sex determination. As Nv-tra2 pRNAi leads to frequent lethality in early developmental stages, maternal provision of Nv-tra2 transcripts is apparently also required for another, non-sex determining function during embryogenesis. In addition, lethality following Nv-tra2 pRNAi appears more pronounced in diploid than in haploid offspring. This diploid lethal effect was also observed following Nv-tra pRNAi, which served as a positive control in our experiments. As diploid embryos from fertilized eggs have a paternal chromosome set in addition to the maternal one, this suggests that either the presence of this paternal chromosome set or the dosage effect resulting from the diploid state is incompatible with the induced male development in N. vitripennis caused by either Nv-tra2 or Nv-tra pRNAi. The role of Nv-tra2 in activating the female sex determination pathway yields more insight into the sex determination mechanism of Nasonia.

OGS2: genome re-annotation of the jewel wasp Nasonia vitripennis

BACKGROUND

Nasonia vitripennis is an emerging insect model system with haplodiploid genetics. It holds a key position within the insect phylogeny for comparative, evolutionary and behavioral genetic studies. The draft genomes for N. vitripennis and two sibling species were published in 2010, yet a considerable amount of transcriptiome data have since been produced thereby enabling improvements to the original (OGS1.2) annotated gene set. We describe and apply the EvidentialGene method used to produce an updated gene set (OGS2). We also carry out comparative analyses showcasing the usefulness of the revised annotated gene set.

RESULTS

The revised annotation (OGS2) now consists of 24,388 genes with supporting evidence, compared to 18,850 for OGS1.2. Improvements include the nearly complete annotation of untranslated regions (UTR) for 97 % of the genes compared to 28 % of genes for OGS1.2. The fraction of RNA-Seq validated introns also grow from 85 to 98 % in this latest gene set. The EST and RNA-Seq expression data provide support for several non-protein coding loci and 7712 alternative transcripts for 4146 genes. Notably, we report 180 alternative transcripts for the gene lola. Nasonia now has among the most complete insect gene set; only 27 conserved single copy orthologs in arthropods are missing from OGS2. Its genome also contains 2.1-fold more duplicated genes and 1.4-fold more single copy genes than the Drosophila melanogaster genome. The Nasonia gene count is larger than those of other sequenced hymenopteran species, owing both to improvements in the genome annotation and to unique genes in the wasp lineage. We identify 1008 genes and 171 gene families that deviate significantly from other hymenopterans in their rates of protein evolution and duplication history, respectively. We also provide an analysis of alternative splicing that reveals that genes with no annotated isoforms are characterized by shorter transcripts, fewer introns, faster protein evolution and higher probabilities of duplication than genes having alternative transcripts.

CONCLUSIONS

Genome-wide expression data greatly improves the annotation of the N. vitripennis genome, by increasing the gene count, reducing the number of missing genes and providing more comprehensive data on splicing and gene structure. The improved gene set identifies lineage-specific genomic features tied to Nasonia's biology, as well as numerous novel genes. OGS2 and its associated search tools are available at http://arthropods.eugenes.org/EvidentialGene/nasonia/ , www.hymenopteragenome.org/nasonia/ and waspAtlas: www.tinyURL.com/waspAtlas . The EvidentialGene pipeline is available at https://sourceforge.net/projects/evidentialgene/ .

An optimized approach to germ-free rearing in the jewel wasp Nasonia

Development of a Nasonia in vitrogerm-free rearing system in 2012 enabled investigation of Nasonia-microbiota interactions and real-time visualization of parasitoid metamorphosis. However, the use of antibiotics, bleach, and fetal bovine serum introduced artifacts relative to conventional rearing of Nasonia. Here, we optimize the germ-free rearing procedure by using filter sterilization in lieu of antibiotics and by removing residual bleach and fetal bovine serum. Comparison of these methods reveals no influence on larval survival or growth, and a 52% improvement in adult production. Additionally, adult males produced in the new germ-free system are similar in size to conventionally reared males. Experimental implications of these changes are discussed.

Global analysis of dorsoventral patterning in the wasp Nasonia reveals extensive incorporation of novelty in a regulatory network

Background

Gene regulatory networks (GRNs) underlie developmental patterning and morphogenetic processes, and changes in the interactions within the underlying GRNs are a major driver of evolutionary processes. In order to make meaningful comparisons that can provide significant insights into the evolution of regulatory networks, homologous networks from multiple taxa must be deeply characterized. One of the most thoroughly characterized GRNs is the dorsoventral (DV) patterning system of the Drosophila melanogaster embryo. We have developed the wasp Nasonia as a comparative DV patterning model because it has shown the convergent evolution of a mode of early embryonic patterning very similar to that of the fly, and it is of interest to know whether the similarity at the gross level also extends to the molecular level.

Results

We used RNAi to dorsalize and ventralize Nasonia embryos, RNAseq to quantify transcriptome-wide expression levels, and differential expression analysis to identify genes whose expression levels change in either RNAi case. This led to the identification of >100 genes differentially expressed and regulated along the DV axis. Only a handful of these genes are shared DV components in both fly and wasp. Many of those unique to Nasonia are cytoskeletal and adhesion molecules, which may be related to the divergent cell and tissue behavior observed at gastrulation. In addition, many transcription factors and signaling components are only DV regulated in Nasonia, likely reflecting the divergent upstream patterning mechanisms involved in producing the conserved pattern of cell fates observed at gastrulation. Finally, several genes that lack Drosophila orthologs show robust and distinct expression patterns. These include genes with vertebrate homologs that have been lost in the fly lineage, genes that are found only among Hymenoptera, and several genes that entered the Nasonia genome through lateral transfer from endosymbiotic bacteria.

Conclusions

Altogether, our results provide insights into how GRNs respond to new functional demands and how they can incorporate novel components.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0285-y) contains supplementary material, which is available to authorized users.

Emerging developmental genetic model systems in holometabolous insects

The number of insect species that are amenable to functional genetic studies is growing rapidly and provides many new research opportunities in developmental and evolutionary biology. The holometabolous insects represent a disproportionate percentage of animal diversity and are thus well positioned to provide model species for a wide variety of developmental processes. Here we discuss emerging holometabolous models, and review some recent breakthroughs. For example, flies and midges were found to use structurally unrelated long-range pattern organizers, butterflies and moths revealed extensive pattern formation during oogenesis, new imaging possibilities in the flour beetle Tribolium castaneum showed how embryos break free of their extraembryonic membranes, and the complex genetics governing interspecies difference in head shape were revealed in Nasonia wasps.

DNA methylation changes induced by long and short photoperiods in Nasonia

Many organisms monitor the annual change in day length and use this information for the timing of their seasonal response. However, the molecular mechanisms underlying photoperiodic timing are largely unknown. The wasp Nasonia vitripennis is an emerging model organism that exhibits a strong photoperiodic response: Short autumnal days experienced by females lead to the induction of developmental arrest (diapause) in their progeny, allowing winter survival of the larvae. How female Nasonia control the developmental trajectory of their offspring is unclear. Here, we took advantage of the recent discovery that DNA methylation is pervasive in Nasonia and tested its role in photoperiodism. We used reduced representation bisulfite sequencing (RRBS) to profile DNA methylation in adult female wasps subjected to different photoperiods and identified substantial differential methylation at the single base level. We also show that knocking down DNA methyltransferase 1a (Dnmt1a), Dnmt3, or blocking DNA methylation pharmacologically, largely disrupts the photoperiodic diapause response of the wasps. To our knowledge, this is the first example for a role of DNA methylation in insect photoperiodic timing.