The Bombyx mori G protein ß subunit 1 (BmGNß1) gene inhibits BmNPV infection.

G protein ß subunit 1 (GNß1) has several functions, including cell growth regulation, the control of second messenger levels, and ion channel switching. Previous transcriptome analyses in our laboratory have shown that BmGNß1 transcription is reduced following infection with Bombyx mori nucleopolyhedrovirus (BmNPV), but it is unknown what role this gene may have in the host response to BmNPV infection. In this study, the BmGNß1 gene was cloned using the RACE method. After BmNPV infection, BmGNß1 was downregulated in Baiyu strains in tissues such as the hemolymph and midgut. Indirect immunofluorescence showed that BmGNß1 was localized to the cytoplasm. We further constructed a BmGNß1-pIZ/V5-His-mCherry overexpression plasmid and designed siRNA to evaluate the role of BmGNß1 in host response to infection. The results showed that BmGNß1 overexpression inhibited BmNPV proliferation, while knockdown of BmGNß1 was correlated with increased BmNPV proliferation. The siRNA-mediated reduction of BmGNß1 was correlated with an increase in BmNPV infection of BmN cells, increased BmNPV vp39 transcription, and reduced survival time of BmNPV-infected B. mori. Overexpression of BmGNß1 in BmN cells was also correlated with apoptosis and a modification in transcript levels of genes involved in host response to BmNPV infection (PI3K, AKT, Bmp53, BmFOXO, Caspase-1, Bmp21, BmPKN and BmCREB), suggesting that BmGNß1 may influence the apoptotic host response of infected B. mori through the PI3K-AKT pathway. This study provides potential targets and theoretical support for breeding BmNPV-resistant silkworm varieties.

Identification of Key Genes Involved in Resistance to Early Stage of BmNPV Infection in Silkworms.

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the most serious pathogens restricting the sustainable development of the sericulture industry. Currently, there is no effective treatment for BmNPV infection in silkworms, and the mechanism underlying BmNPV resistance in silkworms is also not clear. In this study, comparative transcriptome analyses were carried out in midguts of two silkworm varieties, namely BaiyuN, which is a resistance variety, and Baiyu, which is a susceptible variety, at five different time points (i.e., 0, 1, 3, 6, and 9 h) post-BmNPV infection to detect the early-stage transcriptional changes in these silkworms. In total, 1911 and 1577 differentially expressed genes (DEGs) were identified in the Baiyu and BaiyuN varieties, respectively, involving a total of 48 metabolic pathways. Of these pathways, eight were shared by the Baiyu and BaiyuN varieties in response to BmNPV infection. Notably, four genes (i.e., BGIBMGA08815, BGIBMGA003935, BGIBMGA003571, BGIBMGA010059) were upregulated in the Baiyu variety while downregulated in the BaiyuN variety. The inhibited expression of these four genes in the resistant variety highlighted their potential roles in the resistance of early-stage viral replication. Thus, our study provided a new avenue for the further study of the mechanism underlying BmNPV infection in silkworms and the potential treatment of BmNPV infection.

Study on the fusion of sports and medicine in China from 2012 to 2021: A bibliometric analysis via CiteSpace.

Objective: To analyze the research hot spots and frontiers of studies on of the fusion of sports and medicine (FSM) in China in recent decade via CiteSpace. Methods: Relevant publications related FSM published from January 1, 2012 to December 31, 2021 were obtained from the China National Knowledge Infrastructure (CNKI) database. CiteSpace software was used to analyze the amount of publications, institutions and keywords using standard bibliometric indicators. Results: A total of 729 publications on FSM were identified, and 688 qualified records were included in the final analysis. Between 2012 to 2021, the number of publications showed a trend of growth, albeit with certain fluctuations. The authors of these publications were mainly from universities or colleges with sports background. The institution leading the study was the Beijing Sport University (n = 20), the most prolific (n =12) and most-cited (224 times) author was Guo JJ from Capital University of Physical Education Sports. The journal with most publications on FSM was Contemporary Sports Technology (n =74). The analysis of keywords showed that the "FSM" had the highest frequency (n = 269), "integration of sports and medicine" had the strongest citation bursts (4.82), "national fitness" had the highest centrality (0.97) in recent decade, and 15 clusters of keywords were produced by log-likelihood ratio (all silhouette value >0.9). Conclusion: The findings of this bibliometric study analyse the current status and trends in the FSM in China, which may help to identify hot topics, explore new study directions for scholars and policymakers in the future.

High-resolution silkworm pan-genome provides genetic insights into artificial selection and ecological adaptation.

The silkworm Bombyx mori is an important economic insect for producing silk, the "queen of fabrics". The currently available genomes limit the understanding of its genetic diversity and the discovery of valuable alleles for breeding. Here, we deeply re-sequence 1,078 silkworms and assemble long-read genomes for 545 representatives. We construct a high-resolution pan-genome dataset representing almost the entire genomic content in the silkworm. We find that the silkworm population harbors a high density of genomic variants and identify 7308 new genes, 4260 (22%) core genes, and 3,432,266 non-redundant structure variations (SVs). We reveal hundreds of genes and SVs that may contribute to the artificial selection (domestication and breeding) of silkworm. Further, we focus on four genes responsible, respectively, for two economic (silk yield and silk fineness) and two ecologically adaptive traits (egg diapause and aposematic coloration). Taken together, our population-scale genomic resources will promote functional genomics studies and breeding improvement for silkworm.

Metabolic characterization of hemolymph in Bombyx mori varieties after Bombyx mori nucleopolyhedrovirus infection by GC-MS-based metabolite profiling.

The "Huakang 2" silkworm variety, bred by the Sericulture Research Institute of the Chinese Academy of Agricultural Sciences, is highly resistant to Bombyx mori nucleopolyhedrovirus (BmNPV) and effectively solves the issue of frequent Bombyx mori nuclear polyhedrosis in sericultural production. The molecular mechanism of its resistance to BmNPV, however, is still unknown. The purpose of the present study was therefore to identify these anti-BmNPV mechanisms by using metabolomics in combination with transcriptomics after subcutaneous injection of budded virus (BV) with high concentrations of BmNPV from specimens of the Baiyu N variety (which is highly resistant to BmNPV) and the Baiyu variety (which is sensitive to BmNPV). A total of 375 differential metabolites were identified, which mainly included sugars, acids, amines, alcohols, glycosides, and other small molecules. KEGG enrichment analysis and functional clustering of differential metabolites identified possible metabolic pathways, including tyrosine metabolism, oxidative phosphorylation, and alanine, aspartate, and glutamate metabolism. The differentially expressed genes (DEGs) identified by transcriptome analysis were annotated in KEGG. Association analysis showed that the metabolic pathways of different silkworm varieties are affected differently by BmNPV infection, triggering a series of complex physiological and biochemical changes in the organism. In particular, oxidative phosphorylation might be an essential pathway involved in regulation of disease resistance.

Integration of Transcriptomic and Proteomic Analyses Reveals New Insights into the Regulation of Immune Pathways in Midgut of Samia ricini upon SariNPV Infection.

Samia ricini nucleopolyhedrovirus (SariNPV) is one of the main pathogens of S. ricini sericulture and its infection causes severe impacts on economic sericulture development. To explore and reveal the molecular mechanisms of S. ricini in response to SariNPV infection, we employed RNA sequencing (RNA-seq), adopting isobaric tags for relative and absolute quantitation (iTRAQ), and carried out combination analysis of the obtained differentially expressed genes (DEGs) and proteins (DEPs). Through transcriptome sequencing, a total of 2535 DEGs were detected, and with iTRAQ, 434 DEPs with significant expression difference were identified. Through correlation analysis, we found that the expression trends of 116 differentially expressed proteins were the same as those of differentially expressed genes (including 106 up-regulated and 10 down-regulated). Twenty-five key differentially expressed genes (proteins) involved in several signaling and immune related pathways (mainly involving Toll, Imd, Jak-STAT and Wnt signaling pathways, as well as other immune related pathways) were screened through real-time quantitative PCR. Our results, not only provide insights into the pathogenic mechanism of SariNPV infection in ricin silkworm and the immune response mechanism within the host, but also provide a significant contribution for identifying and preventing diseases caused by SariNPV.

Pathogenicity Detection and Genome Analysis of Two Different Geographic Strains of BmNPV.

The pathogenicity of different concentrations of Bombyx mori nuclear polyhedrosis virus- Zhenjiang strain (BmNPV ZJ) and Yunnan strain (BmNPV YN) was assessed in Baiyu larvae. The structures of the two viral strains were observed by negative-staining electron microscopy, and their proliferation was examined by quantitative polymerase chain reaction (qPCR). The genomic sequences of these two viruses were obtained to investigate the differences in their pathogenicity. The lethal concentration 50 (LC50) of BmNPV ZJ against Baiyu larvae was higher than that of BmNPV YN, indicating a relatively more robust pathogenicity in BmNPV YN. Electron microscopic images showed that the edges of BmNPV YN were clearer than those of BmNPV ZJ. The qPCR analysis demonstrated significantly higher relative expressions of immediately early 1 gene (ie-1), p143, vp39, and polyhedrin genes (polh) in BmNPV ZJ than in BmNPV YN at 12-96 h. The complete genomes of BmNPV ZJ and BmNPV YN were, respectively, 135,895 bp and 143,180 bp long, with 141 and 145 coding sequences and 40.93% and 39.71% GC content. Considering the BmNPV ZJ genome as a reference, 893 SNP loci and 132 InDel mutations were observed in the BmNPV YN genome, resulting in 106 differential gene sequences. Among these differential genes, 76 (including 22 hub genes and 35 non-hub genes) possessed amino acid mutations. Thirty genes may have been related to viral genome replication and transcription and five genes may have been associated with the viral oral infection. These results can help in understanding the mechanisms of pathogenicity of different strains of BmNPV in silkworms.

Effect of pyriproxyfen exposure on cocooning and gene expression in the silk gland of Bombyx mori （Linnaeus, 1758）.

Bombyx mori（Linnaeus, 1758） is an important economical insect, and the sericulture is a flourishing industry in many developing countries. Pyriproxyfen, a juvenile hormone pesticide, is often applied to cultivations widely in the world, and its exposure often resulted in silk yield reduction and non-cocooning. However, the effect of pyriproxyfen exposure on cocooning and gene expression level in the silk gland of B. mori has not been studied yet, and this study focused on the above issues. The result indicated that pyriproxyfen exposure can lead to silk gland injury, reduction of silk yield and cocooning rate. Furthermore, the expression levels of silk protein synthesis related genes were down regulated significantly. The same change trends were shown between PI3K/Akt and CncC/Keap1 pathway, which is the expressions of key genes can be elevated by pyriproxyfen exposure. In addition, the activity of detoxification enzymes (P450, GST and CarE) and the expression levels of detoxification genes were elevated after pyriproxyfen exposure, suggesting that detoxification enzymes may play an important role in detoxification of pyriproxyfen in silk gland. These results provided possible clues to the silk gland injury and gene transcriptional level changes in silkworm after pyriproxyfen exposure.

Effects of Pyriproxyfen Exposure on Reproduction and Gene Expressions in Silkworm, Bombyx mori.

The silkworm, Bombyx mori Linnaeus, is an important economic insect and a representative model organism of Lepidoptera, which has been widely used in the study of reproduction and development. The development of the silkworm's reproductive gland is easily affected by many external factors, such as chemical insecticides. After the silkworm larvae were treated with different concentrations of pyriproxyfen, the results showed that the number of eggs and hatching rate of eggs in the silkworm can be reduced by pyriproxyfen, and the concentration effects were displayed. Pyriproxyfen exposure could affect the normal development of the ovary tissue by reducing the number of oocytes and oogonia in the ovaries of silkworm fed with pyriproxyfen. We employed qRT-PCR, to detect the expressions of genes related to ovary development (Vg, Ovo, Otu, Sxl-S and Sxl-L) and hormone regulation (EcR and JHBP2) in silkworm. Our study showed that the transcription levels of Vg, Ovo, Otu, Sxl-S and Sxl-L in the treatment group were lower than those in the control group (6.08%, 61.99%, 83.51%, 99.31% and 71.95%, respectively). The transcription level of ECR was 70.22% for the control group, while that of JHBP2 was upregulated by 3.92-fold. Changes of transcription levels of these genes caused by pyriproxyfen exposure ultimately affect the absorption of nutrients, energy metabolism, ovary development and egg formation of the silkworm, thus leading to reproductive disorders of the silkworm. In general, our study revealed the response of silkworm reproduction to pyriproxyfen exposure and provided a certain reference value for the metabolism of the silkworm to pyriproxyfen.

Effects of pyriproxyfen exposure on immune signaling pathway and transcription of detoxification enzyme genes in fat body of silkworm, Bombyx mori.

Sericulture is a very important and flourishing industry in developing countries. Bombyx mori is a kind of important and well-studied economic insects in the whole world. In China, applying of pyriproxyfen pesticide often resulted in non-cocooning and silk yield reduction. However, the effects of pyriproxyfen exposure on immune signaling pathway in fat body of silkworm has not been reported yet now. In the present study, we found that the growth and development of silkworm were significantly affected by pyriproxyfen exposure and the fat body tissues were injured after treatment. It was also showed that the expressions of key genes of PI3K/Akt and CncC/Keap1 pathway can be elevated at 24-96 h after pyriproxyfen exposure. Furtherly, the relative expression levels of detoxification enzyme genes and the activities of detoxification enzymes were both increased by pyriproxyfen exposure. These results provided comprehensive view of fat body injury and gene expression changes in silkworm after pyriproxyfen exposure.

Metabolic Characterisation of the Midgut of Bombyx mori Varieties after BmNPV Infection Using GC-MS-Based Metabolite Profiling.

Bombyx mori nucleopolyhedrovirus (BmNPV) is a silkworm disease that is especially harmful to cocoon production and seriously restricts sericultural development. Our laboratory successfully cultivated a new highly BmNPV-resistant silkworm variety, Huakang 2; however, its mechanism of BmNPV resistance remains unclear. To understand its resistance mechanism, we conducted a metabolomic and transcriptomic study of the midgut of silkworm varieties, Baiyu N and Baiyu after BmNPV infection. We identified 451 differential metabolites, which were mostly comprised of small molecules, such as saccharides, acids, amines, alcohols, and glycosides. We found that the primary differences in disease resistance between the silkworm varieties are metabolic-pathways, tryptophan metabolism, oxidative phosphorylation, ABC-transporters, beta-alanine metabolism, and phenylalanine metabolism. Combined analysis with transcriptomic data suggested that tryptophan metabolism and oxidative phosphorylation are closely related to the silkworms' BmNPV resistance. We hypothesize that the roles of the two metabolic pathways in the BmNPV resistance mechanism might be the following: Oxidative phosphorylation generates a large amount of adenosine triphosphate (ATP) in response to BmNPV infection to provide silkworms the energy required for establishing BmNPV resistance. Tryptophan metabolism then activates the aryl hydrocarbon receptor (AhR) through the exogenous virus BmNPV, which activates the silkworm's immune system to defeat BmNPV infections.

Effect of Titanium Dioxide Nanoparticles on the Resistance of Silkworm to Cytoplasmic Polyhedrosis Virus in Bombyx mori.

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) is a serious disease harmful to silk industry, which is one of the major sources of financial support for farmers in many developing countries. So far, there is still no good way to prevent or treat this disease. In this study, titanium dioxide nanoparticles (TiO2 NPs) were used to pretreat silkworm larvae, and good results were achieved in improving silkworm immunity and alleviating the damage of cytoplasmic polyhedrosis virus. The results showed that nano-titanium dioxide pretreatment could inhibit the proliferation of BmCPV in the midgut of silkworm, activate JAK/STAT and PI3K-AKT immune signaling pathways, and upregulate the expression of key immune genes, so as to improve the immunity of silkworm and enhance the resistance of silkworm to BmCPV.

Gas Chromatography-Mass Spectrometry Based Midgut Metabolomics Reveals the Metabolic Perturbations under NaF Stress in Bombyx mori.

Fluoride tolerance is an important economic trait in sericulture, especially in some industrial development regions. Analyses of physiological changes involving structural damage to the insect body and molecular analyses of some related genes have focused on this area; however, the changes that occur at the metabolic level of silkworms after eating fluoride-contaminated mulberry leaves remain unclear. Here, metabonomic analysis was conducted using gas chromatography-mass spectrometry (GC-MS) to analyze the changes in midgut tissue after NaF stress using silkworm strains 733xin (susceptible stain) and T6 (strain resistant to fluoride), which were previously reported by our laboratory. Differential metabolomics analysis showed that both T6 and 733xin strains displayed complex responses after exposure to 200 mg/kg NaF. The purine metabolism and arginine and proline metabolic pathways of fluoride-tolerant strains reached significant levels, among which 3'-adenylic acid and hypoxanthine were significantly upregulated, whereas guanine, allantoic acid, xanthine, N-acetyl-L-glutamic acid, and pyruvate were significantly downregulated. These metabolic pathways may be related to the fluoride tolerance mechanism of NaF poisoning and tolerant strains.

Comparative Proteomic Analysis Reveals Immune Competence in Hemolymph of Bombyx mori Pupa Parasitized by Silkworm Maggot Exorista sorbillans.

The silkworm maggot, Exorista sorbillans, is a well-known larval endoparasitoid of the silkworm Bombyx mori that causes considerable damage to the silkworm cocoon crop. To gain insights into the response mechanism of the silkworm at the protein level, we applied a comparative proteomic approach to investigate proteomic differences in the hemolymph of the female silkworm pupae parasitized by E. sorbillans. In total, 50 differentially expressed proteins (DEPs) were successfully identified, of which 36 proteins were upregulated and 14 proteins were downregulated in response to parasitoid infection. These proteins are mainly involved in disease, energy metabolism, signaling pathways, and amino acid metabolism. Eight innate immune proteins were distinctly upregulated to resist maggot parasitism. Apoptosis-related proteins of cathepsin B and 14-3-3 zeta were significantly downregulated in E. sorbillans-parasitized silkworm pupae; their downregulation induces apoptosis. Quantitative PCR was used to further verify gene transcription of five DEPs, and the results are consistent at the transcriptional and proteomic levels. This was the first report on identification of possible proteins from the E. bombycis-parasitized silkworms at the late stage of parasitism, which contributes to furthering our understanding of the response mechanism of silkworms to parasitism and dipteran parasitoid biology.

Transcriptome-wide analysis of the difference of alternative splicing in susceptible and resistant silkworm strains after BmNPV infection.

Novel alternative splicing events were identified from BmNPV-susceptible and -resistant silkworm strains after BmNPV infection using high-throughput RNA-sequencing strategy. In total, 12.82 Gb clean RNA-seq data were generated for the two midgut samples from BmNPV-susceptible and -resistant silkworm strains, and 14.78 Gb clean data for the two fat body samples. The number of alternative splicing events and isoforms in the BmNPV-susceptible silkworm strain was more than that in the BmNPV-resistant silkworm strain. Furthermore, alternative splicing genes uniquely present in BmNPV-resistant silkworm strain were involved in functions about ribosome, whereas, alternative splicing genes uniquely present in BmNPV-susceptible silkworm strain were implicated in functions like DNA helicase activity and signal transduction. Additionally, 33 expressed SR or SR-like proteins were identified, and three genes encoding SR or SR-like proteins (tetratricopeptide repeat protein 14 homolog, ubiquitin carboxyl-terminal hydrolase 32 and zinc finger CCCH domain-containing protein 18) have a higher number of different alternative splicing events between two silkworm strains. The present study suggested BmNPV treatment may have a smaller effect on the mRNA transcription in BmNPV-resistant silkworms than that in BmNPV-susceptible silkworms, and functions of alternative splicing genes are different between the two silkworm strains.

Resistant silkworm strain block viral infection independent of melanization.

Melanization mediated by the prophenoloxidase-activating system (proPO) is an important immune response in invertebrates. However, the role of melanization on viral infection has not been wildly revealed in Bombyx mori (B. mori), silkworm. Here, we investigated the extent of melanization of susceptible (871) and resistant (near-isogenic line 871C) B. mori strains. The result showed that the extent of melanization was significantly higher in the susceptible strain than in the resistant strain. A majority of Serpins were up-regulated in the resistant strain through iTRAQ-based quantitative proteomics, comparing with susceptible strain. Our data further identified that Serpin-5, Serpin-9 and Serpin-19 reduced PO activity, indicating that the menlanization pathway was inhibited in the resistant strain. Moreover, our results indicated that the hemolymph of 871 reduced viral infection in the presence of PTU, indicating that melanization of 871 strain hemolymph blocked viral infection. However, viral infection was significantly suppressed in the hemolymph of 871C strain regardless of the presence of PTU or not, which implied that the resistant strain inhibited viral infection independent of the melanization pathway. Taken together, our findings indicated that the melanization pathway was inhibited in resistant strain. These results expend the analysis of melanization pathway in insects and provide insights into understanding the antiviral mechanism.

The evolutionary road from wild moth to domestic silkworm.

The Silk Road, which derives its name from the trade of silk produced by the domestic silkworm Bombyx mori, was an important episode in the development and interaction of human civilizations. However, the detailed history behind silkworm domestication remains ambiguous, and little is known about the underlying genetics with respect to important aspects of its domestication. Here, we reconstruct the domestication processes and identify selective sweeps by sequencing 137 representative silkworm strains. The results present an evolutionary scenario in which silkworms may have been initially domesticated in China as trimoulting lines, then subjected to independent spreads along the Silk Road that gave rise to the development of most local strains, and further improved for modern silk production in Japan and China, having descended from diverse ancestral sources. We find that genes with key roles in nitrogen and amino acid metabolism may have contributed to the promotion of silk production, and that circadian-related genes are generally selected for their adaptation. We additionally identify associations between several candidate genes and important breeding traits, thereby advancing the applicable value of our resources.

The complete mitochondrial genome of the silkworm, Bombyx mori strain BaiyuN.

Here, we report the complete mitochondrial genome of a Bombyx mori strain BaiyuN, which is identified to be highly resistant to Bombyx mori nucleopolyhedrovirus (BmNPV). Its complete mitochondrial genome is 15,655 bp in length (GenBank accession no. MG797555), consisting of 13 protein-coding genes, 22 tRNA genes, 2 rRNA, and 1 control region (494 bp). The complete mitogenome of the B. mori strain BaiyuN could provide a basic data for further phylogenetics and antivirus research.

Analysis of codon usage bias of envelope glycoprotein genes in nuclear polyhedrosis virus (NPV) and its relation to evolution.

BACKGROUND: Analysis of codon usage bias is an extremely versatile method using in furthering understanding of the genetic and evolutionary paths of species. Codon usage bias of envelope glycoprotein genes in nuclear polyhedrosis virus (NPV) has remained largely unexplored at present. Hence, the codon usage bias of NPV envelope glycoprotein was analyzed here to reveal the genetic and evolutionary relationships between different viral species in baculovirus genus. RESULTS: A total of 9236 codons from 18 different species of NPV of the baculovirus genera were used to perform this analysis. Glycoprotein of NPV exhibits weaker codon usage bias. Neutrality plot analysis and correlation analysis of effective number of codons (ENC) values indicate that natural selection is the main factor influencing codon usage bias, and that the impact of mutation pressure is relatively smaller. Another cluster analysis shows that the kinship or evolutionary relationships of these viral species can be divided into two broad categories despite all of these 18 species are from the same baculovirus genus. CONCLUSIONS: There are many elements that can affect codon bias, such as the composition of amino acids, mutation pressure, natural selection, gene expression level, and etc. In the meantime, cluster analysis also illustrates that codon usage bias of virus envelope glycoprotein can serve as an effective means of evolutionary classification in baculovirus genus.

Transcriptomic Analysis of Resistant and Susceptible Bombyx mori Strains Following BmNPV Infection Provides Insights into the Antiviral Mechanisms.

Purpose. To decipher transcriptomic changes and related genes with potential functions against Bombyx mori nucleopolyhedrovirus infection and to increase the understanding of the enhanced virus resistance of silkworm on the transcriptomic level. Methods. We assembled and annotated transcriptomes of the Qiufeng (susceptible to infection) and QiufengN (resistant to infection) strains and performed comparative analysis in order to decipher transcriptomic changes and related genes with potential functions against BmNPV infection. Results. A total of 78,408 SNPs were identified in the Qiufeng strain of silkworm and 56,786 SNPs were identified in QiufengN strain. Besides, novel AS events were found in these 2 strains. In addition, 1,728 DEGs were identified in the QiufengN strain compared with Qiufeng strain. These DEGs were involved in GO terms related to membrane, metabolism, binding and catalytic activity, cellular processes, and organismal systems. The highest levels of gene representation were found in oxidative phosphorylation, phagosome, TCA cycle, arginine and proline metabolism, and pyruvate metabolism. Additionally, COG analysis indicated that DEGs were involved in "amino acid transport and metabolism" and "carbohydrate transport and metabolism." Conclusion. We identified a series of major pathological changes in silkworm following infection and several functions were related to the antiviral mechanisms of silkworm.