Map-based cloning and functional analysis revealed ABCC2 is responsible for Cry1Ac toxin resistance in Bombyx mori.

Bt toxins are parasporal crystals produced by Bacillus thuringiensis (Bt). They have specific killing activity against various insects and have been widely used to control agricultural pests. However, their widespread use has developed the resistance of many target insects. To maintain the sustainable use of Bt products, the resistance mechanism of insects to Bt toxins must be fully clarified. In this study, Bt-resistant and Bt-susceptible silkworm strains were used to construct genetic populations, and the genetic pattern of silkworm resistance to Cry1Ac toxin was determined. Sequence-tagged site molecular marker technology was used to finely map the resistance gene and to draw a molecular genetic linkage map, and the two closest markers were T1590 and T1581, indicating the resistance gene located in the 155 kb genetic region. After analyzing the sequence of the predicted gene in the genetic region, an ATP binding cassette transporter (ABCC2) was identified as the candidate gene. Molecular modeling and protein-protein docking result showed that a tyrosine insertion in the mutant ABCC2 might be responsible for the interaction between Cry1Ac and ABCC2. Moreover, CRISPR/Cas9-mediated genome editing technology was used to knockout ABCC2 gene. The homozygous mutant ABCC2 silkworm was resistant to Cry1Ac toxin, which indicated ABCC2 is the key gene that controls silkworm resistance to Cry1Ac toxin. The results have laid the foundation for elucidating the molecular resistance mechanism of silkworms to Cry1Ac toxin and could provide a theoretical basis for the biological control of lepidopteran pests.

Differentially expressed genes in the cuticle and hemolymph of the silkworm, Bombyx mori, injected with the fungus Beauveria bassiana.

The most important pathogenic fungus of the silkworm, Bombyx mori L. (Lepidoptera: Bombycidae), is Beauveria bassiana (Balsamo-Crivelli ) Vuillemin (Hypocreales: Clavicipitaceae), which causes significant damage to sericulture production. Therefore, diagnosing fungal disease and developing new control measures are crucial to silk production. To better understand the responsive and interactive mechanisms between the host silkworm and this fungus, variations in silkworm gene expression were investigated using the suppression subtractive hybridization method following the injection of B. bassiana conidia. Two cDNA libraries were constructed, and 140 cDNA clones were isolated. Of the 50 differentially expressed genes identified, 45 (112 clones) were identified in the forward library, and 5 (28 clones) were identified in the reverse library. Expression profiling of six of these genes by quantitative polymerase chain reaction (qPCR) verified that they were induced by the fungal challenge. The present study provides insight into the interaction between lepidopteran insects and pathogenic fungi.

[Molecular mapping of test mapping strain for 18th linkage group recessive genes elp, ch-2 and mln in silkworm (Bombyx mori)].

The ellipsoid egg, the second recessive gene of chocolate larvae, and melanism are controlled by three recessive genes, elp, ch-2, and mln in silkworm, respectively. Their order and genetic distance have been scheduled in established linkage group. Owing to lack of crossing over in females, the reciprocal backcrossed F1(BC1) progenies were bred for linkage analysis using the wild type silkworm strain P50(+elp+ch-2+mln /+elp+ch-2+mln) and W18 with ellipsoid egg, the second recessive gene of chocolate larvae, and melanism (elp ch-2 mln / elp ch-2 mln). In this research, we mapped three mutant genes, elp, ch-2, and mln on the chromosome 18 based on the SSR linkage map and STS markers designed based on silkworm genome sequence. The established linkage group, molecular linkage group, and the physic map of chromosome 18 had been corresponded. The genetic distance for this chromosome in this research was 94.2 cM, and the order of the mutants and molecular markers were consistent with the established silkworm linkage maps and the fine genome sequences. This research will lay important bases for map-based cloning for other mutants on chromosome 18.

CRISPR/Cas9-based knockout reveals that the clock gene timeless is indispensable for regulating circadian behavioral rhythms in Bombyx mori.

Circadian rhythms, which are ubiquitous and adaptive, occur across all species, from microbes to humans, in which they organize and modify behavior and physiology. timeless (tim) is a canonical clock gene. The core composition of the Drosophila melanogaster endogenous circadian clock has been extensively investigated; however, in lepidopteran insects, including Bombyx mori, the mechanism is complicated and little is known regarding the participation of tim in the negative feedback loop responsible for behavioral activities. To arrive at a comprehensive understanding of the role of tim in the B. mori endogenous circadian clock, we exploited the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 gene editing system. We attempted to elucidate the functions of tim in the circadian clock of B. mori using Bmtim mutants. The knockouts affected two circadian behavioral activities: adult emergence and embryo hatching rhythms. Quantitative real-time polymerase chain reaction results confirmed that tim-knockouts induced relative reductions in the expression levels, and thereby the oscillation amplitudes, of Bmper and Bmclk messenger RNAs during both the photophase and scotophase. Additionally, the daily rhythmic expression of Bmdbt was upregulated in the photophase and downregulated in the scotophase in a tim-knockout. Our study reveals that tim is integral to the B. mori circadian clock and may be involved in regulating eclosion and hatching rhythms.

A palmitoyltransferase Approximated gene Bm-app regulates wing development in Bombyx mori.

The silkworm Bombyx mori is an important lepidopteran model insect in which many kinds of natural mutants have been identified. However, molecular mechanisms of most of these mutants remain to be explored. Here we report the identification of a gene Bm-app is responsible for the silkworm minute wing (mw) mutation which exhibits exceedingly small wings during pupal and adult stages. Compared with the wild type silkworm, relative messenger RNA expression of Bm-app is significantly decreased in the u11 mutant strain which shows mw phenotype. A 10 bp insertion in the putative promoter region of the Bm-app gene in mw mutant strain was identified and the dual luciferase assay revealed that this insertion decreased Bm-app promoter activity. Furthermore, clustered regularly interspaced short palindromic repeats/RNA-guided Cas9 nucleases-mediated depletion of the Bm-app induced similar wing defects which appeared in the mw mutant, demonstrating that Bm-app controls wing development in B. mori. Bm-app encodes a palmitoyltransferase and is responsible for the palmitoylation of selected cytoplasmic proteins, indicating that it is required for cell mitosis and growth during wing development. We also discuss the possibility that Bm-app regulates wing development through the Hippo signaling pathway in B. mori.

BmSd gene regulates the silkworm wing size by affecting the Hippo pathway.

Insect wings are developed from the wing disc during metamorphosis. Bombyx mori, a model lepidopteran insect, loses flight ability after long-term domestication from the wild silkworm, Bombyx mandarina. The mw mutant (u11 strain) shows minute wings compared to wild type (e.g., p50 strain) wings. RNA sequencing analysis previously revealed differential Hippo-pathway-related gene expression between the u11 and p50 strains. The Hippo pathway is an evolutionarily conserved signaling cascade that controls organ size during development in animals. In this study, the function of BmSd which has been characterized as one of the Hippo-pathway-related genes was analyzed for silkworm wing development. We found that mats, warts, and hippo expression levels were higher in u11 compared to p50 wing discs. BmSd (scalloped) expression, which encodes a prominent transcriptional partner to Yorkie (Yki), gradually decreased during the wandering stage in u11, but exhibited the opposite expression pattern in p50. When BmSd was knocked down by small interfering RNA during the wandering stage in the p50 strain, 57.9% of the individuals showed minute wings. Additionally, ex, kibra, and wingless expression levels decreased in the BmSd knockdown mutant. Further, BmSd deletion mediated by clustered regularly interspaced short palindromic repeats (CRISPR) / CRISPR-associated protein 9 induced 50% of individuals with minute wings, a phenotype similar to the mw mutant. This result demonstrates that BmSd plays pivotal roles in silkworm wing development. Our results show that the Hippo signaling pathway participates and plays crucial roles in the regulation of silkworm wing development, and our findings provide a basis for further research on B. mori wing development.

Functional analysis of a miRNA-like small RNA derived from Bombyx mori cytoplasmic polyhedrosis virus.

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) is a major pathogen of the economic insect silkworm, Bombyx mori. Virus-encoded microRNAs (miRNAs) have been proven to play important roles in host-pathogen interactions. In this study we identified a BmCPV-derived miRNA-like 21 nt small RNA, BmCPV-miR-1, from the small RNA deep sequencing of BmCPV-infected silkworm larvae by stem-loop quantitative real-time PCR (qPCR) and investigated its functions with qPCR and lentiviral expression systems. Bombyx mori inhibitor of apoptosis protein (BmIAP) gene was predicted by both target prediction software miRanda and Targetscan to be one of its target genes with a binding site for BmCPV-miR-1 at the 5' untranslated region. It was found that the expression of BmCPV-miR-1 and its target gene BmIAP were both up-regulated in BmCPV-infected larvae. At the same time, it was confirmed that BmCPV-miR-1 could up-regulate the expression of BmIAP gene in HEK293T cells with lentiviral expression systems and in BmN cells by transfecting mimics. Furthermore, BmCPV-miR-1 mimics could up-regulate the expression level of BmIAP gene in midgut and fat body in the silkworm. In the midgut of BmCPV-infected larvae, BmCPV-miR-1 mimics could be further up-regulated and inhibitors could lower the virus-mediated expression of BmIAP gene. With the viral genomic RNA segments S1 and S10 as indicators, BmCPV-miR-1 mimics could up-regulate and inhibitors down-regulate their replication in the infected silkworm. These results implied that BmCPV-miR-1 could inhibit cell apoptosis in the infected silkworm through up-regulating BmIAP expression, providing the virus with a better cell circumstance for its replication.

Identification and characterization of two putative microRNAs encoded by Bombyx mori cypovirus.

Viral microRNAs (miRNAs) have been demonstrated to play important roles in virus-host interactions. Some RNA virus-encoded miRNAs have been reported to promote viral replication and may be used as potential drug targets. Bombyx mori cypovirus (BmCPV), an important pathogen of silkworm, is a double-stranded RNA virus frequently causing serious damages in sericulture. Research on miRNA encoded by BmCPV may be useful to elucidate the BmCPV-host interaction and to develop new anti-viral methods. In our previous study, small RNA libraries of the midgut of BmCPV-infected silkworm have been generated by deep sequencing and several BmCPV-encoded putative miRNAs were predicted. In this study, two putative miRNAs encoded by BmCPV were identified and then validated by stem-loop qRT-PCR and northern blot. They are BmCPV-miR-3 encoded by the third genomic RNA segment of BmCPV (478-497bp) and BmCPV-miR-5 encoded by the fifth genomic RNA segment (2481-2500bp), both are 20bp and encoded by ORF regions. miRNA expression could be detected as early as 5h after BmCPV infection, and the expression level of BmCPV-miR-3 is much higher than that of BmCPV-miR-5 in the course of infection. Three potential target genes were predicted in the host genome, two for BmCPV-miR-3 and one for BmCPV-miR-5, but just one in the virus genome for BmCPV-miR-3 only, with the binding sites all in coding regions. Dual luciferase assay and qRT-PCR indicated that BmCPV-miR-3 could down-regulate the expression of both its two target genes, but no regulatory effect by BmCPV-miR-5 on its target gene was detected. In contrast, BmCPV-miR-3 could up-regulate the viral target. This is the first report that an insect double stranded RNA virus may generate miRNAs and the results obtained will benefit the future study of the functions of BmCPV-encoded miRNAs on viral replication and virus-host interaction.

iTRAQ-based quantitative proteomic analysis of midgut in silkworm infected with Bombyx mori cytoplasmic polyhedrosis virus.

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) specifically infects the epithelial cells in the midgut of silkworm and causes them to death, which negatively affects the sericulture industry. In order to determine the midgut response at the protein levels to the virus infection, differential proteomes of the silkworm midgut responsive to BmCPV infection were identified with isobaric tags for relative and absolute quantitation (iTRAQ) labeling followed by liquid chromatography-tandem mass spectrometry (LC-MS/MS). 193, 408, 189 differentially expressed proteins (DEPs) were reliably quantified by iTRAQ analysis in the midgut of BmCPV-infected and control larvae at 24, 48, 72h post infection (hpi) respectively. KEGG enrichment analysis showed that Oxidative phosphorylation, amyotrophic lateral sclerosis, Toll-like receptor signaling pathway, steroid hormone biosynthesis were the significant pathways (Q value≤0.05) both at 24 and 48hpi. qRT-PCR was used to further verify gene transcription of 30 DEPs from iTRAQ, showing that the regulations of 24 genes at the transcript level were consistent with those at the proteomic level. Moreover, the cluster analysis of the three time groups showed that there were seven co-regulated DEPs including BGIBMGA002620-PA, which was a putative p62/sequestosome-1 protein in silkworm. It was upregulated at both the mRNA level and the proteomic level and may play an important role in regulating the autophagy and apoptosis (especially apoptosis) induced by BmCPV infection. This was the first report using an iTRAQ approach to analyze proteomes of the silkworm midgut against BmCPV infection, which contributes to understanding the defense mechanisms of silkworm midgut to virus infection. SIGNIFICANCE: The domesticated silkworm, Bombyx mori, is renowned for silk production as well as being a traditional lepidopteron model insect served as a subject for morphological, genetic, physiological, and developmental studies. Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) specifically infects the epithelial cells in the midgut of silkworm and causes the silkworm to death, which negatively affects the sericulture industry. Studies on insect antiviral immunity and on interactive mechanisms between host cells and BmCPV are in their infancy and remain insufficient. In order to obtain an overall view of silkworm response to BmCPV infection, we performed a proteomic analysis of the midgut of silkworm responses to BmCPV infection by iTRAQ. This was the first report using an iTRAQ approach to analyze proteomes of the silkworm midgut against BmCPV infection, which contributes to understanding the defense mechanisms of silkworm midgut to virus infection.

JAK/STAT signaling pathway-mediated immune response in silkworm (Bombyx mori) challenged by Beauveria bassiana.

Innate immunity was critical in insects defensive system and able to be induced by Janus kinase/signal transducer and activator of transcription cascade transduction (JAK/STAT) signaling pathway. Currently, it had been identified many JAK/STAT signaling pathway-related genes in silkworm, but little function was known on insect innate immunity. To explore the roles of JAK/STAT pathway in antifungal immune response in silkworm (Bombyx mori) against Beauveria bassiana infection, the expression patterns of B. mori C-type lectin 5 (BmCTL5) and genes encoding 6 components of JAK/STAT signaling pathway in silkworm challenged by B. bassiana were analyzed using quantitative real time PCR. Meanwhile the activation of JAK/STAT signaling pathway by various pathogenic micro-organisms and the affect of JAK/STAT signaling pathway inhibitors on antifungal activity in silkworm hemolymph was also detected. Moreover, RNAi assay of BmCTL5 and the affect on expression levels of signaling factors were also analyzed. We found that JAK/STAT pathway could be obviously activated in silkworm challenged with B. bassiana and had no response to bacteria and B. mori cytoplasmic polyhedrosis virus (BmCPV). However, the temporal expression patterns of JAK/STAT signaling pathway related genes were significantly different. B. mori downstream receptor kinase (BmDRK) might be a positive regulator of JAK/STAT signaling pathway in silkworm against B. bassiana infection. Moreover, antifungal activity assay showed that the suppression of JAK/STAT signaling pathway by inhibitors could significantly inhibit the antifungal activity in hemolymph and resulted in increased sensitivity of silkworm to B. bassiana infection, indicating that JAK/STAT signaling pathway might be involved in the synthesis and secretion of antifungal substances. The results of RNAi assays suggested that BmCTL5 might be one pattern recognition receptors for JAK/STAT signaling pathway in silkworm. These findings yield insights for better understand the molecular mechanisms of JAK/STAT signaling pathway in antifungal immune response in silkworm.

dsRNA interference on expression of a RNA-dependent RNA polymerase gene of Bombyx mori cytoplasmic polyhedrosis virus.

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) is one of the major viral pathogens in silkworm. Its infection often results in significant losses to sericulture. Studies have demonstrated that RNAi is one of the important anti-viral mechanisms in organisms. In this study, three dsRNAs targeting the RNA-dependent RNA polymerase (RDRP) gene of BmCPV were designed and synthesized with 2'-F modification to explore their interference effects on BmCPV replication in silkworm larvae. The results showed that injecting dsRNA in the dosage of 4-6 ng per mg body weight into the 5th instar larvae can interfere with the BmCPV-RDRP expression by 93% after virus infection and by 99.9% before virus infection. In addition, the expression of two viral structural protein genes (genome RNA segments 1 and 5) was also decreased with the decrease of RDRP expression, suggesting that RNAi interference of BmCPV-RDRP expression could affect viral replication. The study provides an effective method for investigating virus replication as well as the virus-host interactions in the silkworm larvae using dsRNA.

Cytoplasmic polyhedrosis virus-induced differential gene expression in two silkworm strains of different susceptibility.

Digital gene expression (DGE) was performed to investigate the gene expression profiles of 4008 and p50 silkworm strains at 48 h after oral infection with BmCPV. 3,668,437 clean tags were identified in the BmCPV-infected p50 silkworms and 3,540,790 clean tags in the control p50. By contrast, 4,498,263 clean tags were identified in the BmCPV-infected 4008 silkworms and 4,164,250 clean tags in the control 4008. A total of 691 differentially expressed genes were detected in the infected 4008 DGE library and 185 were detected in the infected p50 DGE library, respectively. The expression profiles identified some important differentially expressed genes involved in signal transduction, enzyme activity and apoptotic changes, some of which were verified using quantitative real-time PCR (qRT-PCR). These results provide important clues on the molecular mechanism of BmCPV invasion and resistance mechanism of silkworms against BmCPV infection.

Cloning and expression analysis of a peptidoglycan recognition protein in silkworm related to virus infection.

In this study, the full-length cDNA of a peptidoglycan recognition protein named BmPGRP-S3 was identified from the silkworm, Bombyx mori by rapid amplification of cDNA ends. It is 807 bp and comprises the following: a 5'-untranslated region (UTR) with a length of 112 bp, a 3'-UTR with a length of 92 bp including a poly-adenylation signal sequence (AATAAA) and a poly(A) tail. The longest open reading frame (ORF) of BmPGRP-S3 is 603 bp and encodes a polypeptide of 200 amino acids with a predicted molecular weight of 22.3 kDa including a PGRP domain. Sequence similarity and phylogenic analysis results indicated that BmPGRP-S3 belongs to the group of insect PGRPs and is closer to BmPGRP-S4 with the highest identity of 68%. Fluorescent quantitative real-time PCR results revealed that the mRNA transcripts of BmPGRP-S3 were presented in all of the tissues, but were highest in the midgut. In the silkworm larvae infected with B. mori cytoplasmic polyhedrosis virus (BmCPV), the relative expression level of BmPGRP-S3 was upregulated. The DNA segment of a mature BmPGRP-S3 peptide was inserted into the expression plasmid pET-28a(+) to construct a recombinant expression plasmid. Western blot results revealed that mature BmPGRP-S3 could be detected in the hemolymph and midgut which were the most important immune tissues in silkworm. All the results suggested that BmPGRP-S3 may play an important role in the immune response of silkworm to BmCPV infection and provided helpful information for further studying the function of BmPGRP-S3 in silkworm.