Genomic and transcriptomic analyses of the subterranean termite Reticulitermes speratus: Gene duplication facilitates social evolution

Genomic and transcriptomic analyses of a social hemipteran provide new insights into insect sociality

The origin of sociality represents one of the most important evolutionary transitions. Insect sociality evolved in some hemipteran aphids, which can produce soldiers and normal nymphs with distinct morphology and behaviour through parthenogenesis. The lack of genomic data resources has hindered the investigations into molecular mechanisms underlying their social evolution. Herein, we generated the first chromosomal-level genome of a social hemipteran (Pseudoregma bambucicola) with highly specialized soldiers and performed comparative genomic and transcriptomic analyses to elucidate the molecular signatures and regulatory mechanisms of caste differentiation. P. bambucicola has a larger known aphid genome of 582.2 Mb with an N50 length of 11.24 Mb, and about 99.6% of the assembly was anchored to six chromosomes with a scaffold N50 of 98.27 Mb. A total of 14,027 protein-coding genes were predicted and 37.33% of the assembly were identified as repeat sequences. The social evolution is accompanied by a variety of changes in genome organization, including expansion of gene families related to transcription factors, transposable elements, as well as species-specific expansions of certain sugar transporters and UGPases involved in carbohydrate metabolism. We also characterized large candidate gene sets linked to caste differentiation and found evidence of expression regulation and positive selection acting on energy metabolism and muscle structure, explaining the soldier-specific traits including morphological and behavioural specialization, developmental arrest and infertility. Overall, this study offers new insights into the molecular basis of social aphids and the evolution of insect sociality and also provides valuable data resources for further comparative and functional studies.

Protein localization and potential function of lipocalin in Reticulitermes speratus queens

To understand the mechanisms underlying social evolution and caste development in social insects, caste-specific organs and genes should be investigated. In the rhinotermitid termite, Reticulitermes speratus, the lipocalin gene RS008881, which encodes a protein transporter, is expressed in the ovarian accessory glands of primary queens. To obtain additional data on its expression and product localization, we conducted real-time quantitative polymerase chain reaction and protein assays using a peptide antibody. Gene expression analysis of the castes revealed that RS008881 was highly expressed in female primary and secondary reproductives. Further analysis of its expression during reproductive caste differentiation showed that its expression levels increased prior to molting into reproductive individuals, even during the winged imago (alates) stage. Western blotting and fluorescent immunohistochemical staining revealed that the RS008881 product was localized in the ovary as well as the eggshells produced by female reproductives. RS008881 may play a significant role in the reproductive biology of R. speratus; protein localization in both the ovary and eggshell suggests multiple functions related to embryo protection and potential pheromone interactions.

Functional redundancy of the three insulin receptors of cockroaches

Gene duplication is a fundamental evolutionary process which provides opportunities to acquire new gene functions. In the case of the insulin receptors (InRs) in cockroaches and close-related insects, two successive duplications determined the occurrence of three InR genes: InR2, InR1 and InR3, the last two forming a sister cluster to InR2. The biological role of each of the gene duplicates and whether they resulted from neofunctionalization or subfunctionalization is still unclear. The analysis of the sequences from different lineages did not detect positive selection as driving the divergence of InR1 and InR3, discarding neofunctionalization, and suggesting that there is no functional divergence between both gene copies. Using the cockroach Blattella germanica as a model, we have determined that BgInR2 is the gene with the highest expression levels in all the tissues analyzed, both in adult females and males, as well as in nymphs and embryos. BgInR3 is second in expression levels while BgInR1 is expressed at lower levels and only in some tissues. The selective depletion by RNAi of each of the three InRs, analyzed in terms of phenotype and fat body transcriptomic profiles, resulted in essentially redundant effects, with a magnitude approximately proportional to the level of expression of the respective InR. Therefore, the results indicate that the InR duplicates likely experienced a subfunctionalization process, by which the three InRs maintained similar functions but contributing to those functions proportionally to their expression levels.

Evolution of the sex-determination gene Doublesex within the termite lineage

The molecular mechanism of sex determination has long been considered conserved in insects. However, recent studies of hemimetabolous insects have challenged this notion. One notable example is termites. In Reticulitermes speratus, a homolog of sex determination gene, Doublesex (RsDsx), exhibits characteristics that are distinct from those of other insects, including sister-group cockroaches. It comprises a single exon, contains only doublesex/mab-3 DNA-binding domain (DM) but lacks a conserved oligomerization domain (OD), and exhibits transcriptional activity only in males. To investigate whether these characteristics are widespread within the termite lineage, we identified Dsx homologs in three different families. The absence of the conserved OD sequences was observed in all termite species examined, whereas the number of exons and expression patterns between sexes varied among families. Particularly, distinctive differences in Dsx were found in species from the Archotermopsidae and Kalotermitidae, both of which have a linear caste developmental pathway. Our findings indicate that diversification of Dsx structure and expression patterns may have contributed to ecological diversification, such as caste developmental pathways, within the termite lineage.

Major changes in domain arrangements are associated with the evolution of termites

Domains as functional protein units and their rearrangements along the phylogeny can shed light on the functional changes of proteomes associated with the evolution of complex traits like eusociality. This complex trait is associated with sterile soldiers and workers, and long-lived, highly fecund reproductives. Unlike in Hymenoptera (ants, bees, and wasps), the evolution of eusociality within Blattodea, where termites evolved from within cockroaches, was accompanied by a reduction in proteome size, raising the question of whether functional novelty was achieved with existing rather than novel proteins. To address this, we investigated the role of domain rearrangements during the evolution of termite eusociality. Analysing domain rearrangements in the proteomes of three solitary cockroaches and five eusocial termites, we inferred more than 5,000 rearrangements over the phylogeny of Blattodea. The 90 novel domain arrangements that emerged at the origin of termites were enriched for several functions related to longevity, such as protein homeostasis, DNA repair, mitochondrial activity, and nutrient sensing. Many domain rearrangements were related to changes in developmental pathways, important for the emergence of novel castes. Along with the elaboration of social complexity, including permanently sterile workers and larger, foraging colonies, we found 110 further domain arrangements with functions related to protein glycosylation and ion transport. We found an enrichment of caste-biased expression and splicing within rearranged genes, highlighting their importance for the evolution of castes. Furthermore, we found increased levels of DNA methylation among rearranged compared to non-rearranged genes suggesting fundamental differences in their regulation. Our findings indicate the importance of domain rearrangements in the generation of functional novelty necessary for termite eusociality to evolve.

Body part-specific development in termite caste differentiation: crosstalk between hormonal actions and developmental toolkit genes

In social insects, interactions among colony members trigger caste differentiation with morphological modifications. During caste differentiation in termites, body parts and caste-specific morphologies are modified during postembryonic development under endocrine controls such as juvenile hormone (JH) and ecdysone. In addition to endocrine factors, developmental toolkit genes such as Hox- and appendage-patterning genes also contribute to the caste-specific body part modifications. These toolkits are thought to provide spatial information for specific morphogenesis. During social evolution, the complex crosstalks between physiological and developmental mechanisms should be established, leading to the sophisticated caste systems. This article reviews recent studies on these mechanisms underlying the termite caste differentiation and addresses implications for the evolution of caste systems in termites.

Pervasive relaxed selection in termite genomes

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Candidate target genes of the male-specific expressed Doublesex in the termite Reticulitermes speratus

Eusocial insects such as termites, ants, bees, and wasps exhibit a reproductive division of labor. The developmental regulation of reproductive organ (ovaries and testes) is crucial for distinguishing between reproductive and sterile castes. The development of reproductive organ in insects is regulated by sex-determination pathways. The sex determination gene Doublesex (Dsx), encoding transcription factors, plays an important role in this pathway. Therefore, clarifying the function of Dsx in the developmental regulation of sexual traits is important to understand the social evolution of eusocial insects. However, no studies have reported the function of Dsx in hemimetabolous eusocial group termites. In this study, we searched for binding sites and candidate target genes of Dsx in species with available genome information as the first step in clarifying the function of Dsx in termites. First, we focused on the Reticulitermes speratus genome and identified 101 candidate target genes of Dsx. Using a similar method, we obtained 112, 39, and 76 candidate Dsx target genes in Reticulitermes lucifugus, Coptotermes formosanus, and Macrotermes natalensis, respectively. Second, we compared the candidate Dsx target genes between species and identified 37 common genes between R. speratus and R. lucifugus. These included several genes probably involved in spermatogenesis and longevity. However, only a few common target genes were identified between R. speratus and the other two species. Finally, Dsx dsRNA injection resulted in the differential expression of several target genes, including piwi-like protein and B-box type zinc finger protein ncl-1 in R. speratus. These results provide valuable resource data for future functional analyses of Dsx in termites.

Genomic signatures of eusocial evolution in insects

The genomes of eusocial insects allow the production and regulation of highly distinct phenotypes, largely independent of genotype. Although rare, eusociality has evolved convergently in at least three insect orders (Hymenoptera, Blattodea and Coleoptera). Despite such disparate origins, eusocial phenotypes show remarkable similarity, exhibiting long-lived reproductives and short-lived sterile workers and soldiers. In this article, we review current knowledge on genomic signatures of eusocial evolution. We confirm that especially an increased regulatory complexity and the adaptive evolution of chemical communication are common to several origins of eusociality. Furthermore, colony life itself can shape genomes of divergent taxa in a similar manner. Future research should be geared towards generating more high-quality genomic resources, especially in hitherto understudied clades, such as ambrosia beetles and termites. The application of more sophisticated tools such as machine learning techniques may allow the detection of more subtle convergent genomic footprints of eusociality.

Uncovering the Chemosensory System of a Subterranean Termite, Odontotermes formosanus (Shiraki) (Isoptera: Termitidae): Revealing the Chemosensory Genes and Gene Expression Patterns

Termites are eusocial insects. Chemical signals between colony members are crucial to the smooth running of colony operations, but little is known about their olfactory system and the roles played by various chemosensory genes in this process. Chemosensory genes are involved in basic olfactory perception in insects. Odontotermes formosanus (Shiraki) is one of the most damaging pests to agricultural crops, forests, and human-made structures. To better understand the olfactory system and the genes involved in olfactory processing in O. formosanus, we produced a transcriptome of worker termites. In this study, we identified 13 OforOBPs, 1 OforCSP, 15 OforORs, 9 OforGRs, and 4 OforSNMPs. Multiple sequence alignments were used in the phylogenetic study, which included data from other termite species and a wide variety of insect species. Moreover, we also investigated the mRNA expression levels using qRT-PCR. The significantly high expression levels of OforCSP1, OforOBP2, OforOR1, and OforSNMP1 suggest that these genes may play important roles in olfactory processing in termite social behavior, including caste differentiation, nestmate and non-nestmate discrimination, and the performance of colony operations among members. Our research establishes a foundation for future molecular-level functional studies of chemosensory genes in O. formosanus, which might lead to the identification of novel targets for termite integrated pest management.

Cell type specific polyploidization in the royal fat body of termite queens

Tissue-specific endopolyploidy is widespread among plants and animals and its role in organ development and function has long been investigated. In insects, the fat body cells of sexually mature females produce substantial amounts of egg yolk precursor proteins (vitellogenins) and exhibit high polyploid levels, which is considered crucial for boosting egg production. Termites are social insects with a reproductive division of labor, and the fat bodies of mature termite queens exhibit higher ploidy levels than those of other females. The fat bodies of mature termite queens are known to be histologically and cytologically specialized in protein synthesis. However, the relationship between such modifications and polyploidization remains unknown. In this study, we investigated the relationship among cell type, queen maturation, and ploidy levels in the fat body of the termite Reticulitermes speratus. We first confirmed that the termite fat body consists of two types of cells, that is, adipocytes, metabolically active cells, and urocytes, urate-storing cells. Our ploidy analysis using flow cytometry has shown that the fat bodies of actively reproducing queens had more polyploid cells than those of newly emerged and pre-reproductive queens, regardless of the queen phenotype (adult or neotenic type). Using image-based analysis, we found that not urocytes, but adipocytes became polyploid during queen differentiation and subsequent sexual maturation. These results suggest that polyploidization in the termite queen fat body is associated with sexual maturation and is regulated in a cell type-specific manner. Our study findings have provided novel insights into the development of insect fat bodies and provide a basis for future studies to understand the functional importance of polyploidy in the fat bodies of termite queens.

Epigenetic regulation and division of labor in social insects

Analogous to multicellular organisms, social insect colonies are characterized by division of labor with queens and workers reflecting germline and soma, respectively. In multicellular organisms, such division is achieved through epigenetic factors regulating cell differentiation during development. Analogously, epigenetic regulation is postulated to regulate caste differences in social insects. We summarize recent findings about the role of epigenetics in social insects, focusing on DNA methylation and histone modifications. We specifically address (i) queen versus worker caste differentiation, (ii) queen versus worker caste differences, and (iii) division of labor among workers. Our review provides an overview of an exciting and controversially discussed field in developmental and molecular biology. It shows that our current understanding about the role of epigenetics in regulating division of labor in social insects is still fragmentary but that refined methods with well-replicated samples and targeted questions offer promising insights into this emerging field of socio-epigenomics.

Evolution of a Cockroach Allergen into the Major Protein of Termite Royal Jelly

Termites live in colonies, and their members belong to different castes that each have their specific role within the termite society. In well-established colonies of higher termites, the only food the founding female, the queen, receives is saliva from workers; such queens can live for many years and produce up to 10,000 eggs per day. In higher termites, worker saliva must thus constitute a complete diet and therein resembles royal jelly produced by the hypopharyngeal glands of honeybee workers that serves as food for their queens; indeed, it might as well be called termite royal jelly. However, whereas the composition of honeybee royal jelly is well established, that of worker termite saliva in higher termites remains largely unknown. In lower termites, cellulose-digesting enzymes constitute the major proteins in worker saliva, but these enzymes are absent in higher termites. Others identified a partial protein sequence of the major saliva protein of a higher termite and identified it as a homolog of a cockroach allergen. Publicly available genome and transcriptome sequences from termites make it possible to study this protein in more detail. The gene coding the termite ortholog was duplicated, and the new paralog was preferentially expressed in the salivary gland. The amino acid sequence of the original allergen lacks the essential amino acids methionine, cysteine and tryptophan, but the salivary paralog incorporated these amino acids, thus allowing it to become more nutritionally balanced. The gene is found in both lower and higher termites, but it is in the latter that the salivary paralog gene got reamplified, facilitating an even higher expression of the allergen. This protein is not expressed in soldiers, and, like the major royal jelly proteins in honeybees, it is expressed in young but not old workers.

Caste-specific expressions and diverse roles of takeout genes in the termite Reticulitermes speratus

Acquisition of novel functions caused by gene duplication may be important for termite social evolution. To clarify this possibility, additional evidence is needed. An important example is takeout, encoding juvenile hormone binding protein. We identified 25 takeouts in the termite Reticulitermes speratus genome. RNA-seq revealed that many genes were highly expressed in specific castes. Two novel paralogs (RsTO1, RsTO2) were tandemly aligned in the same scaffold. Real-time qPCR indicated that RsTO1 and RsTO2 were highly expressed in queens and soldiers, respectively. Moreover, the highest RsTO1 expression was observed in alates during queen formation. These patterns were different from vitellogenins, encoding egg-yolk precursors, which were highly expressed in queens than alates. In situ hybridization showed that RsTO1 mRNA was localized in the alate-frontal gland, indicating that RsTO1 binds with secretions probably used for the defence during swarming flight. In contrast, increased RsTO2 expression was observed approximately 1 week after soldier differentiation. Expression patterns of geranylgeranyl diphosphate synthase, whose product functions in the terpenoid synthesis, were similar to RsTO2 expression. In situ hybridization indicated RsTO2-specific mRNA signals in the soldier-frontal gland. RsTO2 may interact with terpenoids, with a soldier-specific defensive function. It may provide additional evidence for functionalization after gene duplication in termites.

Efficient RNA interference method during caste differentiation with hormone treatment in the termite Reticulitermes speratus (Isoptera: Rhinotermitidae)

Unveiling the proximate mechanism of caste differentiation is crucial for understanding insect social evolution, and gene function analysis is an important tool in this endeavor. The RNA interference (RNAi) technique is useful in termites, but its knockdown effects may differ among species. One of the most important model species in the field of termite sociogenomics is Reticulitermes speratus Kolbe (Isoptera: Rhinotermitidae). Presoldier and worker differentiation of this species can be artificially induced by juvenile hormone and 20-hydroxyecdysone application, respectively. However, appropriate RNAi technique of genes expressed during caste differentiation has never been considered. To clarify this issue, first, we injected nine different volumes of nuclease-free water (NFW, 0-404.8 nL) into workers and found that survival and caste differentiation rates were strongly reduced by the application of the top three largest volumes. Second, we injected double-stranded (ds) RNA of ecdysone receptor homolog (RsEcR) (2.0 µg/151.8 nL NFW) into workers with hormone treatments. The expression levels of RsEcR were significantly reduced at 9 days after dsRNA injection. RsEcR RNAi strongly affected both molting events during presoldier and worker differentiation induced by hormone treatments. The present results highlight the need for caution regarding injection volumes for RNAi experiments using hormone treatments. We suggest that the injection of dsRNA solution (2 µg; approximately 100-200 nL) is suitable for RNAi experiments during caste differentiation induced by hormone application in R. speratus.

Putting hornets on the genomic map

Hornets are the largest of the social wasps, and are important regulators of insect populations in their native ranges. Hornets are also very successful as invasive species, with often devastating economic, ecological and societal effects. Understanding why these wasps are such successful invaders is critical to managing future introductions and minimising impact on native biodiversity. Critical to the management toolkit is a comprehensive genomic resource for these insects. Here we provide the annotated genomes for two hornets, Vespa crabro and Vespa velutina. We compare their genomes with those of other social Hymenoptera, including the northern giant hornet Vespa mandarinia. The three hornet genomes show evidence of selection pressure on genes associated with reproduction, which might facilitate the transition into invasive ranges. Vespa crabro has experienced positive selection on the highest number of genes, including those putatively associated with molecular binding and olfactory systems. Caste-specific brain transcriptomic analysis also revealed 133 differentially expressed genes, some of which are associated with olfactory functions. This report provides a spring-board for advancing our understanding of the evolution and ecology of hornets, and opens up opportunities for using molecular methods in the future management of both native and invasive populations of these over-looked insects.

Gene expression profiles of chemosensory genes of termite soldier and worker antennae

Termites have an elaborate social system that involves cooperation and division of labour among colony members. Although this social system is regulated by chemical signals produced in the colony, it remains unclear how these signals are perceived by other members. Signal transduction is well known to be triggered by the reception of odorant molecules by some binding proteins in the antennae, after which, a signal is transmitted to chemosensory receptors. However, there is insufficient information on the role of chemosensory genes involved in signal transduction in termites. Here, we identified the genes involved in chemosensory reception in the termite Reticulitermes speratus and performed a genome-wide comparative transcriptome analysis of worker and soldier antennae. First, we identified 31 odorant-binding proteins (OBPs), and three chemosensory protein A (CheA) from the genome data. Thereafter, we performed RNA sequencing to compare the expression levels of OBPs, CheAs, and previously identified chemosensory receptor genes between worker and soldier antennae. There were no receptor genes with significant differences in expression between castes. However, the expression levels of three non-receptor odorant-detection/binding proteins (OBP, CheA, and Sensory neuron membrane protein) were significantly different between castes. Real-time qPCR (RT-qPCR) analysis using antennae and other head parts confirmed that these genes were highly expressed in soldier antennae. Finally, independent RT-qPCR analysis showed that the expression patterns of these genes were altered in soldiers from different social contexts. Present results suggest that gene expression levels of some non-receptors are affected by both castes and behavioural interactions among colony members in termites.

Wood feeding and social living: Draft genome of the subterranean termite Reticulitermes lucifugus (Blattodea; Termitoidae)

Termites (Insecta, Blattodea, Termitoidae) are a widespread and diverse group of eusocial insects known for their ability to digest wood matter. Herein, we report the draft genome of the subterranean termite Reticulitermes lucifugus, an economically important species and among the most studied taxa with respect to eusocial organization and mating system. The final assembly (~813 Mb) covered up to 88% of the estimated genome size and, in agreement with the Asexual Queen Succession Mating System, it was found completely homozygous. We predicted 16,349 highly supported gene models and 42% of repetitive DNA content. Transposable elements of R. lucifugus show similar evolutionary dynamics compared to that of other termites, with two main peaks of activity localized at 25% and 8% of Kimura divergence driven by DNA, LINE and SINE elements. Gene family turnover analyses identified multiple instances of gene duplication associated with R. lucifugus diversification, with significant lineage-specific gene family expansions related to development, perception and nutrient metabolism pathways. Finally, we analysed P450 and odourant receptor gene repertoires in detail, highlighting the large diversity and dynamical evolutionary history of these proteins in the R. lucifugus genome. This newly assembled genome will provide a valuable resource for further understanding the molecular basis of termites biology as well as for pest control.

Genome-wide identification and comparative analysis of lipocalin families in Lepidoptera with an emphasis on Bombyx mori

Lipocalins exhibit functional diversity, including roles in retinol transport, invertebrate cryptic coloration, and stress response. However, genome-wide identification and characterization of lipocalin in the insect lineage have not been thoroughly explored. Here, we found that a lineage-specific expansion of the lipocalin genes in Lepidoptera occurred in large part due to tandem duplication events and several lipocalin genes involving insect coloration were expanded more via tandem duplication in butterflies. A comparative analysis of conserved motifs showed both conservation and divergence of lepidopteran lipocalin family protein structures during evolution. We observe dynamic changes in tissue expression preference of paralogs in Bombyx mori, suggesting differential contribution of paralogs to specific organ functions during evolution. Subcellular localization experiments revealed that lipocalins localize to the cytoplasm, nuclear membrane, or nucleus in BmN cells. Moreover, several lipocalin genes exhibited divergent responses to abiotic and biotic stresses, and 1 lipocalin gene was upregulated by 300 fold in B. mori. These results suggest that lipocalins act as signaling components in defense responses by mediating crosstalk between abiotic and biotic stress responses. This study deepens our understanding of the comprehensive characteristics of lipocalins in insects.

King- and queen-specific degradation of uric acid contributes to reproduction in termites

Caste-based reproductive division of labour in social insects is built on asymmetries in resource allocation within colonies. Kings and queens dominantly consume limited resources for reproduction, while non-reproductive castes such as workers and soldiers help reproductive castes. Studying the regulation of such asymmetries in resource allocation is crucial for understanding the maintenance of sociality in insects, although the molecular background is poorly understood. We focused on uric acid, which is reserved and used as a valuable nitrogen source in wood-eating termites. We found that king- and queen-specific degradation of uric acid contributes to reproduction in the subterranean termite Reticulitermes speratus. The urate oxidase gene (RsUAOX), which catalyses the first step of nitrogen recycling from stored uric acid, was highly expressed in mature kings and queens, and upregulated with differentiation into neotenic kings/queens. Suppression of uric acid degradation decreased the number of eggs laid per queen. Uric acid was shown to be provided by workers to reproductive castes. Our results suggest that the capacity to use nitrogen, which is essential for the protein synthesis required for reproduction, maintains colony cohesion expressed as the reproductive monopoly held by kings and queens.

Differential expression of some termite neuropeptides and insulin/IGF-related hormones and their plausible functions in growth, reproduction and caste determination

Background

Insulin-like growth factor (IGF) and other insulin-like peptides (ilps) are important hormones regulating growth and development in animals. Whereas most animals have a single female and male adult phenotype, in some insect species the same genome may lead to different final forms. Perhaps the best known example is the honeybee where females can either develop into queens or workers. More extreme forms of such polyphenism occur in termites, where queens, kings, workers and soldiers coexist. Both juvenile hormone and insulin-like peptides are known to regulate growth and reproduction as well as polyphenism. In termites the role of juvenile hormone in reproduction and the induction of the soldier caste is well known, but the role of IGF and other ilps in these processes remains largely unknown. Here the various termite ilps are identified and hypotheses regarding their functions suggested.

Methods

Genome assemblies and transcriptome short read archives (SRAs) were used to identify insulin-like peptides and neuropeptides in termites and to determine their expression in different species, tissues and castes.

Results and Discussion

Termites have seven different ilps, i.e. gonadulin, IGF and an ortholog of Drosophila insulin-like peptide 7 (dilp7), which are commonly present in insects, and four smaller peptides, that have collectively been called short IGF-related peptides (sirps) and individually atirpin, birpin, cirpin and brovirpin. Gonadulin is lost from the higher termites which have however amplified the brovirpin gene, of which they often have two or three paralogs. Based on differential expression of these genes it seems likely that IGF is a growth hormone and atirpin an autocrine tissue factor that is released when a tissue faces metabolic stress. Birpin seems to be responsible for growth and in the absence of juvenile hormone this may lead to reproductive adults or, when juvenile hormone is present, to soldiers. Brovirpin is expressed both by the brain and the ovary and likely stimulates vitellogenesis, while the function of cirpin is less clear.

Evolution and functionalization of vitellogenin genes in the termite Reticulitermes speratus

Eusociality has been commonly observed in distinct animal lineages. The reproductive division of labor is a particular feature, achieved by the coordination between fertile and sterile castes within the same nest. The sociogenomic approach in social hymenopteran insects indicates that vitellogenin (Vg) has undergone neo-functionalization in sterile castes. Here, to know whether Vgs have distinct roles in nonreproductive castes in termites, we investigated the unique characteristics of Vgs in the rhinotermitid termite Reticulitermes speratus. The four Vgs were identified from R. speratus (RsVg1-4), and RsVg3 sequences were newly identified using the RACE method. Molecular phylogenetic analysis supported the monophyly of the four termite Vgs. Moreover, the termites Vg1-3 and Vg4 were positioned in two different clades. The  dN/dS ratios indicated that the branch leading to the common ancestor of termite Vg4 was under weak purifying selection. Expression analyses among castes (reproductives, workers, and soldiers) and females (nymphs, winged alates, and queens) showed that RsVg1-3 was highly expressed in fertile queens. In contrast, RsVg4 was highly expressed in workers and female nonreproductives (nymphs and winged adults). Localization of RsVg4 messenger RNA was confirmed in the fat body of worker heads and abdomens. These results suggest that Vg genes are functionalized after gene duplication during termite eusocial transition and that Vg4 is involved in nonreproductive roles in termites.

Eusocial Transition in Blattodea: Transposable Elements and Shifts of Gene Expression

(1) Unravelling the molecular basis underlying major evolutionary transitions can shed light on how complex phenotypes arise. The evolution of eusociality, a major evolutionary transition, has been demonstrated to be accompanied by enhanced gene regulation. Numerous pieces of evidence suggest the major impact of transposon insertion on gene regulation and its role in adaptive evolution. Transposons have been shown to be play a role in gene duplication involved in the eusocial transition in termites. However, evidence of the molecular basis underlying the eusocial transition in Blattodea remains scarce. Could transposons have facilitated the eusocial transition in termites through shifts of gene expression? (2) Using available cockroach and termite genomes and transcriptomes, we investigated if transposons insert more frequently in genes with differential expression in queens and workers and if those genes could be linked to specific functions essential for eusocial transition. (3) The insertion rate of transposons differs among differentially expressed genes and displays opposite trends between termites and cockroaches. The functions of termite transposon-rich queen- and worker-biased genes are related to reproduction and ageing and behaviour and gene expression, respectively. (4) Our study provides further evidence on the role of transposons in the evolution of eusociality, potentially through shifts in gene expression.

Comparison of gene expression profiles among caste differentiations in the termite Reticulitermes speratus

Termite castes express specialized phenotypes for their own tasks and are a good example of insect polyphenism. To understand the comprehensive gene expression profiles during caste differentiation, RNA-seq analysis based on the genome data was performed during the worker, presoldier, and nymphoid molts in Reticulitermes speratus. In this species, artificial induction methods for each molt have already been established, and the time scale has been clarified. Three different periods (before the gut purge (GP), during the GP, and after the molt) were discriminated in each molt, and two body parts (head and other body regions) were separately sampled. The results revealed that many differentially expressed genes (head: 2884, body: 2579) were identified in each molt. Based on the independent real-time quantitative PCR analysis, we confirmed the different expression patterns of seven out of eight genes in the presoldier molt. Based on the GO and KEGG enrichment analyses, the expressions of genes related to juvenile hormone titer changes (e.g., JH acid methyltransferase), nutrition status (e.g., Acyl-CoA Delta desaturase), and cell proliferation (e.g., insulin receptor), were shown to specifically fluctuate in each molt. These differences may have a crucial impact on caste differentiation. These data are important resources for future termite sociogenomics.

The Antennal Sensilla and Expression Patterns of Olfactory Genes in the Lower Termite Reticulitermes aculabialis (Isoptera: Rhinotermitidae)

The insect olfactory system plays pivotal roles in insect survival and reproduction through odor detection. Morphological and physiological adaptations are caste-specific and evolved independently in workers, soldiers, and reproductives in termites. However, it is unclear whether the olfactory system is involved in the division of labor in termite colonies. In the present study, the antennal sensilla of alates, workers, soldiers, nymphs, and larvae of the termite Reticulitermes aculabialis Tsai et Hwang ( Isoptera: Rhinotermitidae) were investigated. Transcriptomes were used to detect olfactory genes, and differential expression levels of olfactory genes were confirmed in various castes by qRT-PCR analysis. Nine types of sensilla were identified on the antennae of R. aculabialis, and soldiers possessed all 9 types. In 89,475 assembled unigenes, we found 16 olfactory genes, including 6 chemosensory protein (CSP) and 10 odorant-binding protein (OBP) genes. These OBP genes included 8 general odorant-binding protein genes (GOBPs) and 2 pheromone-binding protein-related protein (PBP) genes. Five CSP genes were more highly expressed in alates than in workers, soldiers, larvae, and nymphs, and the expression levels of CSP6 were significantly higher in nymphs. Seven GOBP and two PBP genes exhibited significantly higher expression levels in alates, and there were no significant differences in the expression levels of GOBP2 among workers, soldiers, alates, and larvae. These results suggest that alates, as primary reproductives, have unique expression patterns of olfactory genes, which play key roles in nuptial flight, mate seeking, and new colony foundation.

Complex regulatory role of DNA methylation in caste- and age-specific expression of a termite

The reproductive castes of eusocial insects are often characterized by extreme lifespans and reproductive output, indicating an absence of the fecundity/longevity trade-off. The role of DNA methylation in the regulation of caste- and age-specific gene expression in eusocial insects is controversial. While some studies find a clear link to caste formation in honeybees and ants, others find no correlation when replication is increased across independent colonies. Although recent studies have identified transcription patterns involved in the maintenance of high reproduction throughout the long lives of queens, the role of DNA methylation in the regulation of these genes is unknown. We carried out a comparative analysis of DNA methylation in the regulation of caste-specific transcription and its importance for the regulation of fertility and longevity in queens of the higher termite Macrotermes natalensis. We found evidence for significant, well-regulated changes in DNA methylation in mature compared to young queens, especially in several genes related to ageing and fecundity in mature queens. We also found a strong link between methylation and caste-specific alternative splicing. This study reveals a complex regulatory role of fat body DNA methylation both in the division of labour in termites, and during the reproductive maturation of queens.

Transcriptomics on Social Interactions in Termites: Effects of Soldier Presence

The organization of social insect colonies requires sophisticated mechanisms to regulate caste composition according to colony demands. In termites, the soldier caste is responsible for the inhibition of soldier differentiation, but the mechanism underlying the regulation of soldier differentiation is still unclear. In this study, we performed transcriptome analyses to identify genes expressed in workers that fluctuated in the presence of soldiers in the subterranean termite Reticulitermes speratus. First, soldier differentiation was artificially induced via juvenile hormone (JH) application, and the inhibitory effects of soldier differentiation on soldier presence were evaluated. Second, transcriptomes were prepared from workers with or without soldiers under JH treatment, and expression analyses were performed to identify differentially expressed genes (DEGs) for each treatment. The expression levels of several DEGs were verified by quantitative real-time PCR. The results indicated that only a small number of DEGs were upregulated by the presence of soldiers. A homology search of DEGs and gene ontology (GO) analysis of the DEGs showed that some genes were responsible for the regulation of hormone levels, social interaction, and response to xenobiotic substances, suggesting that they could be involved in developmental arrest and pheromonal regulation in workers. Moreover, GO analysis indicated that the expression of many genes, including those involved in hormone metabolic processes, fluctuated with JH application. Suppression of soldier differentiation in the presence of soldiers could be accomplished by the expression of a large number of genes required for soldier differentiation.

Termite sociogenomics: evolution and regulation of caste-specific expressed genes

Termite genomes have been sequenced in at least five species from four different families. Genome-based transcriptome analyses have identified large numbers of protein-coding genes with caste-specific expression patterns. These genes include those involved in caste-specific morphologies and roles, for example high fecundity and longevity in reproductives. Some caste-specific expressed genes belong to multi-gene families, and their genetic architecture and expression profiles indicate they have evolved via tandem gene duplication. Candidate regulatory mechanisms of caste-specific expression include epigenetic regulation (e.g. histone modification and non-coding RNA) and diversification of transcription factors and cis-regulatory elements. We review current knowledge in the area of termite sociogenomics, focussing on the evolution and regulation of caste-specific expressed genes, and discuss future research directions.