Search

Abstract Houseflies provide a good experimental model to study the initial evolutionary stages of a primary sex-determining locus because they possess different recently evolved proto-Y chromosomes that contain male-determining loci (M) with the same male-determining gene, Mdmd. We investigate M-loci genomically and cytogenetically revealing distinct molecular architectures among M-loci. M on chromosome V (M V ) has two intact Mdmd copies in a palindrome. M on chromosome III (M III ) has tandem duplications containing 88 Mdmd copies (only one intact) and various repeats, including repeats that are XY-prevalent. M on chromosome II (M II ) and the Y (M Y ) share M III -like architecture, but with fewer repeats. M Y additionally shares M V -specific sequence arrangements. Based on these data and karyograms using two probes, one derives from M III and one Mdmd-specific, we infer evolutionary histories of polymorphic M-loci, which have arisen from unique translocations of Mdmd, embedded in larger DNA fragments, and diverged independently into regions of varying complexity.

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

Adaptation of complex traits depends on standing genetic variation at multiple loci. The allelic variants that have positive fitness effects, however, can differ depending on the genetic background and the selective pressure. Previously, we interrogated the Drosophila melanogaster genome at the population level for polymorphic positions and identified 215 single nucleotide polymorphisms (SNPs) that had significantly changed in frequency after experimental evolution for increased parasitoid resistance. In the current study, we follow up on 11 of these SNPs as putative targets of the experimental selection process (Jalvingh et al., 2014). We study the patterns of genetic variation for these SNPs in several European field populations. Furthermore, we associate the genetic variation of these SNPs to variation in resistance against the parasitoid Asobara tabida, by determining the individual phenotype and SNP genotype for 144 individuals from four Selection lines and four non-selected Control lines and for 400 individuals from 12 Field lines that differ in parasitoid resistance. For the Selection lines we additionally monitored the changes in allele frequencies throughout the five generations of experimental selection. For three genes, mbl (Zn-finger protein), mthl4 (G-protein coupled receptor) and CG17287 (protein-cysteine S-palmitoyltransferase) individual SNP genotypes were significantly associated with resistance level in the Selection and Control lines. Additionally, the minor allele in mbl and mthl4 were consistently and gradually favored throughout the five generations of experimental evolution. However, none of these alleles did appear to be associated to high resistance in the Field lines. We suggest that, within field populations, selection for parasitoid resistance is a gradual process that involves co-adapted gene complexes. Fast artificial selection, however, enforces the sudden cumulating of particular alleles that confer high resistance (genetic sweep). We discuss our findings in the context of local adaptation.

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 (F2), 2148 reciprocally backcrossed females (F3), 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.

Many physiological processes of living organisms show circadian rhythms, governed by an endogenous clock. This clock has a genetic basis and is entrained by external cues, such as light and temperature. Other physiological processes exhibit seasonal rhythms, that are also responsive to light and temperature. We previously reported a natural latitudinal cline of photoperiodic diapause induction in the parasitic wasp Nasonia vitripennis in Europe and a correlated haplotype frequency for the circadian clock gene period (per). To evaluate if this correlation is reflected in circadian behaviour, we investigated the circadian locomotor activity of seven populations from the cline. We found that the proportion of rhythmic males was higher than females in constant darkness, and that mating decreased rhythmicity of both sexes. Only for virgin females, the free running period (τ) increased weakly with latitude. Wasps from the most southern locality had an overall shorter free running rhythm and earlier onset, peak, and offset of activity during the 24 h period, than wasps from the northernmost locality. We evaluated this variation in rhythmicity as a function of period haplotype frequencies in the populations and discussed its functional significance in the context of local adaptation.

An attractive way to improve our understanding of sex determination evolution is to study the underlying mechanisms in closely related species and in a phylogenetic perspective. Hymenopterans are well suited owing to the diverse sex determination mechanisms, including different types of Complementary Sex Determination (CSD) and maternal control sex determination. We investigated different types of CSD in four species within the braconid wasp genus Asobara that exhibit diverse life-history traits. Nine to thirteen generations of inbreeding were monitored for diploid male production, brood size, offspring sex ratio, and pupal mortality as indicators for CSD. In addition, simulation models were developed to compare these observations to predicted patterns for multilocus CSD with up to ten loci. The inbreeding regime did not result in diploid male production, decreased brood sizes, substantially increased offspring sex ratios nor in increased pupal mortality. The simulations further allowed us to reject CSD with up to ten loci, which is a strong refutation of the multilocus CSD model. We discuss how the absence of CSD can be reconciled with the variation in life-history traits among Asobara species, and the ramifications for the phylogenetic distribution of sex determination mechanisms in the Hymenoptera.

Although sex determination is a universal process in sexually reproducing organisms, sex determination pathways are among the most highly variable genetic systems found in nature. Nevertheless, general principles can be identified among the diversity, like the central role of transformer (tra) in insects. When a functional TRA protein is produced in early embryogenesis, the female sex determining route is activated, while prevention of TRA production leads to male development. In dipterans, male development is achieved by prevention of female-specific splicing of tra mRNA, either mediated by X-chromosome dose or masculinizing factors. In Hymenoptera, which have haplodiploid sex determination, complementary sex determination and maternal imprinting have been identified to regulate timely TRA production. In the parasitoid Nasonia, zygotic transformer (Nvtra) expression and splicing is regulated by a combination of maternal provision of Nvtra mRNA and silencing of Nvtra expression in unfertilized eggs. It is unclear, however, if this silencing is directly on the tra locus or whether it is mediated through maternal silencing of a trans-acting factor. Here we show that in Nasonia, female sex determination is dependent on zygotic activation of Nvtra expression by an as yet unknown factor. This factor, which we propose to term womanizer (wom), is maternally silenced during oogenesis to ensure male development in unfertilized eggs. This finding implicates the upstream recruitment of a novel gene in the Nasonia sex determining cascade and supports the notion that sex determining cascades can rapidly change by adding new components on top of existing regulators.

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.

The genetic basis of morphological differences among species is still poorly understood. We investigated the genetic basis of sex-specific differences in wing size between two closely related species of Nasonia by positional cloning a major male-specific locus, wing-size1 (ws1). Male wing size increases by 45% through cell size and cell number changes when the ws1 allele from N. giraulti is backcrossed into a N. vitripennis genetic background. A positional cloning approach was used to fine-scale map the ws1 locus to a 13.5 kilobase region. This region falls between prospero (a transcription factor involved in neurogenesis) and the master sex-determining gene doublesex. It contains the 5'-UTR and cis-regulatory domain of doublesex, and no coding sequence. Wing size reduction correlates with an increase in doublesex expression level that is specific to developing male wings. Our results indicate that non-coding changes are responsible for recent divergence in sex-specific morphology between two closely related species. We have not yet resolved whether wing size evolution at the ws1 locus is caused by regulatory alterations of dsx or prospero, or by another mechanism. This study demonstrates the feasibility of efficient positional cloning of quantitative trait loci (QTL) involved in a broad array of phenotypic differences among Nasonia species.

During the past 50 years the number of badgers (Meles meles) in Denmark has declined by c. 50%. To assess the genetic consequences of the demographic decline, six DNA-microsatellite loci were used to analyse 139 badger tissue-samples, which were collected in 1995–98 from three zones (1, 2 and 3) in Jutland (Denmark). Results from contemporary samples were compared to data obtained from DNA extracted from teeth belonging to 39 badgers collected from zones 1 and 2 in 1957–66. The microsatellites showed a low to moderate polymorphism, and the within area genetic diversity (HE) was relatively low (0.308≤HE≤0.441) as compared to other carnivores and mustelids. Genetic differentiation among the three recent sampling zones was significant (FST=0.072, P<0.0001). The FST values between the two recent (FST=0.092, P<0.0001) and the two historical samples (zones 1 and 2; FST=0.073, P<0.0001) were significant and similar. The HE and FST of the recent and historical zones were not significantly different (P>0.05). Despite the drastic recent decline, Danish badgers did not show apparent signs of genetic bottlenecks in both recent and historical samples. The use of molecular techniques which allowed us to work with small amounts of degraded DNA extracted from old teeth, allowed us to exclude the hypothesis that genetic differentiation between zones 1 and 2 was owing to a barrier, which has appeared in the last 50 years. Furthermore, the low genetic variability and therefore the relatively high genetic differentiation shown by Danish badgers were also shown not to be the result of the recent drastic population decline. [ABSTRACT FROM AUTHOR]