RESUMEN
Analyzing the information-rich content of RNA can help uncover genetic events associated with social insect castes or other social polymorphisms. Here, we exploit a series of cDNA libraries previously derived from whole-body tissue of different castes as well as from three behaviourally distinct populations of the Eastern subterranean termite Reticulitermes flavipes. We found that the number (~0.5 M) of single nucleotide variants (SNVs) was roughly equal between nymph, worker and soldier caste libraries, but dN/dS (ratio of nonsynonymous to synonymous substitutions) analysis suggested that some of these variants confer a caste-specific advantage. Specifically, the dN/dS ratio was high (~4.3) for genes expressed in the defensively specialized soldier caste, relative to genes expressed by other castes (~1.7−1.8) and regardless of the North American population (Toronto, Raleigh, Boston) from which the castes were sampled. The populations, meanwhile, did show a large difference in SNV count but not in the manner expected from known demographic and behavioural differences; the highly invasive unicolonial population from Toronto was not the least diverse and did not show any other unique substitution patterns, suggesting any past bottleneck associated with invasion or with current unicoloniality has become obscured at the RNA level. Our study raises two important hypotheses relevant to termite sociobiology. First, the positive selection (dN/dS > 1) inferred for soldier-biased genes is presumably indirect and of the type mediated through kin selection, and second, the behavioural changes that accompany some social insect urban invasions (i.e., 'unicoloniality') may be detached from the loss-of-diversity expected from invasion bottlenecks.
RESUMEN
The highly organized societies of the Western honey bee Apis mellifera feature a highly reproductive queen at the center of attention and a large cohort of daughters that suppress their own reproduction to help rear more sisters, some of whom become queens themselves. This reproductive altruism is peculiar because in theory it evolves via indirect selection on genes for altruism that are expressed in the sterile workers but not in the reproductive queens. In this study we attempt to situate lists of genes previously implicated in queenright worker sterility into a broader regulatory framework. To do so we use a model bee brain transcriptional regulatory network as a template to infer how sets of genes responsive to ovary-suppressing queen pheromone are functionally interconnected over the model's topology. We predict that genes jointly involved in the regulation of worker sterility should be tightly networked, relative to genes whose functions are unrelated to each other. We find that sets of mapped genes - ranging in size from 17 to 250 - are well dispersed across the network's substructural scaffolds, suggesting that ovary de-activation involves genes that reside within more than one transcriptional regulatory module. For some sets, however, this dispersion is biased into certain areas of the network's substructure. Our analysis identifies the regions enriched for sterility genes and likewise identifies local hub genes that are presumably critical to subnetwork function. Our work offers a glimpse into the gene regulatory context of honey bee worker sterility and uses this context to identify new candidate gene targets for functional analysis. Finally, to the extent that any sterility-related modules identified here have evolved via selection for worker altruism, we can assume that this selection was indirect and of the type specifically invoked by inclusive fitness theory.