Recombination and its impact on the genome of the haplodiploid parasitoid wasp Nasonia

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

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.

Evolution of parasitoid host preference and performance in response to an invasive host acting as evolutionary trap

The invasion of a novel host species can create a mismatch in host choice and offspring survival (performance) when native parasitoids attempt to exploit the invasive host without being able to circumvent its resistance mechanisms. Invasive hosts can therefore act as evolutionary trap reducing parasitoids' fitness and this may eventually lead to their extinction. Yet, escape from the trap can occur when parasitoids evolve behavioral avoidance or a physiological strategy compatible with the trap host, resulting in either host‐range expansion or a complete host‐shift. We developed an individual based model to investigate which conditions promote parasitoids to evolve behavioral preference that matches their performance, including host‐trap avoidance, and which conditions lead to adaptations to the unsuitable hosts. The model was inspired by solitary endo‐parasitoids attacking larval host stages. One important aspect of these conditions was reduced host survival during incompatible interaction, where a failed parasitization attempt by a parasitoid resulted not only in death of her offspring but also in host killing. This non‐reproductive host mortality had a strong influence on the likelihood of establishment of novel host–parasitoid relationship, in some cases constraining adaptation to the trap host species. Moreover, our model revealed that host‐search efficiency and genetic variation in host‐preference play a key role in the likelihood that parasitoids will include the suboptimal host in their host range, or will evolve behavioral avoidance resulting in specialization and host‐range conservation, respectively. Hence, invasive species might change the evolutionary trajectory of native parasitoid species, which is important for predicting biocontrol ability of native parasitoids towards novel hosts.

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.

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.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO₂ detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Lack of parent-of-origin effects in Nasonia jewel wasp: A replication and extension study

In diploid cells, the paternal and maternal alleles are, on average, equally expressed. There are exceptions from this: a small number of genes express the maternal or paternal allele copy exclusively. This phenomenon, known as genomic imprinting, is common among eutherian mammals and some plant species; however, genomic imprinting in species with haplodiploid sex determination is not well characterized. Previous work reported no parent-of-origin effects in the hybrids of closely related haplodiploid Nasonia vitripennis and Nasonia giraulti jewel wasps, suggesting a lack of epigenetic reprogramming during embryogenesis in these species. Here, we replicate the gene expression dataset and observations using different individuals and sequencing technology, as well as reproduce these findings using the previously published RNA sequence data following our data analysis strategy. The major difference from the previous dataset is that they used an introgression strain as one of the parents and we found several loci that resisted introgression in that strain. Our results from both datasets demonstrate a species-of-origin effect, rather than a parent-of-origin effect. We present a reproducible workflow that others may use for replicating the results. Overall, we reproduced the original report of no parent-of-origin effects in the haplodiploid Nasonia using the original data with our new processing and analysis pipeline and replicated these results with our newly generated data.

Next-generation biological control: the need for integrating genetics and genomics

Biological control is widely successful at controlling pests, but effective biocontrol agents are now more difficult to import from countries of origin due to more restrictive international trade laws (the Nagoya Protocol). Coupled with increasing demand, the efficacy of existing and new biocontrol agents needs to be improved with genetic and genomic approaches. Although they have been underutilised in the past, application of genetic and genomic techniques is becoming more feasible from both technological and economic perspectives. We review current methods and provide a framework for using them. First, it is necessary to identify which biocontrol trait to select and in what direction. Next, the genes or markers linked to these traits need be determined, including how to implement this information into a selective breeding program. Choosing a trait can be assisted by modelling to account for the proper agro‐ecological context, and by knowing which traits have sufficiently high heritability values. We provide guidelines for designing genomic strategies in biocontrol programs, which depend on the organism, budget, and desired objective. Genomic approaches start with genome sequencing and assembly. We provide a guide for deciding the most successful sequencing strategy for biocontrol agents. Gene discovery involves quantitative trait loci analyses, transcriptomic and proteomic studies, and gene editing. Improving biocontrol practices includes marker‐assisted selection, genomic selection and microbiome manipulation of biocontrol agents, and monitoring for genetic variation during rearing and post‐release. We conclude by identifying the most promising applications of genetic and genomic methods to improve biological control efficacy.

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.

Chromosome aberrations in pressure-induced triploid Atlantic salmon

BACKGROUND

Triploid organisms have three sets of chromosomes. In Atlantic salmon, hydrostatic pressure treatment of newly fertilized eggs has been extensively used to produce triploids which are functionally sterile due to their unpaired chromosomes. These fish often perform poorly on commercial farms, sometimes without explanation. Inheritance patterns in individuals subjected to pressure treatment have not been investigated in Atlantic salmon thus far. However, work on other species suggests that this treatment can result in aberrant inheritance. We therefore studied this in Atlantic salmon by genotyping 16 polymorphic microsatellites in eyed eggs and juveniles which had been subjected to pressure-induction of triploidy. Communally reared juveniles including fish subjected to pressure-induction of triploidy and their diploid siblings were included as a control.

RESULTS

No diploid offspring were detected in any of the eggs or juveniles which were subjected to hydrostatic pressure; therefore, the induction of triploidy was highly successful. Aberrant inheritance was nevertheless observed in 0.9% of the eggs and 0.9% of the juveniles that had been subjected to pressure treatment. In the communally reared fish, 0.3% of the fish subjected to pressure treatment displayed aberrant inheritance, while their diploid controls displayed 0% aberrant inheritance. Inheritance errors included two eyed eggs lacking maternal DNA across all microsatellites, and, examples in both eggs and juveniles of either the maternal or paternal allele lacking in one of the microsatellites. All individuals displaying chromosome aberrations were otherwise triploid.

CONCLUSIONS

This is the first study to document aberrant inheritance in Atlantic salmon that have been subjected to pressure-induction of triploidy. Our experiments unequivocally demonstrate that even when induction of triploidy is highly successful, this treatment can cause chromosome aberrations in this species. Based upon our novel data, and earlier studies in other organisms, we hypothesize that in batches of Atlantic salmon where low to modest triploid induction rates have been reported, aberrant inheritance is likely to be higher than the rates observed here. Therefore, we tentatively suggest that this could contribute to the unexplained poor performance of triploid salmon that is occasionally reported in commercial aquaculture. These hypotheses require further investigation.

Variation in parasitoidism of Protocalliphora azurea (Diptera: Calliphoridae) by Nasonia vitripennis (Hymenoptera: Pteromalidae) in Spain

Parasitoid wasps may act as hyperparasites and sometimes regulate the populations of their hosts by a top-down dynamic. Nasonia vitripennis (Walker, 1836) is a generalist gregarious parasitoid that parasitizes several host flies, including the blowfly Protocalliphora Hough, 1899 (Diptera, Calliphoridae), which in turn parasitizes bird nestlings. Nonetheless, the ecological factors underlying N. vitripennis prevalence and parasitoidism intensity on its hosts in natural populations are poorly understood. We have studied the prevalence of N. vitripennis in Protocalliphora azurea (Fallén, 1817) puparia parasitizing wild populations of pied flycatcher (Ficedula hypoleuca) and blue tit (Cyanistes caeruleus) birds in two Mediterranean areas in central and southern Spain. We found some evidence that the prevalence of N. vitripennis was higher in moist habitats in southern Spain. A host-dependent effect was found, since the greater the number of P. azurea puparia, the greater the probability and rate of parasitoidism by the wasp. Our results also suggest that N. vitripennis parasitizes more P. azurea puparia in blue tit nests than in pied flycatcher nests as a consequence of a higher load of these flies in the former. Based on the high prevalence of N. vitripennis in P. azurea puparia in nature, we propose that this wasp may regulate blowfly populations, with possible positive effects on the reproduction of both bird species.

Extreme Differences in Recombination Rate between the Genomes of a Solitary and a Social Bee

Social insect genomes exhibit the highest rates of crossing over observed in plants and animals. The evolutionary causes of these extreme rates are unknown. Insight can be gained by comparing recombination rate variation across the genomes of related social and solitary insects. Here, we compare the genomic recombination landscape of the highly social honey bee, Apis mellifera, with the solitary alfalfa leafcutter bee, Megachile rotundata, by analyzing patterns of linkage disequilibrium in population-scale genome sequencing data. We infer that average recombination rates are extremely elevated in A. mellifera compared with M. rotundata. However, our results indicate that similar factors control the distribution of crossovers in the genomes of both species. Recombination rate is significantly reduced in coding regions in both species, with genes inferred to be germline methylated having particularly low rates. Genes with worker-biased patterns of expression in A. mellifera and their orthologs in M. rotundata have higher than average recombination rates in both species, suggesting that selection for higher diversity in genes involved in worker caste functions in social taxa is not the explanation for these elevated rates. Furthermore, we find no evidence that recombination has modulated the efficacy of selection among genes during bee evolution, which does not support the hypothesis that high recombination rates facilitated positive selection for new functions in social insects. Our results indicate that the evolution of sociality in insects likely entailed selection on modifiers that increased recombination rates genome wide, but that the genomic recombination landscape is determined by the same factors.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

Molecular evolutionary trends and feeding ecology diversification in the Hemiptera, anchored by the milkweed bug genome

BACKGROUND

The Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.

RESULTS

The 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.

CONCLUSIONS

With the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

The Genomic Architecture of a Rapid Island Radiation: Recombination Rate Variation, Chromosome Structure, and Genome Assembly of the Hawaiian Cricket Laupala

Phenotypic evolution and speciation depend on recombination in many ways. Within populations, recombination can promote adaptation by bringing together favorable mutations and decoupling beneficial and deleterious alleles. As populations diverge, crossing over can give rise to maladapted recombinants and impede or reverse diversification. Suppressed recombination due to genomic rearrangements, modifier alleles, and intrinsic chromosomal properties may offer a shield against maladaptive gene flow eroding coadapted gene complexes. Both theoretical and empirical results support this relationship. However, little is known about this relationship in the context of behavioral isolation, where coevolving signals and preferences are the major hybridization barrier. Here we examine the genomic architecture of recently diverged, sexually isolated Hawaiian swordtail crickets (Laupala). We assemble a de novo genome and generate three dense linkage maps from interspecies crosses. In line with expectations based on the species' recent divergence and successful interbreeding in the laboratory, the linkage maps are highly collinear and show no evidence for large-scale chromosomal rearrangements. Next, the maps were used to anchor the assembly to pseudomolecules and estimate recombination rates across the genome to test the hypothesis that loci involved in behavioral isolation (song and preference divergence) are in regions of low interspecific recombination. Contrary to our expectations, the genomic region where a male song and female preference QTL colocalize is not associated with particularly low recombination rates. This study provides important novel genomic resources for an emerging evolutionary genetics model system and suggests that trait-preference coevolution is not necessarily facilitated by locally suppressed recombination.

Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics

Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection on linked sites. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant and fungal species. Our quantitative analysis reveals a strong and taxonomically widespread reduction of the crossover rate in the centre of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centres increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their centre, whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossover is initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically widespread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of resulting patterns.

A high-density genetic map reveals variation in recombination rate across the genome of Daphnia magna

BACKGROUND

Recombination rate is an essential parameter for many genetic analyses. Recombination rates are highly variable across species, populations, individuals and different genomic regions. Due to the profound influence that recombination can have on intraspecific diversity and interspecific divergence, characterization of recombination rate variation emerges as a key resource for population genomic studies and emphasises the importance of high-density genetic maps as tools for studying genome biology. Here we present such a high-density genetic map for Daphnia magna, and analyse patterns of recombination rate across the genome.

RESULTS

A F2 intercross panel was genotyped by Restriction-site Associated DNA sequencing to construct the third-generation linkage map of D. magna. The resulting high-density map included 4037 markers covering 813 scaffolds and contigs that sum up to 77 % of the currently available genome draft sequence (v2.4) and 55 % of the estimated genome size (238 Mb). Total genetic length of the map presented here is 1614.5 cM and the genome-wide recombination rate is estimated to 6.78 cM/Mb. Merging genetic and physical information we consistently found that recombination rate estimates are high towards the peripheral parts of the chromosomes, while chromosome centres, harbouring centromeres in D. magna, show very low recombination rate estimates.

CONCLUSIONS

Due to its high-density, the third-generation linkage map for D. magna can be coupled with the draft genome assembly, providing an essential tool for genome investigation in this model organism. Thus, our linkage map can be used for the on-going improvements of the genome assembly, but more importantly, it has enabled us to characterize variation in recombination rate across the genome of D. magna for the first time. These new insights can provide a valuable assistance in future studies of the genome evolution, mapping of quantitative traits and population genetic studies.

Quantitative Trait Locus Analysis of Mating Behavior and Male Sex Pheromones in Nasonia Wasps

A major focus in speciation genetics is to identify the chromosomal regions and genes that reduce hybridization and gene flow. We investigated the genetic architecture of mating behavior in the parasitoid wasp species pair Nasonia giraulti and Nasonia oneida that exhibit strong prezygotic isolation. Behavioral analysis showed that N. oneida females had consistently higher latency times, and broke off the mating sequence more often in the mounting stage when confronted with N. giraulti males compared with males of their own species. N. oneida males produce a lower quantity of the long-range male sex pheromone (4R,5S)-5-hydroxy-4-decanolide (RS-HDL). Crosses between the two species yielded hybrid males with various pheromone quantities, and these males were used in mating trials with females of either species to measure female mate discrimination rates. A quantitative trait locus (QTL) analysis involving 475 recombinant hybrid males (F₂), 2148 reciprocally backcrossed females (F₃), and a linkage map of 52 equally spaced neutral single nucleotide polymorphism (SNP) markers plus SNPs in 40 candidate mating behavior genes revealed four QTL for male pheromone amount, depending on partner species. Our results demonstrate that the RS-HDL pheromone plays a role in the mating system of N. giraulti and N. oneida, but also that additional communication cues are involved in mate choice. No QTL were found for female mate discrimination, which points at a polygenic architecture of female choice with strong environmental influences.

Genetic differentiation within and away from the chromosomal rearrangements characterising hybridising chromosomal races of the western house mouse (Mus musculus domesticus)

The importance of chromosomal rearrangements for speciation can be inferred from studies of genetic exchange between hybridising chromosomal races within species. Reduced fertility or recombination suppression in karyotypic hybrids has the potential to maintain or promote genetic differentiation in genomic regions near rearrangement breakpoints. We studied genetic exchange between two hybridising groups of chromosomal races of house mouse in Upper Valtellina (Lombardy, Italy), using microsatellites. These groups differ by Robertsonian fusions and/or whole-arm reciprocal translocations such that F1 hybrids have a chain-of-five meiotic configuration. Previous studies showed genetic differentiation in two chromosomes in the chain-of-five (10 and 12) close to their centromeres (i.e. the rearrangement breakpoints); we have shown here that the centromeric regions of the other two chromosomes in the chain (2 and 8) are similarly differentiated. The internal chromosomes of the chain (8 and 12) show the greatest differentiation, which may reflect pairing and recombination properties of internal and external elements in a meiotic chain. Importantly, we found that centromeric regions of some non-rearranged chromosomes also showed genetic differentiation between the hybridising groups, indicating a complex interplay between chromosomal rearrangements and other parts of the genome in maintaining or promoting differentiation and potentially driving speciation between chromosomal races.

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.

Causes and consequences of crossing-over evidenced via a high-resolution recombinational landscape of the honey bee

Background

Social hymenoptera, the honey bee (Apis mellifera) in particular, have ultra-high crossover rates and a large degree of intra-genomic variation in crossover rates. Aligned with haploid genomics of males, this makes them a potential model for examining the causes and consequences of crossing over. To address why social insects have such high crossing-over rates and the consequences of this, we constructed a high-resolution recombination atlas by sequencing 55 individuals from three colonies with an average marker density of 314 bp/marker.

Results

We find crossing over to be especially high in proximity to genes upregulated in worker brains, but see no evidence for a coupling with immune-related functioning. We detect only a low rate of non-crossover gene conversion, contrary to current evidence. This is in striking contrast to the ultrahigh crossing-over rate, almost double that previously estimated from lower resolution data. We robustly recover the predicted intragenomic correlations between crossing over and both population level diversity and GC content, which could be best explained as indirect and direct consequences of crossing over, respectively.

Conclusions

Our data are consistent with the view that diversification of worker behavior, but not immune function, is a driver of the high crossing-over rate in bees. While we see both high diversity and high GC content associated with high crossing-over rates, our estimate of the low non-crossover rate demonstrates that high non-crossover rates are not a necessary consequence of high recombination rates.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0566-0) contains supplementary material, which is available to authorized users.

Octopamine-like immunoreactive neurons in the brain and subesophageal ganglion of the parasitic wasps Nasonia vitripennis and N. giraulti

Octopamine is an important neuromodulator in the insect nervous system, influencing memory formation, sensory perception and motor control. In this study, we compare the distribution of octopamine-like immunoreactive neurons in two parasitic wasp species of the Nasonia genus, N. vitripennis and N. giraulti. These two species were previously described as differing in their learning and memory formation, which raised the question as to whether morphological differences in octopaminergic neurons underpinned these variations. Immunohistochemistry in combination with confocal laser scanning microscopy was used to reveal and compare the somata and major projections of the octopaminergic neurons in these wasps. The brains of both species showed similar staining patterns, with six different neuron clusters being identified in the brain and five different clusters in the subesophageal ganglion. Of those clusters found in the subesophageal ganglion, three contained unpaired neurons, whereas the other three consisted in paired neurons. The overall pattern of octopaminergic neurons in both species was similar, with no differences in the numbers or projections of the ventral unpaired median (VUM) neurons, which are known to be involved in memory formation in insects. In one other cluster in the brain, located in-between the optic lobe and the antennal lobe, we detected more neurons in N. vitripennis compared with N. giraulti. Combining our results with findings made previously in other Hymenopteran species, we discuss possible functions and some of the ultimate factors influencing the evolution of the octopaminergic system in the insect brain.

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.

High-resolution linkage map for two honeybee chromosomes: the hotspot quest

Meiotic recombination is a fundamental process ensuring proper disjunction of homologous chromosomes and allele shuffling in successive generations. In many species, this cellular mechanism occurs heterogeneously along chromosomes and mostly concentrates in tiny fragments called recombination hotspots. Specific DNA motifs have been shown to initiate recombination in these hotspots in mammals, fission yeast and drosophila. The aim of this study was to check whether recombination also occurs in a heterogeneous fashion in the highly recombinogenic honeybee genome and whether this heterogeneity can be connected with specific DNA motifs. We completed a previous picture drawn from a routine genetic map built with an average resolution of 93 kb. We focused on the two smallest honeybee chromosomes to increase the resolution and even zoomed at very high resolution (3.6 kb) on a fragment of 300 kb. Recombination rates measured in these fragments were placed in relation with occurrence of 30 previously described motifs through a Poisson regression model. A selection procedure suitable for correlated variables was applied to keep significant motifs. These fine and ultra-fine mappings show that recombination rate is significantly heterogeneous although poorly contrasted between high and low recombination rate, contrarily to most model species. We show that recombination rate is probably associated with the DNA methylation state. Moreover, three motifs (CGCA, GCCGC and CCAAT) are good candidates of signals promoting recombination. Their influence is however moderate, doubling at most the recombination rate. This discovery extends the way to recombination dissection in insects.

Cuticular hydrocarbon divergence in the jewel wasp Nasonia: evolutionary shifts in chemical communication channels?

The evolution and maintenance of intraspecific communication channels constitute a key feature of chemical signalling and sexual communication. However, how divergent chemical communication channels evolve while maintaining their integrity for both sender and receiver is poorly understood. In this study, we compare male and female cuticular hydrocarbon (CHC) profiles in the jewel wasp genus Nasonia, analyse their chemical divergence and investigate their role as species-specific sexual signalling cues. Males and females of all four Nasonia species showed unique, nonoverlapping CHC profiles unambiguously separating them. Surprisingly, male and female phylogenies based on the chemical distances between their CHC profiles differed dramatically, where only male CHC divergence parallels the molecular phylogeny of Nasonia. In particular, N. giraulti female CHC profiles were the most divergent from all other species and very different from its most closely related sibling species N. oneida. Furthermore, although our behavioural assays indicate that female CHC profiles can generally be perceived as sexual cues attracting males in Nasonia, this function has apparently been lost in the highly divergent female N. giraulti CHC profiles. Curiously, N. giraulti males are still attracted to heterospecific, but not to conspecific female CHC profiles. We suggest that this striking discrepancy has been caused by an extensive evolutionary shift in female N. giraulti CHC profiles, which are no longer used as conspecific recognition cues. Our study constitutes the first report of an apparent abandonment of a sexual recognition cue that the receiver did not adapt to.

How hot are drosophila hotspots? examining recombination rate variation and associations with nucleotide diversity, divergence, and maternal age in Drosophila pseudoobscura

Fine scale meiotic recombination maps have uncovered a large amount of variation in crossover rate across the genomes of many species, and such variation in mammalian and yeast genomes is concentrated to <5kb regions of highly elevated recombination rates (10-100x the background rate) called "hotspots." Drosophila exhibit substantial recombination rate heterogeneity across their genome, but evidence for these highly-localized hotspots is lacking. We assayed recombination across a 40Kb region of Drosophila pseudoobscura chromosome 2, with one 20kb interval assayed every 5Kb and the adjacent 20kb interval bisected into 10kb pieces. We found that recombination events across the 40kb stretch were relatively evenly distributed across each of the 5kb and 10kb intervals, rather than concentrated in a single 5kb region. This, in combination with other recent work, indicates that the recombination landscape of Drosophila may differ from the punctate recombination pattern observed in many mammals and yeast. Additionally, we found no correlation of average pairwise nucleotide diversity and divergence with recombination rate across the 20kb intervals, nor any effect of maternal age in weeks on recombination rate in our sample.

Behavioural and genetic analyses of Nasonia shed light on the evolution of sex pheromones

Sex pheromones play a pivotal role in the communication of many sexually reproducing organisms. Accordingly, speciation is often accompanied by pheromone diversification enabling proper mate finding and recognition. Current theory implies that chemical signals are under stabilizing selection by the receivers who thereby maintain the integrity of the signals. How the tremendous diversity of sex pheromones seen today evolved is poorly understood. Here we unravel the genetics of a newly evolved pheromone phenotype in wasps and present results from behavioural experiments indicating how the evolution of a new pheromone component occurred in an established sender-receiver system. We show that male Nasonia vitripennis evolved an additional pheromone compound differing only in its stereochemistry from a pre-existing one. Comparative behavioural studies show that conspecific females responded neutrally to the new pheromone phenotype when it evolved. Genetic mapping and gene knockdown show that a cluster of three closely linked genes accounts for the ability to produce this new pheromone phenotype. Our data suggest that new pheromone compounds can persist in a sender's population, without being selected against by the receiver and without the receiver having a pre-existing preference for the new pheromone phenotype, by initially remaining unperceived. Our results thus contribute valuable new insights into the evolutionary mechanisms underlying the diversification of sex pheromones. Furthermore, they indicate that the genetic basis of new pheromone compounds can be simple, allowing them to persist long enough in a population for receivers to evolve chemosensory adaptations for their exploitation.

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.

The ant genomes have been invaded by several types of mariner transposable elements

To date, only three types of full-length mariner elements have been described in ants, each one in a different genus of the Myrmicinae subfamily: Sinvmar was isolated from various Solenopsis species, Myrmar from Myrmica ruginodis, and Mboumar from Messor bouvieri. In this study, we report the coexistence of three mariner elements (Tnigmar-Si, Tnigmar-Mr, and Tnigmar-Mb) in the genome of a single species, Tapinoma nigerrimum (subfamily Dolichoderinae). Molecular evolutionary analyses of the nucleotide sequence data revealed a general agreement between the evolutionary history of most the elements and the ant species that harbour them, and suggest that they are at the vertical inactivation stage of the so-called Mariner Life Cycle. In contrast, significantly reduced levels of synonymous divergence between Mboumar and Tnigmar-Mb and between Myrmar and Botmar (a mariner element isolated from Bombus terrestris), relative to those observed between their hosts, suggest that these elements arrived to the species that host them by horizontal transfer, long after the species' split. The horizontal transfer events for the two pairs of elements could be roughly dated within the last 2 million years and about 14 million years, respectively. As would be expected under this scenario, the coding sequences of the youngest elements, Tnigmar-Mb and Mboumar, are intact and, thus, potentially functional. Each mariner element has a different chromosomal distribution pattern according to their stage within the Mariner Life Cycle. Finally, a new defective transposable element (Azteca) has also been found inserted into the Tnigmar-Mr sequences showing that the ant genomes have been invaded by at least four different types of mariner elements.

Recombination is associated with the evolution of genome structure and worker behavior in honey bees

The rise of insect societies, marked by the formation of reproductive and sterile castes, represents a major unsolved mystery in evolution. Across several independent origins of sociality, the genomes of social hymenopterans share two peculiar attributes: high recombination and low but heterogeneous GC content. For example, the genome of the honey bee, Apis mellifera, represents a mosaic of GC-poor and GC-rich regions with rates of recombination an order of magnitude higher than in humans. However, it is unclear how heterogeneity in GC content arises, and how it relates to the expression and evolution of worker traits. Using population genetic analyses, we demonstrate a bias in the allele frequency and fixation rate of derived C or G mutations in high-recombination regions, consistent with recombination's causal influence on GC-content evolution via biased gene conversion. We also show that recombination and biased gene conversion actively maintain the heterogeneous GC content of the honey bee genome despite an overall A/T mutation bias. Further, we found that GC-rich genes and intergenic regions have higher levels of genetic diversity and divergence relative to GC-poor regions, also consistent with recombination's causal influence on the rate of molecular evolution. Finally, we found that genes associated with behavior and those with worker-biased expression are found in GC-rich regions of the bee genome and also experience high rates of molecular evolution. Taken together, these findings suggest that recombination acts to maintain a genetically diverse and dynamic part of the genome where genes underlying worker behavior evolve more quickly.

Genetic variability of arrhenotokous and thelytokous Venturia canescens (Hymenoptera)

The ichneumonid wasp Venturia canescens (Hymenoptera) has been studied extensively for foraging behaviour and population dynamics of sexually (arrhenotokous) and parthenogenetically (thelytokous) reproducing individuals. Here we report the development of a set of microsatellite markers for V.canescens and use them to show that arrhenotokous individuals have more genetic variability than thelytokous ones, which are even homozygous for all tested loci. Crosses between arrhenotokous individuals suggested one marker, Vcan071, to be linked with the Complementary Sex Determiner (CSD) locus and one, Vcan109, with the Virus Like Protein (vlp-p40) locus. The genome size of V. canescens was estimated to be 274–279 Mb. We discuss how both reproductive modes can give rise to the observed genetic variability and how the new markers can be used for future genetic studies of V. canescens.

Quantitative trait locus analysis in haplodiploid hymenoptera

This article describes QTL analyses for solitary (Nasonia, a parasitoid wasp) and social hymenopteran species (honeybee and bumblebee). These exemplar QTL analyses determined the genetic basis of morphological, behavioral, and colony level traits. Mapping populations were derived either from lab crosses between highly inbred strains (Nasonia spp.), lab crosses of individuals caught in the field (bumblebees), or offspring from artificially inseminated queens from a managed honeybee population. Using these examples, we demonstrate the importance of a clear understanding of the life history, breeding, and reproductive system of the organism used for a QTL analysis, e.g., haplo-diploidy or reproductive division of labor in social insects. We lead the reader step by step through the process of multiple QTL analyses and describe potential problems and roadblocks (e.g., data handling, statistical, and biological problems) that can obscure or severely impair the results of a QTL study and how to detect and deal with those problems.In particular, we provide a way to empirically estimate the Beavis effect for a larger QTL mapping population and how to estimate a more accurate value of the explained phenotypic variance of each detected QTL using a resampling procedure.

A high recombination rate in eusocial Hymenoptera: evidence from the common wasp Vespula vulgaris

Background

High recombination rates have previously been detected in two groups of eusocial insects; honeybees and ants. In this study we estimate recombination rate in a eusocial wasp Vespula vulgaris that represents a third phylogenetic lineage within eusocial hymenopterans.

Results

A genetic linkage map of V. vulgaris based on 210 markers shows that the total map length is 2129 cM and the recombination rate is 9.7 cM/Mb (or 103 kb/cM). The present estimate in V. vulgaris is somewhat smaller than in the honeybee Apis mellifera and intermediate between the estimates from two ant species (Acromyrmex echinatior, Pogonomyrmex rugosus). Altogether, the estimates from these eusocial species are higher than in any other insect reported so far.

Conlusions

The four species (V. vulgaris, A. mellifera, A. echinatior, P. rugosus) are characterized by advanced eusociality with large colonies, clear queen-worker dimorphism and well developed task specialization. They also have colonies with a single, normally multiply inseminated (polyandrous) queen. Benefits of genotypic diversity within colonies (e.g. through improved task specialization or pathogen and parasite resistance) may have selected for both polyandry and high recombination rate in such advanced eusocial insects.

Caenorhabditis briggsae recombinant inbred line genotypes reveal inter-strain incompatibility and the evolution of recombination

The nematode Caenorhabditis briggsae is an emerging model organism that allows evolutionary comparisons with C. elegans and exploration of its own unique biological attributes. To produce a high-resolution C. briggsae recombination map, recombinant inbred lines were generated from reciprocal crosses between two strains and genotyped at over 1,000 loci. A second set of recombinant inbred lines involving a third strain was also genotyped at lower resolution. The resulting recombination maps exhibit discrete domains of high and low recombination, as in C. elegans, indicating these are a general feature of Caenorhabditis species. The proportion of a chromosome's physical size occupied by the central, low-recombination domain is highly correlated between species. However, the C. briggsae intra-species comparison reveals striking variation in the distribution of recombination between domains. Hybrid lines made with the more divergent pair of strains also exhibit pervasive marker transmission ratio distortion, evidence of selection acting on hybrid genotypes. The strongest effect, on chromosome III, is explained by a developmental delay phenotype exhibited by some hybrid F2 animals. In addition, on chromosomes IV and V, cross direction-specific biases towards one parental genotype suggest the existence of cytonuclear epistatic interactions. These interactions are discussed in relation to surprising mitochondrial genome polymorphism in C. briggsae, evidence that the two strains diverged in allopatry, the potential for local adaptation, and the evolution of Dobzhansky-Muller incompatibilities. The genetic and genomic resources resulting from this work will support future efforts to understand inter-strain divergence as well as facilitate studies of gene function, natural variation, and the evolution of recombination in Caenorhabditis nematodes.

The quantitative genetic basis of sex ratio variation in Nasonia vitripennis: a QTL study

Our understanding of how natural selection should shape sex allocation is perhaps more developed than for any other trait. However, this understanding is not matched by our knowledge of the genetic basis of sex allocation. Here, we examine the genetic basis of sex ratio variation in the parasitoid wasp Nasonia vitripennis, a species well known for its response to local mate competition (LMC). We identified a quantitative trait locus (QTL) for sex ratio on chromosome 2 and three weaker QTL on chromosomes 3 and 5. We tested predictions that genes associated with sex ratio should be pleiotropic for other traits by seeing if sex ratio QTL co-occurred with clutch size QTL. We found one clutch size QTL on chromosome 1, and six weaker QTL across chromosomes 2, 3 and 5, with some overlap to regions associated with sex ratio. The results suggest rather limited scope for pleiotropy between these traits.

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.

Genetic and evolutionary correlates of fine-scale recombination rate variation in Drosophila persimilis

Recombination is fundamental to meiosis in many species and generates variation on which natural selection can act, yet fine-scale linkage maps are cumbersome to construct. We generated a fine-scale map of recombination rates across two major chromosomes in Drosophila persimilis using 181 SNP markers spanning two of five major chromosome arms. Using this map, we report significant fine-scale heterogeneity of local recombination rates. However, we also observed "recombinational neighborhoods," where adjacent intervals had similar recombination rates after excluding regions near the centromere and telomere. We further found significant positive associations of fine-scale recombination rate with repetitive element abundance and a 13-bp sequence motif known to associate with human recombination rates. We noted strong crossover interference extending 5-7 Mb from the initial crossover event. Further, we observed that fine-scale recombination rates in D. persimilis are strongly correlated with those obtained from a comparable study of its sister species, D. pseudoobscura. We documented a significant relationship between recombination rates and intron nucleotide sequence diversity within species, but no relationship between recombination rate and intron divergence between species. These results are consistent with selection models (hitchhiking and background selection) rather than mutagenic recombination models for explaining the relationship of recombination with nucleotide diversity within species. Finally, we found significant correlations between recombination rate and GC content, supporting both GC-biased gene conversion (BGC) models and selection-driven codon bias models. Overall, this genome-enabled map of fine-scale recombination rates allowed us to confirm findings of broader-scale studies and identify multiple novel features that merit further investigation.

The parasitoid wasp Nasonia: an emerging model system with haploid male genetics

Nasonia is a complex of four closely related species of wasps that is rapidly emerging as a model for evolutionary and developmental genetics. It has several features that make it an excellent genetic system, including a short generation time, ease of rearing, interfertile species, visible and molecular markers, and a sequenced genome. The form of sex determination, called "haplodiploidy," makes Nasonia particularly suitable as a genetic tool. Females are diploid and develop from fertilized eggs, whereas males are haploid and develop from unfertilized eggs. This allows geneticists to exploit many of the advantages of haploid genetics in an otherwise complex eukaryotic organism. Nasonia readily inbreeds, permitting production of isogenic lines, and the four species in the genus are interfertile (after removal of the endosymbiont Wolbachia), facilitating movement of genes between the species for efficient positional cloning of quantitative trait loci (QTL). Genome sequencing of the genetic model Nasonia vitripennis and two interfertile species, Nasonia giraulti and Nasonia longicornis, is now completed. This genome project provides a wealth of interspecies polymorphisms (e.g., single nucleotide polymorphisms [SNPs], insertion-deletions [indels], microsatellites) to facilitate positional cloning of genes involved in species differences in behavior, morphology, and development. Advances in the genetics of this system also open a path for improvement of parasitoid insects as agents of pest control.

The genetic basis of interspecies host preference differences in the model parasitoid Nasonia

The genetic basis of host preference has been investigated in only a few species. It is relevant to important questions in evolutionary biology, including sympatric speciation, generalist versus specialist adaptation, and parasite-host co-evolution. Here we show that a major locus strongly influences host preference in Nasonia. Nasonia are parasitic wasps that utilize fly pupae; Nasonia vitripennis is a generalist that parasitizes a diverse set of hosts, whereas Nasonia giraulti specializes in Protocalliphora (bird blowflies). In laboratory choice experiments using Protocalliphora and Sarcophaga (flesh flies), N. vitripennis shows a preference for Sarcophaga, whereas N. giraulti shows a preference for Protocalliphora. Through a series of interspecies crosses, we have introgressed a major locus affecting host preference from N. giraulti into N. vitripennis. The N. giraulti allele is dominant and greatly increases preference for Protocalliphora pupae in the introgression line relative to the recessive N. vitripennis allele. Through the utilization of a Nasonia genotyping microarray, we have identified the introgressed region as 16 Mb of chromosome 4, although a more complete analysis is necessary to determine the exact genetic architecture of host preference in the genus. To our knowledge, this is the first introgression of the host preference of one parasitoid species into another, as well as one of the few cases of introgression of a behavioral gene between species.

Metabolic enzymes associated with xenobiotic and chemosensory responses in Nasonia vitripennis

The numbers of glutathione S-transferase, cytochrome P450 and esterase genes in the genome of the hymenopteran parasitoid Nasonia vitripennis are about twice those found in the genome of another hymenopteran, the honeybee Apis mellifera. Some of the difference is associated with clades of these families implicated in xenobiotic resistance in other insects and some is in clades implicated in hormone and pheromone metabolism. The data support the hypothesis that the eusocial behaviour of the honeybee and the concomitant homeostasis of the nest environment may obviate the need for as many gene/enzyme systems associated with xenobiotic metabolism as are found in other species, including N. vitripennis, that are thought to encounter a wider range of potentially toxic xenobiotics in their diet and habitat.

The distribution of microsatellites in the Nasonia parasitoid wasp genome

Microsatellites are important molecular markers used in numerous genetic contexts. Despite this widespread use, the evolutionary processes governing microsatellite distribution and diversity remain controversial. Here, we present results on the distribution of microsatellites of three species in the parasitic wasp genus Nasonia generated by an in silico data-mining approach. Our results show that the overall microsatellite density in Nasonia is comparable to that of the honey bee, but much higher than in eight non-Hymenopteran arthropods. Across the Nasonia vitripennis genome, microsatellite density varied both within and amongst chromosomes. In contrast to other taxa, dinucleotides are the most abundant repeat type in all four species of Hymenoptera studied. Whether the differences between the Hymenoptera and other taxa are of functional significance remains to be determined.