- Open Access
Reduction of a nymphal instar in a dampwood termite: heterochronic shift in the caste differentiation pathways
© The Author(s) 2019
- Received: 2 February 2019
- Accepted: 16 April 2019
- Published: 16 May 2019
Generally in termites, alates differentiate through multiple nymphal instars which gradually develop wing buds. However, in a dampwood termite, Hodotermopsis sjostedti, alates molt directly from a single nymphal instar with short wing buds. In this study, to examine the mechanism underlying the wing formation during the alate differentiation in H. sjostedti, histological and morphological observations were carried out on the developmental process of wing formation during the nymphal instar, in comparison with those in Zootermopsis nevadensis, which has two nymphal instars. Furthermore, the expression patterns of genes that are thought to be responsible for wing formation, i.e., wing-patterning genes and genes encoding hormone-related factors, were quantified during alate differentiation and compared between the two species.
The results showed that, in H. sjostedti, wings were formed in a complicatedly folded shape, not only inside the wing buds as seen in Z. nevadensis, but also under the dorsal thoracic cuticle, where the wing tips shifted toward the median thoracic part. Accordingly, the wing expansion pattern also differed from that in Z. nevadensis. Furthermore, the results of real-time qRT-PCR on overall expression profiles of wing-patterning genes and hormone-related genes suggest that the single nymphal instar in H. sjostedti well resembles to the second nymphal instar in Z. nevadensis. In particular, significant upregulation of vestigial (vg) and downregulation of Krüppel homolog 1 (Kr-h1) that were observed at the second nymphal instar in Z. nevadensis apparently occurred during the single nymphal instar in H. sjostedti.
The developmental events for wing formation are compacted into a single nymphal instar in H. sjostedti, and as a result, the unique wing formation is seen to compensate for the spatial restriction inside small wing buds, leading to the completion of functional wings.
- Alate differentiation
- Caste differentiation
- Nymphal instar
- Wing formation
- Gene expression
The caste differentiation pathways can be classified into two major patterns depending on the timing of differentiation into reproductive individuals (imagos): linear and forked pathways [3–5]. In the forked pathways, alate development proceeds exclusively via nymphal instars, while worker and soldier development proceed via distinct larval instars. In contrast, in the linear pathways, all individuals go through larval instars without wing buds up to late instars (i.e., pseudergates that perform worker tasks but still possess potential to develop into alates). A fraction of these pseudergates differentiate into alates through nymphal instars with apparent wing buds or wing pads [4, 6, 7]. The forked pathways are seen in Mastotermitidae, Hodotermitidae, and most species of Rhinotermitidae and Termitidae, while the linear pathways are seen in Archotermopsidae, Kalotermitidae, Serritermitidae and a part of Rhinotermitidae [3, 8].
In termites, although the number of moltings that are required for the soldier differentiation is constant among species (2 moltings via presoldier), the number of nymphal instars is diverse among termite lineages [4, 6, 7]. Generally, termites possess multiple nymphal instars through which wing buds are gradually developed, forming complete wings . For example, in Termitidae showing the forked pathways, there are 6 nymphal instars . Species belonging to Archotermopsidae, which are thought to have primitive characteristics, possess 2 or more successive nymphal instars [4, 10].
However, in the Japanese dampwood termite Hodotermopsis sjostedti, which belongs to Archotermopsidae, is distributed in parts of East and Southeast Asia, pseudergates differentiate into alates through only a single nymphal instar with short wing buds (Fig. 1, [11, 12]). Other related species such as Zootermopsis nevadensis (Archotermopsidae) possess 2 nymphal instars, i.e., short and long wing-budded nymphs (Fig. 1, ). In these species, wing buds gradually develop and alate wings form inside the long wing buds (wing pads) at the final nymphal instar. This suggests that the developmental process of alate differentiation in H. sjostedti is distinctive, since wings should be formed inside the small wing buds of the first-instar nymphs.
In this study, therefore, to clarify the mechanism underlying the wing formation during the alate differentiation in H. sjostedti, histological and morphological observations were carried out on the developmental process of wing formation during its nymphal instar, in comparison with those in Z. nevadensis. Furthermore, expression patterns of genes that are thought to be responsible for the wing formation, i.e., wing-patterning genes and genes for hormone-related factors, were quantified during the course of alate development and compared between the two species. The regulatory genes involved in wing formation are known in several insects (e.g., [13–15]). In addition, hormone-related factors involved in molting and metamorphosis have also been reported in many insects (e.g., [16–20]). Altogether, 5 wing-patterning genes and 9 hormone-related genes were identified, and their expression patterns were analyzed during alate differentiation in the two focal species. Results clearly showed that the overall expression patterns at the second nymphal instar in Z. nevadensis apparently shifted to the single nymphal instar in H. sjostedti.
Histological observations during wing formation
In the case of H. sjostedti, at the early stage of nymphs (NE), a single layer of epithelial tissue was observed under the cuticle of short wing buds (Fig. 2m, p). Later, at the mid-stage of nymphs (NM), the epithelial layer started to proliferate, forming a folded structure (Fig. 2n, q). Then, at the late stage of nymphs (NL), the folded structure became thicker and more complicated (Fig. 2o, r). The proliferated epithelial tissues that would form wings could not be contained inside the wing buds. So, as a result, the tissues of folded epithelial layers were seen in the median region of thoracic parts (Fig. 2o).
SEM observations on wing expansion
On the other hand, in H. sjostedti, since folded wings could not be contained inside the short wing buds of nymphs, the folded wings extended into the median region of the thoracic part (Fig. 3d), as seen in the section images (Fig. 2o). After the onset of wing expansion, as seen in Z. nevadensis, wings were expanded from the basal part, so that the folded shape of the apical part of the wing was maintained (Fig. 3e, f). These results showed that, in the case of H. sjostedti, the apical part of the wing did not form at the tip of the wing bud, but rather at a position between the wing-bud tip and the median line of the thorax.
Expression patterns of wing-patterning genes and hormone-related genes
Principal component analysis on gene expression patterns
The differences of wing formation
As shown by the results obtained in this study, in H. sjostedti, the imaginal molt to alates occurs from the first-instar nymphs with short wing buds. Therefore, the developmental sequence of events of its wing formation differs from that of related species that possess multiple nymphal instars, such as Z. nevadensis. Histological and morphological investigations revealed that, in H. sjostedti, the extensive epithelial cell proliferation starts during the short-winged nymphal instars so that the complicatedly folded epithelial tissues are compacted all over the thorax (Fig. 2). As a result, the folding and expansion patterns of wings differ from those in Z. nevadensis (Fig. 3).
Heterochronic shifts in wing-patterning and hormone-related genes
From these patterns of epithelial proliferation and wing expansion, it was hypothesized that the expression of genes for wing formation should occur earlier (in short-wing-budded first-instar nymphs) in H. sjostedti, as the result of deletion of a nymphal stage. As predicted, this study revealed that many of the genes required for wing formation are upregulated in the second nymphal instar of Z. nevadensis, but in the first nymphal instar of H. sjostedti (Figs. 4, 5). Interestingly, the pattern of vg expression was different from other genes (later in development), probably due to the functional differences among genes. The RNA extraction was carried out from whole thoracic segments, so that these expressions can be responsible not only for the wing development, but also for other tissues including flight muscles, which are also required for alate development. In conclusion, it can be said that heterochronic shifts of wing formation have occurred in the lineage of H. sjostedti, so that the developmental events for wing formation that normally happen through two (or more) nymphal instar stages have been compacted within a single nymphal instar.
In endocrine systems in insects, in general, downregulation of the juvenile hormone (JH) pathway is required for the imaginal molt . Therefore, also in termites, a low JH level is known to induce alate differentiation [22–24]. We showed in the present study that the expression level of Kr-h1 was maintained until the early stage of second nymphal instar (N2E) and downregulated at the late stage (N2L) in Z. nevadensis, while it was downregulated at the late stage of nymphs (NL) in H. sjostedti, suggesting that the imaginal molt is apparently induced during the short-wing-budded nymphal instar in H. sjostedti. This situation is similar to the precocious metamorphosis reported in some hemimetabolous insects, in which knockdown of Kr-h1 induced precocious metamorphosis [17, 25, 26]. It is therefore suggested that, during evolution, some adventitious downregulation of JH signaling may have led to a precocious imaginal molt, leading to the lack of second-instar nymphs in H. sjostedti.
Among genes that showed no significance in the expression levels, many genes such as dpp, Scr, BR-C, USP, E75, InR2 and FOXO, showed interesting patterns that can explain the heterochronic shifts; NE and NL in H. sjostedti, respectively, corresponds to N2E and N2L in Z. nevadensis. As for the overall expression patterns, the results of PCA clearly supported the ideas of heterochrony (Fig. 6). All the principal components independent from the species differences supported that N1 and N2E in Z. nevadensis corresponds to NE in H. sjostedti, and N2L in Z. nevadensis corresponds to NL in H. sjostedti (Additional file 4: Table S2).
The developmental constraints of short wing buds
Adaptive significance and intra-colonial conflict
Considering the adaptive significance underlying this unique developmental phenomenon, a candidate factor is thought to be intra-colonial conflict among reproductive candidates. As previously suggested , wing bud mutilation could be related to the intra-colonial conflict over reproduction. Mutilated nymphs often undergo regressive molt (molt of a nymph to a pseudergate) in some termite species. In H. sjostedti, nymphs whose wing buds were mutilated were often found in relatively large and old colonies, and the regressive molt from nymphs to pseudergates was observed in the laboratory rearing experiments . Such wing bud mutilation could be the selective pressure that led to the reduction of second-instar nymphs in H. sjostedti. In some species of Zootermopsis, reproductive soldiers or soldier neotenics were reported, and they are suggested to have evolved as a result of intra- or inter-colonial conflicts [32, 33]. In addition, it is also suggested that the caste differentiation pathways in this family (Archotermopsidae) are more plastic than those in other lineages. In this termite family, therefore, these preadaptive conditions, i.e., intragroup conflicts and developmental plasticity, might have facilitated the appearance of soldier neotenics or the heterochronic reduction of nymphal instars.
In this study, we obtained information about the unique alate differentiation of H. sjostedti from the perspectives of histology and gene expression. Our findings suggested that the developmental events for wing formation are compacted into a single nymphal instar in H. sjostedti, and as a result, a unique wing formation is seen to compensate for the spatial restriction inside small wing buds, leading to the completion of functional wings. It is suggested that there are behavioral or ecological factors, such as intracolony conflicts among the reproductives, underlying the reduction of nymphal instars. In the future, further studies including behavioral observations will be necessary to get insight into the evolutional process in the reduction of nymphal instars.
Colonies of H. sjostedti Holmgren and Z. nevadensis (Hagen) were collected on Yakushima Island, Kagoshima Prefecture, and at Kawanishi City, Hyogo Prefecture, Japan, respectively, in May 2016–2018. Z. nevadensis is known to be a species that was introduced to Japan . Nest logs of both species were maintained in the laboratory as stock colonies at approximately 25 °C under constant darkness, as described previously [11, 12, 35]. For their food, moistened pine lumber was occasionally supplied. In H. sjostedti, since alates are produced in June under natural conditions, many of the colonies sampled in May contain nymphs. In Z. nevadensis, the season for alate production is not so restricted, and the sampled colonies usually contain nymphs throughout most of the year.
Paraffin sections were made for histological observations on the developmental processes of wings in both species. The examined time points were first-instar nymph (N1) and early and late stage of second-instar nymph (N2E and N2L) in Z. nevadensis, and early, mid and late nymphal stage (NE, NM and NL) in H. sjostedti. Dissected thoracic parts of nymphs were fixed in FAA fixative (formaldehyde/acetic acid/ethanol = 6:1:16) and preserved in 70% ethanol. Samples were dehydrated in increasing concentrations of ethanol, then transferred into xylene and finally embedded into paraffin. Serial section (8 µm thick) in sagittal and transverse planes was processed routinely and stained with hematoxylin and eosin according to the methods previously reported . Tissues on slides were observed using a BZ-9000 HS All-in-one fluorescence microscope (KEYENCE, Osaka, Japan).
Scanning electron microscopy
Scanning electron microscopy was carried out to compare the wing bud structures of nymphs and the process of wing expansion during the imaginal molt between Z. nevadensis and H. sjostedti. This process was classified into 3 time points, i.e., prior to imaginal molt (PIM), early stage of imaginal molt (EIM), and late stage of imaginal molt (LIM). Thoracic parts of individuals before and after the imaginal molt were fixed in FAA. After fixation, the thoracic parts were transferred into increasing concentrations of ethanol followed by t-butanol as described in a previous report . In samples at early stages of the imaginal molt, nymphal cuticles were easily removed for observations on newly formed wings. After the transfer into t-butanol, the fixed samples were freeze-dried using a Freeze Dryer ES-2030 (Hitachi Global, Tokyo, Japan), and coated with gold ions with an Ion Sputter E-1010 (Hitachi Global, Tokyo, Japan). Samples were observed with a JSM-5510LV scanning electron microscope (JEOL Ltd., Tokyo, Japan).
Identification of wing-patterning and hormone-related genes
For the comparison of wing development between the two termite species, orthologs of the wing development genes and hormone-related genes were searched in the genome sequence database of Z. nevadensis (http://termitegenome.org; ) and the transcriptome data of H. sjostedti (DDBJ Sequence Read Archive: DRA005483; Sugime et al., .). As a result, 5 genes (ap [apterous], dpp [decapentaplegic], vg [vestigial], Scr [Sex-combs reduced], exd [extradenticle]) were identified as wing development genes and 9 genes (Met [Methoprene-tolerant], Kr-h1 [Krüppel homolog 1], BR-C [Broad-Complex], EcR [Ecdysone Receptor], USP [ultraspiracle], E93, E75, InR2 [Insulin Receptor 2], FOXO [Forkhead box O]) were identified as hormone-related genes. The protein sequences of those candidate genes in Drosophila melanogaster were used as queries for tBLASTn searches against the datasets of the two species. Based on the obtained sequences, primers for real-time qRT-PCR were designed using Primer Express software (ver. 3.0.0, Applied Biosystems, Foster City, CA, USA, Additional file 6: Table S4).
RNA extraction and real-time qRT-PCR
For the investigations on the expression patterns of genes involved in the wing formation, thoracic parts (prothorax, mesothorax and metathorax excluding guts and legs) were dissected from individuals, and immediately frozen in liquid nitrogen. The examined time points were pseudergate (PE), first-instar nymph (N1), early and late stage of second-instar nymph (N2E and N2L) and alate (A) in Z. nevadensis, and pseudergate (PE), early and late nymphal stage (NE and NL) and alate (A) in H. sjostedti. Frozen samples were individually preserved in 1.5 ml polypropylene tubes at − 80 °C until RNA extraction. Total RNA was extracted using RNAiso Plus (Takara Bio, Shiga, Japan) according to the manufacturer’s protocol. After the extraction, samples were treated with DNaseI (Thermo Fisher Scientific, Waltham, MA, USA).
For each sample, 900 ng of total RNA was reverse-transcribed with a High Capacity cDNA Reverse Transcription Kit according to the manufacturer’s instructions (Applied Biosystems, Foster City, CA, USA). Quantifications of the relative amounts of target transcripts were then performed using a Fast SYBR Green Master Mix and the sequence detection system ABI PRISM 7500 (Applied Biosystems, Foster City, CA, USA). For identifying an appropriate endogenous control gene, the suitability of different candidate reference genes (18S rRNA, ribosomal protein L13a [RPL13a] and ribosomal protein 49 [RP49] in Z. nevadensis; 18S rRNA, elongation factor 1 alpha [EF1a], and ribosomal protein 49 [RP49] in H. sjostedti) was evaluated using the software geNorm  and Normfinder . As a result, RPL13a and RP49 were found to be the most appropriate reference genes for Z. nevadensis and H. sjostedti, respectively. Data acquisition and analyses were handled with ABI Prism 7500 software ver. 2.0.4 (Applied Biosystems, Foster City, CA, USA), with the relative standard curve method. For statistical analyses to detect significant differences of expression levels, respectively, in each species, Tukey’s multiple comparisons test (P < 0.05) was performed after one-way ANOVA (P < 0.05), using R 3.5.2 (https://www.r-project.org). The number of samples was n = 5 for all the samples of Z. nevadensis, and n = 3 (PE, A) and n = 7 (NE, NL) for H. sjostedti.
Principal component analysis
To examine the correspondence of nymphal stages between the two focal stages based on gene expression patterns, principal component analysis (PCA) was performed. In the analysis, hypothetical models of nymphal stages were made based on the stages of H. sjostedti, and were compared with the stages in Z. nevadensis. Firstly, the gene expression data that were standardized in the respective species to match the expression ranges among genes, were used for PCA on the correlation matrix to unify the respective gene expression patterns as principal component scores. Thereafter, the respective principal component scores were analyzed by the generalized linear model (GLM) with the Gaussian distribution with identity link. According to the Kaiser’s criterion , 4 principal components were selected (Additional file 1: Fig. S1). Species and nymphal stages were used as the explanatory variables. The significance of the respective coefficients was examined by the Wald test. The best model was chosen in accordance with the Akaike’s information criterion (AIC).
RN and TM designed the study outline. RN and KO performed the experiments and analyzed the obtained data. JS performed the statistical analyses. RN, KO, JS, SK and TM wrote the manuscript. All authors read and approved the manuscript.
We would like to express our gratitude to Run Minoura, Hiroyuki Shimoji, Masaru K. Hojo and Yudai Masuoka for their assistance in field sampling.
The authors declare that they have no competing interests.
Availability of data and materials
Indicated in the text.
Consent for publication
Ethics approval and consent to participate
This work was supported by a Grant-in-Aid for Scientific Research A (No. 25251041) to TM and by Grant-in-Aid for JSPS Fellows to KO (No. 17J06879), from the Ministry of Education, Culture, Sports, Science and Technology of Japan.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Wilson EO. The insect societies. Cambridge: Belknap Press of Harvard University Press; 1971. p. 548.Google Scholar
- Miura T. Developmental regulation of caste-specific characters in social-insect polyphenism. Evol Dev. 2005;7:122–9.View ArticlePubMedGoogle Scholar
- Roisin Y. Diversity and evolution of caste patterns. In: Abe T, Bignell DE, Higashi M, editors. Termites: evolution, sociality, symbioses, ecology. Dordrecht: Kluwer; 2000. p. 95–119.View ArticleGoogle Scholar
- Noirot C. Pathways of caste development in the lower termites. In: Watson JAL, Okot-Kotber BM, Noirot C, editors. Caste differentiation in social insects. New York: Pergamon Press; 1985. p. 41–57.View ArticleGoogle Scholar
- Thompson GJ, Kitade O, Lo N, Crozier RH. Phylogenetic evidence for a single, ancestral origin of a ‘true’ worker caste in termites. J Evol Biol. 2000;13:869–81.View ArticleGoogle Scholar
- Korb J, Hartfelder K. Life history and development: a framework for understanding developmental plasticity in lower termites. Biol Rev. 2008;83:295–313.View ArticlePubMedGoogle Scholar
- Roisin Y, Korb J. Social organisation and the status of workers in termites. In: Bignell DE, Roisin Y, Lo N, editors. Biology of termites: a modern synthesis. Dordrecht: Springer; 2010. p. 133–64.View ArticleGoogle Scholar
- Legendre F, Whiting MF, Grandcolas P. Phylogenetic analyses of termite post-embryonic sequences illuminate caste and developmental pathway evolution. Evol Dev. 2013;15:146–57.View ArticlePubMedGoogle Scholar
- Noirot C. Formation of castes in the higher termites. In: Krishna K, Weesner FM, editors. Biology of termites, vol. 1. New York: Academic; 1969. p. 311–50.View ArticleGoogle Scholar
- Castle GB. The damp-wood termites of the western United States, genus Zootermopsis (formerly, Termopsis). In: Kofoid GA, editor. Termite and termite control. Berkeley: University of California Press; 1934. p. 273–310.Google Scholar
- Miura T, Hirono Y, Machida M, Kitade O, Matsumoto T. Caste developmental system of the Japanese damp-wood termite Hodotermopsis japonica (Isoptera: Termopsidae). Ecol Res. 2000;15:83–92.View ArticleGoogle Scholar
- Miura T, Koshikawa S, Machida M, Matsumoto T. Comparative studies on alate wing formation in two related species of rotten-wood termites: Hodotermopsis sjostedti and Zootermopsis nevadensis (Isoptera, Termopsidae). Insectes Soc. 2004;51:247–52.View ArticleGoogle Scholar
- Weatherbee SD, Halder G, Kim J, Hudson A, Carroll S. Ultrabithorax regulates genes at several levels of the wing-patterning hierarchy to shape the development of the Drosophila haltere. Genes Dev. 1998;12:1474–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Abouheif E, Wray GA. Evolution of the gene network underlying wing polyphenism in ants. Science. 2002;2975579:249–52.View ArticleGoogle Scholar
- Brisson JA, Ishikawa A, Miura T. Wing development genes of the pea aphid and differential gene expression between winged and unwinged morphs. Insect Mol Biol. 2010;19:63–73.View ArticlePubMedGoogle Scholar
- King-Jones K, Thummel CS. Nuclear receptors—a perspective from Drosophila. Nat Rev Genet. 2005;6:311.View ArticlePubMedGoogle Scholar
- Lozano J, Belles X. Conserved repressive function of Krüppel homolog 1 on insect metamorphosis in hemimetabolous and holometabolous species. Sci Rep. 2011;1:163.View ArticlePubMedPubMed CentralGoogle Scholar
- Charles JP, Iwema T, Epa VC, Takaki K, Rynes J, Jindra M. Ligand-binding properties of a juvenile hormone receptor; Methoprene-tolerant. Proc Natl Acad Sci. 2011;108:21128–33.View ArticlePubMedGoogle Scholar
- Belles X, Santos CG. The MEKRE93 (Methoprene tolerant-Krüppel homolog 1-E93) pathway in the regulation of insect metamorphosis, and the homology of the pupal stage. Insect Biochem Mol Biol. 2014;52:60–8.View ArticlePubMedGoogle Scholar
- Kayukawa T, Jouraku A, Ito Y, Shinoda T. Molecular mechanism underlying juvenile hormone-mediated repression of precocious larval–adult metamorphosis. Proc Natl Acad Sci. 2017;114:1057–62.View ArticlePubMedGoogle Scholar
- Goodman WG, Granger NA. In: Gilbert LI, Iatrou K, Gill SS, editors. Comprehensive molecular insect science, vol. 3. Oxford: Elsevier; 2005. p. 319–408.View ArticleGoogle Scholar
- Nijhout HF, Wheeler DE. Juvenile hormone and the physiological basis of insect polymorphisms. Q Rev Biol. 1982;57:109–34.View ArticleGoogle Scholar
- Cornette R, Gotoh H, Koshikawa S, Miura T. Juvenile hormone titers and caste differentiation in the damp-wood termite Hodotermopsis sjostedti (Isoptera, Termopsidae). J Insect Physiol. 2008;54:922–30.View ArticlePubMedGoogle Scholar
- Miura T, Scharf ME. Molecular basis underlying caste differentiation in termites. In: Bignell DE, Roisin Y, Lo N, editors. Biology of termites: a modern synthesis. Dordrecht: Springer; 2011. p. 211–53.Google Scholar
- Konopova B, Smykal V, Jindra M. Common and distinct roles of juvenile hormone signaling genes in metamorphosis of holometabolous and hemimetabolous insects. PLoS ONE. 2011;6:e28728.View ArticlePubMedPubMed CentralGoogle Scholar
- Smykal V, Daimon T, Kayukawa T, Takaki K, Shinoda T, Jindra M. Importance of juvenile hormone signaling arises with competence of insect larvae to metamorphose. Dev Biol. 2014;390:221–30.View ArticlePubMedGoogle Scholar
- Roisin Y. Morphology, development and evolutionary significance of the working stages in the caste system of Prorhinotermes (Insecta, Isoptera). Zoomorphology. 1988;107:339–47.View ArticleGoogle Scholar
- Parmentier D, Roisin Y. Caste morphology and development in Termitogeton nr. planus (Insecta, Isoptera, Rhinotermitidae). J Morphol. 2003;255:69–79.View ArticlePubMedGoogle Scholar
- Roisin Y, Pasteels JM. Replacement of reproductives in Nasutitermes princeps (Desneux) (Isoptera: Termitidae). Behav Ecol Sociobiol. 1986;18:437–42.View ArticleGoogle Scholar
- Roisin Y. Intragroup conflicts and the evolution of sterile castes in termites. Am Nat. 1994;143:751–63.View ArticleGoogle Scholar
- Koshikawa S, Matsumoto T, Miura T. Regressive molt in the Japanese damp-wood termite Hodotermopsis japonica (Isoptera: Termopsidae). Sociobiology. 2001;38:495–500.Google Scholar
- Myles TG. Reproductive soldiers in the Termopsidae (Isoptera). Pan-Pac Entomol. 1986;62:293–9.Google Scholar
- Thorne BL, Breisch NL, Muscedere ML. Evolution of eusociality and the soldier caste in termites: influence of intraspecific competition and accelerated inheritance. Proc Natl Acad Sci. 2003;100:12808–13.View ArticlePubMedGoogle Scholar
- Morimoto K. On the damp-wood termite genus Zootermopsis introduced to Japan. Shiroari. 2000;122:3–8 (in Japanese).Google Scholar
- Miura T, Koshikawa S, Matsumoto T. Winged presoldiers induced by a juvenile hormone analog in Zootermopsis nevadensis: implications for plasticity and evolution of caste differentiation in termites. J Morphol. 2003;257:22–32.View ArticlePubMedGoogle Scholar
- Cornette R, Matsumoto T, Miura T. Histological analysis of fat body development and molting events during soldier differentiation in the damp-wood termite, Hodotermopsis sjostedti (Isoptera, Termopsidae). Zool Sci. 2007;24:1066–74.View ArticlePubMedGoogle Scholar
- Terrapon N, Li C, Robertson HM, Ji L, Meng X, Booth W, Gowin J. Molecular traces of alternative social organization in a termite genome. Nat Commun. 2014;5:3636.View ArticlePubMedGoogle Scholar
- Sugime Y, Oguchi K, Gotoh H, Hayashi Y, Matsunami M, Shigenobu S, Koshikawa S, Miura T. Termite soldier mandibles are elongated by dachshund under hormonal and Hox gene controls. Development 2019;146:dev171942View ArticlePubMedGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:research0034-1.View ArticleGoogle Scholar
- Andersen CL, Jensen JL, Ørntoft TF. Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004;64:5245–50.View ArticlePubMedGoogle Scholar
- Kaiser HF. The application of electronic computers to factor analysis. Educ Psychol Meas. 1960;20:141–51.View ArticleGoogle Scholar