Research article | Open | Published:
Nourishment level affects caste-related gene expression in Polistes wasps
BMC Genomicsvolume 16, Article number: 235 (2015)
Social insects exhibit striking phenotypic plasticity in the form of distinct reproductive (queen) and non-reproductive (worker) castes, which are typically driven by differences in the environment during early development. Nutritional environment and nourishment during development has been shown to be broadly associated with caste determination across social insect taxa such as bees, wasps, and termites. In primitively social insects such as Polistes paper wasps, caste remains flexible throughout adulthood, but there is evidence that nourishment inequalities can bias caste development with workers receiving limited nourishment compared to queens. Dominance and vibrational signaling are behaviors that have also been linked to caste differences in paper wasps, suggesting that a combination of nourishment and social factors may drive caste determination. To better understand the molecular basis of nutritional effects on caste determination, we used RNA-sequencing to investigate the gene expression changes in response to proteinaceous nourishment deprivation in Polistes metricus larvae.
We identified 285 nourishment-responsive transcripts, many of which are related to lipid metabolism and oxidation-reduction activity. Via comparisons to previously identified caste-related genes, we found that nourishment restriction only partially biased wasp gene expression patterns toward worker caste-like traits, which supports the notion that nourishment, in conjunction with social environment, is a determinant of developmental caste bias. In addition, we conducted cross-species comparisons of nourishment-responsive genes, and uncovered largely lineage-specific gene expression changes, suggesting few shared nourishment-responsive genes across taxa.
Overall, the results from this study highlight the complex and multifactorial nature of environmental effects on the gene expression patterns underlying plastic phenotypes.
Phenotypic plasticity provides an important adaptive mechanism by which morphology, physiology, and/or behavior can be adjusted to biotic and abiotic environmental factors including temperature, nutrition, population density, and predator presence . There are a number of striking examples of phenotypic plasticity in insects: dimorphic horn development in dung beetles , seasonal color polyphenism in butterflies , and wing polyphenism in aphids (reviewed in ). One of the best-studied models of insect phenotypic plasticity is reproductive castes in social insects, especially social Hymenoptera (bees, ants, and wasps ). In most social insects, genotypic differences do not account for differences between reproductive (queen) and non-reproductive (worker) castes (reviewed in ). Instead, environmental factors induce differences in hormone titers and gene expression (reviewed in ), leading to the development of queens or workers, which vary in physiology and behavior and for advanced social species, in morphology .
One particular environmental factor, food availability, is especially important for caste polyphenism in social insects. Differential nourishment [7-10] and nutrition-related genes and pathways (e.g. storage proteins such as vitellogenin and hexamerin, insulin/insulin-like signaling (IIS) pathways) have been linked to caste differences in honey bees [11,12], paper wasps [13,14], and termites [15,16]. This suggests that the influence of nutrition on caste formation may be broadly shared across diverse taxa . There have been numerous studies investigating the molecular mechanisms underlying queen-worker caste determination in advanced eusocial species, especially honey bees [18-25]. In honey bees, nutritional differences for larvae fed either royal jelly or worker jelly precede a developmental switch resulting in alternative caste phenotypes . Nutrition is also important for caste differences in primitively eusocial species, which lack morphological castes, but we know much less about the molecular mechanisms that underlie the formation of their more subtle behavioral and physiological castes. In primitively eusocial species such as paper wasps, differential nourishment does not strictly determine caste but can lead to a caste bias, whereby female larvae that are fed larger quantities of food are more likely to be future reproductive queens (called “gynes”) as adults . However, the ultimate caste fate of a female is decided during adulthood. First-brood offspring of an established nest are capable of independent reproduction, but instead they perform allomaternal care (worker behavior) as a response to cues emitted by larvae in the nest . Subsequent social reinforcement of worker behavior often occurs via dominance behaviors by the queen or other workers [29,30].
Primitively eusocial taxa such as Polistes are an informative group for understanding the evolution of eusociality and the origins of castes [31,32]. Across the annual colony cycle of primitively social wasps in the genus Polistes¸ the quantity of larval nourishment changes according to seasonal changes in the adult-to-larva ratio [10,27]. First-brood offspring produced early in the colony cycle have been reared by a single nest-founding queen or few queens, which also perform all colony tasks such as nest building, foraging, egg laying, and brood rearing . At this early time in the colony cycle the adult-to-larva ratio is low [30,33], which leads to low feeding rates and more limited larval nourishment compared to offspring reared by workers and produced later in the colony cycle (future queens or “gynes”). Physiological evidence of nourishment-related differences between workers and gynes collected from naturally-founded colonies in the field include greater fat body stores in gynes [34,35], greater quantities of the storage protein hexamerin 1 in gynes [36,37], and greater quantities of four additional proteins in gynes . Experimental studies show that nourishment inequalities that correspond to early-season and late-season larval development contexts are associated with development of offspring having characteristics of worker and gyne phenotypes, respectively [38,39].
In addition to the established role of nourishment in caste differences in Polistes, social factors, such as dominance behavior [29,31] and maternal influences [40-42] also play a role. Jeanne and Suryanarayanan  propose a hypothesis for caste determination in primitively social wasps that incorporates not only nourishment variability, but also social environmental inputs from maternal care, specifically vibrational signals called antennal drumming. Antennal drumming may be an example of a maternal manipulation  that directs larvae toward a worker developmental trajectory, and this effect may interact with nourishment induced changes in caste phenotype. Thus, nourishment is likely to act in concert with social environmental factors in determining differences in gyne and worker caste development in Polistes.
In this study, we investigated the effect of experimental proteinaceous nourishment deprivation during larval development on caste-related gene expression in a primitively eusocial species, the paper wasp Polistes metricus. We had three main goals. First, we explored the transcriptional responses of wasp larvae to high and low nourishment levels during laboratory rearing using RNA-sequencing. Second, we tested the hypothesis that nourishment level relates to caste-related gene expression; specifically that low nourishment is associated with more worker-like gene expression patterns, and high nourishment is associated with more gyne-like expression patterns. To do this, we compared nourishment differential expression to a set of nearly 800 previously identified genes associated with caste development in field-reared P. metricus , and we did so on multiple levels: individual transcripts, pathways and biological functions. Third, we tested whether the molecular mechanisms underlying the response to nourishment are conserved across taxa by comparing our results to two other nourishment deprivation studies in fruit flies  and dung beetles . Our overall goal was to better understand the extent and nature of the role of nourishment and nourishment-related genes in caste development in a primitively social wasp species.
Differential expression analysis
To assess differential gene expression, we mapped reads to a previously assembled de novo transcriptome for P. metricus; the transcriptome was based on both the sequence data described here in conjunction with additional samples described in a previous study . By comparing expression patterns from head samples from 8 individual wasp larvae under high and low nourishment, we identified 284 P. metricus differentially expressed transcripts (DETs) that differed between low and high nourishment treatments, using DESeq (FDR ≤ 0.05, ), heretofore referred to as “P. metricus nourishment-responsive DETs”. Of these DETs, 207 (72.9%) were upregulated in low compared to high nourishment larvae (Figure 1 is a heat map of the scaled DET read counts across P. metricus nourishment samples. Additional file 1 includes the list of P. metricus nourishment-responsive DETs). Thus, despite the fact that they had less food available, low nourishment led to a majority of genes having higher gene expression, and thus did not simply cause a general shutdown in transcription.
One important feature of these results is the presence of outlier individuals, i.e. one individual low nourishment larva clustered with the high nourishment samples and a high nourishment larva clustered with the low nourishment samples (Figure 1). The variability of gene expression amongst biological samples is not completely surprising; our previous work has found high inter-individual variation in gene expression , effects of lab rearing (JM Jandt, JL Thomson, AC Geffre, AL Toth: Rearing environment may bias social traits: A case study with Polistes wasps, submitted), and variable effects of nourishment level on physiology .
Validation of select expression patterns via comparison to qRT-PCR data
It is important to validate RNA-Seq data using another method such as quantitative reverse transcription polymerase chain reaction (qRT-PCR). However, for this study, we had limited samples and quantities of RNA , so we could not perform qRT-PCR validation on actual samples from this experiment. In lieu of the sample limitations, we instead made a comparison to pre-existing data on nourishment-responsive expression patterns from adult brains of P. metricus under starved vs. ad lib food conditions . Although not ideal because of differences in tissue type and life stage, this comparison still provides a useful point of comparison to validate whether the expression patterns uncovered in this RNA-Seq study are robust.
Of 24 candidate genes, Daugherty et al.  identified 10 genes with differential expression in adult P. metricus brains reared under low and high nourishment conditions using qRT-PCR. We tested for a correlation (Spearman) in the log2 fold changes between the 24 candidate genes from Daugherty et al.  and the orthologous transcripts in this study, i.e. best BLAST hits between the primer sequences from Daugherty et al.  and the P. metricus transcriptome (Additional file 2: Figure S3 is a visual representation of the log2 fold changes for both studies; Additional file 3 lists the log2 fold change for each gene/transcript and the directionality for differentially expressed genes). For some genes, there is more than one best BLAST hit to the P. metricus transcriptome, so all transcripts were used for the correlation analysis. The gene PmTOR is absent from the P. metricus transcriptome, so this gene was removed from the analysis. Although none of the orthologous transcripts are significantly differentially expressed in the RNA-seq study, we nonetheless identified a significant positive correlation in log2 fold changes between these two studies (Spearman ρ: 0.54, p-value = 0.001), which provides support for our observed RNA-Seq results.
One reason we may not have observed statistically significant differential expression of these candidate genes in the RNA-Seq study, despite the strong correlation with the qRT-PCR data, is because of the limited statistical power. Differential expression calls are more stringent in this RNA-Seq study compared to the qRT-PCR analysis because familywise error correction is more severe due to the larger number of genes (over 75,000 compared to 24 genes). Furthermore, the qRT-PCR study has a larger sample size (n = 47) compared to the RNA-Seq study (n = 8), which provides greater statistical power for detecting differential expression.
Comparison to caste-related gene expression
Next, we investigated how P. metricus nourishment-responsive DETs compared to caste-related gene expression. A previous study  compared gene expression in field-collected, early season (worker-destined) larvae to late season (gyne-destined) larvae, and identified 736 caste-related DETs. Both low nourishment (described above) and worker-destined larvae  show a pronounced bias towards upregulated gene expression (72.9% upregulated in low nourishment and 91.7% upregulated in worker-destined larvae). This pattern generally agrees with the prediction that worker-destined and nourishment-deprived larvae have similar transcriptional states.
There was a statistically significant but relatively small overlap (43 transcripts) between the nourishment-responsive and caste-related DETs (Chi-square with Yates’ correction: p < 0.0001). Of these shared DETs, the directionality of gene expression change was unexpected: eight (18.6%) were upregulated in low nourishment larvae compared to high nourishment larvae, whereas 38 (88.4%) were upregulated in worker-destined larvae relative to queen-destined larvae (see Figure 2a for a Venn diagram of the number of unique and overlapping caste and nourishment-responsive DETs, Additional file 2: Figure S1 for a heatmap of the overlapping caste and nourishment-responsive DETs, and Additional file 4 for the list of overlapping DETs). To further explore these data beyond examining an overlap of gene lists, we performed a combined statistical analysis of data from both the nourishment level and caste contrasts, and the results also indicate some overlap in gene expression patterns across the two studies (see Methods, with complete description of the approach and results in Additional file 2 and Additional file 5). Taken together, these results partially support our prediction of gene expression similarity between nourishment-responsive and caste-related gene expression but suggest some unanticipated dissimilarities between them.
Pathway level analysis
Our previous study suggested gene expression similarity across studies is more pronounced on the level of pathways and gene functional categories, rather than specific genes or transcripts . Therefore, we also examined our data at the level of pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) [48,49]. Using best BLAST hits to D. melanogaster, very few transcripts (799 = 1.0% of the transcriptome) were annotated with enzyme codes for the KEGG analysis with Blast2GO  (see Additional file 1 for a list of transcripts annotated with enzyme codes by KEGG pathway). The presence of many transcripts without homology is a shortcoming of the dataset, but is expected and standard for Illumina-based data from non-model species . Only 7% (20) of the nourishment DETs were known members of KEGG pathways, therefore we identified very few (7 of 118) KEGG pathways with at least one DET. We then looked to see if these seven pathways were also related to caste differences . Four of the seven pathways also had caste-related DETs (including glycerolipid metabolism, nitrogen metabolism, and purine metabolism; listed in Additional file 4), but this amount of overlap was not statistically significant (Fisher’s exact test, p-value = 0.12). The majority of the DETs are upregulated in both low nourishment and worker-destined larvae for all four of these shared KEGG pathways (Figure 2b is a Venn diagram of the number of KEGG pathways with nourishment and/or caste DETs), which is in agreement with the prediction that nourishment-responsive biochemical pathways are regulated in a similar direction to what is found between castes in paper wasps.
Gene ontology (GO) enrichment analysis
On the level of gene functional categories, fifty-two GO terms were significantly enriched within the P. metricus nourishment-responsive DETs compared to the remainder of transcriptome (Figure 3; Additional file 1 includes the list of enriched GO terms). These included functions related to lipoprotein metabolism, oxidation reduction activity, and polysaccharide metabolism. More than 30% of GO terms (including terms related to oxidoreductase activity and carbohydrate metabolism) were common to both caste- and nourishment manipulation-related DETs, representing a significant overlap (Chi-square test with Yates’ correction, p < 0.0001; Figure 2b; Additional file 4 lists the shared caste and nourishment enriched GO terms). However, if we examine the direction of differential expression of DETs associated with these GO terms, we do not consistently see the predicted pattern of the same directional bias to both low nourishment larvae and worker-destined larvae (Figure 2b; Additional file 4 lists the direction).
Because nourishment is an important driver of phenotypic plasticity in many species, we were interested in determining whether the molecular mechanisms underlying the response to nourishment identified in paper wasps are conserved in other taxa. After an extensive literature search, we identified two studies from other insects that also examined transcriptional responses to nutritional stress: 1) a microarray study of tissue-specific (fat body and muscle) gene expression derived from fruit fly larvae  and 2) a microarray study of thoracic horn development in female dung beetle pupae . Although not ideal comparisons to our paper wasp dataset because of differences in life stages and sampled tissues, these studies still provide useful preliminary comparisons to begin addressing whether there are any conserved of nourishment-responsive transcripts.
When comparing P. metricus nourishment-responsive transcripts with nourishment-responsive transcripts in D. melanogaster fat body and muscle tissues , we identified small, non-significant overlaps: only 18 and 13 common DETs, respectively (Figure 4a and c; Chi-square tests with Yates’ correction; fat body: p-value = 0.84; muscle: p-value = 0.19). For both tissue types in fruit flies, most transcripts are down-regulated with low nourishment (Additional file 6 lists the common DETs and directionality for each species). Of these shared DETs, the majority are expressed in the same direction in both species: 56% (10 DETs for the fat body dataset; Fisher’s Exact Test, p-value = 0.676) and 54% (7 DETs for the muscle; Fisher’s Exact Test, p-value = 0.730).
At the level of GO categories, there is a statistically significant overlap (six GO terms including carbohydrate metabolic processes [GO:0005975] and oxidation-reduction process [GO:0055114]; see Additional file 6 for the complete list of overlapping GO categories) in the GO terms associated with nourishment-responsive transcripts for P. metricus and D. melanogaster fat bodies (Figure 4b; chi-squared test with Yates’ correction, p-value = 9.557e-11). The consistency in directionality between shared P. metricus and D. melanogaster fat body GO categories (i.e. up-regulation in low nourishment samples for all GO categories, except carbohydrate metabolism with up-regulation in high nourishment samples) further supports common functional changes related to nourishment manipulation in both species. However, this signal is not observed in the comparison between P. metricus nourishment-enriched GO terms and D. melanogaster nourishment-enriched GO terms in the muscle tissue. There is only one common nourishment-enriched GO term (oxidation-reduction process [GO:0055114]; Chi-square test with Yates’ correction, p-value = 0.4235), which is upregulated in the P. metricus low nourishment larvae but upregulated in the high nourishment samples for the D. melanogaster muscle tissue (Figure 4d).
Out of 18,061 homologous transcripts between P. metricus and a dung beetle Onthophagus taurus , we again identified a very small overlap (4 transcripts; not significant; Chi-square test, p-value = 0.57) in the nourishment-responsive transcripts for P. metricus female larvae and nourishment-responsive transcripts in O. taurus female pupal thoracic horns. However, this overlap between paper wasps and dung beetles is again consistent in the directionality of nourishment-responsiveness with shared DETs being upregulated in low nourishment individuals (Figure 4e; Additional file 6). Comparing enriched GO terms between paper wasps and dung beetles, we found only two shared GO categories: aminoglycan metabolic process (GO:0006022) and chitin metabolic process (GO:0006030) (Chi-square test with Yates’ correction, p-value = 0.08; Figure 4f), which suggests few common functional changes related to nourishment in samples from these two species.
In this study, we provide the first genome-wide transcriptional profiling of nourishment response during development in the genus Polistes, a model for understanding the evolution of social castes. We identified 285 nourishment differential expression transcripts (DETs) between larvae raised experimentally on low vs. high nourishment, many of which are associated with lipid metabolism and oxidation-reduction activity. Most (73%) of the Polistes nourishment-responsive transcripts are upregulated in larvae with low nourishment, which is opposite to the pattern observed in other insects including fruit flies  and dung beetles . There were few conserved nourishment-responsive transcripts across species, and these were related to aminoglycan metabolism, carbohydrate metabolism, and oxidation-reduction activity. Among transcripts corresponding to those functions, our data show some cross-species consistency in the direction of expression in response to nourishment deprivation. The overall picture from our preliminary cross-species comparisons is that it is largely different genes that show transcriptional responses to nourishment stress across insect orders. However, it is important to note that transcriptional similarities across these three systems may be underestimated from our analysis because of inconsistencies in the datasets such as analyzing different life stages and tissue types. Further work with directly comparable datasets is needed in order to better understand the extent of conservation of nourishment-responsive gene expression across taxa.
Nourishment inequalities have long been considered the most important environmental determinant of castes in Polistes [8,14,27,39]. If nourishment plays a major role in caste bias in P. metricus, we predicted that worker-destined larvae from natural nests in the field would have similar transcript expression patterns as experimental low nourishment larvae. With respect to directionality of overall transcript expression, this was the case — most DETs were upregulated in worker-destined larvae (92.1%) and in low nourishment larvae (72.9%). This work agrees with previous work in honey bees that also showed more genes with worker-biased expression in larvae, but more genes with queen-biased expression later in development (pupal stage) [18,21]. We also found a small but statistically significant overlap (43 common transcripts) between nourishment-responsive and caste-related DETs in Polistes (Chi-square with Yates’ correction: p < 0.0001; Figure 2; Additional file 2: Figure S1). Focusing on the 43 common DETs, however, shows that the majority of these were upregulated in worker-destined larvae but down-regulated in low nourishment larvae, which is opposite to our prediction and to the pattern of all DETs. This suggests that while nourishment restriction in the laboratory may lead to upregulation of transcript expression in some of the same pathways that are related to worker development, this manipulation did not succeed in causing a full shift to worker-like transcript expression patterns.
There was also a significant overlap in nourishment-responsive and caste-related expression at the GO functional level, with terms related to metabolism, protein binding, and oxidation reduction activity (Figures 2 and 3; Additional file 2: Figure S2; Additional file 4). On the level of pathways, we identified a few shared KEGG pathways between nourishment and caste datasets (Figure 2; Additional file 4), which may be due to the low annotation rate of nourishment-responsive DETs. Although there is limited overlap at the level of KEGG pathways, we did observe the expected pattern of directionality in gene expression, with the majority of DETs upregulated in low nourishment and worker-destined larvae for all pathways.
Overall, our data suggest that nourishment level caused a partial shift in gene expression, with low nourishment individuals being more worker-like and high nourishment individuals being more gyne-like. However, it is pertinent to note that the nourishment manipulation may not have resulted in a strong nutritional stress; i.e., wasps may have compensated for low proteinaceous food availability by consuming more sugar. In a recent study examining adult wasps reared from the same nests as the larvae analyzed here, the nourishment manipulation caused only a partial shift in caste-related physiology . In that study, total lipid and protein hemolymph levels in adults were not affected by low nourishment, whereas adults that had been reared with low nourishment showed greater ovary development after two weeks of being fed while in isolation. This counterintuitive observation corresponds to previous studies showing that adult workers collected in the field have greater ovary development than better-nourished gynes  and substantial evidence that workers are in a physiological state ready to reproduce upon emergence , whereas gynes are in reproductive diapause until after the overwintering period . Taken together, the physiological data  and transcriptomic data (this study) indicate that nourishment alone did not completely shift the developmental trajectory of larvae towards one caste or another.
One important caveat of our study is that the effects of the nourishment treatment on gene expression and physiology may be influenced by laboratory rearing, which may produce conflicting results compared to a natural field setting (JM Jandt, JL Thomson, AC Geffre, AL Toth: Rearing environment may bias social traits: A case study with Polistes wasps, submitted). Lab-reared wasps typically have higher lipid stores perhaps due to both overfeeding and inactivity , and lab-rearing can perturb caste-related gene expression (JM Jandt, JL Thomson, AC Geffre, AL Toth: Rearing environment may bias social traits: A case study with Polistes wasps, submitted). Two genes (inositol oxygenase and Hsp90alpha) known to exhibit caste-specific gene expression are perturbed due to lab rearing, but neither gene was differentially expressed in the current study. In addition, the low and high nourishment samples in our study were reared by single foundresses with unrestricted access to sucrose as opposed to the two feeding levels for caterpillars that they fed to larvae. Diets with a high carbohydrate to protein ratio can cause an increase in lipid levels in the insect fat body [53-55]. The effects of lab-rearing could thus be a factor in some of the non-overlap we observed between nourishment-responsive and caste-related DETs in the current study. Nonetheless, our data do suggest that low nourishment can trigger expression of some of the same genes and pathways that are associated with worker caste development in Polistes metricus.
Our results provide insight into the role of nourishment in the differential gene expression patterns that lead to different castes in Polistes. Our results also support the notion that nourishment inequalities alone cannot explain all caste variability in gene expression. Instead, caste-related gene expression bias is likely to be additionally influenced by social factors such as dominance behavior [29,31], vibrational communication [40-42], and/or epigenetically mediated environmental influences . Further work on gene expression in relation to nourishment, other influences on caste determination in Polistes, and interactions among them can further advance our growing understanding of caste determination in primitively social wasps.
This study provides new data on nourishment-responsive gene expression in the context of caste development for Polistes metricus, a model for studying the evolution of social insect castes. We identified suites of nourishment-responsive transcripts in developing P. metricus larvae. Interestingly, most transcripts were upregulated when larvae experienced proteinaceous nourishment deprivation; thus, reduced food level did not shut down gene expression but instead resulted in active transcription of many genes including several involved in lipid metabolism, carbohydrate metabolism, and oxidation-reduction processes. By comparing to previously reported caste-related gene expression patterns from the same species, we uncovered some similarity in transcripts, pathways, and gene functions related to both nourishment deprivation and worker caste-biased expression. However, many caste-associated genes were not found to be nourishment-responsive, so there are additional factors (likely including many social environmental factors) that influence caste-related gene expression. These results support the notion that nourishment level during development can somewhat bias development into queen or worker caste as adults, but leave room for other factors, and thus underscore the complex and multifactorial nature of caste development.
We collected adult female P. metricus wasps in early March 2009 near Raleigh, NC, as they were exiting the attic of a house in which they had passed the winter in behavioral quiescence and reproductive diapause. Single wasps were placed in cages 30 cm in length, width, and height constructed of clear plastic with plastic screen on the top and two sides. An opening in the top was covered with a piece of hardboard 10 cm square with a nest from a previous season attached to the underside. Nests were trimmed to seven cells ca. 0.5 cm deep, and meconia (larval feces) of the original nest occupants were removed. Each seven-cell nest therefore served as a “starter” nest on which the newly-captured wasps could initiate construction and then expand as their own. Construction paper was the pulp source for nest construction.
Caged wasps were placed on wire racks in a growth chamber 1.2 m × 2.5 m × 2.1 m in the North Carolina State University Phytotron. The chamber contained incandescent lights on a 16 L/8D cycle; fluorescent lights came on 2 h following and went off 2 h before the incandescent lights. Light intensity was 21 micromoles/sec/m2 in the incandescent-only morning and evening periods and 225 micromoles/sec/m2 during the both-lights midday. Chamber temperature was 20°C at night and 30°C during full light, with a gradual ramp-up in the morning and ramp-down in the evening. Each day, we repositioned cages in a rotation pattern to distribute any light or temperature variation in the chamber across all colonies over the course of the experiment.
Each cage was provided with ad libitum sucrose (rock candy) and a water source (15 ml tube plugged with cotton). Both high and low nourishment groups had equal access to sucrose, but differed in the amount of proteinaceous food provided. Proteinaceous nourishment was provided in the form of late 3rd or early 4th-instar Manduca sexta larva ca. 2.0 cm in length. Low nourishment foundresses received one caterpillar every fourth day. After May 22, the feeding schedule for low nourishment foundresses was accelerated to one caterpillar every three days. In the high nourishment treatment, foundresses were provided with caterpillars ad lib, adjusted each day to one caterpillar above the number that had been consumed since the previous day’s provisioning. For this experiment, we collected the largest final instar larva in each nest immediately following spinning of a cocoon by the third larva in the nest to do so. Those three larvae were allowed to develop into adults that then were analyzed for developmental and physiological characteristics. These results were published separately . For our study, four larvae per treatment were flash frozen in liquid nitrogen and stored at -80°C prior to RNA extraction.
RNA-extraction and RNA-sequencing
Polistes possess behaviorally distinct, not morphologically distinct castes, thus we were especially interested in developmental expression changes in larval brains. This led us to choose larval heads as the target tissue for our analysis. Whole heads were used rather than whole brains because of concerns about larval brain dissection quality and low RNA yields. We acknowledge that it is likely that the nourishment manipulation caused additional gene expression changes in other tissues we did not sample. However, we did not use whole bodies for two reasons: 1) our primary interest in brains and behavioral castes, and 2) the fact that brain gene expression differences in social insects can be subtle , meaning that high levels of expression from the fat body and other tissues could have drowned out expression patterns from the brain.
To preserve RNA during dissection, we removed the head region (indicated by dark coloration) of each individual larva, while kept on dry ice, using a sterilized razor blade. From these individual larval heads, we extracted total RNA using an RNeasy Mini Kit (Quigen), which was quality controlled using spectrophotometry (NanoDrop 2000) and a Bioanalyzer (Agilent). The High-Throughput Sequencing and Genotyping Unit of the W.M. Keck Center (University of Illinois at Urbana-Champaign) prepared 16 mRNA Seq libraries (n = 4 per group, high and low nourishment level) derived from two experimental nourishment levels (high and low; this study) and unmanipulated castes (queen- and worker-destined; ) using the “TruSeq RNAseq Sample Prep kit” (Illumina). These prepared libraries were sequenced on a HiSeq 2000 (Illumina) to generate over 1 billion 100 base paired-end reads, which we assembled de novo into a transcriptome (deposited at DDBJ/EMBL/GenBank under the accession GBGV00000000; this study used the first version GBGV01000000). Raw sequence data has been deposited to the National Center for Biotechnology Information’s (NCBI) Short Read Archive (BioProject ID: PRJNA242774, [accession numbers: NCBI:SRX511425, NCBI:SRX511426, NCBI:SRX511427, NCBI:SRX511430, NCBI:SRX511432, NCBI:SRX511433, NCBI:SRX511434, and NCBI:SRX511435]). For a more details about the RNA-Sequencing and the transcriptome assembly, see .
Read mapping, abundance estimation, and differential expression analysis
As described in , we aligned the raw paired-end reads to the reference transcriptome (GBGV01000000) using Bowtie 2 (Version 2.1.0)  with default settings. From these alignments, we quantified transcript abundances for each library using eXpress (Version 1.3.1) , which have been depositied in NCBI’s Gene Expression Omnibus  and are accessible through GEO Series [NCBI:GSE61960] (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE61960). From these raw read counts, we identified differentially expressed transcripts (FDR ≤ 0.05, ) between nourishment manipulation levels (high and low) using the R (Version 3.0.1)  statistical package DESeq (Version 1.12.0)  downloaded from the Bioconductor repository . Raw read counts were normalized by the effective library size, and dispersion factors were determined based on the per-condition method.
Kyoto encyclopedia of genes and genomes (KEGG) and gene ontology (GO) analyses
From functional annotations based on best BLASTx hit (E-value ≤ 1e-3) to Drosophila melanogaster sequences, we used Blast2GO (Version 2.6.5)  to assign enzyme codes, metabolic (KEGG) pathways, and GO terms to P. metricus transcriptomic sequences (see ). Based on these assignments, we identified the KEGG pathways that contain nourishment DETs and assessed enrichment of GO terms (FDR ≤ 0.05; one-tailed, ) between the nourishment DETs compared to the background – the complete P. metricus transcriptome.
Caste and nourishment manipulation comparison
To assess the extent that differential nourishment biased gene expression patterns toward a particular caste phenotype in Polistes paper wasps, we compared molecular signatures of nourishment manipulation and caste determination  at three levels: transcript differential expression, KEGG pathways with DETs, and GO enrichment. At the transcript level, we tested for significance in overlap between nourishment and caste DETs given the number of DETs per each dataset and the total number of P. metricus transcripts by using a Chi-square test with Yates’ correction. For both the nourishment and caste datasets, we identified which KEGG pathways had at least one DET, and we then performed a Fisher exact test to determine whether there is a statistically significant number of shared KEGG pathways with DETs given the number of unique pathways with DETs per each dataset and all pathways without DETs. Finally, we assessed significance in the overlap between nourishment and caste GO enriched categories in relation to the total number of annotated, not enriched, GO terms and unique nourishment/caste enriched GO terms using a Fisher exact test.
Alternative statistical approach to address overlap across nourishment and caste studies
We investigated an alternative approach (beyond comparisons of DET lists) to test for similarity in the transcript expression patterns associated with nourishment restriction and worker-caste development. In some sense, across the two studies we have a replicated experiment with respect to nourishment levels, which utilizes both the lab and field setting. As described in the introduction, worker-destined typically larvae receive more limited nourishment compared to queen-destined larvae [10,27,30,33]. Therefore, we modeled these datasets together while controlling for location (lab or field) to identify nourishment-responsive transcripts. Overall, our results suggest some consistency in the modeling and list comparison approaches. For full details and results using this approach, please see Additional file 2 for Supplemental Methods, Results, Table, and Figures S4-S7 and Additional file 5 for lists of nourishment-responsive transcripts.
Based on our interest in whether nourishment differences induce similar molecular responses in transcript expression and biological processes, we reviewed the literature and searched data repositories (NCBI Gene Expression Omnibus and EMBL ArrayExpress) to identify the most directly comparable studies to the P. metricus nourishment manipulation dataset based on the following criteria: 1) the datasets compared late-stage subadult insects reared on two nourishment levels and 2) the studies used a transcriptome-wide approach, i.e. microarray or RNA-seq. We identified two studies that meet both of these criteria: one in fruit flies (D. melanogaster, ) and the other in dung beetles (Onthophagus Taurus, ). For this analysis, we compared lists of putative orthologous nourishment-biased DETs and nourishment enriched GO terms between paper wasps and fruit flies or dung beetles using Chi-square tests with Yates’ correction. We defined putative orthologous sequences between P. metricus and either D. melanogaster or O. taurus as the best BLAST hit (E-value ≤ 1e-3) between pairs of species. For the fruit fly comparison, we only used the wild-type fruit fly data, and we compared the paper wasp dataset to both tissue types in fruit flies (adipose and muscle). The GO terms listed in the publication of Teleman et al.  were not consistent with the current Gene Ontology database, so we identified fruit fly nourishment enriched using Blast2Go (FDR ≤ 0.05; one-tailed, ). For the dung beetle comparison, we focused our study on a comparison with female thoracic horn dataset, which we felt was the most akin to the paper wasp dataset because samples were derived from tissues in female insects during the pupal stage of development.
Availability of supporting data
Raw sequence data has been deposited to the National Center for Biotechnology Information’s (NCBI) Short Read Archive (BioProject: PRJNA242774, http://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA242774). The transcript abundance data is available in the NCBI’s Gene Expression Omnibus repository (Accession number: GSE61960, http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE61960). All other data sets supporting the results of this article are included within the article and its additional files.
Differentially expressed transcripts
False discovery rate
Kyoto encyclopedia of genes and genomes
National Center for Biotechnology Information
Quantitative reverse transcription polymerase chain reaction
Whitman DW, Agrawal AA. What is phenotypic plasticity and why is it important? In: Whitman DW, Anathakrishnan TN, editors. Phenotypic plasticity of insects. Enfield, NH: Science Publishers; 2009. p. 1–-63.
Emlen DJ. Environmental-control of horn length dimorphism in the beetle onthophagus-acuminatus (Coleoptera, Scarabaeidae). P Roy Soc B-Biol Sci. 1994;256:131–6.
Shapiro AM. Seasonal Polyphenism. In: Hecht MK, Steere WC, Wallace B, editors. Evolutionary Biology. 1st ed. New York: Springer US; 1976. p. 259–333.
Ogawa K, Miura T. Aphid polyphenisms: trans-generational developmental regulation through viviparity. Front Physiol. 2014;5:1.
Wilson EO. The insect societies. Cambridge, Mass: Belknap Press of Harvard University Press; 1971.
Smith CR, Mutti NS, Jasper WC, Naidu A, Smith CD, Gadau J. Patterns of DNA methylation in development, division of labor and hybridization in an ant with genetic caste determination. Plos One. 2012;7:e42433.
Fo H. New observations on the natural history of bees. 3dth ed. Edinburgh: Printed for W. & C. Tait, and Longman, Hurst, Rees, Orme, and Brown, London; 1821.
Hunt JH. Nourishment and the evolution of the social Vespidae. In: Ross KG, Matthews RW, editors. The Social biology of wasps. 1st ed. Ithaca: Comstock Pub. Associates; 1991. p. 426–50.
Hunt JH. Nourishment and social evolution in wasps sensu lato. In: Hunt JH, Nalepa CA, editors. Nourishment and evolution in insect societies. 1st ed. Boulder. New Delhi: Westview Press; Oxford & IBH Pub. Co.; 1994. p. 211–44.
Hunt JH. The evolution of social wasps. Oxford; New York: Oxford University Press; 2007.
Patel A, Fondrk MK, Kaftanoglu O, Emore C, Hunt G, Frederick K, et al. The making of a queen: TOR pathway is a key player in diphenic caste development. Plos One. 2007;2:e509.
Wheeler DE, Buck N, Evans JD. Expression of insulin pathway genes during the period of caste determination in the honey bee, Apis mellifera. Insect Mol Biol. 2006;15:597–602.
Berens AJ, Hunt JH, Toth AL. Comparative transcriptomics of convergent evolution: Different genes but conserved pathways underlie caste phenotypes across lineages of eusocial insects. Mol Biol Evol. 2014 doi:10.1093/molbev/msu1330.
Hunt JH, Wolschin F, Henshaw MT, Newman TC, Toth AL, Amdam GV. Differential gene expression and protein abundance evince ontogenetic bias toward castes in a primitively eusocial wasp. Plos One. 2010;5:e10674.
Scharf ME, Wu-Scharf D, Zhou X, Pittendrigh BR, Bennett GW. Gene expression profiles among immature and adult reproductive castes of the termite Reticulitermes flavipes. Insect Mol Biol. 2005;14:31–44.
Zhou XG, Tarver MR, Scharf ME. Hexamerin-based regulation of juvenile hormone-dependent gene expression underlies phenotypic plasticity in a social insect. Development. 2007;134:601–10.
Toth AL, Robinson GE. Evo-devo and the evolution of social behavior. Trends Genet. 2007;23:334–41.
Barchuk AR, Cristino AS, Kucharski R, Costa LF, Simoes ZLP, Maleszka R. Molecular determinants of caste differentiation in the highly eusocial honeybee Apis mellifera. BMC Dev Biol. 2007;7:70.
Cameron RC, Duncan EJ, Dearden PK. Biased gene expression in early honeybee larval development. BMC Genomics. 2013;14:903.
Chen X, Hu Y, Zheng HQ, Cao LF, Niu DF, Yu DL, et al. Transcriptome comparison between honey bee queen- and worker-destined larvae. Insect Biochem Molec. 2012;42:665–73.
Evans JD, Wheeler DE. Differential gene expression between developing queens and workers in the honey bee, Apis mellifera. Proc Natl Acad Sci U S A. 1999;96:5575–80.
Feldmeyer B, Elsner D, Foitzik S. Gene expression patterns associated with caste and reproductive status in ants: worker-specific genes are more derived than queen-specific ones. Mol Ecol. 2014;23:151–61.
Ometto L, Shoemaker D, Ross KG, Keller L. Evolution of gene expression in fire ants: the effects of developmental stage, caste, and species. Mol Biol Evol. 2011;28:1381–92.
Severson DW, Williamson JL, Aiken JM. Caste-specific transcription in the female honey bee. Insect Biochem. 1989;19:215–20.
Weil T, Korb J, Rehli M. Comparison of queen-specific gene expression in related lower termite species. Mol Biol Evol. 2009;26:1841–50.
Winston ML. The biology of the honey bee. Cambridge, Mass: Harvard University Press; 1987.
Hunt JH, Amdam GV. Bivoltinism as an antecedent to eusociality in the paper wasp genus Polistes. Science. 2005;308:264–7.
Hunt JH. A conceptual model for the origin of worker behaviour and adaptation of eusociality. J Evolution Biol. 2012;25:1–19.
Pardi L. Dominance order in polistes wasps. Physiol Zool. 1948;21:1–13.
Reeve HK. Polistes. In: Ross KG, Matthews RW, editors. The Social biology of wasps. 1st ed. Ithaca: Comstock Pub. Associates; 1991. p. 99–148.
Jandt JM, Tibbetts EA, Toth AL. Polistes paper wasps: a model genus for the study of social dominance hierarchies. Insect Soc. 2014;61:11–27.
Starks PT, Turillazzi S. Polistes paper wasps: emergence of a model genus. Ann Zool Fenn. 2006;43:385–5.
West-Eberhard MJ. The social biology of polistine wasps. Ann Arbor: Museum of Zoology, University of Michigan; 1969.
Eickwort K. Separation of castes of polistes exclamans and notes on its biology (Hym - Vespidae). Insect Soc. 1969;16:67–72.
Strassmann JE, Lee RE, Rojas RR, Baust JG. Caste and sex-differences in cold-hardiness in the social wasps, polistes-annularis and polistes-exclamans (Hymenoptera, Vespidae). Insect Soc. 1984;31:291–301.
Hunt JH, Buck NA, Wheeler DE. Storage proteins in vespid wasps: characterization, developmental pattern, and occurrence in adults. J Insect Physiol. 2003;49:785–94.
Hunt JH, Kensinger BJ, Kossuth JA, Henshaw MT, Norberg K, Wolschin F, et al. A diapause pathway underlies the gyne phenotype in Polistes wasps, revealing an evolutionary route to caste-containing insect societies. Proc Natl Acad Sci U S A. 2007;104:14020–5.
Judd TM, Teal PEA, Hernandez EJ, Choudhury T, Hunt JH. Quantitative differences in nourishment affect caste-related physiology and development in the paper wasp Polistes metricus. Plos One. 2015;10(2):e0116199.
Karsai I, Hunt JH. Food quantity affects traits of offspring in the paper wasp Polistes metricus (Hymenoptera : Vespidae). Environ Entomol. 2002;31:99–106.
Jeanne RL, Suryanarayanan S. A new model for caste development in social wasps. Commun Integr Biol. 2011;4:373–7.
Suryanarayanan S, Hantschel AE, Torres CG, Jeanne RL. Changes in the temporal pattern of antennal drumming behavior across the Polistes fuscatus colony cycle (Hymenoptera, Vespidae). Insect Soc. 2011;58:97–106.
Suryanarayanan S, Hermanson JC, Jeanne RL. A Mechanical signal biases caste development in a social wasp. Curr Biol. 2011;21:231–5.
Alexander RD. The evolution of social behavior. Annu Rev Ecol Syst. 1974;5:325–83.
Teleman AA, Hietakangas V, Sayadian AC, Cohen SM. Nutritional control of protein biosynthetic capacity by insulin via myc in Drosophila. Cell Metab. 2008;7:21–32.
Kijimoto T, Snell-Rood EC, Pespeni MH, Rocha G, Kafadar K, Moczek AP. The nutritionally responsive transcriptome of the polyphenic beetle Onthophagus taurus and the importance of sexual dimorphism and body region. P Roy Soc B-Biol Sci. 2014;281:20142084–4.
Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995;57:289–300.
Daugherty THF, Toth AL, Robinson GE. Nutrition and division of labor: effects on foraging and brain gene expression in the paper wasp Polistes metricus. Mol Ecol. 2011;20:5337–47.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42:D199–205.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21:3674–6.
Ferreira PG, Patalano S, Chauhan R, Ffrench-Constant R, Gabaldon T, Guigo R, et al. Transcriptome analyses of primitively eusocial wasps reveal novel insights into the evolution of sociality and the origin of alternative phenotypes. Genome Biol. 2013;14:R20.
Toth AL, Bilof KBJ, Henshaw MT, Hunt JH, Robinson GE. Lipid stores, ovary development, and brain gene expression in Polistes metricus females. Insect Soc. 2009;56:77–84.
Musselman LP, Fink JL, Narzinski K, Ramachandran PV, Hathiramani SS, Cagan RL, et al. A high-sugar diet produces obesity and insulin resistance in wild-type Drosophila. Dis Model Mech. 2011;4:842–9.
Skorupa DA, Dervisefendic A, Zwiener J, Pletcher SD. Dietary composition specifies consumption, obesity, and lifespan in Drosophila melanogaster. Aging Cell. 2008;7:478–90.
Warbrick-Smith J, Behmer ST, Lee KP, Raubenheimer D, Simpson SJ. Evolving resistance to obesity in an insect. Proc Natl Acad Sci U S A. 2006;103:14045–9.
Weiner SA, Toth AL. Epigenetics in social insects: a new direction for understanding the evolution of castes. Genet Res Int. 2012;2012:609810.
Whitfield CW, Cziko AM, Robinson GE. Gene expression profiles in the brain predict behavior in individual honey bees. Science. 2003;302:296–9.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Roberts A, Pachter L. Streaming fragment assignment for real-time analysis of sequencing experiments. Nat Methods. 2013;10:71–3.
Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.
Team RC. R: A language and environment for statistical computing. In: Book R: A language and environment for statistical computing. 2013.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80.
Talbia Choudhury assisted with the wasp rearing at North Carolina State University. Carole Saravitz and Janet Shurtleff accommodated our work at the North Carolina State University Phytotron, including complimentary use of the growth chamber. Manduca sexta caterpillars were provided by the North Carolina State University Insectary, Beverley Pagura Director. Amy Geffre and Cecile Mercado extracted RNA for sequencing. Members of the Toth laboratory provided insightful comments during the preparation of the manuscript. This work was supported by the National Science Foundation IOS 1146410.
The authors declare that they have no competing interests.
JHH and ALT conceived the study and performed the experimental work. AJB performed the gene expression analyses, comparative studies, and generated the figures. AJB, JHH, and ALT wrote the manuscript. All authors have read and approved the final manuscript.
Polistes metricus differentially expressed transcripts, KEGG pathways, and GO categories. Contains the results of the differential expression analysis on the nourishment data, the list of KEGG pathways with more than one differentially expressed transcript, and the list of nourishment enriched GO categories.
Supplemental methods, results, tables, and figures. Includes supplemental figures for the overlap analysis of nourishment and caste and the results of the validation by comparison to qRT-PCR data. This file also contains the supplemental methods, results, tables, and figures of the additional statistical approach for determining overlap between the nourishment and caste data sets.
Validation by comparison to qRT-PCR data. Contains the log2 fold changes for the homologous sequences between the RNA-seq and qRT-PCR datasets. Sequences that are not differentially expressed are listed as NDE; whereas, the directionality is listed for the differentially expressed transcripts.
Results of nourishment and caste comparison at levels of transcripts, pathways, and biological functions. Lists the common differentially expressed transcripts, KEGG pathways, and GO categories between the nourishment and caste data sets.
Results of additional statistical approach for determining overlap between nourishment and caste data sets. Includes the estimates of nourishment main effects and nourishment by location interaction effects for the transcripts with significant effects.
Nourishment comparison across species. Lists the overlapping nourishment differentially expressed transcripts and GO categories for the comparisons between paper wasps and fruit flies (fat body and muscle tissue) or dung beetles.