TGFβ signaling related genes are involved in hormonal mediation during termite soldier differentiation

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.

Regulation of soldier caste differentiation by microRNAs in Formosan subterranean termite (Coptotermes formosanus Shiraki)

The soldier caste is one of the most distinguished castes inside the termite colony. The mechanism of soldier caste differentiation has mainly been studied at the transcriptional level, but the function of microRNAs (miRNAs) in soldier caste differentiation is seldom studied. In this study, the workers of Coptotermes formosanus Shiraki were treated with methoprene, a juvenile hormone analog which can induce workers to transform into soldiers. The miRNomes of the methoprene-treated workers and the controls were sequenced. Then, the differentially expressed miRNAs (DEmiRs) were corrected with the differentially expressed genes DEGs to construct the DEmiR-DEG regulatory network. Afterwards, the DEmiR-regulated DEGs were subjected to GO enrichment and KEGG enrichment analysis. A total of 1,324 miRNAs were identified, among which 116 miRNAs were screened as DEmiRs between the methoprene-treated group and the control group. A total of 4,433 DEmiR-DEG pairs were obtained. No GO term was recognized as significant in the cellular component, molecular function, or biological process categories. The KEGG enrichment analysis of the DEmiR-regulated DEGs showed that the ribosome biogenesis in eukaryotes and circadian rhythm-fly pathways were enriched. This study demonstrates that DEmiRs and DEGs form a complex network regulating soldier caste differentiation in termites.

Screening of reference genes for microRNA analysis in the study of solider caste differentiation of Formosan subterranean termite Coptotermes formosanus Shiraki

The soldier caste differentiation is a complex process that is governed by the transcriptional regulation and post-transcriptional regulation. microRNAs (miRNAs) are noncoding RNAs that control a wide range of activities. However, their roles in solider caste differentiation are barely studied. RT-qPCR is a powerful tool to study the function of genes. A reference gene is required for normalization for the the relative quantification method. However, no reference gene is available for miRNA quantification in the study of solider caste differentiation of Coptotermes formosanus Shiraki. In this research, in order to screen the suitable reference genes for the study of the roles of miRNAs in solider caste differentiation, the expression levels of 8 candidate miRNA genes were quantified in the head and thorax + abdomen during soldier differentiation. The qPCR data were analyzed using geNorm, NormFinder, BestKeeper, ΔCt method and RefFinder. The normalization effect of the reference genes was evaluated using the let-7-3p. Our study showed that novel-m0649-3p was the most stable reference gene, while U6 was the least stable reference gene. Our study has selected the most stable reference gene, and has paved the way for functional analysis of miRNAs in solider caste differentiation.

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.

Social insect transcriptomics and the molecular basis of caste diversity

Studies of gene expression provide fundamentally important information on the molecular mechanisms underlying variation in phenotype. Recent technological advances have allowed for the robust study of gene expression through analysis of whole transcriptomes. Here, we review current advances in social insect transcriptomics and discuss their implications in understanding phenotypic diversity. Recent transcriptomic studies provide detailed inventories of the genes involved in producing distinct phenotypes in social species. These investigations have identified key genes and networks involved in producing distinct social insect castes. Nevertheless, questions concerning the evolution of gene expression patterns remain. We suggest a path forward for studying gene expression in future studies of biological systems.

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.

The TGF-β Receptor Gene Saxophone Influences Larval-Pupal-Adult Development in Tribolium castaneum

The transforming growth factor-β (TGF-β) superfamily encodes a large group of proteins, including TGF-β isoforms, bone morphogenetic proteins and activins that act through conserved cell-surface receptors and signaling co-receptors. TGF-β signaling in insects controls physiological events, including growth, development, diapause, caste determination and metamorphosis. In this study, we used the red flour beetle, Tribolium castaneum, as a model species to investigate the role of the type I TGF-β receptor, saxophone (Sax), in mediating development. Developmental and tissue-specific expression profiles indicated Sax is constitutively expressed during development with lower expression in 19- and 20-day (6th instar) larvae. RNAi knockdown of Sax in 19-day larvae prolonged developmental duration from larvae to pupae and significantly decreased pupation and adult eclosion in a dose-dependent manner. At 50 ng dsSax/larva, Sax knockdown led to an 84.4% pupation rate and 46.3% adult emergence rate. At 100 ng and 200 ng dsSax/larva, pupation was down to 75.6% and 50%, respectively, with 0% adult emergence following treatments with both doses. These phenotypes were similar to those following knockdowns of 20-hydroxyecdysone (20E) receptor genes, ecdysone receptor (EcR) or ultraspiracle protein (USP). Expression of 20E biosynthesis genes disembodied and spookier, 20E receptor genes EcR and USP, and 20E downstream genes BrC and E75, were suppressed after the Sax knockdown. Topical application of 20E on larvae treated with dsSax partially rescued the dsSax-driven defects. We can infer that the TGF-β receptor gene Sax influences larval-pupal-adult development via 20E signaling in T. castaneum.

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.

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.

Soldier Caste-Specific Protein 1 Is Involved in Soldier Differentiation in Termite Reticulitermes aculabialis

Termite soldiers are a unique caste among social insects, and their differentiation can be induced by Juvenile hormone (JH) from workers through two molts (worker–presoldier–soldier). However, the molecular mechanism underlying the worker-to-soldier transformation in termites is poorly understood. To explore the mechanism of soldier differentiation induced by JH, the gene soldier caste-specific protein 1 (RaSsp1, NCBI accession no: MT861054.1) in R. aculabialis was cloned, and its function was studied. This gene was highly expressed in the soldier caste, and the protein RsSsp1 was similar to the JHBP (JH-binding protein) domain-containing protein by Predict Protein online. In addition, JHIII could be anchored in the hydrophobic cage of RaSsp1 as the epoxide of the JHBP-bound JH according to the protein ligand molecular docking online tool AutoDock. The functional studies indicated that knocking down of the RaSsp1 shorted the presoldier’s head capsule, reduced mandible size, delayed molting time and decreased molting rate (from worker to presoldier) at the beginning of worker gut-purging. Furthermore, knocking down of the RaSsp1 had a more pronounced effect on soldier differentiation (from presoldier to soldier), and manifested in significantly shorter mandibles, rounder head capsules, and lower molting rate (from worker to presoldier) at the beginning of presoldier gut-purging. Correspondingly, the expressions of JH receptor Methoprene-tolerant (Met), the JH-inducible transcription factor Krüppel homolog1 (Kr-h1) and ecdysone signal genes Broad-complex (Br-C) were downregulated when knocking down the RaSsp1 at the above two stages. All these results that RaSsp1 may be involved in soldier differentiation from workers by binding and transporting JH.

Smad on X is vital for larval-pupal transition in a herbivorous ladybird beetle

Insect development is regulated by a combination of juvenile hormone (JH) and 20-hydroxyecdysone (20E). Production of both JH and 20E is regulated by transforming growth factor-β (TGFβ) signaling. TGFβ can be classified into two branches, the Activin and Bone Morphogenetic Protein (BMP) pathways. In Drosophila melanogaster, BMP signaling is critical for JH synthesis, whereas Activin signal is required to generate the large pulse of 20E necessary for entering metamorphosis. However, to which extent the roles of these signals are conserved remains unknown. Here we studied the role of an Activin component Smad on X (Smox) in post-embryonic development in a defoliating ladybird Henosepilachna vigintioctopunctata. RNA interference (RNAi)-aided knockdown of Hvsmox inhibited larval growth, and impaired larval development. All Hvmyo RNAi larvae arrested at the fourth-instar larval stage. Moreover, knockdown of Hvsmox delayed gut and Malpighian tubules remodeling. Furthermore, the expression of a JH biosynthesis gene (Hvjhamt), a JH receptor gene HvMet and a JH response gene HvKr-h1 was greatly enhanced. Conversely, the expression levels of an ecdysteroidogenesis gene (Hvspo), a 20E receptor gene (HvEcR) and six 20E response genes (HvBrC, HvE74, HvE75, HvE93, HvHR3 and HvHR4) were significantly lowered. Knockdown of HvMet partially restored the negative phenotypes in the Hvsmox RNAi beetles. Our results suggest that Smox exerts regulative roles in JH production, ecdysteroidogenesis and organ remodeling, thus contributing to modulate the larva-pupa-adult transformation in H. vigintioctopunctata.

Transcriptomic analysis of the honey bee (Apis mellifera) queen spermathecae reveals genes that may be involved in sperm storage after mating

Honey bee (Apis mellifera) queens have a remarkable organ, the spermatheca, which successfully stores sperm for years after a virgin queen mates. This study uniquely characterized and quantified the transcriptomes of the spermathecae from mated and virgin honey bee queens via RNA sequencing to identify differences in mRNA levels based on a queen's mating status. The transcriptome of drone semen was analyzed for comparison. Samples from three individual bees were independently analyzed for mated queen spermathecae and virgin queen spermathecae, and three pools of semen from ten drones each were collected from three separate colonies. In total, the expression of 11,233 genes was identified in mated queen spermathecae, 10,521 in virgin queen spermathecae, and 10,407 in drone semen. Using a cutoff log2 fold-change value of 2.0, we identified 212 differentially expressed genes between mated and virgin spermathecal queen tissues: 129 (1.4% of total) were up-regulated and 83 (0.9% of total) were down-regulated in mated queen spermathecae. Three genes in mated queen spermathecae, three genes in virgin queen spermathecae and four genes in drone semen that were more highly expressed in those tissues from the RNA sequencing data were further validated by real time quantitative PCR. Among others, expression of Kielin/chordin-like and Trehalase mRNAs was highest in the spermathecae of mated queens compared to virgin queen spermathecae and drone semen. Expression of the mRNA encoding Alpha glucosidase 2 was higher in the spermathecae of virgin queens. Finally, expression of Facilitated trehalose transporter 1 mRNA was greatest in drone semen. This is the first characterization of gene expression in the spermathecae of honey bee queens revealing the alterations in mRNA levels within them after mating. Future studies will extend to other reproductive tissues with the purpose of relating levels of specific mRNAs to the functional competence of honey bee queens and the colonies they head.

Diversity of Insect Sesquiterpenoid Regulation

Insects are arguably the most successful group of animals in the world in terms of both species numbers and diverse habitats. The sesquiterpenoids juvenile hormone, methyl farnesoate, and farnesoic acid are well known to regulate metamorphosis, reproduction, sexual dimorphism, eusociality, and defense in insects. Nevertheless, different insects have evolved with different sesquiterpenoid biosynthetic pathway as well as products. On the other hand, non-coding RNAs such as microRNAs have been implicated in regulation of many important biological processes, and have recently been explored in the regulation of sesquiterpenoid production. In this review, we summarize the latest findings on the diversity of sesquiterpenoids reported in different groups of insects, as well as the recent advancements in the understanding of regulation of sesquiterpenoid production by microRNAs.

The making of the defensive caste: Physiology, development, and evolution of the soldier differentiation in termites

Termites (Blattodea, Termitoidea, or Isoptera) constitute one of the major lineages of eusocial insects. In termite societies, multiple types of functional individuals, that is, castes, perform divisions of labors to coordinate social behaviors. Among other castes, the soldier caste is distinctive since it is sterile and exclusively specialized into defensive behavior with largely modified morphological features. Therefore, many of the previous studies have been focused on soldiers, in terms of ecology, behavior, and evolution as well as developmental and physiological mechanisms. This article overviews the accumulation of studies especially focusing on the developmental and physiological mechanisms underlying the soldier differentiation in termites. Furthermore, the evolutionary trajectories that have led the acquisition of soldier caste and have diversified the soldier characteristics in association with the social evolution are discussed.

Methoprene-Induced Genes in Workers of Formosan Subterranean Termites (Coptotermes formosanus Shiraki)

Termites have a distinct polyphenism controlled by concise hormonal and molecular mechanisms. Workers undergo double molts to transform into soldiers (worker–presoldier–soldier). Juvenile hormone analogs, such as methoprene, can induce workers to transform into presoldiers. However, the molecular mechanism underlying the worker-to-presoldier transformation in Coptotermes formosanus Shiraki is still not clear. We sequenced the transcriptome of workers four days after they had fed on methoprene-treated filter paper and control group workers, which fed on acetone-treated filter paper. The transcriptome of C. formosanus was assembled using the de novo assembly method. Expression levels of unigenes in the methoprene-treated group and the control group were compared. The differentially expressed genes were further analyzed by Gene Ontology (GO) term enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis. Tetrapyrrole binding, oxidoreductase activity, and metal ion binding were the only three enriched GO terms. Juvenile hormone synthesis was the first ranked enriched pathway. Carbohydrate, amino acid, and lipid metabolism pathways were also enriched. These three pathways may be related to fat body development, which is critical for presoldier formation. Our results have demonstrated the significance of JH synthesis pathways, and pathways related to fat body development in the artificial induction of presoldiers.

Transcriptomic and Functional Analyses of Phenotypic Plasticity in a Higher Termite, Macrotermes barneyi Light

Eusocial termites have a complex caste system, which leads to the division of labor. Previous studies offered some insight into the caste differentiation in lower termites; however, few studies were focusing on the molecular mechanisms of higher termites with sophisticated societies. Comparative transcriptomic analyses of five immature castes of a higher termite, Macrotermes barneyi Light, suggest that phenotypic plasticity is modulated by an array of transcriptional changes, including differentially expressed genes (e.g., caste-biased genes Vtg and TnC), co-expression networks (e.g., genes associated with nymph reproduction), and alternative splicing (e.g., events related to muscle development in presoldiers). Transcriptional (RT-PCR and RT-qPCR) and functional (in vivo RNAi) validation studies reveal multiple molecular mechanisms contributing to the phenotypic plasticity in eusocial termites. Molecular mechanisms governing the phenotypic plasticity in M. barneyi could be a rule rather than an exception in the evolution of sociality.

Histone modifying genes are involved in the molting period during soldier differentiation in Zootermopsis nevadensis

Caste differentiation in eusocial insects is an outstanding example of phenotypic plasticity. Recent studies indicate that epigenetic regulation, including DNA methylation and histone modification, play a role in the morphological and behavioral polyphenism observed in the caste differentiation of hymenopteran insects. The role of epigenetic regulation in termite caste differentiation, however, is still obscure. In this study, we performed a functional analysis of epigenetic-related genes during soldier differentiation in Zootermopsis nevadensis, for which the entire genome sequence is available. In an incipient colony of this species, the oldest 3rd instar larva (No. 1 larva) always differentiates into a presoldier (intermediate stage of soldier), and the next-oldest 3rd instar larva (No. 2 larva) molts into a 4th instar (which functions as a worker). First, we detected seven epigenetic-related genes with significantly increased expression levels in the soldier-destined No. 1 larvae using RNA-seq data. Second, RNA interference (RNAi) of these seven genes was performed in the No. 1 larvae. RNAi of three histone modifying genes extended the presoldier molting period. Furthermore, these RNAi treatments reduced the expression levels of genes involved in juvenile hormone (JH) synthesis, binding and signaling. These results indicate that epigenetic-related genes do not directly affect termite soldier differentiation; nonetheless, some histone modifying genes have an effect on molting periods, possibly due to the regulation of JH action during soldier differentiation.

Transcriptomic changes during caste development through social interactions in the termite Zootermopsis nevadensis

One of the most striking examples of phenotypic plasticity is the different phenotypes (i.e., castes) within a same nest of social insects. Castes usually derive from a single genotype initially by receiving social cues among individuals during development. Specific gene expression changes may be involved in caste differentiation, and thus, the regulatory mechanism of these changes should be clarified in order to understand social maintenance and evolution. The damp-wood termite Zootermopsis nevadensis is one of the most important model termite species, due to not only the availability of genomic, transcriptomic, and epigenomic information but also evidence that soldier- and worker-destined individuals can be identified in natural conditions. Given that the nutritional intakes via social interactions are crucial for caste differentiation in this species, there is a possibility that transcriptomic changes are influenced by the nutritional difference among these individuals. Here, whole body RNA-seq analysis of 3rd-instar larvae with biological replications and Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses were conducted. We found the drastic expression differences during caste developments between soldier- and worker-destined individuals. The results indicated that there are several key signaling pathways responsible for caste formations, which are involved in developments and social interactions. Particularly, the nutritional sensitive signaling was upregulated in soldier-destined individuals, while some metabolic pathways were identified in worker-destined individuals. These bioinformatic data obtained should be utilized to examine the molecular mechanisms of caste determination in social insects.

A Crucial Caste Regulation Gene Detected by Comparing Termites and Sister Group Cockroaches

Sterile castes are a defining criterion of eusociality; investigating their evolutionary origins can critically advance theory. In termites, the soldier caste is regarded as the first acquired permanently sterile caste. Previous studies showed that juvenile hormone (JH) is the primary factor inducing soldier differentiation, and treatment of workers with artificial JH can generate presoldier differentiation. It follows that a shift from a typical hemimetabolous JH response might be required for soldier formation during the course of termite evolution within the cockroach clade. To address this possibility, analysis of the role of JH and its signaling pathway was performed in the termite Zootermopsis nevadensis and compared with the wood roach Cryptocercus punctulatus, a member of the sister group of termites. Treatment with a JH analog (JHA) induced a nymphal molt in C. punctulatus RNA interference (RNAi) of JH receptor Methoprene tolerant (Met) was then performed, and it inhibited the presoldier molt in Z. nevadensis and the nymphal molt in C. punctulatus Knockdown of Met in both species inhibited expression of 20-hydroxyecdysone (20E; the active form of ecdysone) synthesis genes. However, in Z. nevadensis, several 20E signaling genes were specifically inhibited by Met RNAi. Consequently, RNAi of these genes were performed in JHA-treated termite individuals. Knockdown of 20E signaling and nuclear receptor gene, Hormone receptor 39 (HR39/FTZ-F1β) resulted in newly molted individuals with normal worker phenotypes. This is the first report of the JH-Met signaling feature in termites and Cryptocercus JH-dependent molting activation is shared by both taxa and mediation between JH receptor and 20E signalings for soldier morphogenesis is specific to termites.