Fine-scale mapping of the Nasonia genome to chromosomes using a high-density genotyping microarray

Genetic and genomic architecture of species-specific cuticular hydrocarbon variation in parasitoid wasps

Cuticular hydrocarbons (CHCs) serve two fundamental functions in insects: protection against desiccation and chemical signalling. How the interaction of genes shapes CHC profiles, which are essential for insect survival, adaptation and reproductive success, is still poorly understood. Here we investigate the genetic and genomic basis of CHC biosynthesis and variation in parasitoid wasps of the genus Nasonia. We mapped 91 quantitative trait loci (QTL) explaining the variation of a total of 43 CHCs in F₂ hybrid males from interspecific crosses between three Nasonia species. To identify candidate genes, we localized orthologues of CHC biosynthesis-related genes in the Nasonia genomes. We discovered multiple genomic regions where the location of QTL coincides with the location of CHC biosynthesis-related candidate genes. Most conspicuously, on a region close to the centromere of chromosome 1, multiple CHC biosynthesis-related candidate genes co-localize with several QTL explaining variation in methyl-branched alkanes. The genetic underpinnings behind this compound class are not well understood so far, despite their high potential for encoding chemical information as well as their prevalence in hymenopteran CHC profiles. Our study considerably extends our knowledge on the genetic architecture governing this important compound class, establishing a model for methyl-branched alkane genetics in the Hymenoptera in general.

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

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

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

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

Assemblies of the genomes of parasitic wasps using meta-assembly and scaffolding with genetic linkage

Safe, effective biological-control introductions against invasive pests depend on narrowly host-specific natural enemies with the ability to adapt to a changing environment. As part of a project on the genetic architectures of these traits, we assembled and annotated the genomes of two aphid parasitoids, Aphelinus atriplicis and Aphelinus certus. We report here several assemblies of A. atriplicis made with Illumina and PacBio data, which we combined into a meta-assembly. We scaffolded the meta-assembly with markers from a genetic map of hybrids between A. atriplicis and A. certus. We used this genetic-linkage scaffolded (GLS) assembly of A. atriplicis to scaffold a de novo assembly of A. certus. The de novo assemblies of A. atriplicis differed in contiguity, and the meta-assembly of these assemblies was more contiguous than the best de novo assembly. Scaffolding with genetic-linkage data allowed chromosomal-level assembly of the A. atriplicis genome and scaffolding a de novo assembly of A. certus with this GLS assembly, greatly increased the contiguity of the A. certus assembly to the point where it was also at the chromosomal-level. However, completeness of the A. atriplicis assembly, as measured by percent complete, single-copy BUSCO hymenopteran genes, varied little among de novo assemblies and was not increased by meta-assembly or genetic scaffolding. Furthermore, the greater contiguity of the meta-assembly and GLS assembly had little or no effect on the numbers of genes identified, the proportions with homologs or functional annotations. Increased contiguity of the A. certus assembly provided modest improvement in assembly completeness, as measured by percent complete, single-copy BUSCO hymenopteran genes. The total genic sequence increased, and while the number of genes declined, gene length increased, which together suggest greater accuracy of gene models. More contiguous assemblies provide uses other than gene annotation, for example, identifying the genes associated with quantitative trait loci and understanding of chromosomal rearrangements associated with speciation.

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.

A chromosome-level genome assembly of the parasitoid wasp Pteromalus puparum

Parasitoid wasps represent a large proportion of hymenopteran species. They have complex evolutionary histories and are important biocontrol agents. To advance parasitoid research, a combination of Illumina short-read, PacBio long-read and Hi-C scaffolding technologies was used to develop a high-quality chromosome-level genome assembly for Pteromalus puparum, which is an important pupal endoparasitoid of caterpillar pests. The chromosome-level assembly has aided in studies of venom and detoxification genes. The assembled genome size is 338 Mb with a contig N50 of 38.7 kb and a scaffold N50 of 1.16 Mb. Hi-C analysis assembled scaffolds onto five chromosomes and raised the scaffold N50 to 65.8 Mb, with more than 96% of assembled bases located on chromosomes. Gene annotation was assisted by RNA sequencing for the two sexes and four different life stages. Analysis detected 98% of the BUSCO (Benchmarking Universal Single-Copy Orthologs) gene set, supporting a high-quality assembly and annotation. In total, 40.1% (135.6 Mb) of the assembly is composed of repetitive sequences, and 14,946 protein-coding genes were identified. Although venom genes play important roles in parasitoid biology, their spatial distribution on chromosomes was poorly understood. Mapping has revealed venom gene tandem arrays for serine proteases, pancreatic lipase-related proteins and kynurenine-oxoglutarate transaminases, which have amplified in the P. puparum lineage after divergence from its common ancestor with Nasonia vitripennis. In addition, there is a large expansion of P450 genes in P. puparum. These examples illustrate how chromosome-level genome assembly can provide a valuable resource for molecular, evolutionary and biocontrol studies of parasitoid wasps.

Genome elimination mediated by gene expression from a selfish chromosome

Numerous plants and animals harbor selfish B chromosomes that "drive" or transmit themselves at super-Mendelian frequencies, despite long-term fitness costs to the organism. Currently, it is unknown how B chromosome drive is mediated, and whether B-gene expression plays a role. We used modern sequencing technologies to analyze the fine-scale sequence composition and expression of paternal sex ratio (PSR), a B chromosome in the jewel wasp Nasonia vitripennis. PSR causes female-to-male conversion by destroying the sperm's hereditary material in young embryos to drive. Using RNA interference, we demonstrate that testis-specific expression of a PSR-linked gene, named haploidizer, facilitates this genome elimination-and-sex conversion effect. haploidizer encodes a putative protein with a DNA binding domain, suggesting a functional link with the sperm-derived chromatin.

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.

Trait mapping in diverse arthropods by bulked segregant analysis

Bulked segregant analysis (BSA) is a cross-based method for genetic mapping in sexually reproducing organisms. The method's use of bulked (pooled) samples markedly reduces the genotyping effort associated with traditional linkage mapping studies. Further, it can be applied to species with life histories or physical attributes (as for micro-insects) that render genetic mapping with other methods impractical. Recent studies in both insects and mites have revealed that advanced BSA experimental designs can resolve causal loci to narrow genomic intervals, facilitating follow-up investigations. As high-quality genomes become more widely available, BSA methods are poised to become an increasingly important tool for the rapid mapping of both monogenic and polygenic traits in diverse arthropod species.

Expansion of the fatty acyl reductase gene family shaped pheromone communication in Hymenoptera

Fatty acyl reductases (FARs) are involved in the biosynthesis of fatty alcohols that serve a range of biological roles. Insects typically harbor numerous FAR gene family members. While some FARs are involved in pheromone biosynthesis, the biological significance of the large number of FARs in insect genomes remains unclear.

Using bumble bee (Bombini) FAR expression analysis and functional characterization, hymenopteran FAR gene tree reconstruction, and inspection of transposable elements (TEs) in the genomic environment of FARs, we uncovered a massive expansion of the FAR gene family in Hymenoptera, presumably facilitated by TEs. The expansion occurred in the common ancestor of bumble bees and stingless bees (Meliponini). We found that bumble bee FARs from the expanded FAR-A ortholog group contribute to the species-specific pheromone composition. Our results indicate that expansion and functional diversification of the FAR gene family played a key role in the evolution of pheromone communication in Hymenoptera.

The genomic architecture and molecular evolution of ant odorant receptors

The massive expansions of odorant receptor (OR) genes in ant genomes are notable examples of rapid genome evolution and adaptive gene duplication. However, the molecular mechanisms leading to gene family expansion remain poorly understood, partly because available ant genomes are fragmentary. Here, we present a highly contiguous, chromosome-level assembly of the clonal raider ant genome, revealing the largest known OR repertoire in an insect. While most ant ORs originate via local tandem duplication, we also observe several cases of dispersed duplication followed by tandem duplication in the most rapidly evolving OR clades. We found that areas of unusually high transposable element density (TE islands) were depauperate in ORs in the clonal raider ant, and found no evidence for retrotransposition of ORs. However, OR loci were enriched for transposons relative to the genome as a whole, potentially facilitating tandem duplication by unequal crossing over. We also found that ant OR genes are highly AT-rich compared to other genes. In contrast, in flies, OR genes are dispersed and largely isolated within the genome, and we find that fly ORs are not AT-rich. The genomic architecture and composition of ant ORs thus show convergence with the unrelated vertebrate ORs rather than the related fly ORs. This might be related to the greater gene numbers and/or potential similarities in gene regulation between ants and vertebrates as compared to flies.

The Maternal Effect Gene Wds Controls Wolbachia Titer in Nasonia

Maternal transmission of intracellular microbes is pivotal in establishing long-term, intimate symbioses. For germline microbes that exert negative reproductive effects on their hosts, selection can theoretically favor the spread of host genes that counteract the microbe's harmful effects. Here, we leverage a major difference in bacterial (Wolbachia pipientis) titers between closely related wasp species with forward genetic, transcriptomic, and cytological approaches to map two quantitative trait loci that suppress bacterial titers via a maternal effect. Fine mapping and knockdown experiments identify the gene Wolbachia density suppressor (Wds), which dominantly suppresses bacterial transmission from mother to embryo. Wds evolved by lineage-specific non-synonymous changes driven by positive selection. Collectively, our findings demonstrate that a genetically simple change arose by positive Darwinian selection in less than a million years to regulate maternally transmitted bacteria via a dominant, maternal effect gene.

Functional characterization of the transcriptional regulatory elements of three highly expressed constitutive genes in the jewel wasp, Nasonia vitripennis

The jewel wasp, Nasonia vitripennis Ashmead (Hymenoptera: Pteromalidae), is an easily reared parasitoid that is providing an ever increasingly malleable model for examining the biology and genetics of Hymenoptera. Utilizing genomic and transcriptome resources, 5' upstream transcriptional regulatory sequences (TREs) from three highly expressed genes were identified and cloned. Criteria for TRE selection included the presence of an adjacent gene 5' of the translation initiation site. One gene was methylated whereas the other two were nonmethylated. Each TRE, heat-shock protein 70 (hsp70), activator of 90 kDa hsp ATPase protein 1 (hsp90A), and lipid storage droplet surface-binding protein 1 (lsdp) was linked with enhanced green fluorescent protein (EGFP) coding sequence and cloned into both pDP9e somatic and piggyBac germline transformation vectors. EGFP expression patterns under control of each TRE were compared with patterns of DsRed fluorescence produced from the transformation vector cassette. Functional activity of each TRE was observed in cultured Spodoptera frugiperda 9 (Sf9) cells and Drosophila melanogaster as well as in N. vitripennis embryos demonstrating that all three sequences had functional transcriptional regulatory activity in three different insect orders. Identification and functional characterization of these three TREs will provide critical and necessary resources for further genetic analyses of N. vitripennis, Hymenoptera and other insects.

Generation of heritable germline mutations in the jewel wasp Nasonia vitripennis using CRISPR/Cas9

The revolutionary RNA-guided endonuclease CRISPR/Cas9 system has proven to be a powerful tool for gene editing in a plethora of organisms. Here, utilizing this system we developed an efficient protocol for the generation of heritable germline mutations in the parasitoid jewel wasp, Nasonia vitripennis, a rising insect model organism for the study of evolution, development of axis pattern formation, venom production, haplo-diploid sex determination, and host–symbiont interactions. To establish CRISPR-directed gene editing in N. vitripennis, we targeted a conserved eye pigmentation gene cinnabar, generating several independent heritable germline mutations in this gene. Briefly, to generate these mutants, we developed a protocol to efficiently collect N. vitripennis eggs from a parasitized flesh fly pupa, Sarcophaga bullata, inject these eggs with Cas9/guide RNA mixtures, and transfer injected eggs back into the host to continue development. We also describe a flow for screening mutants and establishing stable mutant strains through genetic crosses. Overall, our results demonstrate that the CRISPR/Cas9 system is a powerful tool for genome manipulation in N. vitripennis, with strong potential for expansion to target critical genes, thus allowing for the investigation of several important biological phenomena in this organism.

Disentangling a Holobiont - Recent Advances and Perspectives in Nasonia Wasps

The parasitoid wasp genus Nasonia (Hymenoptera: Chalcidoidea) is a well-established model organism for insect development, evolutionary genetics, speciation, and symbiosis. The host-microbiota assemblage which constitutes the Nasonia holobiont (a host together with all of its associated microbes) consists of viruses, two heritable bacterial symbionts and a bacterial community dominated in abundance by a few taxa in the gut. In the wild, all four Nasonia species are systematically infected with the obligate intracellular bacterium Wolbachia and can additionally be co-infected with Arsenophonus nasoniae. These two reproductive parasites have different transmission modes and host manipulations (cytoplasmic incompatibility vs. male-killing, respectively). Pioneering studies on Wolbachia in Nasonia demonstrated that closely related Nasonia species harbor multiple and mutually incompatible Wolbachia strains, resulting in strong symbiont-mediated reproductive barriers that evolved early in the speciation process. Moreover, research on host-symbiont interactions and speciation has recently broadened from its historical focus on heritable symbionts to the entire microbial community. In this context, each Nasonia species hosts a distinguishable community of gut bacteria that experiences a temporal succession during host development and members of this bacterial community cause strong hybrid lethality during larval development. In this review, we present the Nasonia species complex as a model system to experimentally investigate questions regarding: (i) the impact of different microbes, including (but not limited to) heritable endosymbionts, on the extended phenotype of the holobiont, (ii) the establishment and regulation of a species-specific microbiota, (iii) the role of the microbiota in speciation, and (iv) the resilience and adaptability of the microbiota in wild populations subjected to different environmental pressures. We discuss the potential for easy microbiota manipulations in Nasonia as a promising experimental approach to address these fundamental aspects.

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/ .

Dissection of the complex genetic basis of craniofacial anomalies using haploid genetics and interspecies hybrids in Nasonia wasps

The animal head is a complex structure where numerous sensory, structural and alimentary structures are concentrated and integrated, and its ontogeny requires precise and delicate interactions among genes, cells, and tissues. Thus, it is perhaps unsurprising that craniofacial abnormalities are among the most common birth defects in people, or that these defects have a complex genetic basis involving interactions among multiple loci. Developmental processes that depend on such epistatic interactions become exponentially more difficult to study in diploid organisms as the number of genes involved increases. Here, we present hybrid haploid males of the wasp species pair Nasonia vitripennis and Nasonia giraulti, which have distinct male head morphologies, as a genetic model of craniofacial development that possesses the genetic advantages of haploidy, along with many powerful genomic tools. Viable, fertile hybrids can be made between the species, and quantitative trail loci related to shape differences have been identified. In addition, a subset of hybrid males show head abnormalities, including clefting at the midline and asymmetries. Crucially, epistatic interactions among multiple loci underlie several developmental differences and defects observed in the F2 hybrid males. Furthermore, we demonstrate an introgression of a chromosomal region from N. giraulti into N. vitripennis that shows an abnormality in relative eye size, which maps to a region containing a major QTL for this trait. Therefore, the genetic sources of head morphology can, in principle, be identified by positional cloning. Thus, Nasonia is well positioned to be a uniquely powerful model invertebrate system with which to probe both development and complex genetics of craniofacial patterning and defects.

Allele-Specific Transcriptome and Methylome Analysis Reveals Stable Inheritance and Cis-Regulation of DNA Methylation in Nasonia

Gene expression divergence between closely related species could be attributed to both cis- and trans- DNA sequence changes during evolution, but it is unclear how the evolutionary dynamics of epigenetic marks are regulated. In eutherian mammals, biparental DNA methylation marks are erased and reset during gametogenesis, resulting in paternal or maternal imprints, which lead to genomic imprinting. Whether DNA methylation reprogramming exists in insects is not known. Wasps of the genus Nasonia are non-social parasitoids that are emerging as a model for studies of epigenetic processes in insects. In this study, we quantified allele-specific expression and methylation genome-wide in Nasonia vitripennis and Nasonia giraulti and their reciprocal F1 hybrids. No parent-of-origin effect in allelic expression was found for >8,000 covered genes, suggesting a lack of genomic imprinting in adult Nasonia. As we expected, both significant cis- and trans- effects are responsible for the expression divergence between N. vitripennis and N. giraulti. Surprisingly, all 178 differentially methylated genes are also differentially methylated between the two alleles in F1 hybrid offspring, recapitulating the parental methylation status with nearly 100% fidelity, indicating the presence of strong cis-elements driving the target of gene body methylation. In addition, we discovered that total and allele-specific expression are positively correlated with allele-specific methylation in a subset of the differentially methylated genes. The 100% cis-regulation in F1 hybrids suggests the methylation machinery is conserved and DNA methylation is targeted by cis features in Nasonia. The lack of genomic imprinting and parent-of-origin differentially methylated regions in Nasonia, together with the stable inheritance of methylation status between generations, suggests either a cis-regulatory motif for methylation at the DNA level or highly stable inheritance of an epigenetic signal in Nasonia.

RNA-sequencing and whole-genome bisulfite sequencing in the hybrid offspring of two Nasonia parasitoid wasp species revealed strong cis-regulation of methylation and allele-specific expression. No gene was found to display genomic imprinting.

The relationship between methylation of genomic DNA and expression of the genes that it encodes—and how this relationship changes during evolution—has been widely studied in mammals, but remains less well understood for insects. Here we analyze the expressed mRNA transcripts and genomic DNA methylation of the hybrid offspring of a pair of Nasonia parasitoid wasp species, producing a wealth of information about the regulation of gene expression. We find that variation in DNA sequence impacts expression on the same strand (called “cis-regulation”), and that cytosine methylation state is also associated in cis with the regulatory consequences of this base alteration. We show that these wasp species lack differential expression dependent on parent-of-origin (called “genomic imprinting”), and that in the hybrids the alleles retain the methylation status of the parental species in a strong cis-regulated fashion. Transcript abundances were also largely driven in a cis-regulated manner, consistent with a correlation between methylation status and expression levels. Despite the many differences between Nasonia and mammals in the impact of genomic DNA methylation, in both groups the use of methylated cytosine has been co-opted in ways that help tune gene expression.

The expanding genetic toolbox of the wasp Nasonia vitripennis and its relatives

The parasitoid wasp Nasonia represents a genus of four species that is emerging as a powerful genetic model system that has made and will continue to make important contributions to our understanding of evolutionary biology, development, ecology, and behavior. Particularly powerful are the haplodiploid genetics of the system, which allow some of the advantages of microbial genetics to be applied to a complex multicellular eukaryote. In addition, fertile, viable hybrids can be made among the four species in the genus. This makes Nasonia exceptionally well suited for evolutionary genetics approaches, especially when combined with its haploid genetics and tractability in the laboratory. These features are complemented by an expanding array of genomic, transcriptomic, and functional resources, the application of which has already made Nasonia an important model system in such emerging fields as evolutionary developmental biology and microbiomics. This article describes the genetic and genomic advantages of Nasonia wasps and the resources available for their genetic analysis.

Identification of Genes Uniquely Expressed in the Germ-Line Tissues of the Jewel Wasp Nasonia vitripennis

The jewel wasp Nasonia vitripennis is a rising model organism for the study of haplo-diploid reproduction characteristic of hymenopteran insects, which include all wasps, bees, and ants. We performed transcriptional profiling of the ovary, the female soma, and the male soma of N. vitripennis to complement a previously existing transcriptome of the wasp testis. These data were deposited into an open-access genome browser for visualization of transcripts relative to their gene models. We used these data to identify the assemblies of genes uniquely expressed in the germ-line tissues. We found that 156 protein-coding genes are expressed exclusively in the wasp testis compared with only 22 in the ovary. Of the testis-specific genes, eight are candidates for male-specific DNA packaging proteins known as protamines. We found very similar expression patterns of centrosome associated genes in the testis and ovary, arguing that de novo centrosome formation, a key process for development of unfertilized eggs into males, likely does not rely on large-scale transcriptional differences between these tissues. In contrast, a number of meiosis-related genes show a bias toward testis-specific expression, despite the lack of true meiosis in N. vitripennis males. These patterns may reflect an unexpected complexity of male gamete production in the haploid males of this organism. Broadly, these data add to the growing number of genomic and genetic tools available in N. vitripennis for addressing important biological questions in this rising insect model organism.

Genetic and epigenetic architecture of sex-biased expression in the jewel wasps Nasonia vitripennis and giraulti

There is extraordinary diversity in sexual dimorphism (SD) among animals, but little is known about its epigenetic basis. To study the epigenetic architecture of SD in a haplodiploid system, we performed RNA-seq and whole-genome bisulfite sequencing of adult females and males from two closely related parasitoid wasps, Nasonia vitripennis and Nasonia giraulti. More than 75% of expressed genes displayed significantly sex-biased expression. As a consequence, expression profiles are more similar between species within each sex than between sexes within each species. Furthermore, extremely male- and female-biased genes are enriched for totally different functional categories: male-biased genes for key enzymes in sex-pheromone synthesis and female-biased genes for genes involved in epigenetic regulation of gene expression. Remarkably, just 70 highly expressed, extremely male-biased genes account for 10% of all transcripts in adult males. Unlike expression profiles, DNA methylomes are highly similar between sexes within species, with no consistent sex differences in methylation found. Therefore, methylation changes cannot explain the extensive level of sex-biased gene expression observed. Female-biased genes have smaller sequence divergence between species, higher conservation to other hymenopterans, and a broader expression range across development. Overall, female-biased genes have been recruited from genes with more conserved and broadly expressing "house-keeping" functions, whereas male-biased genes are more recently evolved and are predominately testis specific. In summary, Nasonia accomplish a striking degree of sex-biased expression without sex chromosomes or epigenetic differences in methylation. We propose that methylation provides a general signal for constitutive gene expression, whereas other sex-specific signals cause sex-biased gene expression.

Introgression study reveals two quantitative trait loci involved in interspecific variation in memory retention among Nasonia wasp species

Genes involved in the process of memory formation have been studied intensively in model organisms; however, little is known about the mechanisms that are responsible for natural variation in memory dynamics. There is substantial variation in memory retention among closely related species in the parasitic wasp genus Nasonia. After a single olfactory conditioning trial, N. vitripennis consolidates long-term memory that lasts at least 6 days. Memory of the closely related species N. giraulti is present at 24 h but is lost within 2 days after a single trial. The genetic basis of this interspecific difference in memory retention was studied in a backcrossing experiment in which the phenotype of N. giraulti was selected for in the background of N. vitripennis for up to five generations. A genotyping microarray revealed five regions that were retained in wasps with decreased memory retention. Independent introgressions of individual candidate regions were created using linked molecular markers and tested for memory retention. One region on chromosome 1 (spanning ∼5.8 cM) and another on chromosome 5 (spanning ∼25.6 cM) resulted in decreased memory after 72 h, without affecting 24-h-memory retention. This phenotype was observed in both heterozygous and homozygous individuals. Transcription factor CCAAT/enhancer-binding protein and a dopamine receptor, both with a known function in memory formation, are within these genomic regions and are candidates for the regulation of memory retention. Concluding, this study demonstrates a powerful approach to study variation in memory retention and provides a basis for future research on its genetic basis.

Development of a Nasonia vitripennis outbred laboratory population for genetic analysis

The parasitoid wasp genus Nasonia has rapidly become a genetic model system for developmental and evolutionary biology. The release of its genome sequence led to the development of high-resolution genomic tools, for both interspecific and intraspecific research, which has resulted in great advances in understanding Nasonia biology. To further advance the utility of Nasonia vitripennis as a genetic model system and to be able to fully exploit the advantages of its fully sequenced and annotated genome, we developed a genetically variable and well-characterized experimental population. In this study, we describe the establishment of the genetically diverse HVRx laboratory population from strains collected from the field in the Netherlands. We established a maintenance method that retains genetic variation over generations of culturing in the laboratory. As a characterization of its genetic composition, we provide data on the standing genetic variation and estimate the effective population size (N(e)) by microsatellite analysis. A genome-wide description of polymorphism is provided through pooled resequencing, which yielded 417,331 high-quality SNPs spanning all five Nasonia chromosomes. The HVRx population and its characterization are freely available as a community resource for investigators seeking to elucidate the genetic basis of complex trait variation using the Nasonia model system.