Genomics and proteomics combined analysis revealed the toxicity response of silkworm Bombyx mori to the environmental pathogen Bacillus cereus ZJ-4.

Bacterial contamination has caused a major public health problem worldwide. Bacillus cereus is a conditional environmental pathogenic bacteria that can cause food poisoning. Whether environmental pathogens can cause widespread transmission in the insect kingdom is unclear. In this study, a Bacillus cereus ZJ-4 was isolated from the hospital environment of Zhenjiang City, Jiangsu Province, China. It was fatal by injection into the silkworm hemolymph. To investigated the potential toxic factors of ZJ-4 and clarified the toxicity response mechanism of silkworm by the ZJ-4 infection. Then, the whole genome of ZJ-4 was sequenced, and the immune mechanism of silkworm fat body to ZJ-4 pathogen was studied by HE pathological section and proteomics. Bacterial genome sequencing indicated that ZJ-4 had 352 drug resistance genes and 6 virulence genes. After 36 h of subcutaneous puncture with ZJ-4 suspension, the pathological changes were obviously found in HE pathological sections of fat body tissue. Comparative proteomic results indicated that differentially expressed proteins are mainly involved in stress reactions, biological regulation, and innate immunity. The qRT-PCR analysis showed that the expressions of ß-GRP, Spaetzle, MyD88, Tube and Dorsal genes in Toll pathway were up-regulated, while Pell and Cactus genes were down-regulated; in the antimicrobial peptide pathway, Glv2, Lzm, Mor, and Leb3 genes were up-regulated, while attacin1 and defensin genes were down-regulated; Sod gene was up-regulated, while Cat gene was down-regulated in the antioxidant pathway; Ldh, Sdh, and Mdh genes were down-regulated in glucose metabolism pathway. These results indicated that ZJ-4 can damage the innate immune pathway of silkworm, and also affect the normal immune function of fat body cells.

Increased expression of Suppressor of cytokine signaling 2 (BmSOCS2) is correlated with suppression of Bombyx mori nucleopolyhedrovirus replication in silkworm larval tissues and cells.

The resistance of silkworm to infection by Bombyx mori nuclear polyhedrosis virus (BmNPV) is a main focus of sericultural research. Previously, a BmNPV-resistant strain, NB, was identified among a collection of Chinese silkworm strains in our lab. To better understand the molecular mechanism of NB strain resistance, the patterns of host immune response gene transcription in resistant (NB) and susceptible (306) strains were examined. Quantative real-time PCR (qRT-PCR) revealed that multiple insect innate immune signaling pathways (Toll, Imd and JAK/STAT) were strongly activated upon infection with BmNPV. Notably, Suppressor of cytokine signaling 2 (BmSOCS2) mRNA expression was significantly up-regulated in midgut tissues of the resistant NB strain, suggesting that the BmSOCS2 gene product may be involved in host immune defense against BmNPV infection. A significant inhibition of BmNPV replication was also observed in BmN cells transfected with a vector encoding BmSOCS2. The results suggest that BmSOCS2 is a key gene involved in the resistance of the NB silkworm strain to BmNPV infection.

High-Quality de novo Chromosome-Level Genome Assembly of a Single Bombyx mori With BmNPV Resistance by a Combination of PacBio Long-Read Sequencing, Illumina Short-Read Sequencing, and Hi-C Sequencing.

The reference genomes of Bombyx mori (B. mori), Silkworm Knowledge-based database (SilkDB) and SilkBase, have served as the gold standard for nearly two decades. Their use has fundamentally shaped model organisms and accelerated relevant studies on lepidoptera. However, the current reference genomes of B. mori do not accurately represent the full set of genes for any single strain. As new genome-wide sequencing technologies have emerged and the cost of high-throughput sequencing technology has fallen, it is now possible for standard laboratories to perform full-genome assembly for specific strains. Here we present a high-quality de novo chromosome-level genome assembly of a single B. mori with nuclear polyhedrosis virus (BmNPV) resistance through the integration of PacBio long-read sequencing, Illumina short-read sequencing, and Hi-C sequencing. In addition, regular bioinformatics analyses, such as gene family, phylogenetic, and divergence analyses, were performed. The sample was from our unique B. mori species (NB), which has strong inborn resistance to BmNPV. Our genome assembly showed good collinearity with SilkDB and SilkBase and particular regions. To the best of our knowledge, this is the first genome assembly with BmNPV resistance, which should be a more accurate insect model for resistance studies.

Transcriptomics and proteomics-based analysis of heterosis on main economic traits of silkworm, Bombyx mori.

The application of silkworm hybrids have promoted the innovation and development of agricultural technology, but the mechanism of heterosis in silkworm has not been explained clearly. In this study, the heterosis of silkworm in the aspects of body weight, silk gland and cocoon weight was investigated by means of silkworm hybridization and multi-omics approaches, including transcriptome and proteome. The results showed that heterosis of silkworm body weight, silk gland and cocoon weight was overdominant, but only part of genes and proteins were overdominant, and most of genes and proteins were non-additive. Combined analysis obtained six up-regulated genes and four down-regulated genes that were consistent both in transcriptome and proteome. Gene functional enrichment analysis indicated that most up-regulated genes and proteins were mostly related to metabolism, which led to accelerated metabolism and protein synthesis and contributing to improved heterosis. The up-regulation of 6-phosphate glucose dehydrogenase (G6PDH), phosphatidylethanolamine-binding protein (PEBP) and sHSP20.4, which are involved in metabolism, might be related to silk gland heterosis. SIGNIFICANCE: A combination of transcriptomic and proteomic analysis was used to understand the molecular mechanism of silkworm heterosis. We found that the phenotypic traits of silkworm are overdominant, while the analysis of transcriptome and proteome showed that only part of genes and proteins were overdominant, and most of genes and proteins were non-additive. Some of the genes had unique expression in F1, which was speculated that genes under heterozygous condition may result in rearrangement and cause metabolic changes in the hybrids. Those both up-regulated in transcriptomic and proteomic analysis were found to be involved in various metabolic processes, so as to accelerate metabolism and protein synthesis, thus exhibiting heterosis.