Differential Selection on Caste-Associated Genes in a Subterranean Termite

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.