Microarray analysis of the gene expression profile in the midgut of silkworm infected with cytoplasmic polyhedrosis virus

The mechanisms and factors that induce trained immunity in arthropods and mollusks

Besides dividing the organism's immune system into adaptive and innate immunity, it has long been thought that only adaptive immunity can establish immune memory. However, many studies have shown that innate immunity can also build immunological memory through epigenetic reprogramming and modifications to resist pathogens' reinfection, known as trained immunity. This paper reviews the role of mitochondrial metabolism and epigenetic modifications and describes the molecular foundation in the trained immunity of arthropods and mollusks. Mitochondrial metabolism and epigenetic modifications complement each other and play a key role in trained immunity.

The Intestinal Immune Defense System in Insects

Over a long period of evolution, insects have developed unique intestinal defenses against invasion by foreign microorganisms, including physical defenses and immune responses. The physical defenses of the insect gut consist mainly of the peritrophic matrix (PM) and mucus layer, which are the first barriers to pathogens. Gut microbes also prevent the colonization of pathogens. Importantly, the immune-deficiency (Imd) pathways produce antimicrobial peptides to eliminate pathogens; mechanisms related to reactive oxygen species are another important pathway for insect intestinal immunity. The janus kinase/STAT signaling pathway is involved in intestinal immunity by producing bactericidal substances and regulating tissue repair. Melanization can produce many bactericidal active substances into the intestine; meanwhile, there are multiple responses in the intestine to fight against viral and parasitic infections. Furthermore, intestinal stem cells (ISCs) are also indispensable in intestinal immunity. Only the coordinated combination of the intestinal immune defense system and intestinal tissue renewal can effectively defend against pathogenic microorganisms.

Identification of a Novel Densovirus in Aphid, and Uncovering the Possible Antiviral Process During Its Infection

Densoviruses (DVs) are single-stranded DNA viruses and exclusively happen in invertebrates. Most of DVs reported in insects are pathogenic to their native hosts, however, no pathogenic effect of them has been examined in vertebrates. Hence, DVs are the potential agents used in pest managements. Aphids are the primary vectors of plant viruses. In this study, we identified a novel DV in Chinese Sitobion miscanthi population, provisionally named “Sitobion miscanthi densovirus” (SmDV). Taxonomically, SmDV belongs to genus Hemiambidensovirus. In S. miscanthi, SmDV is hosted in diverse cells and can be horizontally transmitted via wheat feeding. Subject to SmDV, aphids activate their intrinsic antiviral autophagy pathway. Grouped with ascorbate and aldarate metabolism, chlorophyll metabolism, p450 related drug metabolism, and retinoid metabolism, aphids form a complex immune network response to the infection of SmDV. Obviously, it works as elder aphids still alive even they contain the highest examined concentration of SmDV. This study provides a foundation for the identifications of novel DVs, and further improves the understanding of the molecular interactions between insects and DVs.

BmCPV-Derived Circular DNA vcDNA-S7 Mediated by Bombyx mori Reverse Transcriptase (RT) Regulates BmCPV Infection

Circular DNAs derived from single-stranded RNA viruses play important roles in counteracting viral infection. However, whether double-stranded RNA viruses generate functional circular DNAs is still unknown. Using circDNA sequencing, divergent PCR, DNA in situ hybridization and rolling circular amplification, we presently confirmed that in silkworm, Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), a double-stranded RNA virus belonging to cypovirus, is prone to produce a BmCPV-derived circular DNA termed as vcDNA-S7. We have also found that vcDNA-S7 formation is mediated by endogenous reverse transcriptase (RT), and the proliferation of BmCPV can be inhibited by vcDNA-S7 in vitro and in vivo. Moreover, we have discovered that the silkworm RNAi immune pathway is activated by vcDNA-S7, while viral small interfering RNAs (vsiRNAs) derived from transcribed RNA by vcDNA-S7 can be detected by small RNA deep sequencing. These results suggest that BmCPV-derived vcDNA-S7, mediated by RT, can serve as a template for the biogenesis of antiviral siRNAs, which may lead to the repression of BmCPV infection. To our knowledge, this is the first demonstration that a circular DNA, produced by double stranded RNA viruses, is capable of regulating virus infection.

Transcriptional responses of Daphnis nerii larval midgut to oral infection by Daphnis nerii cypovirus-23

BACKGROUND

Daphnis nerii cypovirus-23 (DnCPV-23) is a new type of cypovirus and has a lethal effect on the oleander hawk moth, Daphnis nerii which feeds on leave of Oleander and Catharanthus et al. After DnCPV-23 infection, the change of Daphnis nerii responses has not been reported.

METHODS

To better understand the pathogenic mechanism of DnCPV-23 infection, 3rd-instar Daphnis nerii larvae were orally infected with DnCPV-23 occlusion bodies and the transcriptional responses of the Daphnis nerii midgut were analyzed 72 h post-infection using RNA-seq.

RESULTS

The results showed that 1979 differentially expressed Daphnis nerii transcripts in the infected midgut had been identified. KEGG analysis showed that protein digestion and absorption, Toll and Imd signaling pathway were down-regulated. Based on the result, we speculated that food digestion and absorption in insect midgut might be impaired after virus infection. In addition, the down-regulation of the immune response may make D. nerii more susceptible to bacterial infections. Glycerophospholipid metabolism and xenobiotics metabolism were up-regulated. These two types of pathways may affect the viral replication and xenobiotic detoxification of insect, respectively.

CONCLUSION

These results may facilitate a better understanding of the changes in Daphnis nerii metabolism during cypovirus infection and serve as a basis for future research on the molecular mechanism of DnCPV-23 invasion.

Vegetative Insecticidal Protein (Vip): A Potential Contender From Bacillus thuringiensis for Efficient Management of Various Detrimental Agricultural Pests

Bacillus thuringiensis (Bt) bacterium is found in various ecological habitats, and has natural entomo-pesticidal properties, due to the production of crystalline and soluble proteins during different growth phases. In addition to Cry and Cyt proteins, this bacterium also produces Vegetative insecticidal protein (Vip) during its vegetative growth phase, which is considered an excellent toxic candidate because of the difference in sequence homology and receptor sites from Cry proteins. Vip proteins are referred as second-generation insecticidal proteins, which can be used either alone or in complementarity with Cry proteins for the management of various detrimental pests. Among these Vip proteins, Vip1 and Vip2 act as binary toxins and have toxicity toward pests belonging to Hemiptera and Coleoptera orders, whereas the most important Vip3 proteins have insecticidal activity against Lepidopteran pests. These Vip3 proteins are similar to Cry proteins in terms of toxicity potential against susceptible insects. They are reported to be toxic toward pests, which can’t be controlled with Cry proteins. The Vip3 proteins have been successfully pyramided along with Cry proteins in transgenic rice, corn, and cotton to combat resistant pest populations. This review provides detailed information about the history and importance of Vip proteins, their types, structure, newly identified specific receptors, and action mechanism of this specific class of proteins. Various studies conducted on Vip proteins all over the world and the current status have been discussed. This review will give insights into the significance of Vip proteins as alternative promising candidate toxic proteins from Bt for the management of pests in most sustainable manner.

Host Response against Virus Infection in an Insect: Bidensovirus Infection Effect on Silkworm (Bombyx mori)

Silk cocoons obtained from silkworms are the primary source of commercial silk, making the silkworm an economically important insect. However, the silk industry suffers significant losses due to various virus infections. Bombyx mori bidensovirus (BmBDV) is one of the pathogens that cause flacherie disease in silkworms. Most silkworm strains die after BmBDV infection. However, certain silkworm strains show resistance to the virus, which is determined by a single recessive gene, nsd-2. The +^nsd-2 gene (allele of nsd-2; the susceptibility gene) encodes a putative amino acid transporter expressed only in the insect’s midgut, where BmBDV can infect, suggesting that this membrane protein may function as a receptor for BmBDV. Interestingly, the expression analysis revealed no changes in the +^nsd-2 gene expression levels in virus-uninfected silkworms, whereas the gene expression drastically decreased in the virus-infected silkworm. This condition indicates that the host factor’s expression, the putative virus receptor, is affected by BmBDV infection. It has recently been reported that the expression levels of some host genes encoding cuticle, antioxidant, and immune response-related proteins were significantly regulated by BmBDV infection. In this review, we discuss the host response against virus infection based on our knowledge and long-term research experience in this field.

Heat shock protein 19.9 (Hsp19.9) from Bombyx mori is involved in host protection against viral infection

Adverse environmental conditions cause serious economic losses in sericulture; Bombyx mori nucleopolyhedrovirus (BmNPV) is the primary biotic stress and high temperature is the major abiotic stress in this industry. B. mori heat shock protein 19.9 (Bmhsp19.9) overexpression was previously demonstrated to protect transgenic silkworm H19.9 against extreme temperature. This study analyzed the role of Bmhsp19.9 in H19.9A and H19.9B silkworm lines and BmE cells infected with BmNPV at regular and high temperatures. qPCR results showed that Bmhsp19.9 expression was upregulated in BmE cells and silkworm after BmNPV challenge. Bmhsp19.9 overexpression significantly inhibited BmNPV proliferation in BmE cells. The viral DNA content was significantly decreased in transgenic H19.9 silkworm compared to the control. These results suggested that Bmhsp19.9 was involved in antiviral immunity against BmNPV. Furthermore, Bmhsp19.9 overexpression protected BmE cells against BmNPV under high temperature shock. This indicates that Bmhsp19.9 is a promising candidate for improving silkworm resistance to biotic and abiotic stresses, thereby reducing sericulture losses.

Comprehensive analysis of lncRNA-mRNA regulatory network in BmNPV infected cells treated with Hsp90 inhibitor

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the main pathogens that seriously affect the sustainable development of sericulture industry. Inhibition of Hsp90 by Hsp90 inhibitor, geldanamycin (GA) significantly suppresses BmNPV proliferation in Bombyx mori, while the functional mechanism is not clear. LncRNA has been widely reported to play an important role in immune responses and host-virus interactions in mammalian. However, related research has been rarely reported on silkworm. In this study, firstly, we confirmed the decrease of BmNPV ORF75 protein in the BmNPV-infected BmN cells treated with GA. Next, by using a genome-wide transcriptome analysis, we compared the lncRNA and mRNA expression profiles in BmNPV infected BmN cells treated with or without GA and identified a total of 282 differentially expressed lncRNAs (DElncRNAs) and 523 DEmRNAs. KEGG pathway analysis revealed DEmRNA were mainly involved in ubiquitin mediated proteolysis, spliceosome, RNA transport and oxidative phosphorylation. Further, we selected 27 immune-related DEmRNAs, which displayed the similar changes of expression patterns on both protein level and transcript level to construct DElncRNA-DEmRNA network. In addition, based on the DElncRNA-bmo-miR-278-3p-BmHSC70 regulatory network, we explored the potential function of several lncRNAs as sponges to inhibit the regulatory effect of bmo-278-3p on Bombyx mori heat shock protein cognate 70 (BmHSC70). Our finding suggests that lncRNAs play a role in the regulation of BmNPV proliferation by Hsp90.

Antiviral defense against Cypovirus 1 (Reoviridae) infection in the silkworm, Bombyx mori

Recent years have shown a large increase in studies of infection of the silkworm (Bombyx mori) with Cypovirus 1 (previously designated as B. mori cytoplasmic polyhedrosis virus), that causes serious damage in sericulture. Cypovirus 1 has a single-layered capsid that encapsulates a segmented double-strand RNA (dsRNA) genome which are attractive features for the establishment of a biotechnological platform for the production of specialized gene silencing agents, either as recombinant viruses or as viral-like particles with nonreplicative dsRNA cargo. For both combatting viral disease and application of Cypovirus-based pest control, however, a better understanding is needed of the innate immune response caused by Cypovirus infection of the midgut of lepidopteran larvae. Studies of deep sequencing of viral small RNAs have indicated the importance of the RNA interference pathway in the control of Cypovirus infection although many functional aspects still need to be elucidated and conclusive evidence is lacking. A considerable number of transcriptome studies were carried out that revealed a complex response that hitherto remains uncharacterized because of a dearth in functional studies. Also, the uptake mechanism of Cypovirus by the midgut cells remains unclarified because of contrasting mechanisms revealed by electron microscopy and functional studies. The field will benefit from an increase in functional studies that will depend on transgenic silkworm technology and reverse genetics systems for Cypovirus 1.

Identification and Expression Analysis of Four Small Heat Shock Protein Genes in Cigarette Beetle, Lasioderma serricorne (Fabricius)

Small heat shock proteins (sHsps) are molecular chaperones that play crucial roles in the stress adaption of insects. In this study, we identified and characterized four sHsp genes (LsHsp19.4, 20.2, 20.3, and 22.2) from the cigarette beetle, Lasioderma serricorne (Fabricius). The four cDNAs encoded proteins of 169, 180, 181, and 194 amino acids with molecular weights of 19.4, 20.2, 20.3, and 22.2 kDa, respectively. The four LsHsp sequences possessed a typical sHsp domain structure. Quantitative real-time PCR analyses revealed that LsHsp19.4 and 20.3 transcripts were most abundant in pupae, whereas the transcript levels of LsHsp20.2 and 22.2 were highest in adults. Transcripts of three LsHsp genes were highly expressed in the larval fat body, whereas LsHsp20.2 displayed an extremely high expression level in the gut. Expression of the four LsHsp genes was dramatically upregulated in larvae exposed to 20-hydroxyecdysone. The majority of the LsHsp genes were significantly upregulated in response to heat and cold treatments, while LsHsp19.4 was insensitive to cold stress. The four genes were upregulated when challenged by immune triggers (peptidoglycan isolated from Staphylococcus aureus and from Escherichia coli 0111:B4). Exposure to CO₂ increased LsHsp20.2 and 20.3 transcript levels, but the LsHsp19.4 transcript level declined. The results suggest that different LsHsp genes play important and distinct regulatory roles in L. serricorne development and in response to diverse stresses.

Applying microarray-based technique to study and analyze silkworm (Bombyx mori) transcriptomic response to long-term high iron diet

To investigate the biological processes affected by long-term iron supplementation, newly hatched silkworms were exposed to high iron mulberry diet (10 and 100 ppm) and its effect on silkworm transcriptom was determined. The results showed that the silkworm was responsive to iron by increasing iron concentration and ferritin levels in the hemolymph and by regulating the expression of many other genes. A total of 523 and 326 differentially expressed genes were identified in 10 and 100 ppm Fe group compared to the control, respectively. Of these genes, 249 were shared between in both the 10 ppm and 100 ppm Fe group, including 152 up-regulated and 97 down-regulated genes. These shared genes included 19 known Fe regulated, 24 immune-related, 12 serine proteases and serine proteases homologs, 41 cuticular and cuticle genes. Ten genes (carboxypeptidases A, serine protease homologs 85, fibrohexamerin/P25, transferrin, sex-specific storage-protein 2, fungal protease inhibitor F, insect intestinal mucin, peptidoglycan recognition protein B, cuticle protein CPH45, unknown gene) were involved in the regulation of iron overload responses.

Molecular and biochemical responses in the midgut of the silkworm, Bombyx mori, infected with Nosema bombycis

BACKGROUND

Microsporidia are a group of eukaryotic intracellular parasites that infect almost all vertebrates and invertebrates. However, there is little information available of how microsporidia obtain nutrients and energy from host cells. The purpose of this study was to investigate the energy and material requirements of Nosema bombycis for the invasion procedure through analyzing the global variation of the gene expression, protein abundance, fatty acids level and ATP flux induced by the microsporidia N. bombycis infection in the midgut of the silkworm Bombyx mori.

METHODS

A suppression subtractive hybridization (SSH) and quantitative real-time PCR (qPCR) analysis were performed to identify the genes upregulated in the midgut of B. mori 48 h following N. bombycis infection. Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used to annotate and summarize the differentially expressed genes, according to the categories 'molecular function', 'cellular component' and 'biological process'. To evaluate the nutrition material and energy costs in B.mori infected by N. bombycis, biochemical analysis was performed to determine the variation of protein abundance, fatty acid levels and ATP flux with or without the microsporidia N. bombycis infection in the midgut of the silkworm B. mori.

RESULTS

A total of 744 clones were obtained, 288 clones were randomly selected for sequencing, and 110 unigenes were generated. Amongst these, 49.21%, 30.16% and 14.29% genes were involved in 19 molecular functions, 19 biological processes and nine cellular components, respectively. A total of 11 oxidative phosphorylation- and eight proton-coupled ATP synthesis-related genes were upregulated. Seven protein degradation-, three fat degradation-related genes were upregulated, and no genes related to the de novo synthesis of amino acids and fatty acids were significantly upregulated. The data from the biochemical analysis showed the contents of total protein and ATP of B. mori midgut tissues decreased significantly, whereas the fatty acid content did not significantly change after four days of N. bombycis infection. Microsporidia N. bombycis infection upregulated the expression level of genes involved in host ATP synthesis, protein and fat degradation, which eventually causes the obvious decline of protein content and ATP synthesis in the host midgut, whereas the fatty acids content did not change significantly.

CONCLUSIONS

This study suggested to some extent that N. bombycis invasion can activate the host protein degradation and accelerate the production of host ATP. Microsporidia of N. bombycis show preference for proteins rather than fatty acids from the host to ensure the material preparation required by their parasitic life-cycle. Requirements of N. bombycis for energy were also mainly dependent on the host ATP production. This study provides a new data that may help our understanding of the molecular mechanisms of obtaining energy and nutrients from the host by the microsporidium N. bombycis.

Discovery of anti-viral molecules and their vital functions in Bombyx mori

The silkworm Bombyx mori (B. mori), a lepidopteran model organism, has become an important model for molecular biology researches with its genome completely sequenced. Silkworms confront different types of virus diseases, mainly including those caused by Bombyx mori nucleopolyhedrovirus (BmNPV), Bombyx mori densovirus type 1 (BmDNV-1), Bombyx mori bidesovirus (BmBDV) which was termed as Bombyx mori densovirus type 2 (BmDNV-2) or Bombyx mori parvo-like virus (BmPLV) before in sericulture. B. mori offers excellent models to study the molecular mechanisms of insect innate immune responses to viruses. A variety of molecules and pathways have been identified to be involved in the immune responses in the silkworm to viruses, such as the antimicrobial peptides, prophenoloxidase-activating system, apoptosis, ROS, small RNA and related molecules. Here in this review, we summarize the current research advances in molecules involved in silkworm anti-virus pathways. Moreover, taking BmNPV as an example, we proposed a schematic model of molecules and pathways involved in silkworm immune responses against virus infection. We hope this review can facilitate further study of antiviral mechanisms in silkworm, and provide a reference for virus diseases in other organisms.

DNA methylation in silkworm genome may provide insights into epigenetic regulation of response to Bombyx mori cypovirus infection

DNA methylation is an important epigenetic modification that regulates a wide range of biological processes including immune response. However, information on the epigenetics-mediated immune mechanisms in insects is limited. Therefore, in this study, we examined transcriptomes and DNA methylomes in the fat body and midgut tissues of silkworm, Bombyx mori with or without B. mori cytoplasmic polyhedrosis virus (BmCPV) infection. The transcriptional profile and the genomic DNA methylation patterns in the midgut and fat body were tissue-specific and dynamically altered after BmCPV challenge. KEGG pathway analysis revealed that differentially methylated genes (DMGs) could be involved in pathways of RNA transport, RNA degradation, nucleotide excision repair, DNA replication, etc. 27 genes were shown to have both differential expression and differential methylation in the midgut and fat body of infected larvae, respectively, indicating that the BmCPV infection-induced expression changes of these genes could be mediated by variations in DNA methylation. BS-PCR validated the hypomethylation of G2/M phase-specific E3 ubiquitin-protein ligase-like gene in the BmCPV infected midgut. These results demonstrated that epigenetic regulation may play roles in host-virus interaction in silkworm and would be potential value for further studies on mechanism of BmCPV epithelial-specific infection and epigenetic regulation in the silkworm.

Molecular analysis of Culex quinquefasciatus larvae responses to Lysinibacillus sphaericus Bin toxin

Lysinibacillus sphaericus produces the mosquito larvicidal binary toxin consisting of BinA and BinB, which are both required for toxicity against Culex and Anopheles larvae. The molecular mechanisms behind Bin toxin-induced damage remain unexplored. We used whole-genome microarray-based transcriptome analysis to better understand how Culex larvae respond to Bin toxin treatment at the molecular level. Our analyses of Culex quinquefasciatus larvae transcriptome changes at 6, 12, and 18 h after Bin toxin treatment revealed a wide range of transcript signatures, including genes linked to the cytoskeleton, metabolism, immunity, and cellular stress, with a greater number of down-regulated genes than up-regulated genes. Bin toxin appears to mainly repress the expression of genes involved in metabolism, the mitochondrial electron transport chain, and the protein transporter of the outer/inner mitochondrial membrane. The induced genes encode proteins linked to mitochondrial-mediated apoptosis and cellular detoxification including autophagic processes and lysosomal compartments. This study is, to our knowledge, the first microarray analysis of Bin toxin-induced transcriptional responses in Culex larvae, providing a basis for an in-depth understanding of the molecular nature of Bin toxin-induced damage.

Transcriptome Analysis of the Midgut of the Chinese Oak Silkworm Antheraea pernyi Infected with Antheraea pernyi Nucleopolyhedrovirus

The Antheraea pernyi nucleopolyhedrovirus (ApNPV) is an exclusive pathogen of A. pernyi. The intense interactions between ApNPV and A. pernyi cause a series of physiological and pathological changes to A. pernyi. However, no detailed report exists regarding the molecular mechanisms underlying the interactions between ApNPV and A. pernyi. In this study, four cDNA libraries of the A. pernyi midgut, including two ApNPV-infected groups and two control groups, were constructed for transcriptomic analysis to provide new clues regarding the molecular mechanisms that underlie these interactions. The transcriptome of the A. pernyi midgut was de novo assembled using the Trinity platform because of the lack of a genome resource for A. pernyi. Compared with the controls, a total of 5,172 differentially expressed genes (DEGs) were identified, including 2,183 up-regulated and 2,989 down-regulated candidates, of which 2,965 and 911 DEGs were classified into different GO categories and KEGG pathways, respectively. The DEGs involved in A. pernyi innate immunity were classified into several categories, including heat-shock proteins, apoptosis-related proteins, serpins, serine proteases and cytochrome P450s. Our results suggested that these genes were related to the immune response of the A. pernyi midgut to ApNPV infection via their essential roles in regulating a variety of physiological processes. Our results may serve as a basis for future research not only on the molecular mechanisms of ApNPV invasion but also on the anti-ApNPV mechanism of A. pernyi.

Analysis of gene expression in the midgut of Bombyx mori during the larval molting stage

Background

Insects can be models for understanding human intestinal infection and pathology. Molting, a special period during which the old insect cuticle is shed and a new one is produced, is crucial for insect development. Holometabolous insects may experience several larva-to-larva moltings to become larger, a pupal molt and adult eclosion to become adults. During the larval molts, they stop feeding and become quiescent. Although the molting larvae become quiescent, it is not known if changes in microbiome, physiology, development and immunity of midguts occur.

Results

Transcriptome analysis indicated that functions such as metabolism, digestion, and transport may become reduced due to the downregulated expression of many associated genes. During the molting stage, midguts harbor less microflora and DNA synthesis decreases. Both ecdysone and juvenile hormone in the larval midgut likely degrade after entering the larva-to-larva molting stage. However, at 12 h after ecdysis, the feeding larvae of 5th instars that were injected with 20-hydroxyecdysone entered a molting-like stage, during which changes in midgut morphology, DNA synthesis, gene expression, and microflora exhibited the same patterns as observed in the actual molting state.

Conclusion

This study is important for understanding insect midgut physiology, development and immunity during a special development stage when no food is ingested. Although the molting larva becomes immobile and quiescent, we demonstrate that numerous changes occur in midgut morphology, physiology, metabolism and microbiome during this period.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3162-8) contains supplementary material, which is available to authorized users.

Transcriptome analysis of interactions between silkworm and cytoplasmic polyhedrosis virus

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) specifically infects silkworm midgut (MG) and multiplication occurs mainly in posterior midgut (PM). In this study, MG and fat body (FB) were extracted at 0, 3, 24, and 72 h after BmCPV infection. The total sequence reads of each sample were more than 1510000, and the mapping ratio exceeded 95.3%. Upregulated transcripts increased in MG during the infection process. Gene ontology (GO) categories showed that antioxidants were all upregulated in FB but not in MG. BGI001299, BGI014434, BGI012068, and BGI009201 were MG-specific genes with transmembrane transport function, the expression of which were induced by BmCPV. BGI001299, BGI014434, and BGI012068 expressed in entire MG and may be involved in BmCPV invasion. BGI009201 expressed only in PM and may be necessary for BmCPV proliferation. BmPGRP-S2 and BGI012452 (a putative serine protease) were induced by BmCPV and may be involved in immune defense against BmCPV. The expression level of BmCPV S1, S2, S3, S6, and S7 was high and there was no expression of S9 in MG 72 h, implying that the expression time of structural protein coding genes is earlier. These results provide insights into the mechanism of BmCPV infection and host defense.

Lack of Connection Between Midgut Cell Autophagy Gene Expression and BmCPV Infection in the Midgut of Bombyx mori

Autophagy is associated with multiple biological processes and has protective and defensive functions with respect to immunity, inflammation, and resistance to microbial infection. In this experiment, we wished to investigate whether autophagy is a factor in the midgut cell response of Bombyx mori to infection by the B. mori cytoplasmic polyhedrosis virus (BmCPV). Our results indicated that the expression of three autophagy-related genes (BmAtg8, BmAtg5, and BmAtg7) in the midgut did not change greatly after BmCPV infection in B. mori. Basal ATG8/ATG8PE protein expression was detected in different B. mori tissues by using western blot analysis. Immunohistochemistry showed that the ATG8/ATG8PE proteins were located mainly in the cytoplasm. ATG8/ATG8PE protein levels decreased at 12 and 16 h after BmCPV infection. Our results indicate that autophagy responded slightly to BmCPV infection, but could not prevent the invasion and replication of the virus.

Proteomics analysis of digestive juice from silkworm during Bombyx mori nucleopolyhedrovirus infection

Previous studies have analyzed the midgut transcriptome and proteome after challenge with Bombyx mori nucleopolyhedrovirus (BmNPV), however little information is available on the digestive juice proteome after BmNPV challenge. This study investigated BmNPV infection-induced protein changes in the digestive juice of silkworms using shotgun proteomics and MS sequencing. From the digestive juice of normal third-day, fifth-instar silkworm larvae, 75 proteins were identified, 44 of which were unknown; from larvae 6 h after inoculation with BmNPV, 106 proteins were identified, of which 39 were unknown. After BmNPV challenge, more secreted proteins appeared that had antiviral and digestive features. GO annotation analysis clustered most proteins in the lumen into catalytic, binding, and metabolic processes. Numerous proteins were reported to have BmNPV interactions. Hsp70 protein cognate, lipase-1, and chlorophyllide A-binding protein precursor were upregulated significantly after BmNPV challenge. Levels of trypsin-like serine protease, beta-1,3-glucanase, catalase, and serine protease transcripts decreased or were not significantly change after BmNPV challenge. Taken together, these findings provided insights into the interaction between host and BmNPV and revealed potential functions of digestive juice after per os BmNPV infection.

The gene expression profile of resistant and susceptible Bombyx mori strains reveals cypovirus-associated variations in host gene transcript levels

High-throughput paired-end RNA sequencing (RNA-Seq) was performed to investigate the gene expression profile of a susceptible Bombyx mori strain, Lan5, and a resistant B. mori strain, Ou17, which were both orally infected with B. mori cypovirus (BmCPV) in the midgut. There were 330 and 218 up-regulated genes, while there were 147 and 260 down-regulated genes in the Lan5 and Ou17 strains, respectively. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment for differentially expressed genes (DEGs) were carried out. Moreover, gene interaction network (STRING) analyses were performed to analyze the relationships among the shared DEGs. Some of these genes were related and formed a large network, in which the genes for B. mori cuticular protein RR-2 motif 123 (BmCPR123) and the gene for B. mori DNA replication licensing factor Mcm2-like (BmMCM2) were key genes among the common up-regulated DEGs, whereas the gene for B. mori heat shock protein 20.1 (Bmhsp20.1) was the central gene among the shared down-regulated DEGs between Lan5 vs Lan5-CPV and Ou17 vs Ou17-CPV. These findings established a comprehensive database of genes that are differentially expressed in response to BmCPV infection between silkworm strains that differed in resistance to BmCPV and implied that these DEGs might be involved in B. mori immune responses against BmCPV infection.

Transcriptome analysis of Bombyx mori larval midgut during persistent and pathogenic cytoplasmic polyhedrosis virus infection

Many insects can be persistently infected with viruses but do not show any obvious adverse effects with respect to physiology, development or reproduction. Here, Bombyx mori strain Daizo, persistently infected with cytoplasmic polyhedrosis virus (BmCPV), was used to study the host's transcriptional response after pathogenic infection with the same virus in midgut tissue of larvae persistently and pathogenically infected as 2nd and 4th instars. Next generation sequencing revealed that from 13,769 expressed genes, 167 were upregulated and 141 downregulated in both larval instars following pathogenic infection. Several genes that could possibly be involved in B. mori immune response against BmCPV or that may be induced by the virus in order to increase infectivity were identified, whereas classification of differentially expressed transcripts (confirmed by qRT-PCR) resulted in gene categories related to physical barriers, immune responses, proteolytic/metabolic enzymes, heat-shock proteins, hormonal signaling and uncharacterized proteins. Comparison of our data with the available literature (pathogenic infection of persistently vs. non-persistently infected larvae) unveiled various similarities of response in both cases, which suggests that pre-existing persistent infection does not affect in a major way the transcriptome response against pathogenic infection. To investigate the possible host's RNAi response against BmCPV challenge, the differential expression of RNAi-related genes and the accumulation of viral small RNAs (vsRNAs) were studied. During pathogenic infection, siRNA-like traces like the 2-fold up-regulation of the core RNAi genes Ago-2 and Dcr-2 as well as a peak of 20 nt small RNAs were observed. Interestingly, vsRNAs of the same size were detected at lower rates in persistently infected larvae. Collectively, our data provide an initial assessment of the relative significance of persistent infection of silkworm larvae on the host response following pathogenic infection with CPV, while they also highlight the relative importance of RNAi as an antiviral mechanism.

Differentially expressed genes in the fat body of Bombyx mori in response to phoxim insecticide

The silkworm, Bombyx mori, is an economically important insect. However, poisoning of silkworms by organophosphate pesticides causes tremendous loss to the sericulture. The fat body is the major tissue involved in detoxification and produces antimicrobial peptides and regulates hormones. In this study, a microarray system comprising 22,987 oligonucluotide 70-mer probes was employed to examine differentially expressed genes in the fat body of B. mori exposed to phoxim insecticide. The results showed that a total of 774 genes were differentially expressed upon phoxim exposure, including 500 up-regulated genes and 274 down-regulated genes. The expression levels of eight detoxification-related genes were up-regulated upon phoxim exposure, including six cytochrome P450s and two glutathione-S-transferases. It was firstly found that eight antimicrobial peptide genes were down-regulated, which might provide important references for studying the larvae of B. mori become more susceptible to microbial infections after phoxim treatment. In addition, we firstly detected the expression level of metamorphosis-related genes after phoxim exposure, which may lead to impacted reproduction. Our results may facilitate the overall understanding of the molecular mechanism of multiple pathways following exposure to phoxim insecticide in the fat body of B. mori.

Comparative transcriptome analysis of the lichen-forming fungus Endocarpon pusillum elucidates its drought adaptation mechanisms

The lichen-forming fungus was isolated from the desert lichen Endocarpon pusillum that is extremely drought resistant. To understand the molecular mechanisms of drought resistance in the fungus, we employed RNA-seq and quantitative real-time PCR to compare and characterize the differentially expressed genes in pure culture at two different water levels and with that in desiccated lichen. The comparative transcriptome analysis indicated that a total of 1781 genes were differentially expressed between samples cultured under normal and PEG-induced drought stress conditions. Similar to those in drought resistance plants and non-lichenized fungi, the common drought-resistant mechanisms were differentially expressed in E. pusillum. However, the expression change of genes involved in osmotic regulation in E. pusillum is different, which might be the evidence for the feature of drought adaptation. Interestingly, different from other organisms, some genes involved in drought adaption mechanisms showed significantly different expression patterns between the presence and absence of drought stress in E. pusillum. The expression of 23 candidate stress responsive genes was further confirmed by quantitative real-time PCR using dehydrated E. pusillum lichen thalli. This study provides a valuable resource for future research on lichen-forming fungi and shall facilitate future functional studies of the specific genes related to drought resistance.

sHsp22.6, an intronless small heat shock protein gene, is involved in stress defence and development in Apis cerana cerana

Small heat shock proteins (sHSPs) play an important role in protecting against stress-induced cell damage and fundamental physiological processes. In this study, we identified an intronless sHsp gene from Apis cerana cerana (AccsHsp22.6). The open reading frame of AccsHsp22.6 was 585 bp and encoded a 194 amino acid protein. Furthermore, a 2064 bp 5'-flanking region was isolated, and potential transcription factor binding sites associated with development and stress response were identified. Quantitative PCR and western blot analyses demonstrated that AccsHsp22.6 was detected at higher levels in the midgut than in other tissues tested, and it is highly expressed during the shift to different development stages. Moreover, AccsHsp22.6 was significantly up-regulated by abiotic and biotic stresses, such as 4 °C, 16 °C, 42 °C, cyhalothrin, pyridaben, H2O2, UV, CdCl2, 20-hydroxyecdysone and Ascosphaera apis treatments. However, AccsHsp22.6 was slightly repressed by other stresses, including 25 °C, phoxim, paraquat and HgCl2 treatments. The recombinant AccsHSP22.6 also exhibited significant temperature tolerance, antioxidation and molecular chaperone activity. In addition, we found that knockdown of AccsHsp22.6 by RNA interference remarkably reduced temperature tolerance in A. cerana cerana. Taken together, these results suggest that AccsHsp22.6 plays an essential role in the development stages and defence against cellular stress.

Transfection of BmCPV genomic dsRNA in silkmoth-derived Bm5 cells: stability and interactions with the core RNAi machinery

While several studies have been conducted to investigate the stability of dsRNA in the extracellular medium (hemolymph, gut content, saliva), little is known regarding the persistence of dsRNA once it has been introduced into the cell. Here, we investigate the stability of Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) genomic dsRNA fragments after transfection into Bombyx-derived Bm5 cells. Using RT-PCR as a detection method, we found that dsRNA could persist for long periods (up to 8 days) in the intracellular environment. While the BmCPV genomic dsRNA was processed by the RNAi machinery, its presence had no effects on other RNAi processes, such as the silencing of a luciferase reporter by dsLuc. We also found that transfection of BmCPV genomic dsRNA could not establish a viral infection in the Bm5 cells, even when co-transfections were carried out with dsRNAs targeting Dicer and Argonaute genes, suggesting that the neutralization by RNAi does not play a role in the establishment of an in vitro culture system. The mechanism of the dsRNA stability in Bm5 cells is discussed, as well as the implications for the establishment for an in vitro culture system for BmCPV.

Transcriptome analysis of silkworm, Bombyx mori, during early response to Beauveria bassiana challenges

Host–pathogen interactions are complex processes and it is a central challenge to reveal these interactions. Fungal infection of silkworm, Bombyx mori, may induce a variety of responsive reaction. However, little is known about the molecular mechanism of silkworm immune response against the fungal infection. To obtain an overview of the interaction between silkworm and an entomopathogenic fungus Beauveria bassiana, Digital Gene Expression profiling, a tag based high-throughput transcriptome sequencing method, was employed to screen and identify differentially expressed genes (DEGs, FDR≤0.001, ∣log₂ratio∣≥1) of silkworm larvae during early response against B. bassiana infection. Total 1430 DEGs including 960 up-regulated and 470 down-regulated ones were identified, of which 627 DEGs can be classified into GO categories by Gene Ontology (GO) analysis. KEGG pathways analysis of these DEGs suggested that many biological processes, such as defense and response, signal transduction, phagocytosis, regulation of gene expression, RNA splicing, biosynthesis and metabolism, protein transport etc. were involved in the interaction between the silkworm and B. bassiana. A number of differentially expressed fungal genes were also identified by mapping the sequencing tags to B. bassiana genome. These results provided new insights to the molecular mechanism of silkworm immune response to B. bassiana infection.

Comparison of RNA isolation methods from insect larvae

Isolating RNA from insects is becoming increasingly important in molecular entomology. Four methods including three commercial kits RNeasy Mini Kit (Qiagen), SV Total RNA isolation system (Promega), TRIzol reagent (Invitrogen), and a cetyl trimethylammonium bromide (CTAB)-based method were compared regarding their ability to isolate RNA from whole-body larvae of Thaumatotibia leucotreta (Meyrick), Thanatophilus micans (F.), Plutella xylostella (L.), and Tenebrio molitor (L.). A difference was observed among the four methods regarding RNA quality but not quantity. However, RNA quality and quantity obtained was not dependent on the insect species. The CTAB-based method produced low-quality RNA and the Trizol reagent produced partially degraded RNA, whereas the RNeasy Mini Kit and SV Total RNA isolation system produced RNA of consistently high quality. However, after reverse transcription to cDNA, RNA produced using all four extraction methods could be used to successfully amplify a 708 bp fragment of the cytochrome oxidase I gene. Of the four methods, the SV Total RNA isolation system showed the least amount of DNA contamination with the highest RNA integrity number and is thus recommended for stringent applications where high-quality RNA is required. This is the first comparison of RNA isolation methods among different insect species and the first to compare RNA isolation methods in insects in the last 20 years.

Comprehensive analysis of gene expression profiles of the beet armyworm Spodoptera exigua larvae challenged with Bacillus thuringiensis Vip3Aa toxin

Host-pathogen interactions result in complex relationship, many aspects of which are not completely understood. Vip proteins, which are Bacillus thuringensis (Bt) insecticidal toxins produced during the vegetative stage, are selectively effective against specific insect pests. This new group of Bt proteins represents an interesting alternative to the classical Bt Cry toxins because current data suggests that they do not share the same mode of action. We have designed and developed a genome-wide microarray for the beet armyworm Spodoptera exigua, a serious lepidopteran pest of many agricultural crops, and used it to better understand how lepidopteran larvae respond to the treatment with the insecticidal protein Vip3Aa. With this approach, the goal of our study was to evaluate the changes in gene expression levels caused by treatment with sublethal doses of Vip3Aa (causing 99% growth inhibition) at 8 and 24 h after feeding. Results indicated that the toxin provoked a wide transcriptional response, with 19% of the microarray unigenes responding significantly to treatment. The number of up- and down-regulated unigenes was very similar. The number of genes whose expression was regulated at 8 h was similar to the number of genes whose expression was regulated after 24 h of treatment. The up-regulated sequences were enriched for genes involved in innate immune response and in pathogen response such as antimicrobial peptides (AMPs) and repat genes. The down-regulated sequences were mainly unigenes with homology to genes involved in metabolism. Genes related to the mode of action of Bt Cry proteins were found, in general, to be slightly overexpressed. The present study is the first genome-wide analysis of the response of lepidopteran insects to Vip3Aa intoxication. An insight into the molecular mechanisms and components related to Vip intoxication will allow designing of more effective management strategies for pest control.

The possible impact of persistent virus infection on the function of the RNAi machinery in insects: a hypothesis

RNAi experiments in insects are characterized by great variability in efficiency; for instance beetles and locusts are very amenable to dsRNA-mediated gene silencing, while other insect groups, most notably lepidopterans, are more refractory to RNAi. Several factors can be forwarded that could affect the efficiency of RNAi, such as the composition and function of the intracellular RNAi machinery, the mechanism of dsRNA uptake, the presence of dsRNA- and siRNA-degrading enzymes and non-specific activation of the innate immune response. In this essay, we investigate the evidence whether persistent infection with RNA viruses could be a major factor that affects the response to exogenous dsRNA in insects. The occurrence of RNA viruses in different insect groups will be discussed, as well as several mechanisms by which viruses could interfere with the process of RNAi. Finally, the impact of RNA virus infection on the design of dsRNA-based insect control strategies will be considered.

Characteristics of phoxim-exposed gene transcription in the silk gland of silkworms

Silkworm (Bombyx mori), a model Lepidoptera insect, is an important economic insect. Its silk gland is the important organ for silk protein synthesis and secretion. Phoxim exposure causes deficient cocooning of silkworm and has become one of the major negative factors for the silk industry. To study the impact of phoxim exposure on silk gland, using gene chip technology, we examined differentially expressed genes in silk gland after silkworms were exposed to phoxim (4.0μg/mL) for 24h. Functional annotation, classification and KEGG signaling pathway analysis were performed. The results showed that out of 3206 genes detected in silk gland after phoxim exposure, 270 were differentially expressed significantly, including 249 up-regulated genes and 21 down-regulated genes. These differentially expressed genes related to apoptosis, detoxification and protein degradation were selected. Using qRT-PCR, the expression levels of 9 genes involved in apoptosis, detoxification and protein degradation were validated. In addition, the expression profiles of three related fibroin synthesis genes (Fib-H, Fib-L and P25) were analyzed. Our results showed that phoxim exposure induced apoptosis of silk gland cells and inhibition of fibroin synthesis. This may be the cause of deficient silkworm cocooning.

Involvement of microRNAs in infection of silkworm with bombyx mori cytoplasmic polyhedrosis virus (BmCPV)

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) is one of the most important pathogens of silkworm. MicroRNAs (miRNAs) have been demonstrated to play key roles in regulating host-pathogen interaction. However, there are limited reports on the miRNAs expression profiles during insect pathogen challenges. In this study, four small RNA libraries from BmCPV-infected midgut of silkworm at 72 h post-inoculation and 96 h post-inoculation and their corresponding control midguts were constructed and deep sequenced. A total of 316 known miRNAs (including miRNA*) and 90 novel miRNAs were identified. Fifty-eight miRNAs displayed significant differential expression between the infected and normal midgut (P value < = 0.01 and fold change > = 2.0 or < = 0.5), among which ten differentially expressed miRNA were validated by qRT-PCR method. Further bioinformatics analysis of predicted target genes of differentially expressed miRNAs showed that the miRNA targets were involved in stimulus and immune system process in silkworm.

Microarray analysis of gene expression profile in resistant and susceptible Bombyx mori strains reveals resistance-related genes to nucleopolyhedrovirus

To investigate the molecular mechanism of silkworm resistance to BmNPV infection, we constructed a near-isogenic line (BC8) with BmNPV resistance using highly resistant (NB) and highly susceptible parental strains (306). We investigated variations in the gene expression in the midguts of BmNPV-infected BC8 and 306 at 12 h pi using the microarray. 92 differentially expressed genes were identified. Real-time qPCR analysis confirmed that 10 genes were significantly up-regulated or down-regulated in the midguts of BC8 and NB compared to 306. To our knowledge, we first defined the role of the amino acid transporter and 26S proteasome in insect antiviral. However, serine protease was not completely consistent with data of reported previously in insect antiviral. The role of the 5 genes (Bm123, Bm122, COP β', aquaporin, glycoside hydrolases) was also demonstrated in insect antiviral. Our results provided new insights into the molecular mechanism of the Bombyx mori immune response against BmNPV infection.

Transcriptional characteristics of gene expression in the midgut of domestic silkworms (Bombyx mori) exposed to phoxim

Silkworm (Bombyx mori) is not only an economically important insect but also a model system for lepidoptera. As a vital organ of digestion and nutrient absorption, the midgut of insects also serves as the first physiological barrier to chemical pesticides. In this study, microarray was performed to profile the gene expression changes in the midgut of silkworms exposed to phoxim. After 24h of phoxim exposure (4.0μg/mL), 266 genes displayed at least 2.0-fold changes in expression levels. Among them, 192 genes were up-regulated, and 74 genes were down-regulated. The most significant changes were 14.88-fold up-regulation and 23.36-fold down-regulation. According to gene ontology annotation and pathway analysis, differentially expressed genes were mainly classified into different groups based on their potential involvements in detoxification, immunne response, stress response, energy metabolism and transport. Particularly, the transcription levels of detoxification-related genes were up-regulated, such as cytochrome P450s, esterases and glutathione-S-transferase (GST), indicating increased detoxification activity in the midgut. Our study provides new insights into the molecular mechanism of pesticide metabolism in the midgut of insects, which may promote the development of highly efficient insecticides.

cDNA cloning and characterization of LASP1 from silkworm, Bombyx mori, involved in cytoplasmic polyhedrosis virus infection

Full-length cDNA of a LIM and SH3 contained protein 1 (named BmLASP1) was identified from the silkworm, Bombyx mori, for the first time by rapid amplification of cDNA ends. The full-length cDNA of BmLASP1 is 2094 bp, consisting of a 5'-terminal untranslated region (UTR) of 117 bp, and a 3'-UTR of 610 bp with two poly-adenylation signal sequence AATAAA and a poly (A) tail. The BmLASP1 cDNA encodes a polypeptide comprising 455 amino acids, including a LIM domain, two nebulin domains and an SH3 domain. The theoretical isoelectric point is 7.07 and the predicted molecular weight is 51.8 kDa. BmLASP1 has no signal peptide but three potential N-glycosylation sites. Sequence similarity and phylogenic analyses indicated that BmLASP1 belonged to the group of insect LASP1 with a longer linker region which is different from vertebrate LASP1. The LASP1 in silkworm contained eight exons in its coding regions, and the last exon-intron boundary was conserved the same as in mammalian and Ciona intestinalis LASP1 genes. By fluorescent quantitative real-time polymerase chain reaction, the mRNA transcripts of BmLASP1 were mainly detected in the gonad, head, and spiracle, and slightly in the silk gland, vasa mucosa, midgut, fat body, and hemocytes. After silkworm larvae were infected by B. mori cytoplasmic polyhedrosis virus (BmCPV), the relative expression level of BmLASP1 was down-regulated in the midgut. This result suggested that BmLASP1 may play an important role in the response of silkworm to BmCPV infection.

Molecular cloning and characterization of Hsp27.6: the first reported small heat shock protein from Apis cerana cerana

Small heat shock proteins (sHSPs) play an important role in the cellular defense of prokaryotic and eukaryotic organisms against a variety of internal and external stressors. In this study, a cDNA clone encoding a member of the α-crystallin/sHSP family, termed AccHsp27.6, was isolated from Apis cerana cerana. The full-length cDNA is 1,014 bp in length and contains a 708-bp open reading frame encoding a protein of 236 amino acids with a calculated molecular weight of 27.6 kDa and an isoelectric point of 7.53. Seven putative heat shock elements and three NF-κB binding sites were present in the 5'-flanking region, suggesting a possible function in immunity. A semi-quantitative RT-PCR analysis indicated that AccHsp27.6 was expressed in all tested tissues and at different developmental stages. Furthermore, expression of the AccHsp27.6 transcript was induced by exposure to heat shock, H(2)O(2), a number of different chemicals (including SO(2), formaldehyde, alcohol, acetone, chloroform, and the pesticides phoxime and acetamiprid), and the microbes Staphylococcus aureus and Micrococcus luteus. In contrast, the mRNA expression could be repressed by CO(2), the pesticides pyriproxyfen and cyhalothrin, and the microbes Bacillus subtilis and Pseudomonas aeruginosa. Notably, the recombinant AccHsp27.6 protein exhibited significant in vitro molecular chaperone activity and antimicrobial activity. Taken together, these results suggest that AccHsp27.6 might play an important role in the response to abiotic and biotic stresses and in immune reactions.

Transcriptional profiling of midgut immunity response and degeneration in the wandering silkworm, Bombyx mori

BACKGROUND

Lepidoptera insects have a novel development process comprising several metamorphic stages during their life cycle compared with vertebrate animals. Unlike most Lepidoptera insects that live on nectar during the adult stage, the Bombyx mori silkworm adults do not eat anything and die after egg-laying. In addition, the midguts of Lepidoptera insects produce antimicrobial proteins during the wandering stage when the larval tissues undergo numerous changes. The exact mechanisms responsible for these phenomena remain unclear.

PRINCIPAL FINDINGS

We used the silkworm as a model and performed genome-wide transcriptional profiling of the midgut between the feeding stage and the wandering stage. Many genes concerned with metabolism, digestion, and ion and small molecule transportation were down-regulated during the wandering stage, indicating that the wandering stage midgut loses its normal functions. Microarray profiling, qRT-PCR and western blot proved the production of antimicrobial proteins (peptides) in the midgut during the wandering stage. Different genes of the immune deficiency (Imd) pathway were up-regulated during the wandering stage. However, some key genes belonging to the Toll pathway showed no change in their transcription levels. Unlike butterfly (Pachliopta aristolochiae), the midgut of silkworm moth has a layer of cells, indicating that the development of midgut since the wandering stage is not usual. Cell division in the midgut was observed only for a short time during the wandering stage. However, there was extensive cell apoptosis before pupation. The imbalance of cell division and apoptosis probably drives the continuous degeneration of the midgut in the silkworm since the wandering stage.

CONCLUSIONS

This study provided an insight into the mechanism of the degeneration of the silkworm midgut and the production of innate immunity-related proteins during the wandering stage. The imbalance of cell division and apoptosis induces irreversible degeneration of the midgut. The Imd pathway probably regulates the production of antimicrobial peptides in the midgut during the wandering stage.