Transcriptome analyses of primitively eusocial wasps reveal novel insights into the evolution of sociality and the origin of alternative phenotypes
© Ferreira et al.; licensee BioMed Central Ltd. 2013
Received: 17 December 2012
Accepted: 26 February 2013
Published: 26 February 2013
Understanding how alternative phenotypes arise from the same genome is a major challenge in modern biology. Eusociality in insects requires the evolution of two alternative phenotypes - workers, who sacrifice personal reproduction, and queens, who realize that reproduction. Extensive work on honeybees and ants has revealed the molecular basis of derived queen and worker phenotypes in highly eusocial lineages, but we lack equivalent deep-level analyses of wasps and of primitively eusocial species, the latter of which can reveal how phenotypic decoupling first occurs in the early stages of eusocial evolution.
We sequenced 20 Gbp of transcriptomes derived from brains of different behavioral castes of the primitively eusocial tropical paper wasp Polistes canadensis. Surprisingly, 75% of the 2,442 genes differentially expressed between phenotypes were novel, having no significant homology with described sequences. Moreover, 90% of these novel genes were significantly upregulated in workers relative to queens. Differential expression of novel genes in the early stages of sociality may be important in facilitating the evolution of worker behavioral complexity in eusocial evolution. We also found surprisingly low correlation in the identity and direction of expression of differentially expressed genes across similar phenotypes in different social lineages, supporting the idea that social evolution in different lineages requires substantial de novo rewiring of molecular pathways.
These genomic resources for aculeate wasps and first transcriptome-wide insights into the origin of castes bring us closer to a more general understanding of eusocial evolution and how phenotypic diversity arises from the same genome.
Phenotypic plasticity is a fundamental biological process that allows organisms to adapt to changes in their environment [1, 2]. Examples of plastic phenotypes include insect castes [3, 4], horn-polyphenic beetles  and sex differences , where they play a crucial role in shaping the ecology and evolution of species and ecosystems. Understanding how alternative phenotypes arise from the same genome is one of the most challenging questions in modern biology [1, 2].
The most impressive examples of phenotypic plasticity are found in the eusocial Hymenoptera (ants, some bees and some aculeate wasps), where specialized reproductive or non-reproductive phenotypes (castes) play a pivotal role in the ecological success of these insects . Eusocial insect castes usually arise through differential expression of shared genes [8, 9], but studies to date have been largely restricted to the highly eusocial bees and ants, where caste phenotypes are evolutionarily derived, being greatly modified from their ancestral state (for example, with morphological adaptations that are established early in development). Thus, some of the molecular processes underlying the origins of castes are likely to differ substantially from those underlying highly eusocial species [10, 11]. Understanding how castes first arise is best studied in primitively eusocial species, but little is known about their molecular evolution, particularly in the aculeate wasps [12–15]. Thus, there is a need for comprehensive gene-level studies on primitively eusocial wasps.
In the evolution of eusociality, behavioral and physiological traits of the non-social ancestor become decoupled into complementary queen and worker phenotypes , but we understand little about the molecular patterns and processes by which this originates and to what extent molecular ground plans are re-organized. These questions are of general biological importance in understanding the evolution of phenotypic diversity [1, 2]. Empirical assessment of how gene transcription is decoupled into alternative phenotypes has been difficult to obtain because it requires unbiased quantification of gene expression (of both known and novel (undescribed) genes) across whole genomes . Microarray analyses of the highly eusocial honeybee (Apis mellifera) and fire ant (Solenopsis invicta) suggest marked functional decoupling of molecular processes in queens and workers, involving large numbers of genes [18–21]. Little is known about the patterns and processes of transcriptional decoupling for castes in primitively eusocial insects [12, 14]. A simple prediction is that patterns of transcriptional and functional decoupling should be less pronounced than in highly eusocial species, because phenotypes are in the first stages of merely behavioral decoupling [9, 12].
Processes of phenotypic decoupling at the molecular level may occur via several mechanisms [5, 22]. Conserved toolkits of ancestral molecular processes may be redeployed to regulate similar alternative phenotypes across species and lineages [23–26]. In eusocial insects, conserved toolkit genes appear to be differentially expressed in castes across species, including bees and wasps [12, 27]. An alternative, untested prediction is that at least some of the molecular processes underlying social phenotypes, in different social lineages and among different levels of social complexity, will differ substantially. This is because extensive molecular and developmental re-wiring or new gene evolution may be required for the loss of phenotypic plasticity and elaboration of phenotype-specific traits [10, 28]. There may also be rapid evolution of genes associated with the emergence of alternative phenotypes [6, 17, 29], and there is evidence for this in the highly eusocial honeybees and fire ants [30, 31]. Finally, the importance of novel genes (that is, previously undescribed genes that lack any detectable protein-coding homologues from existing sequence data ) in the evolution of phenotypic innovations has recently emerged in animals as varied as yeast , hydra, reptiles [34–36], and also honeybees . To date the eusocial insect literature has largely focused on conserved molecular pathways for social behavior [12, 26, 27, 38, 39], and the relative roles of the different mechanisms of phenotypic decoupling are unknown.
Here, we generate a genome-wide catalogue of aculeate wasp genes expressed in adult females to provide a resource for genomic analyses. We then conduct the first RNA-seq analyses of caste-biased expression in a wasp and primitively eusocial insect to test the above hypotheses on the patterns and processes of molecular decoupling of alternative phenotypes at the early stages in the evolution of eusocial behavior. These datasets allow a first simultaneous assessment of the role of conserved genetic toolkits, novel genes and gene evolution when social behavior first evolves. We also use our data to re-assess the phylogenetic relationships between the three aculeate Hymenoptera subfamilies (bees, wasps and ants), which is vital information for interpreting common processes and lineage-specific novelties underlying the evolution of alternative phenotypes. Recently, there has been much debate over the relationships of aculeate bees, wasps and ants, challenging the morphology-based view that Vespoidea (ants plus aculeate wasps) are monophyletic [46, 47]. Genome sequences (or similar catalogue of genes) for aculeate wasps will help resolve this controversy and provide a basis for interpreting comparative data on social evolution.
Results and discussion
Gene assembly for an aculeate wasp
We first generated a reference assembly of the genes expressed by P. canadensis adult female phenotypes (queens, workers, foundresses and callows) (Figure 1; Section 1 in Additional file 1) by sequencing normalized 454 libraries from phenotype-specific pools (37 individuals in total; Section 2 in Additional file 1). Using transcriptome rather than genome sequencing data to provide a gene set is a powerful and accessible approach for initiating genomic analyses on non-model organisms. Library normalization meant that we were able to capture transcripts at all levels of the expression spectrum. The longer read lengths provided by 454 technology and the high coverage of these datasets increase the likelihood of full transcript assembly. The 454 reads were assembled into a reference transcriptome set of 26,284 isogroups - henceforth referred to as 'genes' for simplicity (Section 2 in Additional file 1; Additional file 2). This gene set included 98% of the core eukaryotic genes with equal coverage for each phenotype (Section 3 in Additional file 1). Additionally, sequences were 100% identical to existing overlapping Sanger sequences for this species (n = 36 genes) . Genes were equally spread across honeybee chromosomes and included 94% of publicly available ESTs for the temperate paper wasp Polistes metricus (n = 422 genes ; Section 3 in Additional file 1). Although this assembly is restricted to the genes expressed in adult females, it provides an important step forward in generating a comprehensive genome-wide sequence resource for aculeate wasps.
Number of best BLAST hits between the genes of sequenced genomes for aculeate hymenopterans
The characteristics of isotigs lacking homology - namely shorter ORFs, lower protein coding potential, lower GC content and lower expression values  - suggest that a large proportion of novel isotigs correspond to non-protein coding transcripts. But low sequence conservation across species, lack of genomic information for aculeate wasps and the typical tissue specificity of the non-coding RNAs  may also explain the novelty of these genes.
We also verified that the abundance of novel genes was not an artifact resulting from assembly of transcript fragments rather than complete transcript sequences (Section 6 in Additional file 1). Our analyses suggested that less than 10% of genes were affected by this, and thus do not alter significantly our main analyses (below). However, without a genome sequence we cannot exclude the possibility that some genes may represent incomplete assemblies.
Paraphyly of Vespoidea
We further tested the phylogenetic relationships using a Bayesian framework. In contrast, this gave greater support to the classical view of aculeate wasps and ants being monophyletic (Figure 3b; Section 10 in Additional file 1). The reason for the discrepancy in the results of these two methods is likely to be due to the wasp clade being the most poorly represented in our analysis (with only one species), and moreover, it is the clade with the longest branch. Our dataset offers a first chance to examine subfamily relationships across large numbers of genes but, clearly, phylogenomic data on more species of aculeate wasps are required to determine whether the term 'Vespoidea' should be dropped, or reclassified as paraphyletic. This finding would have important general implications for our understanding of eusociality as it would suggest that bees and ants shared an aculeate wasp-like common ancestor, that ants are wingless wasps , and that bees are wasps that lost predacious behaviors. Primitively eusocial wasps may therefore hold some important clues to the early stages of eusociality, if they are the basal descendants of non-aculeate Hymenoptera.
Transcriptional and functional decoupling at the origin of alternative phenotypes
We next conducted RNA-seq analyses of phenotype-biased expression to provide a first genome-wide look at molecular decoupling of alternative phenotypes at the early stages of eusocial evolution. We quantified genome-wide transcription in individuals (n = 2 to 12 per phenotype) from each of the four adult female phenotypes by sequencing at least 3 Gbp of Illumina short-reads of brain cDNA per phenotype (Additional files 6, 7, 8 and 9). These samples appeared to capture most transcription across the genome since 99.8% of the 26,284 genes in the reference transcriptome were detected in our pooled Illumina dataset (Section 7 in Additional file 1).
Isotig and isogroup differential expression statistics for the four adult female phenotypes in P. canadensis
Up-regulated (q > 0.6) isogroups (isotigs)
Isogroups with recognized homologs
Isogroups with annotated homologs
Down-regulated (q > 0.6) isogroups (isotigs)
Isogroups with recognized homologs
Isogroups with annotated homologs
Isotig and isogroup differential expression in the comparison of only queen and worker castes
Up-regulated (q > 0.6) isogroups (isotigs)
Isogroups with recognized homologs
Isogroups with annotated homologs
Putative lncRNAs (total identified = 3,824)
Interestingly, foundresses and callows accounted for the majority (89%) of down-regulated genes. Some of these could be attributed directly to a corresponding up-regulation in workers: specifically, 55% and 75%, respectively, of the down-regulated genes in foundresses and callows. But, the remaining foundress and callow down-regulated genes were not upregulated in other castes. This is also reflected in the GO classification (Additional file 5), where none of the 22 GO terms over-represented among the down-regulated genes in foundresses and only 5 of the 23 terms over-represented among the down-regulated genes in callows were up-regulated in workers. Although preliminary, these results raise some important questions: they suggest that callows are born neither worker- nor queen-like, but perhaps instead as an undifferentiated state, as expected in a species where all females retain full reproductive totipotency; down-regulation in callows may also reflect maturation processes that take place in the first days after emergence as an adult. Foundresses are thought to be hopeful reproductives, co-founding a nest in the hope of becoming the dominant egg-layer (queen); our data suggest that foundresses are quite different from established queens, even though they fit the queen phenotype in all other respects (that is, they are active egglayers, rarely absent from the nest, and are behaviorally dominant). Taken together, these results highlight how the definition of a phenotype in this species is more complex than a simple distinction between queens and workers.
Using the 268 genes differentially expressed between queens and workers that could be annotated (Table 3), we then looked for any signs of functional specialization between these castes. Of the GO terms corresponding to these genes (Additional file 4), 2.4% (53 out of 2,242) were significantly enriched between queens and workers (Additional file 5; 40 terms (1.7%) over-represented in queens, 13 terms (0.57%) over-represented in workers). Although many of our genes did not have any functional information, there are some interesting first insights from these analyses. The putative lack of substantial functional specialization between castes in Polistes contrasts with the marked functional specialization in highly eusocial species, such as honeybees . This may reflect selection on females to retain reproductive totipotency in primitive societies, like P. canadensis, in order that they can exploit alternative reproductive strategies , and minimize the risks of specialization associated with small colony size . General terms overrepresented in workers include cellular protein modification processes, membrane and ion transport, which may suggest elevated cellular activity in worker brains. Among the over-represented terms in queens, several candidates associated with cell structure and cytoskeleton were found, such as homologs of myosin chains, actin isoforms and troponin, which may be associated with synaptic plasticity and memory . However, because so few of the genes up-regulated in workers could be annotated, and because of the lack of significant levels of enrichment between castes, these inferred functional differences between castes should be taken with caution and require further analyses once some functionality can be assigned to the unknown genes.
Role of novel genes in the evolution of alternative phenotypes
We explored the hypothesis that novel genes (those lacking homologs in NR databases; see above) may play an important role in the early stages of phenotypic decoupling in P. canadensis. The importance of novel genes was recently highlighted in organisms as varied as Hydra to reptiles for their contribution to phenotypic diversity, through generating both behavioral and morphological variation [34–36]. Of the 2,442 up-regulated caste-biased genes (Table 2), 24.4% (n = 596 genes) had known homologues, indicating that both novel and described genes are important in P. canadensis social behavior. Interestingly, 81.5% of worker-biased genes lacked homologs, compared to only 17.1% of queen-biased genes (Table 2). Foundresses and callows show similar levels of homology to queens. Lack of homology cannot be attributed to fragmentation or other anomalies arising from using transcriptome data (Section 6 in Additional file 1) and they are also valid for the stricter comparison of queens versus workers (Table 3). These results are exactly in line with Johnson and Tsutsui  and Barchuk et al. , who found an over-representation of novel (taxonomically restricted) genes up-regulated in honeybee workers (adults  and brood ), and suggested that these novel genes may facilitate expansion of worker behavioral complexity. Our data provide the first suggestion that novel gene transcription may be important at the early stages of worker evolution, confirming the potentially broad role of novel genes in eusocial evolution  and adding to the emerging general importance of novel genes in the evolution of phenotypic diversity [34–36]. It is important to note that microarray methods (as used in the honeybee studies) are largely biased towards conserved genes , as sequence conservation facilitates gene prediction. Thus, RNA-seq methods as used in our study offer new potential to explore the role of putative novel genes in phenotypic evolution.
Our ability to interpret what these novel genes are is limited without genome and proteome analyses, which is beyond the scope of this study. But, we explored the hypothesis that many of the caste-biased novel genes were in fact non-coding RNAs, as identified from our whole transcriptome analysis (above). Interestingly, we found a significant difference in the proportion of novel transcripts that were putatively lncRNAs, with an over-representation in workers relative to queens (Table 3). lncRNAs are thought to play a role in gene regulation , but also may be evidence of de novo gene origin . Future work will determine the significance of these processes in worker, rather than queen, evolution.
Conserved molecular toolkits associated with alternative phenotypes
Social lineages evolve from non-social ancestors . Alternative phenotypes (social insect castes) may therefore evolve through the decoupling of conserved sets of genes that regulated changes in hormone titers underlying the provisioning and reproductive phases of the non-social ancestor . A prevailing hypothesis is that changes in the patterns of expression of these ancestral ground-plan genes could account for caste evolution across social taxa [26, 28]. Using a candidate gene approach, sets of such 'toolkit' genes and molecular processes have been shown to have conserved roles in caste regulation across bees and wasps [12, 27, 38]. We also found some evidence of this in P. canadensis: vitellogenin (four isogroups), insulin (one isogroup) and major royal jelly protein (two isogroups) were queen-biased; juvenile hormone (two isogroups) and methyltransferases (two isogroups) were worker-biased. However, other toolkit genes (hexamerin, Malvolio, Amfor, P450) were present in our dataset, but were not identified as significantly caste-biased (Additional file 4) [12, 27, 65, 66].
As a genome-wide, unbiased representation of genes expressed by alternative phenotypes in wasps, our datasets also offer the opportunity to test to what extent a conserved molecular toolkit in general explains alternative social phenotypes across eusocial taxa [12, 27]. Interestingly, we found little correlation in the general identity and direction of expression of genes underlying castes in other eusocial insects, as identified from microarray analyses. For example, only 6.5% of honeybee (A. mellifera) caste-biased genes  were caste-biased in P. canadensis and there was no consistency in the direction of expression. Of 50 cDNAs predictive of honeybee worker behavior , only 8 were found in P. canadensis and they were not significantly worker biased (Additional file 5). There was no correlation in direction and identity of gene expression in castes of the highly eusocial fire ant S. invicta (r ~ 0) . Finally, only 22 of the 423 caste-biased genes of the temperate wasp P. metricus  were also caste-biased in P. canadensis, and they showed a low correlation in the direction of phenotype-specific expression between species (Pearson correlation of 0.25 for worker, 0.14 for queen, 0.24 for foundress and 0.23 between gyne and callow; Additional file 5). This was unexpected since these two species are close relatives, share a recent common ancestor, and have similar social and behavioral traits. Differences in ecology of tropical and temperate Polistes, and possible secondary adaptation to temperate living , may contribute to these differences. These analyses suggest that the social toolkit represents a small portion of the genome involved in phenotype decoupling of P. canadensis, and that there is much lineage and ecology-specific transcription underlying caste differences. This supports the hypothesis that the genes underlying castes in different social lineages will not necessarily be conserved, because extensive molecular and developmental re-wiring may be required in the evolution of caste commitment and eusociality [10, 11].
Insights into gene evolution in alternative phenotypes
Rapid evolution of genes has been associated with alternative phenotypes in a range of animals [6, 17, 29]. Caste-biased genes evolve rapidly in the honeybee and fire ant [30, 31], possibly because fast-evolving genes are preferentially recruited into caste-biased gene expression . Interestingly, we found no evidence of this in P. canadensis (Section 10 in Additional file 1; Additional file 5), suggesting that at the early stages of phenotypic decoupling genes are maintained through purifying selection. This agrees with social evolutionary theory that predicts that females retain the ability to exploit alternative reproductive opportunities (for example, to switch from worker to queen) to maximize individual-level inclusive fitness at the early stages of sociality . Selection may be relaxed after workers commit to a non-reproductive pathway, as found in highly eusocial species [10, 28], resulting in gene evolution, new genes and new gene networks [31, 67–69]. Analyses of incipient phenotypic decoupling in other animals will help determine whether accelerated gene evolution only occurs after key innovations or fundamental shifts in life-history traits.
We next looked for evidence of processes that may have been lost or gained during phenotypic decoupling in eusocial evolution (Section 10 in Additional file 1; Additional file 5). We identified 431 caste-biased genes from across eusocial insect species that are also conserved in the 19 insect species genomes available. From this conserved set of genes, 12 GO terms were over-represented in highly eusocial insects (see Section 10 in Additional file 1 for species list) relative to P. canadensis, and 217 terms were under-represented (Section 10 in Additional file 1; Additional file 5). Under-represented groups (which include regulatory and signaling processes) may be less important in complex eusociality and hence not maintained by stabilizing selection. However, we cannot exclude the possibility that these differences are wasp-specific. Over-represented terms may be important in complex eusocial behavior . These included processes previously highlighted for putative roles in caste determination (for example, oxidation reduction, metabolic processes) , suggesting they reflect the level of eusociality rather than differences in ecology or lineage.
Our RNAseq analyses of phenotypic decoupling in P. canadensis provide some intriguing novel insights into the evolution of alternative phenotypes, and the early stages of caste evolution and sociality. First, they highlight the importance of novel genes in phenotypic diversity, emphasizing the need to expand existing genomic datasets beyond established model organisms to include a wider range of taxonomic groups. Secondly, they add important insights to the current emerging picture that much of the molecular changes that accompany social evolution occur predominantly in the worker rather than queen caste, with an over-representation of novel caste-biased genes in workers. These features may be important in facilitating the evolution of behavioral complexity in the worker caste. We show that these patterns are evident even at the early stages of sociality, where castes retain plasticity. Further, contrary to highly eusocial species, genes involved in caste differentiation do not appear to be subject to relaxed selection at the early stages of sociality. Diversity in the molecular regulation of castes across social lineages is expected if many pleiotropic genes involved with sexual conflict are lost in the early stages of social evolution, when ancestral monogamy is required to generate the conditions for the evolution of worker behavior . Subsequent evolution of queen and worker behaviors, therefore, may take very different molecular pathways in different lineages, requiring new gene networks to evolve independently . Further RNA-seq studies on other species (for example, ) and also genome sequencing to explore the potential role of alternative splicing (for example, ) will help determine the complementary roles of conserved and novel molecular processes in shaping social and other polyphenisms, bringing us closer to understanding how genomes give rise to phenotypic diversity in general.
Materials and methods
P. canadensis wasps of known behavioral repertoires were collected from wild populations in Panama, in July 2009 (Punta Galeta, Colon). All wasps were collected directly off their nests with forceps around midday during the active periods (that is, sunny weather) and preserved immediately in RNAlater (Ambion, Invitrogen, Applied Biosystems), and stored at -20°C until analysis (Section 1 in Additional file 1).
Transcriptome sequencing, assembly and analyses
454 sequencing of pooled samples of 37 wasps across phenotypes (5 to 18 individuals per phenotype; 2.1 million reads, 80% brain, 10% abdomen, and 10% antennae) was used to generate a reference transcriptome (Section 2 in Additional file 1). Newbler v2.3 was used to generate the final assembled gene set (Table S1 in Additional file 1). Transcripts were annotated using GO categories assigned using BLASTx of GenBank NR databases with a conservative e-value threshold of 10-5, and Blast2Go was used to assess enrichment of GO terms among phenotypes (Section 4 in Additional file 1). Illumina sequencing of 14 biological replicates (>377 million reads) across 5 lanes was conducted to quantify differential gene expression, expressed as RPKM (reads per kilobase per million) values (Section 8 in Additional file 1). We trialed a number of methods for identifying differentially expressed genes and settled on a novel non-parametric method (NOISeq ) This method infers the noise distribution from the data and performs pairwise comparison of the samples to identify differentially expressed genes. A variety of probability thresholds were tested (Section 8 in Additional file 1). For the GO analysis we used a q-value >0.6 that represents a 50% chance that the gene is differentially expressed rather than not differentially expressed.
Protein sequences were aligned using MUSCLE , with default parameters. This is a multiple sequence aligner that includes an iterative alignment refinement phase to overcome known pitfalls of the progressive alignment strategy. Subsequently, poorly aligned regions of the alignment were trimmed with trimAl v1.3.  to remove columns with gaps in more than 30% of the sequences. A maximum-likelihood analysis was conducted on the concatenated alignment containing 33,506 sites using PhyML v3.0 . Out of a total of five evolutionary models (LG, JTT, WAG, VT, BLOSUM62), the general replacement model LG (after 'Lee and Gascuel') was found to be the best fitting model using the AIC criterion . In all cases four categories of evolutionary rates were used, estimating the gamma shape parameter and the proportion of invariable sites from the data. Branch supports were obtained using an approximate likelihood ratio test as implemented in PhyML ('Minimum of SH and chi-square' option). The resulting topology was compared to an alternative topology placing Polistes as a sister group to the ants. To do so, the phylogeny was re-computed, fixing the monophyly of Polistes and ants, but allowing the rest of the topology to be optimized. Support for the two topologies were compared using a Shimodaira-Hasegawa test, as implemented in CONSEL . A Bayesian analysis was conducted as implemented in PhyloBayes , using the default CAT model and running two independent MCMC runs during 300,000 generations, and sampling every 100 generations. Consensus trees were built after removing the first 20% sampled trees and using a majority consensus rule.
Raw sequence data are available at the European Read Archive (accession number ERP001342). The Transcriptome Shotgun Assembly project has been deposited at DDBJ/EMBL/GenBank under the accession GAFR01000001-GAFR01045087. The version described in this paper is the first version, GAFR01000000. All data and datasets can also be accessed at .
long non-coding RNA
open reading frame
reads per kilobase per million.
We thank Centre for Genomic Regulation (CRG) (Heinz Himmelbauer and Francisco Câmara), GenePool and the Babraham Institute (Wolf Reik and Simon Andrews) for support and access to their sequencing facilities. We thank RCUK, CONSOLIDER (CSD2007-0005), the Spanish Ministry of Science (BIO2006-03380, BFU2009-09168), GTEX (NIH Project 1RO1MH090941-01), and a Wellcome Trust Senior Investigator Award to Wolf Reik. SS was supported by a UK Research Council (RCUK) fellowship, and PGF by a FCT postdoctoral grant (SFRH/BPD/42003/2007) from the Portuguese Ministry of Science and Technology. We thank the Smithsonian Tropical Research Institute (particularly William Wcislo, Mary Jane West-Eberhard, Jorge Morales and staff at the Galeta Field Station) for hosting the fieldwork, and Stephanie Dreier and Thibault Lengronne for assistance in the field. The fieldwork was conducted under ANAM Permit No. SE/A-33-09.
- West-Eberhard MJ: Developmental Plasticity and Evolution. 2003, New York: Oxford University PressGoogle Scholar
- Nijhout HF: Development and evolution of adaptive polyphenisms. Evol Devel. 2003, 5: 9-18. 10.1046/j.1525-142X.2003.03003.x.View ArticleGoogle Scholar
- Wilson EO: The evolution of caste systems in social insects. Proc Amer Phil Soc. 1979, 123: 204-211.Google Scholar
- Holldobler B, Wilson E: The Ants. 1990, Berlin: SpringerView ArticleGoogle Scholar
- Moczek AP: Phenotypic plasticity and diversity in insects. Phil Trans R Soc Lond B. 2010, 365: 593-603. 10.1098/rstb.2009.0263.View ArticleGoogle Scholar
- Ellegren H, Parsch J: The evolution of sex-biased genes and sex-biased gene expression. Nat Rev Genet. 2007, 8: 689-698. 10.1038/nrg2167.PubMedView ArticleGoogle Scholar
- Oster GF, Wilson EO: Caste and Ecology in the Social Insects. 1978, Princeton: Princeton University PressGoogle Scholar
- Smith CR, Toth AL, Suarez AV, Robinson GE: Genetic and genomic analyses of the division of labour in insect societies. Nature Rev Genet. 2008, 9: 735-748.PubMedView ArticleGoogle Scholar
- Sumner S: Determining the molecular basis of sociality in insects: prospects, progress and potential in sociogenomics. Ann Zool Fenn. 2006, 43: 423-442.Google Scholar
- Boomsma JJ: Lifetime monogamy and the evolution of eusociality. Phil Trans R Soc Lond B. 2009, 364: 3191-3207. 10.1098/rstb.2009.0101.View ArticleGoogle Scholar
- Linksvayer TA, Wade MJ: Genes with social effects are expected to harbor more sequence variation within and between species. Evolution. 2009, 63: 1685-1696. 10.1111/j.1558-5646.2009.00670.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Toth AL, Varala K, Henshaw MT, Rodriguez-Zas SL, Hudson ME, Robinson GE: Brain transcriptomic analysis in paper wasps identifies genes associated with behavior across social insect lineages. Proc R Soc Lond B. 2010, 277: 2139-2148. 10.1098/rspb.2010.0090.View ArticleGoogle Scholar
- Toth AL, Varala K, Newman TC, Miguez FE, Hutchison SK, Willoughby DA, Simons JF, Egholm M, Hunt JH, Hudson ME, Robinson GE: Wasp gene expression supports an evolutionary link between maternal behavior and eusociality. Science. 2007, 318: 441-444. 10.1126/science.1146647.PubMedView ArticleGoogle Scholar
- Sumner S, Pereboom JJM, Jordan WC: Differential gene expression and phenotypic plasticity in behavioral castes of the primitively eusocial wasp, Polistes canadensis. Proc R Soc Lond B. 2006, 273: 19-26. 10.1098/rspb.2005.3291.View ArticleGoogle Scholar
- Hoffman EA, Goodisman MAD: Gene expression and the evolution of phenotypic diversity in social wasps. BMC Biol. 2007, 5: 23-10.1186/1741-7007-5-23.PubMedPubMed CentralView ArticleGoogle Scholar
- West-Eberhard MJ: Wasp societies as microcosms for the study of development and evolution. Natural History and Evolution of Paper-Wasps. Edited by: Turillazzi S, West-Eberhard MJ. 1996, Oxford Science Publications, 290-317.Google Scholar
- Snell-Rood EC, Cash A, Han MV, Kijimoto T, Andrews J, Moczek AP: Developmental decoupling of alternative phenotypes: insights from the transcriptomes of horn-polyphenic beetles. Evolution. 2011, 65: 231-245. 10.1111/j.1558-5646.2010.01106.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Whitfield CW, Cziko AM, Robinson GE: Gene expression profiles in the brain predict behavior in individual honey bees. Science. 2003, 302: 296-299. 10.1126/science.1086807.PubMedView ArticleGoogle Scholar
- Grozinger CM, Fan Y, Hoover SE, Winston ML: Genome-wide analysis reveals differences in brain gene expression patterns associated with caste and reproductive status in honeybees (Apis mellifera). Mol Ecol. 2007, 16: 4837-4848. 10.1111/j.1365-294X.2007.03545.x.PubMedView ArticleGoogle Scholar
- Alaux C, Sinha S, Hasadsri L, Hunt GJ, Guzmán-Novoa E, DeGrandi-Hoffman G, Uribe-Rubio JL, Southey BR, Rodriguez-Zas S, Robinson GE: Honey bee aggression supports a link between gene regulation and behavioral evolution. Proc Natl Acad Sci USA. 2009, 106: 15400-15405. 10.1073/pnas.0907043106.PubMedPubMed CentralView ArticleGoogle Scholar
- 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-1392. 10.1093/molbev/msq322.PubMedView ArticleGoogle Scholar
- Hughes AL: Evolution of adaptive phenotypic traits without positive Darwinian selection. Heredity. 2012, 108: 347-353. 10.1038/hdy.2011.97.PubMedPubMed CentralView ArticleGoogle Scholar
- Brakefield P, Gates J, Keys D, Kesbeke F, Wijngaarden P, Monteiro A, French V, Carroll S: Development, plasticity and evolution of butterfly eyespot patterns. Nature. 1996, 384: 236-242. 10.1038/384236a0.PubMedView ArticleGoogle Scholar
- Abouheif E, Wray GA: Evolution of the gene network underlying wing polyphenism in ants. Science. 2002, 297: 249-252. 10.1126/science.1071468.PubMedView ArticleGoogle Scholar
- Moczek AP, Nagy LM: Diverse developmental mechanisms contribute to different levels of diversity in horned beetles. Evol Dev. 2005, 7: 175-185. 10.1111/j.1525-142X.2005.05020.x.PubMedView ArticleGoogle Scholar
- Toth AL, Robinson GE: Evo-devo and the evolution of social behavior. Trends Genet. 2007, 23: 334-341. 10.1016/j.tig.2007.05.001.PubMedView ArticleGoogle Scholar
- Amdam GV, Norberg K, Fondrk MK, Page RE: Reproductive ground plan may mediate colony-level selection effects on individual foraging behavior in honey bees. Proc Natl Acad Sci USA. 2004, 101: 11350-11355. 10.1073/pnas.0403073101.PubMedPubMed CentralView ArticleGoogle Scholar
- Johnson BR, Linksvayer TA: Deconstructing the superorganism: social physiology, groundplans, and sociogenomics. Quart Rev Biol. 2010, 85: 57-79. 10.1086/650290.PubMedView ArticleGoogle Scholar
- Suzuki Y, Nijhout HF: Evolution of a polyphenism by genetic accommodation. Science. 2006, 311: 650-652. 10.1126/science.1118888.PubMedView ArticleGoogle Scholar
- Hunt BG, Ometto L, Wurm Y, Shoemaker DD, Soojin VY, Keller L, Goodisman MAD: Relaxed selection is a precursor to the evolution of phenotypic plasticity. Proc Natl Acad Sci USA. 2011, 108: 15936e15941-View ArticleGoogle Scholar
- Hunt BG, Wyder S, Elango N, Werren JH, Zdobnov EM, Yi SV, Goodisman MAD: Sociality is linked to rates of protein evolution in a highly social insect. Mol Biol Evol. 2010, 27: 497-500. 10.1093/molbev/msp225.PubMedView ArticleGoogle Scholar
- Ding Y, ZHou Q, Wang W: Origins of new genes and evolution of their novel functions. Ann Rev Ecol Syst. 2012, 43: 345-363. 10.1146/annurev-ecolsys-110411-160513.View ArticleGoogle Scholar
- Carvunis AR, Rolland T, Wapinski I, Calderwood MA, Yildirim MA, Hidalgo A, Barbette J, Santhanam B, Brar GA, Simonis N, Charloteaux B, Weissman JS, Regev A, Thierry-mieg N, Cusick ME, Vidal M: Proto-genes and de novo gene birth. Nature. 2012, 487: 370-374. 10.1038/nature11184.PubMedPubMed CentralView ArticleGoogle Scholar
- Dai H, Chen Y, Chen S, Mao Q, Kennedy D, Landback P, Eyre-Walker A, Du W, Long M: The evolution of courtship behaviors through the origination of a new gene in Drosophila. Proc Natl Acad Sci USA. 2008, 105: 7478-7483. 10.1073/pnas.0800693105.PubMedPubMed CentralView ArticleGoogle Scholar
- Khalturin K, Anton-Erxleben F, Sassmann S, Wittlieb , Hemmrich G, Bosch TCG: A novel gene family controls species-specific morphological traits in hydra. PLoS Biol. 2008, 6: e278-10.1371/journal.pbio.0060278.PubMedPubMed CentralView ArticleGoogle Scholar
- Fry BG, Roelants K, Winter K, Hodgson WC, Griesman L, Kwok HF, Scanlon D, Karas J, Shaw C, Wong L, Norman JA: Novel venom proteins produced by differential domain-expression strategies in beaded lizards and gila monsters (genus Heloderma). Mol Biol Evol. 2010, 27: 395-407. 10.1093/molbev/msp251.PubMedView ArticleGoogle Scholar
- Johnson BR, Tsutsui ND: Taxonomically restricted genes are associated with the evolution of sociality in the honey bee. BMC Genomics. 2011, 12: 164-10.1186/1471-2164-12-164.PubMedPubMed CentralView ArticleGoogle Scholar
- Woodard SH, Fischman BJ, Venkat A, Hudson ME, Varala K, Cameron SA, Clark AG, Robinson GE: Genes involved in convergent evolution of eusociality in bees. Proc Natl Acad Sci USA. 2011, 108: 7472-7477. 10.1073/pnas.1103457108.PubMedPubMed CentralView ArticleGoogle Scholar
- Fischman BJ, Woodard SH, Robinson GE: Molecular evolutionary analyses of insect societies. Proc Natl Acad Sci USA. 2011, 108: 10847-10854. 10.1073/pnas.1100301108.PubMedPubMed CentralView ArticleGoogle Scholar
- Turillazzi S, West-Eberhard MJ: Natural History and Evolution of Paper Wasps. 1996, Oxford: Oxford University PressGoogle Scholar
- Hunt JH, Amdam GV: Bivoltinism as an antecedent to eusociality in the paper wasp genus Polistes. Science. 2005, 308: 264-267. 10.1126/science.1109724.PubMedPubMed CentralView ArticleGoogle Scholar
- Carpenter JM: Phylogeny and biogeography of Polistes. Natural History and Evolution of Paper Wasps. Edited by: Turillazzi S, West-Eberhard MJ. 1996, Oxford: Oxford University Press, 18-57.Google Scholar
- Leadbeater E, Carruthers JM, Green JP, Rosser NS, Field JP: Nest inheritance is the missing source of direct fitness in a primitively eusocial insect. Science. 2011, 333: 874-876. 10.1126/science.1205140.PubMedView ArticleGoogle Scholar
- Hughes CR, Beck MO, Strassmann JE: Queen succession in the social wasp, Polistes annularis. Ethology. 1987, 76: 124-132.View ArticleGoogle Scholar
- Sumner S, Kelstrup H, Fanelli D: Reproductive constraints, direct fitness and indirect fitness explain helping behavior in the primitively eusocial wasp, Polistes canadensis. Proc R Soc Lond B. 2010, 277: 1721-1728. 10.1098/rspb.2009.2289.View ArticleGoogle Scholar
- Savard J, Tautz D, Richards S, Weinstock GM, Gibbs RA, Werren JH, Tettelin H, Lercher MJ: Phylogenomic analysis reveals bees and wasps (Hymenoptera) at the base of the radiation of Holometabolous insects. Genome Res. 2006, 16: 1334-1338. 10.1101/gr.5204306.PubMedPubMed CentralView ArticleGoogle Scholar
- Heraty J, Ronquist F, Carpenter JM, Hawks D, Schulmeister S, Dowling AP, Murray D, Munro J, Wheeler WC, Schiff N, Sharkey M: Evolution of the hymenopteran megaradiation. Mol Phylo Evol. 2011, 60: 73-88. 10.1016/j.ympev.2011.04.003.View ArticleGoogle Scholar
- Smith CR, Smith CD, Robertson HM, Helmkampf M, Zimin A, Yandell M, Holt C, Hu H, Abouheif E, Benton R, Cash E, Croset V, Currie CR, Elhaik E, Elsik CG, Fave M-J, Fernandes V, Gibson JD, Graur D, Gronenberg W, Grubbs KJ, Hagen DE, Viniegra ASI, Johnson BR, Johnson RM, Khila A, Kim JW, Mathis KA, Munoz-Torres MC, Murphy MC, et al: Draft genome of the red harvester ant Pogonomyrmex barbatus. Proc Natl Acad Sci USA. 2011, 108: 5667-5672. 10.1073/pnas.1007901108.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith CD, Zimin A, Holt C, Abouheif E, Benton R, Cash E, Croset V, Currie CR, Elhaik E, Elsik CG, Fave M-J, Fernandes V, Gadau , Gibson JD, Graur D, Grubbs KJ, Hagen DE, Helmkampf M, Holley J-A, Hu H, Viniegra ASI, Johnson BR, Johnson RM, Khila A, Kim JW, Laird J, Mathis KA, Moeller JA, Munoz-Torres MC, Murphy MC, et al: Draft genome of the globally widespread and invasive Argentine ant (Linepithema humile). Proc Natl Acad Sci USA. 2011, 108: 5673-5678. 10.1073/pnas.1008617108.PubMedPubMed CentralView ArticleGoogle Scholar
- Eddy SR: Accelerated profile HMM searches. PloS Comput Biol. 2011, 7: e1002195-10.1371/journal.pcbi.1002195.PubMedPubMed CentralView ArticleGoogle Scholar
- Arrial RT, Togawa RC, Brigido MdeM: Screening non-coding RNAs in transcriptomes from neglected species using PORTRAIT: case study of the pathogenic fungus Paracoccidioides brasiliensis. BMC Bioinformat. 2009, 10: 239-10.1186/1471-2105-10-239.View ArticleGoogle Scholar
- Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, Guernec G, Martin D, Merkel A, Knowles DG, Lagarde J, Veeravalli L, Ruan X, Ruan Y, Lassmann T, Carninci P, Brown JB, Lipovich L, Gonzalez JM, Thomas M, Davis CA, Shiekhattar R, Gingeras TR, Hubbard TJ, Notredame C, Harrow J, Guigó R, et al: The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012, 22: 1775-1789. 10.1101/gr.132159.111.PubMedPubMed CentralView ArticleGoogle Scholar
- Davis R, Baldauf S, Mayhew P: The origins of species richness in the Hymenoptera: insights from a family-level supertree. BMC Evol Biol. 2011, 10: 109-View ArticleGoogle Scholar
- Pilgrim EK, von Dohlen CD, Pitts JP: Molecular phylogenetics of Vespoidea indicate paraphyly of the superfamily and novel relationships of its component families and subfamilies. Zool Script. 2008, 37: 539-560. 10.1111/j.1463-6409.2008.00340.x.View ArticleGoogle Scholar
- Peters R, Meyer B, Krogmann L, Borner J, Meusemann K, Schutte K, Niehuis O, Misof B: The taming of an impossible child: a standardized all-in approach to the phylogeny of Hymenoptera using public database sequences. BMC Biol. 2011, 9: 55-10.1186/1741-7007-9-55.PubMedPubMed CentralView ArticleGoogle Scholar
- Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, Sorensen MV, Haddock SHD, Schmidt-Rhaesa A, Okusu A, Kristensen RM, Wheeler WC, Martindale MQ, Giribet G, et al: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature. 2008, 452: 745-749. 10.1038/nature06614.PubMedView ArticleGoogle Scholar
- Huerta-Cepas J, Marcet-Houben M, Pignatelli M, Moya A, Gabaldon T: The pea aphid phylome: a complete catalogue of evolutionary histories and arthropod orthology and paralogy relationships for Acyrthosiphon pisum genes. Insect Mol Biol. 2010, 19 (Suppl 2): 13-21.PubMedView ArticleGoogle Scholar
- Wilson EO, Holldobler B: The rise of the ants: A phylogenetic and ecological explanation. Proc Natl Acad Sci USA. 2005, 102: 7411-7414. 10.1073/pnas.0502264102.PubMedPubMed CentralView ArticleGoogle Scholar
- Tarazona S, García-Alcalde F, Dopazo J, Ferrer A, Conesa A: Differential expression in RNA-seq: A matter of depth. Genome Res. 2012, 21: 2213-2223.View ArticleGoogle Scholar
- Bourke AFG: Principles of Social Evolution. 2011, Oxford: Oxford University PressView ArticleGoogle Scholar
- Rex CS, Gavin CF, Rubio MD, Kramar EA, Chen LY, Jia Y, Huganir RL, Muzyczka N, Gall CM, Miller CA, Lynch G, Rumbaugh G: Myosin IIb regulates actin dynamics during synaptic plasticity and memory formation. Neuron. 2010, 67: 603-617. 10.1016/j.neuron.2010.07.016.PubMedPubMed CentralView ArticleGoogle Scholar
- 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-10.1186/1471-213X-7-70.PubMedPubMed CentralView ArticleGoogle Scholar
- Guttman M, Amit I, Garber M, French C, Lin MF, Feldser D, Huarte M, Zuk O, Carey BW, Cassady JP, Cabili MN, Jaenisch R, Mikkelsen TS, Jacks T, Hacohen N, Bernstein BE, Kellis M, Regev A, Rinn JL, Lander ES: Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 2009, 458: 223-227. 10.1038/nature07672.PubMedPubMed CentralView ArticleGoogle Scholar
- Hughes WOH, Oldroyd BP, Beekman M, Ratnieks FLW: Ancestral monogamy shows kin selection is key to the evolution of eusociality. Science. 2008, 320: 1213-1216. 10.1126/science.1156108.PubMedView ArticleGoogle Scholar
- Ament SA, Corona M, Pollock HS, Robinson GE: Insulin signaling is involved in the regulation of worker division of labor in honey bee colonies. Proc Natl Acad Sci USA. 2008, 105: 4226-4231. 10.1073/pnas.0800630105.PubMedPubMed CentralView ArticleGoogle Scholar
- Giray T, Giovanetti M, West-Eberhard MJ: Juvenile hormone, reproduction, and worker behavior in the neotropical social wasp Polistes canadensis. Proc Natl Acad Sci USA. 2005, 102: 3330-3335. 10.1073/pnas.0409560102.PubMedPubMed CentralView ArticleGoogle Scholar
- Gadagkar R: The evolution of caste polymorphism in social insects: genetic release followed by diversifying evolution. J Genet. 1997, 76: 167-179. 10.1007/BF02932215.View ArticleGoogle Scholar
- Robinson GE, Ben-Shahar Y: Social behavior and comparative genomics: new genes or new gene regulation?. Genes Brain Behav. 2002, 1: 197-203. 10.1034/j.1601-183X.2002.10401.x.PubMedView ArticleGoogle Scholar
- Linksvayer TA, Wade MJ: The evolutionary origin and elaboration of sociality in the aculeate hymenoptera: Maternal effects, sib-social effects, and heterochrony. Quart Rev Biol. 2005, 80: 317-336. 10.1086/432266.PubMedView ArticleGoogle Scholar
- Colgan TC, Carolan JC, Bridgett SJ, Sumner S, Blaxter ML, Brown MJF: Polyphenism in social insects: Insights from a transcriptome-wide analysis of gene expression in the life stages of the key pollinator, Bombus terrestris. BMC Genomics. 2011, 12: 623-10.1186/1471-2164-12-623.PubMedPubMed CentralView ArticleGoogle Scholar
- Lyko F, Foret S, Kucharski R, Wolf S, Falckenhayn C, Maleszka R: The honey bee epigenomes: differential methylation of brain DNA in queens and workers. PLoS Biol. 2010, 8: e1000506-10.1371/journal.pbio.1000506.PubMedPubMed CentralView ArticleGoogle Scholar
- Edgar R: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.PubMedPubMed CentralView ArticleGoogle Scholar
- Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T: trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009, 25: 1972-1973. 10.1093/bioinformatics/btp348.PubMedPubMed CentralView ArticleGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59: 307-321. 10.1093/sysbio/syq010.PubMedView ArticleGoogle Scholar
- Le SQ, Gascuel O: An improved general amino acid replacement matrix. Molec Biol Evol. 2008, 25: 1307-1320. 10.1093/molbev/msn067.PubMedView ArticleGoogle Scholar
- Shimodaira H, Hasegawa M: CONSEL: for assessing the confidence of phylogenetic tree selection. Bioinformatics. 2001, 17: 1246-1247. 10.1093/bioinformatics/17.12.1246.PubMedView ArticleGoogle Scholar
- Lartillot N, Philippe H: Bayesian phylogenetic software based on mixture models. Molec Biol Evol. 2004, 21: 1095-1109. 10.1093/molbev/msh112.PubMedView ArticleGoogle Scholar
- Data and Datasets. [http://genome.crg.es/~pferreira/pcandata/pcan.htm]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.