Correction
23 Oct 2013: Pannebakker BA, Trivedi U, Blaxter ML, Watt R, Shuker DM (2013) Correction: The Transcriptomic Basis of Oviposition Behaviour in the Parasitoid Wasp Nasonia vitripennis. PLOS ONE 8(10): 10.1371/annotation/820ede3e-65d1-404a-9ddd-4db415ea760e. https://doi.org/10.1371/annotation/820ede3e-65d1-404a-9ddd-4db415ea760e View correction
Figures
Abstract
Linking behavioural phenotypes to their underlying genotypes is crucial for uncovering the mechanisms that underpin behaviour and for understanding the origins and maintenance of genetic variation in behaviour. Recently, interest has begun to focus on the transcriptome as a route for identifying genes and gene pathways associated with behaviour. For many behavioural traits studied at the phenotypic level, we have little or no idea of where to start searching for “candidate” genes: the transcriptome provides such a starting point. Here we consider transcriptomic changes associated with oviposition in the parasitoid wasp Nasonia vitripennis. Oviposition is a key behaviour for parasitoids, as females are faced with a variety of decisions that will impact offspring fitness. These include choosing between hosts of differing quality, as well as making decisions regarding clutch size and offspring sex ratio. We compared the whole-body transcriptomes of resting or ovipositing female Nasonia using a “DeepSAGE” gene expression approach on the Illumina sequencing platform. We identified 332 tags that were significantly differentially expressed between the two treatments, with 77% of the changes associated with greater expression in resting females. Oviposition therefore appears to focus gene expression away from a number of physiological processes, with gene ontologies suggesting that aspects of metabolism may be down-regulated during egg-laying. Nine of the most abundant differentially expressed tags showed greater expression in ovipositing females though, including the genes purity-of-essence (associated with behavioural phenotypes in Drosophila) and glucose dehydrogenase (GLD). The GLD protein has been implicated in sperm storage and release in Drosophila and so provides a possible candidate for the control of sex allocation by female Nasonia during oviposition. Oviposition in Nasonia therefore clearly modifies the transcriptome, providing a starting point for the genetic dissection of oviposition.
Citation: Pannebakker BA, Trivedi U, Blaxter MA, Watt R, Shuker DM (2013) The Transcriptomic Basis of Oviposition Behaviour in the Parasitoid Wasp Nasonia vitripennis. PLoS ONE 8(7): e68608. https://doi.org/10.1371/journal.pone.0068608
Editor: Hiroaki Matsunami, Duke University, United States of America
Received: April 29, 2013; Accepted: May 29, 2013; Published: July 19, 2013
Copyright: © 2013 Pannebakker et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The authors are very grateful to the Natural Environment Research Council (NERC: http://www.nerc.ac.uk/) for funding this work. BAP was also funded by the Netherlands Organization for Scientific Research (NWO: http://www.nwo.nl/. Veni grant no. 863.08.008) and the Netherlands Genomics Initiative (NGI: http://www.genomics.nl/. Horizon Breakthrough and Zenith grants no. 935.19.006 and 935.11.041). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
A key challenge facing biologists is to link the genotype to the phenotype [1]. For behaviour this is a particular challenge, as behaviours may be complex, highly environment-dependent phenotypes [2]–[4]. The rationale for trying to understand the genotype-phenotype map is either mechanistic (how are organisms made?) or functional (what genes generate the variation that evolution acts upon?). Traditional approaches to this problem have either worked from the bottom-up, exploring the phenotypic consequences of mutants to elucidate the pathways that generate phenotypes up the hierarchy of biological organisation (from cells, to tissues, to whole organisms), or instead they have worked from the top-down, seeking to describe whole organism variation in terms of genes or genomic regions. Both approaches have proved successful (for textbook treatments see [5]–[8]). However, it has become clear that, despite this success, we need to bridge the gap between these two approaches more completely if we are to map DNA sequences onto phenotypic variation, and thus link the molecular and phenotypic evolution of behaviour [9], [10].
Part of the problem has been that early optimism regarding the genetic basis of behaviour has been confounded by the number (and to some extent the identity) of genes associated with behavioural phenotypes [11]. In some fortunate cases, one or a few known genes do appear to provide a good understanding of the genetic basis of a given behaviour, and naturally-occurring variation in those genes has been identified. Perhaps two of the best examples involve the rover-sitter polymorphism in Drosophila melanogaster larvae that influences foraging strategies and length of foraging trails (associated with the foraging (for) gene, now known to be a cGMP-dependent protein kinase [12], and also the mechanism underlying circadian behavioural rhythms, again in Drosophila (associated with genes such as period and timeless [13]). However, these apparently simple single- or few-gene examples are rare and perhaps more complex than originally thought. It is becoming increasingly clear that many hundreds of sequences, coding and non-coding, may act in concert to influence behavioural phenotypes, and that individual mutations may vary in their phenotypic effects depending upon the genetic background (for instance outside of the standard laboratory genetic backgrounds of model organisms [11]). Indeed, despite the success of quantitative trait loci (QTL) studies in identifying genomic regions associated with many traits [14], including complex behaviour (such as sex allocation [15]), more recent work has confirmed that QTL or genome-wide association studies (GWAS) can give an unrealistically simple view of genetic architecture (for critiques see [9], [16]). This is exemplified by the “missing heritability” problem [17]. Even when a gene is found to be associated with a particular behaviour it may be a highly upstream gene with many pleiotropic effects (e.g. fruitless in Drosophila melanogaster [18]). Mutations of such genes may yield behavioural pathologies, but it remains an open question as to whether such loci, or their regulatory sequences, typically contribute to the segregating variation in natural populations of interest to evolutionary biologists.
Recent developments in microarray and sequencing technologies have offered an alternative approach that may help fill this gap: the large-scale analysis of the genes actually expressed in cells or tissues, or during key periods of development or behaviour [3], [4], [19]–[23]. These transcriptomic approaches offer a means of directly exploring which genes are associated with complex phenotypes such as behaviour, and experiments can reveal changes in gene regulation associated with changes in those phenotypes. For instance, recent work has revealed some of the genes and genetic networks associated with changes in female behaviour after mating in species of Drosophila [24]–[30], the Mediterranean fruit fly Ceratitis capitata [31], the honey bee Apis mellifera [32], [33], and the parasitoid wasp Nasonia vitripennis (R. Watt, U. Trivedi, T.J. Park, M. Blaxter and D.M. Shuker, unpublished data). These data may help identify physiological and neurological mechanisms associated with behaviours, and thus “candidate genes”, which is important as many of the behaviours studied by behavioural ecologists may be far removed from the often simpler behaviours that are studied in behaviour genetics laboratories. However, an animal’s behavioural repertoire is also influenced by the existing neural architecture and its physiology, such as titres of hormones and other signalling molecules, and so it remains an open question as to whether changes in gene expression contemporaneous to the performing of a behavioural act can meaningfully identify genes associated with that behaviour.
Here we consider whether a transcriptomic approach can reveal changes in gene expression during oviposition in the parasitoid wasp Nasonia vitripennis. Oviposition is a fundamental physiological process for females that need to invest limited resources into the production and maturation of eggs, which could otherwise have been invested in other somatic functions [34]. Furthermore, oviposition is a crucial reproductive decision that females of many species have to make. In terms of insects such as parasitoid wasps, oviposition is a key determinant of larval success and is the only form of parental care a female will provide for her offspring [35]. For Nasonia, females have to find and choose suitable hosts (blowfly pupae of the appropriate species, quality and developmental stage). They then have to decide how many eggs to lay (Nasonia are gregarious, laying multiple eggs per host puparium), and what sex ratio to produce. There is substantial evidence from Nasonia that females use an array of cues when finding, choosing and utilising hosts for egg-laying [36]–[40]. However, our understanding of the genetics underlying these important processes and decisions remains very limited [15], [34], [41], [42].
Importantly, a fully sequenced genome is available for Nasonia vitripennis and two of its con-geners (N. longicornis and N. giraulti), providing a starting point for genomic analyses of complex traits [43]. We have taken advantage of the availability of whole-genome information and associated gene annotations to explore what genes are expressed by female Nasonia vitripennis during oviposition. Insight into the genetic background of parasitoid oviposition will be crucial to understand how evolution acts upon a this fundamental physiological and behavoural trait. To the best of our knowledge, this is the first description of the oviposition transcriptome for any parastioid wasp.
Materials and Methods
Study organism
Nasonia vitripennis (Hymenoptera: Chalcidoidea) is a generalist gregarious wasp that parasitises large dipteran pupae (including Sarcophagidae and Calliphoridae [43], [44]). Depending on the host species, females oviposit between 20–50 eggs on an individual host, with male offspring emerging just before females (after approximately 14 days at 25°C). Females are synovigenic, i.e. they are born with a limited number of mature eggs in their ovarioles, but can produce and mature further eggs provided they have protein resources available [45]. Males have small wings and are unable to fly, remaining close to the emergence patch where they compete with each other for matings with emerging females. Females are fully winged and typically mate only once before dispersing to find new hosts.
The wasps used in this experiment were from the AsymC strain, which was isolated in 1986 by curing the wild-type strain LabII of Wolbachia [46]. Wasps have been maintained on Calliphora vomitoria or C. vicina hosts at 25°C, 16L: 8D light conditions since this time, allowing AsymC to become highly inbred. The AsymC line was sequenced and annotated in the Nasonia genome project [43], allowing for the direct mapping of our transcriptomic data on to the available genomic resources.
Oviposition experiment
In order to control for possible host and genotype effects, we isolated a single 2-day old mated AsymC female in a glass vial and provided her with a single host to produce F1 daughters. Eight 2-day old mated F1 females were subsequently provided with three hosts to produce the F2 test females. We randomly selected one host from each F1 female, and isolated 16 2-day old mated F2 test females in glass tubes, of which eight were randomly allocated to the oviposition treatment and eight to the resting treatment. We provided the test females with a single host for 24 hours as pre-treatment to facilitate egg development [45]. We then discarded the pre-treatment hosts and gave each female a piece of chromatography paper soaked in honey solution for a further 24 hours.
For the experiment, we transferred the females to 1.5 mL Eppendorf tubes that contained a single host for the oviposition experiment, or were empty for the resting treatment. After 60 mins, females were flash-frozen in liquid nitrogen and stored on dry-ice until the addition of RNAlater-ICE (Ambion, Austin, TX, USA) after which they were transferred to −20oC. All females in the oviposition treatment were observed to have commenced ovipositing. We pooled the F2 test females from each F1 mother according to treatment, generating a total of eight pooled samples per treatment (consisting of 8 females per pool) for RNA isolation and sequencing (i.e. N = 16 for the whole experiment).
RNA isolation
RNA was isolated using Qiagen RNeasy Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s protocol. We further purified our samples using Turbo DNase (Ambion, Austin, TX, USA) and tested the integrity and concentration of our samples using a NanoDrop spectrophotometer (Nanodrop Technologies, Wilmington, DE, USA).
Digital gene expression tag profiling (DGE)
Tag library preparation was done using the DGE-Tag Profiling for NlaIII Sample Prep kit (Illumina, San Diego, CA, USA) according to the manufacturer’s protocol. In short, mRNA was first captured from the total RNA by magnetic oligo(dT)beads. mRNA was then converted into cDNA and bead-bound cDNA was subsequently digested with NlaIII. Fragments other than the 3' cDNA fragments attached the beads were washed away, and a NlaIII-adapter containing a MmeI recognition site was ligated to the 5' end of the bead-bound cDNA. MmeI cuts to 21bp downstream from its recognition site, creating 21bp tags that start with the NlaIII recognition site, CATG. A second adapter was then ligated at the MmeI site, allowing for PCR amplification, followed by attachment of the enriched tags to the surface of the sequencing flowcell in which they were sequenced by synthesis using the Illumina Genome Analyzer I system at GenePool, the University of Edinburgh's in-house sequencing facility. The raw data (tag sequences and counts) were deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE43352 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE43352).
DGE tag mapping and filtering
Reads were first aligned against the Nasonia vitripennis genome (Nvit1.0, ftp://ftp.hgsc.bcm.tmc.edu/pub/data/Nvitripennis/fasta/Nvit_1.0/linearized_sequence) using MAQ 0.6.8 [47], allowing for a 2bp mismatch between the 21bp tag and reference sequence. Only tags with a phred mapping quality threshold of 30 were used in further analysis, corresponding to one wrongly aligned sequence in 1000. We then aligned the remaining tags to the Nasonia Official Gene Set (OGS) v.1.2 (http://hymenopteragenome.org/Nasonia/) using the same mapping settings. Tags that failed to map to the Nasonia genome or the OGS were classed as null and not considered in further analysis. Whilst DGE tag profiling produces both sense and antisense tags, the majority of antisense tags produced in our set-up are most likely the result of technical artefacts [48]. We therefore deemed it appropriate to remove all antisense tags from our analysis. Transcripts present at extremely low abundance are a common source of noise in these kinds of studies and likely include technical artefacts. We therefore also excluded tags that had 15 or less reads across all 16 samples (i.e. a mean of less than one read per replicate). In terms of the dataset, we initially had 74,173,884 reads comprising 215,130 tags. Following removal of anti-sense and null sequences we had 32,036,887 reads from 88,627 tags (Table S1). Removing the very low abundance transcripts left us with 31,807,132 reads from 30,334 tags for our analysis (i.e. the last step removed 65.8% of tags but only 0.7% of reads).
Statistical analysis
We tested for differential expression between ovipositing and resting females using DESeq v2.10 [49], run with R version 2.15.0 [50]. Briefly, DESeq employs a negative binomial error distribution to model transcript abundance in high-throughput datasets. With the comparatively high level of replication (for this kind of study) used in this experiment (N = 8 replicates per treatment), we used empirical estimates of the deviance of tag counts for each tag (using the option “gene-est-only”). DESeq estimates the significance of differential expression for each tag, and then corrects for the multiple comparisons using the false discovery rate (FDR) correction of Benjamini & Hochberg [51]. We considered tags to be differentially expressed if they had an adjusted p-value after FDR correction of p<0.1. We also report which of these tags that have greater than 4-fold log2 gene expression differences (another common measure of differential expression).
Annotation and Gene Ontology
We annotated our tags using the Nasonia Official Gene Set version 2 (OGS2, January 2012, http://arthropods.eugenes.org/genes2/Nasonia/). Briefly, annotation for the tags aligned to the genome was done on the basis of their genomic position, with tags that were first aligned to the OGS1.2 converted to OGS2 based on gene equivalence. In total 25,427 out of 30,344 (83.8%) tags were annotated with the OGS2. Each tag was then assigned gene ontology (GO) terms using the Nasonia GO annotation for the OGS2 generated as part of an effort by the Nasonia community (http://hymenopteragenome.org/Nasonia/?q=evidential_gene_data). Singular enrichment analysis of the GO terms was performed in agriGO (http://bioinfo.cau.edu.cn/agriGO [52]) using chi-square tests to compare differentially expressed genes to the Nasonia GO annotation.
Results
Oviposition in N. vitripennis females is associated with changes in gene expression at the whole-body level. From a total of 30,334 tags, 322 showed significant differential expression at a FDR of 0.1 (unadjusted p-values: all p<0.0012), with 209 significant at a FDR of 0.05 (Table S2). These differentially expressed tags represent 84,326 reads out of a total of 31,807,132, or 0.27% of the total transcriptome. Oviposition appears to be associated with a focusing of gene expression: 73 tags (22.7%) show significantly greater expression in ovipositing females, whilst 249 (77.3%) show greater expression in resting females (Figure 1). In addition, 9 of the 12 differentially expressed tags with total counts across all replicates in excess of 1000 are associated with greater expression in ovipositing females (see also Table 1). Forty-four tags (13.7% of our differentially expressed tags) had a greater than 4-fold log2 change in expression between treatments (Table 2). Six of these tags showed higher expression in ovipositing females, whilst the other 38 tags showed higher expression in resting females.
Total transcript abundance by treatment for the 322 tags differentially expressed between resting and ovipositing Nasonia vitripennis females (summed across N = 8 biological replicates per treatment). Approximately 77% of these tags show higher expression levels in resting females. The dotted line indicates a 1:1 relationship. All tags with total abundance across both treatments of <16 counts were excluded from the analysis.
The most abundant tag of all mapped to the gene elongation factor 1-alpha 1 (with an average abundance across all replicates of approximately 60080 reads) which was reassuringly not differentially expressed between our treatments (this gene, with its universal role in translation, is typically used a control for quantitative-PCR studies of differential expression of individual loci for this reason). In total, 17 tags had average abundances in excess of 10,000 reads per replicate. Our most abundant differentially expressed tag mapped to the glucose dehydrogenase gene, and had an average abundance of 787 reads per replicate. Forty differentially expressed tags were unique to one treatment or another (32 unique to resting females and 8 unique to ovipositing females), but all of these tags had total counts between 16 (our minimum threshold abundance) and 37.
From the total tag pool of 30,334 tags, 25,427 (83.8%) were annotated, covering 9,182 unique genes in the Nasonia OGS2 (37.4% out of total 24,525 genes in OGS2). 4,907 tags (16.2%) did not match any genes in the OGS2. Of the 322 tags showing differentially expression, 286 (88.8%) were annotated, while 36 (11.2%) did not match any OGS2 genes. Of the annotated differentially expressed tags, 163 (57.0%) could be assigned a Gene Ontology (GO) classification representing 38 unique GO terms for the up-regulated tags, and 134 unique GO terms for the down-regulated tags (Figures S1 and S2). Singular Enrichment Analysis of the associated GO terms revealed an overrepresentation of 5 GO terms in the down-regulated tags, all of which related to metabolic processes (Table 3).
Discussion
Our results show that oviposition is associated with changes in whole-body patterns of gene expression in female Nasonia vitripennis wasps. Approximately three-quarters of the changes involve greater expression in resting females (or put another way, down-regulation in ovipositing females), and our enrichment analysis suggests that these down-regulated genes are more likely to be involved in various metabolic processes than expected by chance (Table 3). We therefore suggest that as female wasps move from resting to ovipositing, aspects of their metabolism are down-regulated, focusing gene expression on other processes. Similar gene expression focusing has been reported in Drosophila responding to environmental stress, where metabolism related genes are down-regulated and only few stimulus specific stress genes are up-regulated [53]. However, from the outset we note that the changes in gene expression we have uncovered represent less than half of one percent of the total (whole-body) transcriptome, so that at least at this level of resolution many genes are unaffected by the switch in behaviour. At such an early stage in the development of the field of behavioural transcriptomics it is unclear how common a finding this will be. Of course, many physiological functions have to carry on regardless of what an animal is doing, and all of our most highly expressed tags were not differentially expressed (indeed, the top 362 tags in terms of average counts were not differentially expressed).
Alongside the key result that oviposition is associated with broad-scale changes in gene expression, with metabolic processes affected more than expected by chance, to what extent have we been able to identify likely candidate genes that influence oviposition behaviour? Nine of our most abundant differentially expressed tags showed greater levels of expression in ovipositing females (Table 1). These included the gene purity-of-essence (poe), which is also known as pushover (push) in Drosophila melanogaster [54]. The purity-of-essence protein is an evolutionarily conserved, large membrane protein containing two zinc finger domains (a structural protein motif associated with binding DNA, RNA, or proteins [55]) that influences behaviour and synaptic transmission in D. melanogaster. Two mutants of poe/push in Drosophila both cause increased nervous excitability and reduced motor function [54] and mutants also influence peripheral nerve morphology [55]. The gene potentially has other functions (perhaps associated with synaptic transmission), with mutants being associated with sterility in male Drosophila [54], and an incomplete push/poe protein has also been identified as a calmodulin-binding protein expressed in fly photoreceptors [56]. This is perhaps our best candidate for a gene that might influence the motor control associated with handling hosts and/or the resulting drilling into the host puparium, followed by the mechanical control of egg laying.
Also in this group of up-regulated genes is glucose dehydrogenase, the protein of which is associated with successful sperm uptake and release in Drosophila melanogaster [57]. This was the most abundant gene that was differentially expressed in our experiment. Whilst in our experiment the resting females were mated and so also had sperm present in their spermatheca, our finding is suggestive that glucose dehydrogenase is needed during oviposition for the successful release of sperm. As female Nasonia need to control sperm release very precisely in order to allocate sex (with fertilised ova developing into daughters, and unfertilised ova developing into sons), variation in the expression of the glucose dehydrogenase gene may contribute to phenotypic variation in sex ratio and constrain a female’s ability to adaptively allocate sex as predicted by Hamilton’s theory of Local Mate Competition (LMC [58]; see also [59]–[61]). As the focus of much of the behavioural work associated with oviposition in Nasonia is associated with sex allocation, this will be an extremely interesting gene to consider in more detail, and we are currently exploring changes in gene expression in Nasonia vitripennis females when we experimentally manipulate the cues females use to make their sex allocation decisions.
Another encouraging candidate is found amongst the 44 of our tags that exhibited greater than 4-fold log2 differential expression (13.7% of our DE tags). Six of these showed greater levels of expression in ovipositing females, including a neurocalcin homolog. The neurocalcin protein is a neuronal calcium-binding protein that may be involved in neurotransmitter release and the regulation of neural function [62]. In vertebrates, neurocalcin may be important for regulating sexual dimorphism of the neural song system in birds [63] and sexually dimorphic patterns of expression have also been shown in rodents [64]. Whilst the possible functions of neurocalcin-like proteins are poorly known in insects, a neurocalcin is expressed in Drosophila neurones [65] and as such this is again the kind of gene we would expect to show changes in the pattern of expression during a behaviour in which females had to either neurally process information and/or perform a complex set of motor patterns.
Otherwise, the links to behaviour are less immediate. For instance, we found that the genes Scavenger receptor class B member 1 (SR-B1) and acyl-CoA delta-11-desaturase were up-regulated in ovipositing females. These genes code for proteins associated in lipid transport and fatty acid metabolism [66], [67], and so they may influence how lipids are mobilised during oviposition either in terms of releasing energy reserves for the behaviours themselves or in terms of energy reserves needed for future oogenesis to replace the eggs being deposited. The gene bcl-2-related ovarian killer (BOK) was also up-regulated during oviposition. BOK homologues regulate programmed cell death (apoptosis) during Drosophila oogenesis [68], which perhaps suggests that apoptosis is occurring during oviposition in Nasonia. Other up-regulated genes include a ubiquinone biosynthesis gene, which again is associated with energy production via the electron transport chain [69]. All these latter cases are perhaps better interpreted as changes that result from the physiological mechanisms involved in oviposition (for instance in terms of the organismal and life-history changes associated with oviposition and the energy utilised during egg-laying), rather than changes that initiate the oviposition behaviours themselves. This fits with the predominant pattern in our data mentioned above, namely the reduction in expression of genes associated with a range of metabolic processes when females are ovipositing. As such, we might be picking up more of the genetic basis of the physiological and life-history consequences of oviposition, and in particular the change in energy-use and metabolism associated with reproductive behaviour, than we are picking up the genetic signals of the behaviour itself. In the very simple experiment undertaken here, disentangling cause and effect for genes and behaviour is clearly difficult, especially with our limited understanding of how cellular processes interact with whole-organism behaviour.
In fact, one of the crucial limitations of behavioural transcriptomics at the moment is that our annotations of genes are associated with molecular functions at a cellular or tissue level that may seem a long way from the regulation or control of a focal behaviour (unless we hit a key developmental gene or signalling pathway). This limitation has recently been articulated by Pavey et al. [70], who suggest that new gene ontologies with explicit ecological (or in our case behavioural) functions might be a useful way to generalise the patterns of associations between molecular and whole-organism phenotypes. For example, it might be that pathways associated with the control of nutrition (such as the insulin-like signalling pathway) end up being implicated in numerous behaviours and life-history decisions because of the fundamental importance of energy acquisition and allocation. These patterns would be easier to recognise if we had more explicit ecological or behavioural designations for genes, built up through studies such as the one we present here. That said, even this will only be a start, as functional analyses in which genes are knocked-down will remain crucial to real progress. Fortunately, techniques such as RNA-interference (RNAi), whereby gene transcripts of focal genes can be targeted and rendered non-functional through cleavage (i.e. transcriptional silencing), are becomingly increasingly available in non-model organisms, including Nasonia, crickets, beetles and bugs (e.g. [71]–[77] ).
To conclude, oviposition modifies the pattern of gene expression in female Nasonia vitripennis, suggestive of a down-regulation of aspects of metabolism during oviposition behaviour. At the moment we know rather little about the genetic basis of oviposition behaviour in Nasonia, or indeed in any other insect. What we do know so far is mostly in terms of genetic variation in oviposition decisions, both in Nasonia and elsewhere (e.g. [15], [42], [78], [79] ). Here we have been able to begin the process of identifying putative candidate genes associated with oviposition, including identifying a protein known to influence sperm usage in flies (glucose dehydrogenase). However, more generally we hope that our results encourage more animal behaviour researchers to begin to explore how patterns of gene expression are associated with their own behaviours of interest, using the sequencing technologies now readily available. Indeed, currently perhaps the most valuable aspect of studies such as ours is in terms of generating hypotheses about the mechanisms underlying the behaviour we study, hypotheses that can and should be tested comparatively across as many species as possible.
Supporting Information
Figure S1.
Proportion of Gene Ontology (GO) terms for biological process (a), molecular function (b) and cellular component (c) in up-regulated tag set (white bars) and in Official Gene Set 2 (black bars) in ovipositing versus resting Nasonia vitripennis females.
https://doi.org/10.1371/journal.pone.0068608.s001
(PDF)
Figure S2.
Proportion of Gene Ontology (GO) terms for biological process (a), molecular function (b) and cellular component (c) in down-regulated tag set (white bars) and in Official Gene Set 2 (black bars) in ovipositing versus resting Nasonia vitripennis females. Significantly enriched terms are indicated with an asterisk.
https://doi.org/10.1371/journal.pone.0068608.s002
(PDF)
Table S1.
Valid tags used for analysis of digital gene expression between ovipositing and resting Nasonia vitripennis females. Table shows tag identifier (TagID), location of tag in N. vitripennis genome (Nvit 1.0) or Official Gene Set v1.2 (Location) and tag sequence (TagSequence).
https://doi.org/10.1371/journal.pone.0068608.s003
(XLSX)
Table S2.
Differentially expressed tags in ovipositing versus resting Nasonia vitripennis females. Table shows tag identifier (TagID), tag sequence (TagSequence), mean counts, averaged over all samples from both conditions (baseMean), mean counts from ovipositing females (baseMeanA), mean counts from resting females (baseMeanB), fold change from ovipositing versus resting females (foldChange), the logarithm (to basis 2) of the fold change from ovipositing versus resting females (log2FoldChange), P-value for the statistical signicance of the change in expression (P.value) and P-value adjusted for multiple testing with the Benjamini-Hochberg procedure (FDR.adjusted.P.value).
https://doi.org/10.1371/journal.pone.0068608.s004
(XLSX)
Acknowledgments
Our thinking about behavioural transcriptomics has been influenced by many in the community, but we are particularly grateful to Mike Ritchie, John Colbourne and Jack Werren for discussions and inspiration. We also would like to thank Dave Wheeler for early access to the Nasonia annotation data files, Nicki Cook for comments on an early version of the manuscript, and Bertanne Visser for feedback on Nasonia metabolic pathways.
Author Contributions
Conceived and designed the experiments: BAP MAB DMS. Performed the experiments: BAP RW DMS. Analyzed the data: BAP UT DMS. Contributed reagents/materials/analysis tools: BAP UT MAB DMS. Wrote the paper: BAP DMS.
References
- 1. Houle D (2010) Colloquium papers: Numbering the hairs on our heads: the shared challenge and promise of phenomics. Proc Natl Acad Sci USA 107 Suppl: 1793–1799
- 2. Sokolowski MB (2001) Drosophila: genetics meets behaviour. Nat Rev Genet 2: 879–890.
- 3. Boake CRB, Arnold SJ, Breden F, Meffert LM, Ritchie MG, et al. (2002) Genetic tools for studying adaptation and the evolution of behavior. Am Nat 160 Suppl: S143–159
- 4. Anholt RRH, Mackay TFC (2004) Quantitative genetic analyses of complex behaviours in Drosophila. Nat Rev Genet 5: 838–849.
- 5.
Falconer DS, Mackay TFC (1996) Introduction to Quantitative Genetics. Fourth edi. Longmans Green, Harlow, Essex, UK.
- 6.
Lynch M, Walsh B (1998) Genetics and Analysis of Quantitative Traits. Sinauer Associates, Sunderland.
- 7.
Anholt RRH, Mackay TFC (2009) Principles of Behavioral Genetics. London UK: Academic Press.
- 8.
Griffiths AJF, Wessler SR, Carroll S, Doebley JF (2011) Introduction to Genetic Analysis. New York: W.H. Freeman & Co.
- 9. Erickson DL, Fenster CB, Stenøien HK, Price D (2004) Quantitative trait locus analyses and the study of evolutionary process. Mol Ecol 13: 2505–2522.
- 10. Springer SA, Crespi BJ, Swanson WJ (2011) Beyond the phenotypic gambit: molecular behavioural ecology and the evolution of genetic architecture. Mol Ecol 20: 2240–2257.
- 11.
Flint J, Greenspan RJ, Kendler KS (2010) How genes influence behavior. Oxford UK: Oxford University Press.
- 12. Osborne KA, Robichon A, Burgess E, Butland S, Shaw R, et al. (1997) Natural behavior polymorphism due to a cGMP-dependent protein kinase of Drosophila. Science 277: 834–836.
- 13. Williams JA, Sehgal A (2001) Molecular components of the circadian system in Drosophila. Annu Rev Physiol 63: 729–755.
- 14. Mackay TFC (2001) The genetic architecture of quantitative traits. Annu Rev Genet 35: 303–339.
- 15. Pannebakker BA, Watt R, Knott SA, West SA, Shuker DM (2011) The quantitative genetic basis of sex ratio variation in Nasonia vitripennis: a QTL study. J Evol Biol 24: 12–22.
- 16. Rockman MV (2012) The QTN program and the alleles that matter for evolution: all that’s gold does not glitter. Evolution 66: 1–17.
- 17. Yang J, Benyamin B, McEvoy BP, Gordon S, Henders AK, et al. (2010) Common SNPs explain a large proportion of the heritability for human height. Nat Genet 42: 565–569.
- 18. Greenspan RJ, Ferveur JF (2000) Courtship in Drosophila. Annu Rev Genet 34: 205–232.
- 19. Gibson G (2002) Microarrays in ecology and evolution: a preview. Mol Ecol11: 17–24.
- 20. Ranz JM, Machado CA (2006) Uncovering evolutionary patterns of gene expression using microarrays. Trends Ecol Evol 21: 29–37.
- 21. Hoheisel JD (2006) Microarray technology: beyond transcript profiling and genotype analysis. Nat Rev Genet 7: 200–210.
- 22. Rokas A, Abbot P (2009) Harnessing genomics for evolutionary insights. Trends Ecol Evol 24: 192–200.
- 23. Stapley J, Reger J, Feulner PGD, Smadja C, Galindo J, et al. (2010) Adaptation Genomics: the next generation. Trends Ecol Evol 25: 705–712.
- 24. Lawniczak MKN, Begun DJ (2004) A genome-wide analysis of courting and mating responses in Drosophila melanogaster females. Genome 47: 900–910.
- 25. McGraw LA, Gibson G, Clark AG, Wolfner MF, Carolina N (2004) Genes regulated by mating, sperm, or seminal proteins in mated female Drosophila melanogaster. Curr Biol 14: 1509–1514.
- 26. McGraw LA, Clark AG, Wolfner MF (2008) Post-mating gene expression profiles of female Drosophila melanogaster in response to time and to four male accessory gland proteins. Genetics 179: 1395–1408.
- 27. Mack PD, Kapelnikov A, Heifetz Y, Bender M (2006) Mating-responsive genes in reproductive tissues of female Drosophila melanogaster. Proc Natl Acad Sci USA 103: 10358–10363.
- 28. Innocenti P, Morrow EH (2009) Immunogenic males: a genome-wide analysis of reproduction and the cost of mating in Drosophila melanogaster females. J Evol Biol 22: 964–973.
- 29. Bono JM, Matzkin LM, Kelleher ES, Markow TA (2011) Postmating transcriptional changes in reproductive tracts of con- and heterospecifically mated Drosophila mojavensis females. Proc Natl Acad Sci USA 108: 7878–7883.
- 30. Immonen E, Ritchie MG (2012) The genomic response to courtship song stimulation in female Drosophila melanogaster. Proc Biol Sci. 279 1359–1365.
- 31. Gomulski LM, Dimopoulos G, Xi Z, Scolari F, Gabrieli P, et al. (2012) Transcriptome profiling of sexual maturation and mating in the Mediterranean fruit fly, Ceratitis capitata. PloS one 7: e30857.
- 32. Kocher SD, Ayroles JF, Stone EA, Grozinger CM (2010) Individual variation in pheromone response correlates with reproductive traits and brain gene expression in worker honey bees. PloS one 5: e9116.
- 33. Kocher SD, Richard FJ, Tarpy DR, Grozinger CM (2008) Genomic analysis of post-mating changes in the honey bee queen (Apis mellifera). BMC Genomics 9: 232.
- 34. Jervis M, Ellers J, Harvey JA (2008) Resource acquisition, allocation, and utilization in parasitoid reproductive strategies. Annu Rev Entomol. 53: 361–385.
- 35.
Godfray HCJ (1994) Parasitoids: behavioral and evolutionary ecology. Princeton University Press.
- 36. Werren JH (1983) Brood size and sex ratio regulation in the parasitic wasp Nasonia vitripennis (Walker) (Hymenoptera: Pteromalidae). Neth J Zool 34: 123–143.
- 37. Flanagan KE, West SA, Godfray HCJ (1998) Local mate competition, variable fecundity and information use in a parasitoid. Anim Behav 56: 191–198.
- 38. Shuker DM, West SA (2004) Information constraints and the precision of adaptation: sex ratio manipulation in wasps. Proc Natl Acad Sci USA 101: 10363–10367.
- 39. Rivero A, West SA (2005) The costs and benefits of host feeding in parasitoids. Anim Behav 69: 1293–1301.
- 40. Schurmann D, Collatz J, Hagenbucher S, Ruther J, Steidle J (2009) Olfactory host finding, intermediate memory and its potential ecological adaptation in Nasonia vitripennis. Naturwissenschaften 96: 383–391.
- 41. Pannebakker BA, Halligan DL, Reynolds KT, Shuker DM, Ballantyne GA, et al. (2008) Effects of spontaneous mutation accumulation on sex ratio traits in a parasitoid wasp. Evolution 62: 1921–1935.
- 42. Desjardins CA, Perfectti F, Bartos JD, Enders LS, Werren JH (2010) The genetic basis of interspecies host preference differences in the model parasitoid Nasonia. Heredity 104: 270–277.
- 43. Werren JH, Richards S, Desjardins CA, Niehuis O, Gadau J, et al. (2010) Functional and evolutionary insights from the genomes of three parasitoid Nasonia species. Science 327: 343–348.
- 44. Whiting AR (1967) The biology of the parasitic wasp Mormoniella vitripennis (Walker). Q Rev Biol 42: 333–406.
- 45. Rivero A, West SA (2005) The costs and benefits of host feeding in parasitoids. Anim Behav 69: 1293–1301.
- 46. Perrot-Minnot MJ, Werren JH (1999) Wolbachia infection and incompatibility dynamics in experimental selection lines. J Evol Biol 12: 272–282.
- 47. Li H, Ruan J, Durbin R (2008) Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res 18: 1851–1858.
- 48.
Watt R (2012) The Transcriptomic Basis of Post-Mating Responses in Females of the Parasitic Wasp Nasonia vitripennis University of Edinburgh.
- 49. Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biol 11: R106.
- 50.
R Development Core Team (2012) R: A language and enviroment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0. Available: http://www.R-project.org/.
- 51. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc B 57: 289–300.
- 52. Du Z, Zhou X, Ling Y, Zhang Z, Su Z (2010) agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res 38: W64–70.
- 53. Sørensen JG, Nielsen MM, Kruhøffer M, Justesen J, Loeschcke V (2005) Full genome gene expression analysis of the heat stress response in Drosophila melanogaster. Cell Stress Chaperon 10: 312–328.
- 54. Richards S, Hillman T, Stern M (1996) Mutations in the Drosophila pushover gene confer increased neuronal excitability and spontaneous synaptic vesicle fusion. Genetics 142: 1215–1223.
- 55. Yager J, Richards S, Hekmat-Scafe DS, Hurd DD, Sundaresan V, et al. (2001) Control of Drosophila perineurial glial growth by interacting neurotransmitter-mediated signaling pathways. Proc Natl Acad Sci USA 98: 10445–10450.
- 56. Xu XZS, Wes PD, Chen H, Li HS, Yu M, et al. (1998) Retinal targets for calmodulin include proteins implicated in synaptic transmission. J Biol Chem 273: 31297–31307.
- 57. Iida K, Cavener D (2004) Glucose dehydrogenase is required for normal sperm storage and utilization in female Drosophila melanogaster. J Exp Biol 207: 675–681.
- 58. Hamilton WD (1967) Extraordinary sex ratios. Science 156: 477–488.
- 59. Werren JH (1980) Sex ratio adaptations to local mate competition in a parasitic wasp. Science 208: 1157–1159.
- 60. Werren JH (1983) Sex ratio evolution under local mate competition in a parasitic wasp. Evolution 37: 116–124.
- 61. Shuker DM, Pen I, Duncan AB, Reece SE, West SA (2005) Sex ratios under asymmetrical local mate competition: theory and a test with parasitoid wasps. Am Nat 166: 301–316.
- 62. Burgoyne RD (2007) Neuronal calcium sensor proteins: generating diversity in neuronal Ca2+ signalling. Nat Rev Neurosci 8: 182–193.
- 63. Long P, Veney SL (2010) Neurocalcin protein labeling reveals a dimorphism within the developing zebra finch brain. Brain Res 1315: 11–18.
- 64. Murias AR, Weruaga E, Recio JS, Airado C, Díaz D, et al. (2007) Distribution of neurocalcin-containing neurons reveals sexual dimorphism in the mouse olfactory bulb. Chem Senses 32: 673–680.
- 65. Faurobert E, Chen CK, Hurley JB, Teng DHF (1996) Drosophila neurocalcin, a fatty acylated, Ca2+-binding protein that associates with membranes and inhibits in vitro phosphorylation of bovine rhodopsin. J Biol Chem 271: 10256–10262.
- 66. Acton S, Rigotti A, Landschulz K, Xu S, Hobbs H, et al. (1996) Identification of scavenger receptor SR-BI as a high density lipoprotein receptor. Science 271: 460–461.
- 67. Behrouzian B, Buist PH (2002) Fatty acid desaturation: variations on an oxidative theme. Curr Opin Chem Biol 6: 577–582.
- 68. Tanner EA, Blute TA, Brachmann CB, McCall K (2011) Bcl-2 proteins and autophagy regulate mitochondrial dynamics during programmed cell death in the Drosophila ovary. Development 138: 327–338.
- 69. Jonassen T, Proft M, Randez-Gil F, Schultz JR, Marbois BN, et al. (1998) Yeast Clk-1 homologue (Coq7/Cat5) is a mitochondrial protein in coenzyme Q synthesis. J Biol Chem 273: 3351–3357.
- 70. Pavey SA, Bernatchez L, Aubin-Horth N, Landry CR (2012) What is needed for next-generation ecological and evolutionary genomics? Trends Ecol Evol 27: 673–678.
- 71. Lynch JA, Brent AE, Leaf DS, Pultz MA, Desplan C (2006) Localized maternal orthodenticle patterns anterior and posterior in the long germ wasp Nasonia. Nature 439: 728–732.
- 72. Lynch JAJA, Desplan C (2006) A method for parental RNA interference in the wasp Nasonia vitripennis. Nat Protoc 1: 486–494.
- 73. Hrycaj S, Mihajlovic M, Mahfooz N, Couso JP, Popadić A (2008) RNAi analysis of nubbin embryonic functions in a hemimetabolous insect, Oncopeltus fasciatus. Evol Dev 10: 705–716.
- 74. Marshall JL, Huestis DL, Hiromasa Y, Wheeler S, Oppert C, et al. (2009) Identification, RNAi knockdown, and functional analysis of an ejaculate protein that mediates a postmating, prezygotic phenotype in a cricket. PloS one 4: e7537.
- 75. Aspiras AC, Smith FW, Angelini DR (2011) Sex-specific gene interactions in the patterning of insect genitalia. Dev Biol 360: 369–380.
- 76. Zwier MV, Verhulst EC, Zwahlen RD, Beukeboom LW, Van de Zande L (2012) DNA methylation plays a crucial role during early Nasonia development. Insect Mol Biol 21: 129–138.
- 77. Konuma T, Morooka N, Nagasawa H, Nagata S (2012) Knockdown of the adipokinetic hormone receptor increases feeding frequency in the two-spotted cricket Gryllus bimaculatus. Endocrinology 153: 3111–3122.
- 78. Nylin S, Nygren G, Windig J, Janz N, Bergstrom A (2005) Genetics of host-plant preference in the comma butterfly Polygonia c-album (Nymphalidae), and evolutionary implications. Biol J Linn Soc Lond 84: 755–765.
- 79. Forister ML, Ehmer AG, Futuyma DJ (2007) The genetic architecture of a niche: variation and covariation in host use traits in the Colorado potato beetle. J Evol Biol 20: 985–996.