Nasonia Genome Research Articles

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

BackgroundGene 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.ResultsWe 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.ConclusionsAltogether, our results provide insights into how GRNs respond to new functional demands and how they can incorporate novel components.Electronic supplementary materialThe online version of this article (doi:10.1186/s12915-016-0285-y) contains supplementary material, which is available to authorized users.

Fine-Scale Mapping of the Nasonia Genome to Chromosomes Using a High-Density Genotyping Microarray

Nasonia, a genus of four closely related parasitoid insect species, is a model system for genetic research. Their haplodiploid genetics (haploid males and diploid females) and interfertile species are advantageous for the genetic analysis of complex traits and the genetic basis of species differences. A fine-scale genomic map is an important tool for advancing genetic studies in this system. We developed and used a hybrid genotyping microarray to generate a high-resolution genetic map that covers 79% of the sequenced genome of Nasonia vitripennis. The microarray is based on differential hybridization of species-specific oligos between N. vitripennis and Nasonia giraulti at more than 20,000 markers spanning the Nasonia genome. The map places 729 scaffolds onto the five linkage groups of Nasonia, including locating many smaller scaffolds that would be difficult to map by other means. The microarray was used to characterize 26 segmental introgression lines containing chromosomal regions from one species in the genetic background of another. These segmental introgression lines have been used for rapid screening and mapping of quantitative trait loci involved in species differences. Finally, the microarray is extended to bulk-segregant analysis and genotyping of other Nasonia species combinations. These resources should further expand the usefulness of Nasonia for studies of the genetic basis and architecture of complex traits and speciation.

Comparative Analyses of DNA Methylation and Sequence Evolution Using Nasonia Genomes

The functional and evolutionary significance of DNA methylation in insect genomes remains to be resolved. Nasonia is well situated for comparative analyses of DNA methylation and genome evolution, since the genomes of a moderately distant outgroup species as well as closely related sibling species are available. Using direct sequencing of bisulfite-converted DNA, we uncovered a substantial level of DNA methylation in 17 of 18 Nasonia vitripennis genes and a strong correlation between methylation level and CpG depletion. Notably, in the sex-determining locus transformer, the exon that is alternatively spliced between the sexes is heavily methylated in both males and females, whereas other exons are only sparsely methylated. Orthologous genes of the honeybee and Nasonia show highly similar relative levels of CpG depletion, despite ~190 My divergence. Densely and sparsely methylated genes in these species also exhibit similar functional enrichments. We found that the degree of CpG depletion is negatively correlated with substitution rates between closely related Nasonia species for synonymous, nonsynonymous, and intron sites. This suggests that mutation rates increase with decreasing levels of germ line methylation. Thus, DNA methylation is prevalent in the Nasonia genome, may participate in regulatory processes such as sex determination and alternative splicing, and is correlated with several aspects of genome and sequence evolution.

Genetics of cuticular hydrocarbon differences between males of the parasitoid wasps Nasonia giraulti and Nasonia vitripennis

Many insects rely on cuticular hydrocarbons (CHCs) as major recognition signals between individuals. Previous research on the genetics of CHCs has focused on Drosophila in which the roles of three desaturases and one elongase were highlighted. Comparable studies in other insect taxa have not been conducted so far. Here, we explore the genetics of CHCs in hybrids of the jewel wasps Nasonia giraulti and Nasonia vitripennis. We analyzed the CHC profiles of pure strain and of F(2) hybrid males using gas chromatography coupled with mass spectrometry and distinguished 54 peaks, of which we identified 52 as straight-chain, monounsaturated, or methyl-branched CHCs. The latter compound class proved to be particularly abundant and diverse in Nasonia. Quantitative trait locus (QTL) analysis suggests fixed genetic differences between the two strains in 42 of the 54 studied traits, making Nasonia a promising genetic model for identifying genes involved in CHC biosynthesis. QTL for methyl-branched CHCs partly clustered in genomic regions with high recombination rate: a possible indication for pleiotropic genes that control their biosynthesis, which is largely unexplored so far. Finally, we identified and mapped genes in the Nasonia genome with high similarity to genes that have been implicated in alkene biosynthesis in Drosophila and discuss those that match in their position with predicted QTL for alkenes.

Genomics and Peptidomics of Neuropeptides and Protein Hormones Present in the Parasitic Wasp Nasonia vitripennis

Neuropeptides and protein hormones constitute a very important group of signaling molecules, regulating central physiological processes such as reproduction, development, and behavior. Using a bioinformatics approach, we screened the recently sequenced genome of the parasitic wasp, Nasonia vitripennis, for the presence of these signaling molecules and annotated 30 precursor genes encoding 51 different mature neuropeptides or protein hormones. Twenty-four of the predicted mature Nasonia neuropeptides could be experimentally confirmed by mass spectrometry. We also discovered a completely novel neuropeptide gene in Nasonia, coding for peptides containing the C-terminal sequence RYamide. This gene has orthologs in nearly all arthropods with a sequenced genome, and its expression in mosquitoes was confirmed by mass spectrometry. No precursor could be identified for N-terminally extended FMRFamides, even though their putative G protein coupled receptor (GPCR) is present in the Nasonia genome. Neither the precursor nor the putative receptor could be identified for allatostatin-B, capa, the glycoprotein hormones GPA2/GPB5, kinin, proctolin, sex peptide, and sulfakinin, arguing that these signaling systems are truly absent in the wasp. Also, antidiuretic factors, allatotropin, and NPLP-like precursors are missing in Nasonia, but here the receptors have not been identified in any insect, so far. Nasonia (Hymenoptera) has the lowest number of neuropeptide precursor genes compared to Drosophila melanogaster, Aedes aegypti (both Diptera), Bombyx mori (Lepidoptera), Tribolium castaneum (Coleoptera), Apis mellifera (Hymenoptera), and Acyrthosiphon pisum (Hemiptera). This lower number of neuropeptide genes might be related to Nasonia's parasitic life.

A comparison of recombination frequencies in intraspecific versus interspecific mapping populations of Nasonia

We present the first intraspecific linkage map for Nasonia vitripennis based on molecular markers. The map consists of 36 new microsatellite markers, extracted from the Nasonia genome sequence, and spans 515 cM. The five inferred linkage groups correspond to the five chromosomes of Nasonia. Comparison of recombination frequencies of the marker intervals spread over the whole genome (N=33 marker intervals) between the intraspecific N. vitripennis map and an interspecific N. vitripennis x N. giraulti map revealed a slightly higher (1.8%) recombination frequency in the intraspecific cross. We further considered an N. vitripennis x N. longicornis map with 29 microsatellite markers spanning 430 cM. Recombination frequencies in the two interspecific crosses differed neither between reciprocal crosses nor between mapping populations of embryos and adults. No major chromosomal rearrangements were found for the analyzed genomic segments. The observed differential F(2) hybrid male mortality has no significant effect on the genome-wide recombination frequency in Nasonia. We conclude that interspecific crosses between the different Nasonia species, a hallmark of Nasonia genetics, are generally suitable for mapping quantitative and qualitative trait loci for species differences.

Recombination and Its Impact on the Genome of the Haplodiploid Parasitoid Wasp Nasonia

Homologous meiotic recombination occurs in most sexually reproducing organisms, yet its evolutionary advantages are elusive. Previous research explored recombination in the honeybee, a eusocial hymenopteran with an exceptionally high genome-wide recombination rate. A comparable study in a non-social member of the Hymenoptera that would disentangle the impact of sociality from Hymenoptera-specific features such as haplodiploidy on the evolution of the high genome-wide recombination rate in social Hymenoptera is missing. Utilizing single-nucleotide polymorphisms (SNPs) between two Nasonia parasitoid wasp genomes, we developed a SNP genotyping microarray to infer a high-density linkage map for Nasonia. The map comprises 1,255 markers with an average distance of 0.3 cM. The mapped markers enabled us to arrange 265 scaffolds of the Nasonia genome assembly 1.0 on the linkage map, representing 63.6% of the assembled N. vitripennis genome. We estimated a genome-wide recombination rate of 1.4–1.5 cM/Mb for Nasonia, which is less than one tenth of the rate reported for the honeybee. The local recombination rate in Nasonia is positively correlated with the distance to the center of the linkage groups, GC content, and the proportion of simple repeats. In contrast to the honeybee genome, gene density in the parasitoid wasp genome is positively associated with the recombination rate; regions of low recombination are characterized by fewer genes with larger introns and by a greater distance between genes. Finally, we found that genes in regions of the genome with a low recombination frequency tend to have a higher ratio of non-synonymous to synonymous substitutions, likely due to the accumulation of slightly deleterious non-synonymous substitutions. These findings are consistent with the hypothesis that recombination reduces interference between linked sites and thereby facilitates adaptive evolution and the purging of deleterious mutations. Our results imply that the genomes of haplodiploid and of diploid higher eukaryotes do not differ systematically in their recombination rates and associated parameters.

Comparative analysis of nuclear tRNA genes of Nasonia vitripennis and other arthropods, and relationships to codon usage bias

Using bioinformatics methods, we identified a total of 221 and 199 tRNA genes in the nuclear genomes of Nasonia vitripennis and honey bee (Apis mellifera), respectively. We performed comparative analyses of Nasonia tRNA genes with honey bee and other selected insects to understand genomic distribution, sequence evolution and relationship of tRNA copy number with codon usage patterns. Many tRNA genes are located physically close to each other in the form of small clusters in the Nasonia genome. However, the number of clusters and the tRNA genes that form such clusters vary from species to species. In particular, the Ala-, Pro-, Tyr- and His-tRNA genes tend to accumulate in clusters in Nasonia but not in honey bee, whereas the bee contains a long cluster of 15 tRNA genes (of which 13 are Gln-tRNAs) that is absent in Nasonia. Though tRNA genes are highly conserved, contrasting patterns of nucleotide diversity are observed among the arm and loop regions of tRNAs between Nasonia and honey bee. Also, the sequence convergence between the reconstructed ancestral tRNAs and the present day tRNAs suggests a common ancestral origin of Nasonia and honey bee tRNAs. Furthermore, we also present evidence that the copy number of isoacceptor tRNAs (those having a different anticodon but charge the same amino acid) is correlated with codon usage patterns of highly expressed genes in Nasonia.

Maintenance of multiple lineages of R1 and R2 retrotransposable elements in the ribosomal RNA gene loci of Nasonia

Sequencing reads from the Nasonia genome project were used to study the ribosomal RNA gene loci and the retrotransposons R1 and R2 that insert specifically into the 28S genes. Five highly divergent R1 and five highly divergent R2 families were identified in the three sequenced species, as well as a non-autonomous element that appears to use the retrotransposition machinery of R1. A duplication of the R1 target site within the spacer region of the rDNA units was also found to be extensively utilized by R1 elements. We document numerous instances where the R1 and R2 families appropriated parts of the retrotransposition machinery of other lineages and speculate that this enables rapid adaptation and the maintenance of multiple R1 and R2 families.

Species-Diagnostic Single-Nucleotide Polymorphism and Sequence-Tagged Site Markers for the Parasitic Wasp Genus <I>Nasonia</I> (Hymenoptera: Pteromalidae)

Wasps of the genus Nasonia are important biological control agents of house flies and related filth flies, which are major vectors of human pathogens. Species of Nasonia (Hymenoptera: Pteromalidae) are not easily differentiated from one another by morphological characters, and molecular markers for their reliable identification have been missing so far. Here, we report eight single-nucleotide polymorphism and three sequence-tagged site markers derived from expressed sequenced tag libraries for the two closely related and regionally sympatric species N. giraulti and N. vitripennis. We studied variation of these markers in natural populations of the two species, and we mapped them in the Nasonia genome. The markers are species-diagnostic and evenly spread over all five chromosomes. They are ideal for rapid species identification and hybrid recognition, and they can be used to map economically relevant quantitative trait loci in the Nasonia genome.

Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements.

As an accompanying manuscript to the release of the honey bee genome, we report the entire sequence of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) ribosomal RNA (rRNA)-encoding gene sequences (rDNA) and related internally and externally transcribed spacer regions of Apis mellifera (Insecta: Hymenoptera: Apocrita). Additionally, we predict secondary structures for the mature rRNA molecules based on comparative sequence analyses with other arthropod taxa and reference to recently published crystal structures of the ribosome. In general, the structures of honey bee rRNAs are in agreement with previously predicted rRNA models from other arthropods in core regions of the rRNA, with little additional expansion in non-conserved regions. Our multiple sequence alignments are made available on several public databases and provide a preliminary establishment of a global structural model of all rRNAs from the insects. Additionally, we provide conserved stretches of sequences flanking the rDNA cistrons that comprise the externally transcribed spacer regions (ETS) and part of the intergenic spacer region (IGS), including several repetitive motifs. Finally, we report the occurrence of retrotransposition in the nuclear large subunit rDNA, as R2 elements are present in the usual insertion points found in other arthropods. Interestingly, functional R1 elements usually present in the genomes of insects were not detected in the honey bee rRNA genes. The reverse transcriptase products of the R2 elements are deduced from their putative open reading frames and structurally aligned with those from another hymenopteran insect, the jewel wasp Nasonia (Pteromalidae). Stretches of conserved amino acids shared between Apis and Nasonia are illustrated and serve as potential sites for primer design, as target amplicons within these R2 elements may serve as novel phylogenetic markers for Hymenoptera. Given the impending completion of the sequencing of the Nasonia genome, we expect our report eventually to shed light on the evolution of the hymenopteran genome within higher insects, particularly regarding the relative maintenance of conserved rDNA genes, related variable spacer regions and retrotransposable elements.