Caste‐biases in gene expression are specific to developmental stage in the ant Formica exsecta

Understanding how a single genome creates and maintains distinct phenotypes is a central goal in evolutionary biology. Social insects are a striking example of co‐opted genetic backgrounds giving rise to dramatically different phenotypes, such as queen and worker castes. A conserved set of molecular pathways, previously envisioned as a set of ‘toolkit’ genes, has been hypothesized to underlie queen and worker phenotypes in independently evolved social insect lineages. Here, we investigated the toolkit from a developmental point of view, using RNA‐Seq to compare caste‐biased gene expression patterns across three life stages (pupae, emerging adult and old adult) and two female castes (queens and workers) in the ant Formica exsecta. We found that the number of genes with caste‐biased expression increases dramatically from pupal to old adult stages. This result suggests that phenotypic differences between queens and workers at the pupal stage may derive from a relatively low number of caste‐biased genes, compared to higher number of genes required to maintain caste differences at the adult stage. Gene expression patterns were more similar among castes within developmental stages than within castes despite the extensive phenotypic differences between queens and workers. Caste‐biased expression was highly variable among life stages at the level of single genes, but more consistent when gene functions (gene ontology terms) were investigated. Finally, we found that a large part of putative toolkit genes were caste‐biased at least in some life stages in F. exsecta, and the caste‐biases, but not their direction, were more often shared between F. exsecta and other ant species than between F. exsecta and bees. Our results indicate that gene expression should be examined across several developmental stages to fully reveal the genetic basis of polyphenisms.     

Making a queen: an epigenetic analysis of the robustness of the honeybee (Apis mellifera) queen developmental pathway

Specialized castes are considered a key reason for the evolutionary and ecological success of the social insect lifestyle. The most essential caste distinction is between the fertile queen and the sterile workers. Honeybee (Apis mellifera) workers and queens are not genetically distinct, rather these different phenotypes are the result of epigenetically regulated divergent developmental pathways. This is an important phenomenon in understanding the evolution of social insect societies. Here, we studied the genomic regulation of the worker and queen developmental pathways, and the robustness of the pathways by transplanting eggs or young larvae to queen cells. Queens could be successfully reared from worker larvae transplanted up to 3 days age, but queens reared from older worker larvae had decreased queen body size and weight compared with queens from transplanted eggs. Gene expression analysis showed that queens raised from worker larvae differed from queens raised from eggs in the expression of genes involved in the immune system, caste differentiation, body development and longevity. DNA methylation levels were also higher in 3‐day‐old queen larvae raised from worker larvae compared with that raised from transplanted eggs identifying a possible mechanism stabilizing the two developmental paths. We propose that environmental (nutrition and space) changes induced by the commercial rearing practice result in a suboptimal queen phenotype via epigenetic processes, which may potentially contribute to the evolution of queen–worker dimorphism. This also has potentially contributed to the global increase in honeybee colony failure rates.     

Evolution of gene expression in fire ants: the effects of developmental stage, caste, and species

Ants provide remarkable examples of equivalent genotypes developing into divergent and discrete phenotypes. Diploid eggs can develop either into queens, which specialize in reproduction, or workers, which participate in cooperative tasks such as building the nest, collecting food, and rearing the young. In contrast, the differentiation between males and females generally depends upon whether eggs are fertilized, with fertilized (diploid) eggs giving rise to females and unfertilized (haploid) eggs giving rise to males. To obtain a comprehensive picture of the relative contributions of gender (sex), caste, developmental stage, and species divergence to gene expression evolution, we investigated gene expression patterns in pupal and adult queens, workers, and males of two species of fire ants, Solenopsis invicta and S. richteri. Microarray hybridizations revealed that variation in gene expression profiles is influenced more by developmental stage than by caste membership, sex, or species identity. The second major contributor to variation in gene expression was the combination of sex and caste. Although workers and queens share equivalent diploid nuclear genomes, they have highly distinctive patterns of gene expression in both the pupal and the adult stages, as might be expected given their extraordinary level of phenotypic differentiation. Overall, the difference in the proportion of differentially expressed genes was greater between workers and males than between workers and queens or queens and males, consistent with the fact that workers and males share neither gender nor reproductive capability. Moreover, between-species comparisons revealed that the greatest difference in gene expression patterns occurred in adult workers, a finding consistent with the fact that adult workers most directly experience the distinct external environments characterizing the different habitats occupied by the two species. Thus, much of the evolution of gene expression in ants may occur in the worker caste, despite the fact that these individuals are largely or completely sterile. Analyses of gene expression evolution revealed a combination of positive selection and relaxation of stabilizing selection as important factors driving the evolution of such genes.     

Differential gene expression and phenotypic plasticity in behavioural castes of the primitively eusocial wasp, Polistes canadensis

Understanding how a single genome can produce a variety of different phenotypes is of fundamental importance in evolutionary and developmental biology. One of the most striking examples of phenotypic plasticity is the female caste system found in eusocial insects, where variation in reproductive (queens) and non-reproductive (workers) phenotypes results in a broad spectrum of caste types, ranging from behavioural through to morphological castes. Recent advances in genomic techniques allow novel comparisons on the nature of caste phenotypes to be made at the level of the genes in organisms for which there is little genome information, facilitating new approaches in studying social evolution and behaviour. Using the paper wasp Polistes canadensis as a model system, we investigated for the first time how behavioural castes in primitively eusocial insect societies are associated with differential expression of shared genes. We found that queens and newly emerged females express gene expression patterns that are distinct from each other whilst workers generally expressed intermediate patterns, as predicted by Polistes biology. We compared caste-associated genes in P. canadensis with those expressed in adult queens and workers of more advanced eusocial societies, which represent four independent origins of eusociality. Nine genes were conserved across the four taxa, although their patterns of expression and putative functions varied. Thus, we identify several genes that are putatively of evolutionary importance in the molecular biology that underlies a number of caste systems of independent evolutionary origin.     

Gene duplication and the evolution of phenotypic diversity in insect societies

Gene duplication is an important evolutionary process thought to facilitate the evolution of phenotypic diversity. We investigated if gene duplication was associated with the evolution of phenotypic differences in a highly social insect, the honeybee Apis mellifera. We hypothesized that the genetic redundancy provided by gene duplication could promote the evolution of social and sexual phenotypes associated with advanced societies. We found a positive correlation between sociality and rate of gene duplications across the Apoidea, indicating that gene duplication may be associated with sociality. We also discovered that genes showing biased expression between A. mellifera alternative phenotypes tended to be found more frequently than expected among duplicated genes than singletons. Moreover, duplicated genes had higher levels of caste‐, sex‐, behavior‐, and tissue‐biased expression compared to singletons, as expected if gene duplication facilitated phenotypic differentiation. We also found that duplicated genes were maintained in the A. mellifera genome through the processes of conservation, neofunctionalization, and specialization, but not subfunctionalization. Overall, we conclude that gene duplication may have facilitated the evolution of social and sexual phenotypes, as well as tissue differentiation. Thus this study further supports the idea that gene duplication allows species to evolve an increased range of phenotypic diversity.     

A simple genetic basis for complex social behaviour mediates widespread gene expression differences

A remarkable social polymorphism is controlled by a single Mendelian factor in the fire ant Solenopsis invicta. A genomic element marked by the gene Gp‐9 determines whether workers tolerate one or many fertile queens in their colony. Gp‐9 was recently shown to be part of a supergene with two nonrecombining variants, SB and Sb. SB/SB and SB/Sb queens differ in how they initiate new colonies, and in many physiological traits, for example odour and maturation rate. To understand how a single genetic element can affect all these traits, we used a microarray to compare gene expression patterns between SB/SB and SB/Sb queens of three different age classes: 1‐day‐old unmated queens, 11‐day‐old unmated queens and mated, fully reproductive queens collected from mature field colonies. The number of genes that were differentially expressed between SB/SB and SB/Sb queens of the same age class was smallest in 1‐day‐old queens, maximal in 11‐day‐old queens and intermediate in reproductive queens. Gene ontology analysis showed that SB/SB queens upregulate reproductive genes faster than SB/Sb queens. For all age classes, genes inside the supergene were overrepresented among the differentially expressed genes. Consistent with the hypothesized greater number of transposons in the Sb supergene, 13 transposon genes were upregulated in SB/Sb queens. Viral genes were also upregulated in SB/Sb mature queens, consistent with the known greater parasite load in colonies headed by SB/Sb queens compared with colonies headed by SB/SB queens. Eighteen differentially expressed genes between reproductive queens were involved in chemical signalling. Our results suggest that many genes in the supergene are involved in regulating social organization and queen phenotypes in fire ants.     

Transcriptomic analysis of epigenetic modification genes in the termite Reticulitermes speratus

Eusocial insects display a caste system in which different castes are morpho-logically and physiologically specialized for different tasks.Recent studies have revealed that epigenetic modifications,including DNA methylation and histone modification,me-diate caste determination and differentiation,longevity,and polyethism in eusocial insects.Although there has been a growing interest in the relationship between epigenetic mech-anisms and phenotypic plasticity in termites,there is ltte information about differential expression levels among castes and expression sites for these genes in termites.Here we show royaltissuc-specific expression of epigenetic modification genes in the termite Reticulitermes speratus.Using RNA-seq,we identified 74 genes,including three DNA methyltransferases,seven sirtuins,48 Trithorax group proteins,and 16 Polycomb group proteins.Among these genes,15 showed king-specific expression,and 52 showed age-dependent differential expression in kings and queens.Quantitative real-time PCR revealed that DNA methyltransferase 3 is expressed specifically in the king's testis and fat body,whereas some histone modification genes are remarkably expressed in the king's testis and queen's ovary.These findings imply that epigenetic modification plays important roles in the gamete production process in termite kings and queens.     

Relaxed genetic constraint is ancestral to the evolution of phenotypic plasticity

Phenotypic plasticity--the capacity of a single genotype to produce different phenotypes in response to varying environmental conditions--is widespread. Yet, whether, and how, plasticity impacts evolutionary diversification is unclear. According to a widely discussed hypothesis, plasticity promotes rapid evolution because genes expressed differentially across different environments (i.e., genes with "biased" expression) experience relaxed genetic constraint and thereby accumulate variation faster than do genes with unbiased expression. Indeed, empirical studies confirm that biased genes evolve faster than unbiased genes in the same genome. An alternative hypothesis holds, however, that the relaxed constraint and faster evolutionary rates of biased genes may be a precondition for, rather than a consequence of, plasticity's evolution. Here, we evaluated these alternative hypotheses by characterizing evolutionary rates of biased and unbiased genes in two species of frogs that exhibit a striking form of phenotypic plasticity. We also characterized orthologs of these genes in four species of frogs that had diverged from the two plastic species before the plasticity evolved. We found that the faster evolutionary rates of biased genes predated the evolution of the plasticity. Furthermore, biased genes showed greater expression variance than did unbiased genes, suggesting that they may be more dispensable. Phenotypic plasticity may therefore evolve when dispensable genes are co-opted for novel function in environmentally induced phenotypes. Thus, relaxed genetic constraint may be a cause--not a consequence--of the evolution of phenotypic plasticity, and thereby contribute to the evolution of novel traits.     

Soldier‐biased gene expression in a subterranean termite implies functional specialization of the defensive caste

In a termite colony, reproduction is typically monopolized by a small number of sexuals that are supported by reproductively altruistic soldiers and workers. We expect caste differentiation to be associated with clear‐cut differences in gene expression, and for these differences to reflect caste function and development. Here, we use RNA‐Sequencing to compare the gene expression profiles of sexual nymphs and two non‐reproductive helper castes (i.e., workers and soldiers) of the Eastern subterranean termite Reticulitermes flavipes. We found that of n = 93 genes that are strictly expressed as a function of caste, a majority (78%) show a soldier‐specific pattern. This conspicuous soldier‐bias in genome‐wide expression suggests that this defensively specialized caste is functionally well‐differentiated from both the reproductive and the other non‐reproductive caste of this species, despite a shared developmental program with workers. Gene ontology analysis supports the notion of functional specialization by soldiers, as soldier‐biased gene sets are enriched for novel biological processes. Whether this pattern reflects ancient or more recent bouts of selection for caste novelty at the gene‐regulatory level is not known, but because soldiers are sterile and thus have no direct fitness, any selection for novelty must have been mediated indirectly, through reproducing relatives.     

Gene expression patterns associated with caste and reproductive status in ants: worker‐specific genes are more derived than queen‐specific ones

Variation in gene expression leads to phenotypic diversity and plays a central role in caste differentiation of eusocial insect species. In social Hymenoptera, females with the same genetic background can develop into queens or workers, which are characterized by divergent morphologies, behaviours and lifespan. Moreover, many social insects exhibit behaviourally distinct worker castes, such as brood‐tenders and foragers. Researchers have just started to explore which genes are differentially expressed to achieve this remarkable phenotypic plasticity. Although the queen is normally the only reproductive individual in the nest, following her removal, young brood‐tending workers often develop ovaries and start to reproduce. Here, we make use of this ability in the ant Temnothorax longispinosus and compare gene expression patterns in the queens and three worker castes along a reproductive gradient. We found the largest expression differences between the queen and the worker castes (~2500 genes) and the smallest differences between infertile brood‐tenders and foragers (~300 genes). The expression profile of fertile workers is more worker‐like, but to a certain extent intermediate between the queen and the infertile worker castes. In contrast to the queen, a high number of differentially expressed genes in the worker castes are of unknown function, pointing to the derived status of hymenopteran workers within insects.     

Genome-wide analysis reveals differences in brain gene expression patterns associated with caste and reproductive status in honey bees (Apis mellifera)

A key characteristic of eusocial species is reproductive division of labour. Honey bee colonies typically have a single reproductive queen and thousands of sterile workers. Adult queens differ dramatically from workers in anatomy, physiology, behaviour and lifespan. Young female workers can activate their ovaries and initiate egg laying; these 'reproductive' workers differ from sterile workers in anatomy, physiology, and behaviour. These differences, however, are on a much smaller scale than those observed between the queen and worker castes. Here, we use microarrays to monitor expression patterns of several thousand genes in the brains of same-aged virgin queens, sterile workers, and reproductive workers. We found large differences in expression between queens and both worker groups (~2000 genes), and much smaller differences between sterile and reproductive workers (221 genes). The expression patterns of these 221 genes in reproductive workers are more queen-like, and may represent a core group of genes associated with reproductive physiology. Furthermore, queens and reproductive workers preferentially up-regulate genes associated with the nurse bee behavioural state, which supports the hypothesis of an evolutionary link between worker division of labour and molecular pathways related to reproduction. Finally, several functional groups of genes associated with longevity in other species are significantly up-regulated in queens. Identifying the genes that underlie the differences between queens, sterile workers, and reproductive workers will allow us to begin to characterize the molecular mechanisms underlying the evolution of social behaviour and large-scale remodelling of gene networks associated with polyphenisms.     

mRNA expression and DNA methylation in three key genes involved in caste differentiation in female honeybees(Apis mellifera)

In honeybee(Apis mellifera)colonies,queens and workers are alternative forms of the adult female honeybee that develop from genetically identical zygotes but that depend on differential nourishment.Queens and workers display distinct morphologies,anatomies and behavior,better known as caste differentiation.Despite some basic insights,the exact mechanism responsible for this phenomenon,especially at the molecular level,remains unclear although some progress has been achieved.In this study,we examined mRNA levels of the TOR(target of rapamycin)and Dnmt3(DNA methyltransferase 3)genes,closely related to caste differentiation in honeybees.We also investigated mRNA expression of the S6K(similar to RPS6-p70-protein kinase)gene linked closely to organismal growth and development in queen and worker larvae(1-day and 3-day old).Last,we investigated the methylation status of these three genes in corresponding castes.We found no difference in mRNA expression for the three genes between 1st instar queen and worker larvae;however,3rd instar queen larvae had a higher level of TOR mRNA than worker larvae.Methylation levels of all three genes were lower in queen larvae than worker larvae but the differences were not statistically significant.These findings provide basic data for broadening our understanding of caste differentiation in female honeybees.     

Patterns of DNA methylation in development, division of labor and hybridization in an ant with genetic caste determination

Background

DNA methylation is a common regulator of gene expression, including acting as a regulator of developmental events and behavioral changes in adults. Using the unique system of genetic caste determination in Pogonomyrmex barbatus, we were able to document changes in DNA methylation during development, and also across both ancient and contemporary hybridization events.

Methodology/Principal Findings

Sodium bisulfite sequencing demonstrated in vivo methylation of symmetric CG dinucleotides in P. barbatus. We also found methylation of non-CpG sequences. This validated two bioinformatics methods for predicting gene methylation, the bias in observed to expected ratio of CpG dinucleotides and the density of CpG/TpG single nucleotide polymorphisms (SNP). Frequencies of genomic DNA methylation were determined for different developmental stages and castes using ms-AFLP assays. The genetic caste determination system (GCD) is probably the product of an ancestral hybridization event between P. barbatus and P. rugosus. Two lineages obligately co-occur within a GCD population, and queens are derived from intra-lineage matings whereas workers are produced from inter-lineage matings. Relative DNA methylation levels of queens and workers from GCD lineages (contemporary hybrids) were not significantly different until adulthood. Virgin queens had significantly higher relative levels of DNA methylation compared to workers. Worker DNA methylation did not vary among developmental stages within each lineage, but was significantly different between the currently hybridizing lineages. Finally, workers of the two genetic caste determination lineages had half as many methylated cytosines as workers from the putative parental species, which have environmental caste determination.

Conclusions/Significance

These results suggest that DNA methylation may be a conserved regulatory mechanism moderating division of labor in both bees and ants. Current and historic hybridization appear to have altered genomic methylation levels suggesting a possible link between changes in overall DNA methylation and the origin and regulation of genetic caste determination in P. barbatus.     

Expression profiles during honeybee caste determination

Background

Depending on their larval environment, female honeybees develop into either queens or workers. As in other polyphenisms, this developmental switch depends not on genomic differences between queens and workers but on the differential expression of entire suites of genes involved with larval fate. As such, this and other polyphenic systems can provide a novel tool for understanding how genomes and environmental conditions interact to produce different developmental trajectories. Here we use gene-expression profiles during honeybee caste determination to present the first genomic view of polyphenic development.

Results

Larvae raised as queens or workers differed greatly in their gene-expression patterns. Workers remained more faithful than queens to the expression profiles of younger, bipotential, larvae. Queens appeared to both downregulate many of the genes expressed by bipotential larvae and turn on a distinct set of caste-related genes. Queens overexpressed several metabolic enzymes, workers showed increased expression of a member of the cytochrome P450 family, hexameric storage proteins and dihydrodiol dehydrogenase, and young larvae overexpressed two putative heat-shock proteins (70 and 90 kDa), and several proteins related to RNA processing and translation.

Conclusions

Large differences in gene expression between queens and workers indicate that social insect castes have faced strong directional selection pressures. Overexpression of metabolic enzymes by queen-destined larvae appears to reflect the enhanced growth rate of queens during late larval development. Many of the differently expressed genes we identified have been tied to metabolic rates and cellular responses to hormones, a result consistent with known physiological differences between queen and worker larvae.