De novo transcriptome analysis for examination of the nutrition metabolic system related to the evolutionary process through which stick insects gain the ability of flight (Phasmatodea).

OBJECTIVE: Insects are the most evolutionarily successful groups of organisms, and this success is largely due to their flight ability. Interestingly, some stick insects have lost their flight ability despite having wings. To elucidate the shift from wingless to flying forms during insect evolution, we compared the nutritional metabolism system among flight-winged, flightless-winged, and flightless-wingless stick insect groups. RESULTS: Here, we report RNA sequencing of midgut transcriptome of Entoria okinawaensis, a prominent Japanese flightless-wingless stick insect, and the comparative analysis of its transcriptome in publicly available midgut transcriptomes obtained from seven stick insect species. A gene enrichment analysis for differentially expressed genes, including those obtained from winged vs wingless and flight vs flightless genes comparisons, revealed that carbohydrate metabolic process-related genes were highly expressed in the winged stick insect group. We also found that the expression of the mitochondrial enolase superfamily member 1 transcript was significantly higher in the winged stick insect group than in the wingless stick insect group. Our findings could indicate that carbohydrate metabolic processes are related to the evolutionary process through which stick insects gain the ability of flight.

Construction of TUATinsecta database that integrated plant and insect database for screening phytophagous insect metabolic products with medicinal potential.

Phytophagous insect larvae feed on plants containing secondary metabolic products with biological activity against other predatory organisms. Phytophagous insects can use their specialised metabolic systems to covert these secondary metabolic products into compounds with therapeutic properties useful to mankind. Some Asians drink tea decoctions made from phytophagous insect frass which is believed to be effective against inflammatory diseases. However, insects that can convert plant-derived secondary metabolic products into useful human therapeutic agents remain poorly studied. Here, we constructed the TUATinsecta database by integrating publicly plant/insect datasets for the purpose of selecting insect species. Using TUAT-insecta we selected the Asian swallowtail butterfly, Papilio xuthus larvae fed on several species of Rutaceous plants and examined whether the plant-derived secondary metabolites, especially those present in frass, were chemically altered or not. We extracted metabolic products from frass using three organic solvents with different polarities, and evaluated solvent fractions for their cytotoxic effects against several human cell lines. We found that chloroform frass extracts from P. xuthus larvae fed on Poncirus trifoliata leaves contained significant cytotoxic activity. Our findings demonstrate that screening of insect species using the 'TUATinsecta' database provides an important pipeline for discovering novel therapeutic agents that might be useful for mankind.

Analysis of molecular mechanism for acceleration of polyembryony using gene functional annotation pipeline in Copidosoma floridanum.

BACKGROUND: Polyembryony is defined as the formation of several embryos from a single egg. This phenomenon can occur in humans, armadillo, and some endoparasitoid insects. However, the mechanism underlying polyembryogenesis in animals remains to be elucidated. The polyembryonic parasitoid wasp Copidosoma floridanum oviposits its egg into an egg of the host insect; eventually, over 2000 individuals will arise from one egg. Previously, we reported that polyembryogenesis is enhanced when the juvenile hormone (JH) added to the culture medium in the embryo culture. Hence, in the present study, we performed RNA sequencing (RNA-Seq) analysis to investigate the molecular mechanisms controlling polyembryogenesis of C. floridanum. Functional annotation of genes is not fully available for C.floridanum; however, whole genome assembly has been archived. Hence, we constructed a pipeline for gene functional annotation in C. floridanum and performed molecular network analysis. We analyzed differentially expressed genes between control and JH-treated molura after 48 h of culture, then used the tblastx program to assign whole C. floridanum transcripts to human gene. RESULTS: We obtained 11,117 transcripts in the JH treatment group and identified 217 differentially expressed genes compared with the control group. As a result, 76% of C. floridanum transcripts were assigned to human genes. Gene enrichment analysis revealed genes associated with platelet degranulation, fatty acid biosynthesis, cell morphogenesis in the differentiation and integrin signaling pathways were fluctuated following JH treatment. Furthermore, Cytoscape analysis revealed a molecular interaction that was possibly associated with polyembryogenesis . CONCLUSIONS: We have constructed a pipeline for gene functional annotation of C. floridanum, and identified transcripts with high similarity to human genes during early embryo developmental. Additionally, this study reveals new molecular interactions associated with polyembryogenesis; these interactions could indicate the molecular mechanisms underlying polyembryony. Our results highlight the potential utility of molecular interaction analysis in human twins.

Superoxide dismutase down-regulation and the oxidative stress is required to initiate pupation in Bombyx mori.

Perhaps, oxidative stress progresses pupation in some Lepidopteran insects; however, the reasons for this remain obscure. In our previous study, we clarified Bombyx mori SOD1 (BmSOD1) and B. mori SOD2 (BmSOD2) proteins respond in common to ultraviolet irradiation (UV) oxidative stress and metamorphosis. This result strongly suggested pupation initiates by oxidative stress and might mediate by down-regulation of expression of BmSOD1 and BmSOD2 proteins. Thus, we examined about these relationships in B. mori in this study. In the microarray data reanalysis, we found the Notch signaling pathways as the common pathways in pupation and UV oxidative stress in B. mori. Also, we showed a molting hormone, 20-hydroxyecdysone, leads not only generation of superoxide but also downregulation of the expression of BmSOD proteins during pupation in B. mori. Our findings can contribute to a deeper understanding of how biological defense systems work against environmental oxidative stress.

Mass killing by female soldier larvae is adaptive for the killed male larvae in a polyembryonic wasp.

Self-sacrifice is very rare among organisms. Here, we report a new and astonishing case of adaptive self-sacrifice in a polyembryonic parasitic wasp, Copidosoma floridanum. This wasp is unique in terms of its larval cloning and soldier larvae. Male clone larvae have been found to be killed by female soldier larvae, which suggests intersexual conflict between male and female larvae. However, we show here that mass killing is adaptive to all the killed males as well as the female soldiers that have conducted the killing because the killing increases their indirect fitness by promoting the reproduction of their clone sibs. We construct a simple model that shows that the optimal number of surviving males for both male and female larvae is very small but not zero. We then compare this prediction with the field data. These data agree quite well with the model predictions, showing an optimal killing rate of approximately 94-98% of the males in a mixed brood. The underlying mechanism of this mass kill is almost identical to the local competition for mates that occurs in other wasp species. The maternal control of the sex ratio during oviposition, which is well known in other hymenopterans, is impossible in this polyembryonic wasp. Thus, this mass kill is necessary to maximize the fitness of the female killers and male victims, which can be seen as an analogy of programmed cell death in multicellular organisms.

Comparative analysis of seven types of superoxide dismutases for their ability to respond to oxidative stress in Bombyx mori.

Insects are well adapted to changing environmental conditions. They have unique systems for eliminating reactive oxygen species (ROS). Superoxide dismutase (SOD) is a key enzyme that plays a primary role in removing ROS. Bombyx mori is a lepidopteran insect, whose body size is larger than the model insect Drosophila melanogaster, which enabled us to more easily examine gene expression at the tissue level. We searched B. mori SOD (BmSOD) genes using genome database, and we analyzed their function under different type of oxidative stress. Consequently, we identified four new types of BmSODs in addition to the three types already known. Two of the seven types had a unique domain architecture that has not been discovered previously in the SOD family, and they were expressed in different tissues and developmental stages. Furthermore, these BmSODs responded differently to several kinds of stressors. Our results showed that the seven types of BmSODs are likely to play different roles in B. mori; therefore, B. mori could be used to distinguish the functions of each SOD for resistance to oxidative stress that changes with the environmental conditions.

Construction of a simple evaluation system for the intestinal absorption of an orally administered medicine using Bombyx mori larvae.

Human intestinal absorption is estimated using a human colon carcinoma cell line (Caco-2) cells from human colorectal adenocarcinoma, intestinal perfusion, or a mammalian model. These current evaluation systems are limited in their ability to estimate human intestinal absorption. In addition, in vivo evaluation systems using laboratory animals such as mice and rats entail animal ethics problems, and it is difficult to screen compounds on a large scale at the drug discovery stage. Thus, we propose the use of Bombyx mori larvae for evaluation of intestinal absorption of compounds as an alternative system in this study. First, to compare the characteristics among Caco-2 cells, human intestine, and B. mori larval midgut, we analyzed their RNA-seq data, and we found 26 drug transporters common to humans and B. mori. Next, we quantitatively developed an oral administration technique in B. mori and established a method using silkworm B. mori larvae that can easily estimate the intestinal permeability of compounds. Consequently, we could determine the dose and technique for oral administration in B. mori larvae. We also developed a B. mori model to evaluate the intestinal permeability of orally administered. Our constructed evaluation system will be useful for evaluating intestinal permeability in medical drug development.

Effects of pre-overwintering conditions on eupyrene and apyrene spermatogenesis after overwintering in Polygonia c-aureum (Lepidoptera: Nymphalidae).

Sperm polymorphism is widely known in invertebrates. In insects, Lepidoptera has two types of sperm: nucleated eupyrene (fertile) sperm and anucleated apyrene (unfertile) sperm. These sperm types are produced during post-embryogenesis, and eupyrene spermatogenesis precedes apyrene spermatogenesis. During overwintering, spermatogenesis stops and a portion of undifferentiated-stage spermatocytes degenerate. After overwintering, spermatogenesis restarts with unaffected spermatogonia. However, how new spermatozoa arise in the adult testes after overwintering is not known in Lepidoptera. In this study, we investigated the spermatogenesis events in the nymphalid butterfly Polygonia c-aureum after overwintering under various environmental conditions. Our results indicate that both eupyrene and apyrene spermatogenesis restart at any stopping stage and sperm of these types are regenerated in no particular order after adult insect overwintering. This suggests that the spermatogenesis occurring after overwintering proceeds without embryogenetic restrictions related to the developmental sequence.

Relevant principal factors affecting the reproducibility of insect primary culture.

The primary culture of insect cells often suffers from problems with poor reproducibility in the quality of the final cell preparations. The cellular composition of the explants (cell number and cell types), surgical methods (surgical duration and surgical isolation), and physiological and genetic differences between donors may be critical factors affecting the reproducibility of culture. However, little is known about where biological variation (interindividual differences between donors) ends and technical variation (variance in replication of culture conditions) begins. In this study, we cultured larval fat bodies from the Japanese rhinoceros beetle, Allomyrina dichotoma, and evaluated, using linear mixed models, the effect of interindividual variation between donors on the reproducibility of the culture. We also performed transcriptome analysis of the hemocyte-like cells mainly seen in the cultures using RNA sequencing and ultrastructural analyses of hemocytes using a transmission electron microscope, revealing that the cultured cells have many characteristics of insect hemocytes.

Identification of functional enolase genes of the silkworm Bombyx mori from public databases with a combination of dry and wet bench processes.

BACKGROUND: Various insect species have been added to genomic databases over the years. Thus, researchers can easily obtain online genomic information on invertebrates and insects. However, many incorrectly annotated genes are included in these databases, which can prevent the correct interpretation of subsequent functional analyses. To address this problem, we used a combination of dry and wet bench processes to select functional genes from public databases. RESULTS: Enolase is an important glycolytic enzyme in all organisms. We used a combination of dry and wet bench processes to identify functional enolases in the silkworm Bombyx mori (BmEno). First, we detected five annotated enolases from public databases using a Hidden Markov Model (HMM) search, and then through cDNA cloning, Northern blotting, and RNA-seq analysis, we revealed three functional enolases in B. mori: BmEno1, BmEno2, and BmEnoC. BmEno1 contained a conserved key amino acid residue for metal binding and substrate binding in other species. However, BmEno2 and BmEnoC showed a change in this key amino acid. Phylogenetic analysis showed that BmEno2 and BmEnoC were distinct from BmEno1 and other enolases, and were distributed only in lepidopteran clusters. BmEno1 was expressed in all of the tissues used in our study. In contrast, BmEno2 was mainly expressed in the testis with some expression in the ovary and suboesophageal ganglion. BmEnoC was weakly expressed in the testis. Quantitative RT-PCR showed that the mRNA expression of BmEno2 and BmEnoC correlated with testis development; thus, BmEno2 and BmEnoC may be related to lepidopteran-specific spermiogenesis. CONCLUSIONS: We identified and characterized three functional enolases from public databases with a combination of dry and wet bench processes in the silkworm B. mori. In addition, we determined that BmEno2 and BmEnoC had species-specific functions. Our strategy could be helpful for the detection of minor genes and functional genes in non-model organisms from public databases.

Sexual complementarity between host humoral toxicity and soldier caste in a polyembryonic wasp.

Defense against enemies is a type of natural selection considered fundamentally equivalent between the sexes. In reality, however, whether males and females differ in defense strategy is unknown. Multiparasitism necessarily leads to the problem of defense for a parasite (parasitoid). The polyembryonic parasitic wasp Copidosoma floridanum is famous for its larval soldiers' ability to kill other parasites. This wasp also exhibits sexual differences not only with regard to the competitive ability of the soldier caste but also with regard to host immune enhancement. Female soldiers are more aggressive than male soldiers, and their numbers increase upon invasion of the host by other parasites. In this report, in vivo and in vitro competition assays were used to test whether females have a toxic humoral factor; if so, then its strength was compared with that of males. We found that females have a toxic factor that is much weaker than that of males. Our results imply sexual complementarity between host humoral toxicity and larval soldiers. We discuss how this sexual complementarity guarantees adaptive advantages for both males and females despite the one-sided killing of male reproductives by larval female soldiers in a mixed-sex brood.

Expression of the fructose receptor BmGr9 and its involvement in the promotion of feeding, suggested by its co-expression with neuropeptide F1 in Bombyx mori.

Insect gustatory receptors (Grs) are members of a large family of proteins with seven transmembrane domains that provide insects with the ability to detect chemical signals critical for feeding, mating, and oviposition. To date, 69 Bombyx mori Grs (BmGrs) genes have been identified via genome studies. BmGr9 has been shown to respond specifically to fructose and to function as a ligand-gated ion channel selectively activated by fructose. However, the sites where this Gr are expressed remain unclear. We demonstrated using reverse transcription (RT)-PCR that BmGr9 is widely expressed in the central nervous system (CNS), as well as oral sensory organs. Additionally, immunohistochemistry was performed using anti-BmGr9 antiserum to show that BmGr9 is expressed in cells of the oral sensory organs, including the maxillary galea, maxillary palps, labrum, and labium, as well as in putative neurosecretory cells of the CNS. Furthermore, double immunohistochemical analysis showed that most BmGr9-expressing cells co-localized with putative neuropeptide F1-expressing cells in the brain, suggesting that BmGr9 is involved in the promotion of feeding behaviors. In addition, a portion of BmGr9-expressing cells in the brain co-localized with cells expressing BmGr6, a molecule of the sugar receptor clade, suggesting that sugars other than fructose are involved in the regulation of feeding behaviors in B. mori larvae.

Comparison between the Amount of Environmental Change and the Amount of Transcriptome Change.

Cells must coordinate adjustments in genome expression to accommodate changes in their environment. We hypothesized that the amount of transcriptome change is proportional to the amount of environmental change. To capture the effects of environmental changes on the transcriptome, we compared transcriptome diversities (defined as the Shannon entropy of frequency distribution) of silkworm fat-body tissues cultured with several concentrations of phenobarbital. Although there was no proportional relationship, we did identify a drug concentration "tipping point" between 0.25 and 1.0 mM. Cells cultured in media containing lower drug concentrations than the tipping point showed uniformly high transcriptome diversities, while those cultured at higher drug concentrations than the tipping point showed uniformly low transcriptome diversities. The plasticity of transcriptome diversity was corroborated by cultivations of fat bodies in MGM-450 insect medium without phenobarbital and in 0.25 mM phenobarbital-supplemented MGM-450 insect medium after previous cultivation (cultivation for 80 hours in MGM-450 insect medium without phenobarbital, followed by cultivation for 10 hours in 1.0 mM phenobarbital-supplemented MGM-450 insect medium). Interestingly, the transcriptome diversities of cells cultured in media containing 0.25 mM phenobarbital after previous cultivation (cultivation for 80 hours in MGM-450 insect medium without phenobarbital, followed by cultivation for 10 hours in 1.0 mM phenobarbital-supplemented MGM-450 insect medium) were different from cells cultured in media containing 0.25 mM phenobarbital after previous cultivation (cultivation for 80 hours in MGM-450 insect medium without phenobarbital). This hysteretic phenomenon of transcriptome diversities indicates multi-stability of the genome expression system. Cellular memories were recorded in genome expression networks as in DNA/histone modifications.

The role of the effector caspases drICE and dcp-1 for cell death and corpse clearance in the developing optic lobe in Drosophila.

In the developing Drosophila optic lobe, cell death occurs via apoptosis and in a distinctive spatio-temporal pattern of dying cell clusters. We analyzed the role of effector caspases drICE and dcp-1 in optic lobe cell death and subsequent corpse clearance using mutants. Neurons in many clusters required either drICE or dcp-1 and each one is sufficient. This suggests that drICE and dcp-1 function in cell death redundantly. However, dying neurons in a few clusters strictly required drICE but not dcp-1, but required drICE and dcp-1 when drICE activity was reduced via hypomorphic mutation. In addition, analysis of the mutants suggests an important role of effecter caspases in corpse clearance. In both null and hypomorphic drICE mutants, greater number of TUNEL-positive cells were observed than in wild type, and many TUNEL-positive cells remained until later stages. Lysotracker staining showed that there was a defect in corpse clearance in these mutants. All the results suggested that drICE plays an important role in activating corpse clearance in dying cells, and that an additional function of effector caspases is required for the activation of corpse clearance as well as that for carrying out cell death.

Establishment and characterization of two new cell lines from the mosquito Armigeres subalbatus (Coquillett) (Diptera: Culicidae).

Armigeres subalbatus (Coquillett) is a medically important mosquito and a model species for immunology research. We successfully established two cell lines from the neonate larvae of A. subalbatus using two different media. To our knowledge, this is the first report of an established Armigeres mosquito cell line. The cell lines, designated as Ar-3 and Ar-13, consist of adherent and diploid cells with compact colonies. Both these cell lines grow slowly after passage at a split ratio of 1:5 and a population doubling time of 2.7 and 3.0 d, respectively. Random amplified polymorphic DNA polymerase chain reaction (RAPD-PCR) was used to confirm that these lines correspond to the species of origin and are clearly distinct from seven other insect cell lines. Furthermore, reverse-transcription PCR was used to demonstrate that the Ar-3 cell line is susceptible to the Japanese encephalitis virus and two insect flaviviruses associated with Culex and Aedes mosquitoes but relatively insensitive to dengue virus. These data indicate that the newly established cell lines are cellular models of A. subalbatus as well as beneficial tools for the propagation of viruses associated with the Armigeres mosquito.

Structure and development of male pheromone gland of longicorn beetles and its phylogenetic relationships within the tribe Clytini.

The male sex pheromone of the longicorn beetle, Xylotrechus pyrrhoderus pyrrhoderus Bates (Cerambycidae: Tribe Clytini) plays an important role in attracting females. This pheromone is produced by the pheromone gland located in the prothorax. However, the detailed structure and underlying developmental process of this gland are still unknown. We investigated the gland structure by using histological analysis and confirmed that the gland consists of the following parts: gland cell mass, a unique spherical space in the cuticle layer, and ductules connecting the gland cells with the spherical space and conducting canals to the outer opening. The gland structure first appeared male-specific in the late pupal stage, during which the epidermal cells began depositing the exocuticle; the development of the gland was completed after adult emergence. Furthermore, we verified the structural equivalents of the X. p. pyrrhoderus male pheromone gland in 11 species of 2 tribes, Clytini and Anaglyptini. The glands of these insects could be classified into four types on the basis of the absence or presence of the spherical space and the division of the gland cell mass layer. Most noteworthy, all the species with the spherical space and division-type gland were restricted to the Xylotrechus clade, as inferred from the molecular phylogenetic analysis. These results suggest that Clytini and Anaglyptini species share a fundamental process of male pheromone gland development, and that the Japanese Xylotrechus species might have established their current status by developing distinct structural features in the male pheromone gland.

Gene expression of protein-coding and non-coding RNAs related to polyembryogenesis in the parasitic wasp, Copidosoma floridanum.

Polyembryony is a unique form of development in which many embryos are clonally produced from a single egg. Polyembryony is known to occur in many animals, but the underlying genetic mechanism responsible is unknown. In a parasitic wasp, Copidosoma floridanum, polyembryogenesis is initiated during the formation and division of the morula. In the present study, cDNA libraries were constructed from embryos at the cleavage and subsequent primary morula stages, times when polyembryogenesis is likely to be controlled genetically. Of 182 and 263 cDNA clones isolated from these embryos, 38% and 70%, respectively, were very similar to protein-coding genes obtained from BLAST analysis and 55 and 65 clones, respectively, were stage-specific. In our libraries we also detected a high frequency of long non-coding RNA. Some of these showed stage-specific expression patterns in reverse transcription quantitative polymerase chain reaction (RT-qPCR) analysis. The stage-specificity of expression implies that these protein-coding and non-coding genes are related to polyembryogenesis in C. floridanum. The non-coding genes are not similar to any known non-coding RNAs and so are good candidates as regulators of polyembryogenesis.

Carbonyl reduction in the biosynthesis of a male sex pheromone secreted by the grape borer Xylotrechus pyrrhoderus.

Males of the cerambycid beetle Xylotrechus pyrrhoderus release a mixture of (S)-2-hydroxy-3-octanone [(S)-1] and (2S,3S)-2,3-octanediol [(2S,3S)-2] as a sex pheromone that attracts conspecific females. The chemical structures of these pheromone components include a common motif and are assumed to be biosynthetically related. Here, we show that deuterated (S)-1, applied on the cuticle of a pronotal pheromone gland, was converted into (2S,3S)-2, that included deuterium atoms, but a reverse conversion did not take place. These results reveal a carbonyl reductase to be active in the pheromone gland, and that the ketol is a biosynthetic precursor of the diol. Males did not produce (R)-1; however, deuterated (R)-1 was converted into (2R,3R)-2, indicating an attack of the enzyme from the opposite side of the hydroxyl group at the 2-position. Furthermore, to understand the substrate specificity of the enzyme, racemates of 2-hydroxy-3-hexanone and 2-hydroxy-3-decanone were synthesized and applied to the gland. Their conversion into the corresponding diols suggests that the enzyme reduces the carbonyl group at the 3-position, regardless of the chain length.

Sex differences in the protection of host immune systems by a polyembryonic parasitoid.

Endoparasitoids have the ability to evade the cellular immune responses of a host and to create an environment suitable for survival of their progeny within a host. Generally, the host immune system is suppressed by endoparasitoids. However, polyembryonic endoparasitoids appear to invade their hosts using molecular mimicry rather than immune system suppression. It is not known how the host immune system is modified by polyembryonic endoparasitoids. Using haemocyte counts and measurement of cellular immune responses, we evaluated modification of the host immune system after separate infestations by a polyembryonic parasitoid (Copidosoma floridanum) and another parasitoid (Glyptapanteles pallipes) and by both together (multi-parasitism). We found that the polyembryonic parasitoid maintains and enhances the host immune system, whereas the other parasitoid strongly suppresses the immune system. Multi-parasitization analysis revealed that C. floridanum cancelled the immune suppression by G. pallipes and strengthened the host immunity. This enhancement was much stronger with male than with female C. floridanum.