There are tens of thousands of RNAs in all organisms. Because of the existence of transcriptomics, a large number of RNAs have been found continuously, and this number constantly rises, no matter in which species (Mallory and Shkumatava, 2015). These RNAs play various vital functions in organisms, including immunity, growth and development, behavioral influence, etc. (Mallory and Shkumatava, 2015). Many differences have also been found in the silkworm transcriptome, which are also related to viral infection and resistance (Wang et al., 2015). Among them, the distinction between Dazao-infected and uninfected brains showed 742 differentially expressed genes (DEGs), including 218 upregulated genes and 524 downregulated genes (Wang et al., 2015). Among these differentially expressed genes, 2 genes were significantly related to circadian regulation, 9 genes were related to synaptic transmission, and 3 genes were tightly involved in serum signaling pathways, all of which are speculated to be closely associated with BmNPV infection in Bombyx mori (Wang et al., 2015; Abbas et al., 2020). Because irregular movements occur in silkworms in the late stages of infection, these irregular movements might amplify the spread of the polyhedrosis virus among silkworms again (Wang et al., 2015). The differential expression of these brain genes might be due to the impact of the virus, thus indirectly causing their behavioral disorders.

BLAST was used to find differentially expressed genes from previous transcriptomic analysis, BmβGRP4, which was found to have an essential inhibitory effect on the replication of BmNPV (Wang J. et al., 2021). RT-qPCR results between different tissues showed that BmβGRP4 was highly expressed in essential tissues such as the head, skin, midgut, hemolymph, and fat body. At the same time, larvae also showed higher expression compared to adult insects (Wang J. et al., 2021). The infection experiments of P50 showed that BmβGRP4 had a significant down-regulated expression when BmNPV infect and replicate in vivo, while the RNAi of BmβGRP4 showed a dramatic increase of viral VP39 (The main nucleocapsid protein of BmNPV)at 24–72 hpi (Wang J. et al., 2021). However, overexpression of BmβGRP4 in BmN cells showed the opposite result, with a significant reduction in the viral replication (Wang J. et al., 2021). The main reason was that BmβGRP4 could induce upregulation of several apoptosis-related genes, including BmApaf1, BmDredd, BmCaspaseNC, BmICE, and BmCaspase1 to promote virus-induced apoptosis (Wang J. et al., 2021). Previous studies have identified various genes, mainly in the midgut or midgut juice of the silkworm. Still, very few studies have examined the reason for these differences between different tissues. However, studies have shown that the virus replicates in the midgut through the periplasm into the silkworm (Jiang and Xia, 2014).

Nevertheless, numerous studies have concluded that the infection of BmNPV is systemic and affects the whole host body to different degrees, whether it is the brain, midgut, or hemolymph. Therefore, experiments on the entire body needed to be extensively studied to explore the correlation between the differences. In 2006, it was shown that a total of 78, 408 SNPs were identified from the fat body of BmNPV-infected Qiufeng (susceptible strain) compared with the infected Qiufeng N (resistant strain) (Li et al., 2016). In contrast, 56, 786 SNPs were identified in Qiufeng N. Significant nucleic acid sequence differences existed between the two strains (Li et al., 2016). Despite the sequence differences, there were 1,728 DEGs, most of which showed down-regulation in Qiufeng N (Li et al., 2016). Most of the GO enrichment analysis of these DEGs focused on the membrane, metabolism, binding and catalytic activity, cellular processes, and organismal systems. In contrast, KEGG concentrated on oxidative phosphorylation, phagosome, and TCA cycle (Li et al., 2016). These results were similar to an earlier analysis by Xue et al. (2012). Although Xu et al. experimented in the midgut of the silkworm, both experiments involved multiple metabolic pathways and oxidative phosphorylation (Xue et al., 2012; Li et al., 2016). While another analysis of different strains, YeA (resistant strain) and YeB (susceptible strain), showed that after 48 h, 7 potential DEGs between the two groups contained BmAtlatin-n, BmFerHCH, BmTHY, Bmseroin, Bmseroin2, BmNHR96, and BmSINAL10. qPCR verified the expression of these DEGs in the midgut and revealed that these genes were differentially expressed between the two after 48 h (Wang X. Y. et al., 2021). All were thought to be potentially involved in the resistant immune response of the silkworm (Wang X. Y. et al., 2021). Although there are commonalities between these studies, further studies are still needed to determine the function of these individual genes.

Resistant genes in domestic silkworms are mainly paternally biased and dominantly inherited. Therefore, resistant genes may have a strong association with sex in silkworms. In a comparative analysis of the highly resistant strain NB and the highly susceptible strain 306 of different sexes, He et al. found that for different sexes of different strains, the number of DEGs either upregulated or downregulated was larger in male NB_V than in female NB_V after infection, while in 306_V, the trend showed the opposite. Eighteen DEGs were uniquely expressed only in NB_V which included 14 upregulated DEGs. The result is consistent with the current view that dominant genes mainly control resistance to BmNPV. Using classical genetics experiments on Cydia pomonella’s baculovirus resistance experiment, Asser-Kaiser proved that a dominant gene controls the resistance gene on the Z chromosome. After 10 years of research, Professor Chen’s team at Jiangsu University found that the BmNPV-resistant NB strain turned out to be a heterozygous population containing three genotypes (++, +−, −). Still, after more than 20 generations of isolation and purification, the laboratory finally obtained a completely homozygous antiviral strain NB (++). Combined with genetic experiments, the resistance of BmNPV in NB might be controlled by a single dominant gene (Hu et al., 2023). These results suggest that a single dominant gene may also control resistance genes in silkworms.

The transcriptomics of different strains and tissues, and even between different sexes, have demonstrated that a large number of DEGs produced in silkworms after infection were inextricably linked to the virus, and these DEGs are concentrated in immune, membrane function and metabolic pathways among different strains in different tissues. Although the large number of candidate DEGs screened by transcriptome may play an important role, there needs to be a more precise explanation for the exact differences in these DEGs. Further detailed experiments are still needed to verify them.

After being infected by BmNPV, the number and expression of proteins in silkworms will significantly change. The protein expression of different silkworms was also different (Liu et al., 2010). Previous studies have also identified multiple resistance-related proteins in the larval midgut (Ponnuvel et al., 2003; Pandian et al., 2008). Proteins might have an essential role in anti-BmNPV. Ten years earlier, a common proteomics analysis of NB ♀, ♂ and 306 ♀, ♂, their direct cross-group (NB ♀, 306 ♂, F1 hybrid) and reciprocal cross-group (306 ♀, NB ♂, F1 hybrid), revealed that 53 and 21 uniquely expressed differential proteins in the two F1 offspring, respectively (Qin et al., 2012). These results also suggest that resistance genes may have a genetic effect and that breeding resistant strains can focus on male hybrids. Because the transcriptome results are the same, the inheritance of resistance genes was more pronounced in male-resistant strains NB (He et al., 2021). Comparative proteome results of both antagonistic and perceptual strains were also shown that only caspase-1 and serine protease were expressed in resistant silkworms but not in the susceptible strain or reciprocal cross group, a result confirmed by WB (Qin et al., 2012). Label-free quantitative proteomic analysis of midgut fluid from infected A35 and P50, the candidate BmTA (trypsin, alkaline A) was significantly upregulated in A35 after infection (Zhang S.-Z. et al., 2020). Studies found that a specific concentration of recombinant BmTA was associated with the virus co-incubation and had an inhibitory effect on the subsequent replication of ODVs and BVs. The pIZT/V5-His-mCherry vector was used to overexpress BmTA in BmN cells. Compared with the control group, the viral VP39 was significantly lower than that of the control group (Zhang S.-Z. et al., 2020).

BC9, a resistant near-isogenic silkworm line generated from resistant line NB and susceptible line 306, the proteomics of its fat bodies after infection showed that the proteome differed after infection. Two of the final selected four proteins involved in energy metabolism, their sequences shared high similarity with phosphoglycerate kinase (PGK) and arginine kinase (AK), respectively (Liu et al., 2010). The remaining protein shared sequence similarity with Pol polyprotein and a novel unidentified protein (Liu et al., 2010). DEPs were mainly involved in energy metabolism since replication of BmNPV requires membrane fusion and DNA replication and translation, all of which require large amounts of energy (Xue et al., 2012). DEPs associated with energy metabolism might then be used by viral replication. Comparative subcellular protein analysis based on the midgut of BC9 and P50 showed that 87 DEPs were screened from them by excluding genetic background (BC9 – vs. P50 –) and individual differences (P50 + vs. P50 –) (Yao et al., 2003; Wang X.-Y. et al., 2017). Energy and protein metabolism, signaling pathways, disease, and transport were predominantly enriched. In particular, protein changes in microsomes were obvious (Wang X.-Y. et al., 2017). The iTRAQ quantitative proteomic showed that between the BC9 and P50, 793 DEPs were found (Yu et al., 2017). After excluding the genetic background and individual differences, 84 total DEPs, 15 cytoskeletons, 12 apoptosis, 8 ubiquitinated, 14 immune, 6 cell signaling regulation, 5 endocytosis-acting, 11 translational, and 14 endopeptidase-associated proteins were selected, all of which are thought to be important in the silkworm (Yu et al., 2017). The fact is that during the infected phase, energy-related, membrane-associated, and molecular functions are significantly different in the host, suggesting the presence of important biological functions.

In recent years, studies of the heat shock protein family have shown the exact resistant function to the virus, but it also played a positive role in the viral replication (Wu et al., 2019b). Wu et al. (2019b) demonstrated that stable replication of BmNPV was affected by inhibition of the HSP90 expression (Shang et al., 2020). Of the 195 DEPs screened in BmN cells inhibited by GA (HSP90 inhibitor), 136 were upregulated DEPs. Expression of immune-related proteins, cellular DNA repair-related proteins, and zinc finger proteins was significantly increased, while kinases were significantly decreased (Wu et al., 2019b). The PPI network showed that 4 heat shock proteins (Hsp12, Hsp20, Hsp23.7, and Hsp70), 5 immune-related proteins, 7 protein kinases, and 10 zinc finger proteins were inextricably linked to HSP90 (Wu et al., 2019b). Further, HSP90 experiments showed its exact positive function to BmNPV (Katsuma, 2021). Other HSP90-inhibited experiments by 17-AAG in BmN cells detected a significant reduction of the viral particles. In contrast, the opposite trend was shown in the control group (treated with a low dose of inhibitor), demonstrating the apparent role of HSP90 in promoting the BmNPV replication (Katsuma, 2021). WB experiments were performed to investigate the specific effect of HSP90 on viral replication. Although HSP90 had no significant effect on the very early gene, it did have significant effects on DBP, BRO, and POLH (early, delayed early, and very late stages of BmNPV, respectively) (Katsuma, 2021).

The metabolic environment of the silkworm, such as the midgut and fat body, changed dramatically after infection. Studies have shown a dynamic equilibrium between the gut and microbial communities in which microorganisms regulate different biological processes (Ha et al., 2009; Bae et al., 2010; Lee and Hase, 2014). These changes were inextricably linked to the infection of BmNPV. LC–MS/MS analysis indicated that in BmE cells at different infection stages (3 ~ 72 hpi), metabolite pathways were concentrated in ABC transporter, amino acid transport RNA, and purine metabolism signaling pathways (Hua et al., 2021). Cluster analysis was gathered in 5-pyridoxolactone, GPC, 2-OH-Ade, γ-Glu-Cys, and hydroxybutyramide, which all showed a constant upregulation in these periods (Hua et al., 2021). At the same time, only 5-pyridoxolactone strongly inhibited BmNPV replication in the phenotype experiments (Hua et al., 2021). It is shown that metabolites may have a complex resistance mechanism to BmNPV.

Comparing the midgut metabolites of infected and uninfected (Jingsong X Haoyue) silkworms, there were 13 differentially upregulated metabolites, including aspartic acid, asparagine, cysteine, isoleucine, lysine, methionine, phenylalanine, proline, threonine, tyrosine, valine, leucine, and serine (Shi et al., 2021). The authors pointed out that, except for cysteine, isoleucine, proline, and serine, the other nine amino acids were essential for the growth and development of the silkworm (Shi et al., 2021).

The hemolymph metabolites of the YeA (a resistant strain) and YeB (a susceptible strain) found that the expression of trehalose in the YeA continued to rise after infection. Also, the metabolites related to TCA metabolism were actively expressed, and branched-chain amino acids (BCAA) produced irregular fluctuations on the timeline (Wang X. Y. et al., 2019). In another experiment, the microenvironment of the hemolymph, midgut, fat body, and other organs was drastically altered after infection. A total of 296, 215, and 108 DEMs were screened by KEGG and HMDB annotation, and all these changes were thought to be inextricably linked to the BmNPV infection (Wang G. et al., 2021). The above three tissues’ metabolites enrichment results mainly focused on the metabolism of chitin, pyruvate, and glyoxylate, multiple metabolisms of carbohydrates and fatty acids, also with a variety of amino acids such as tryptophan, tyrosine, and phenylalanine (Wang G. et al., 2021). Only tryptophan has been shown to aid in the transformation of indole-3-aldehyde (IAId), indole-3-acetic acid (IAA), and other ligands of the silkworm aryl hydrocarbon receptor (AhR) in the intestines, which can activate a variety of immune cells to be beneficial for immune homeostasis, thereby producing resistance to BmNPV (Gao et al., 2018; Qian et al., 2020). Through metabolomic studies, it was ascertained that metabolites played an important role in silkworm disease resistance; however, problems still need to be solved, including the multiple pathways involved in metabolites and the proportion of key and secondary metabolites. These problems demonstrate the difficulty of studying metabolites in the silkworm; thus, further experiments on disease-resistant mechanisms are required.

BmNPV-resistant genes are present in domestic silkworms and have different roles in different sexes and strains before and after infection (Lu et al., 2018). These differences found by omics such as DEGs found in the transcriptomics involved multiple regulatory pathways, mainly focused on metabolism and membrane transport (He et al., 2021; Katsuma, 2021; Wang X. Y. et al., 2021). Numerous experiments have also shown that these DEGs had resistance to the virus, or among the HSP family screened by proteomics, HSP90 had a facilitative effect on viral replication, while HSP is also regulated by a variety of ncRNAs (Katsuma, 2021; Zhao Z.-M. et al., 2022), other proteins such as serine proteins, VTPases involved in energy metabolism were all shown to differ in the domestic silkworm (Qin et al., 2012; Zhang S.-Z. et al., 2020), suggested a complex mechanism for the role of proteins; moreover, the strong inhibitory effect of 5-pyridolactone on BmNPV (Hua et al., 2021). Also, tryptophan and 9 other amino acids were all proved to be essential to the growth and development of the silkworm and were discovered through metabolomics (Shi et al., 2021).

As two classes of short-stranded regulatory noncoding RNAs, also known as miRNAs, the study of piRNA and siRNA in BmNPV resistance is equally important. piRNA can use transposon cleavage to influence the expression of target genes, targeting the degradation of target transcripts through PIWI proteins interacting with it. In addition, piRNA can also cope with viral infection (Nie et al., 2013; Xia et al., 2022). Feng et al. identified a large number of piRNAs from the location transposon region in the fat bodies and midgut of the DaZao strain infected with BmNPV. These piRNAs had significant differential expression before and after BmNPV infection; compared with normal tissues, 322 piRNAs were significantly upregulated, and 276 were significantly downregulated in the fat bodies after BmNPV infection. In the midgut, 129 piRNAs were significantly upregulated, and 117 were downregulated considerably, indicating that piRNA has a specific role in silkworms (Feng et al., 2021). In the latest study, the replication of AcMNPV in lepidopteran insect cells sf9 requires the expression of piRNA. Xia et al. found that PIWI-like-1 and PIWI-like-2-3 can inhibit the replication of AcMNPV. At the same time, PIWI-like-4-5 can promote the replication of the virus (Xia et al., 2022).

siRNA is a class of double-strand noncoding RNAs that mainly function as foreign RNA and are, therefore, active in studying viral genes (Zhao et al., 2021). Long-term evolution of viruses has also resulted in the development of RNAi that can evade host immune responses. This includes siRNA, a type of exogenous dsRNA, which is first cleaved and bound to Dcr-2 and to form an RNA-induced silencing complex (RISC) with Ago2. Then the siRNA is cut off a strand, and the remaining vsiRNA plays a corresponding RNA degradation (Zhao et al., 2021). Dcr-2 expression was elevated in silkworms infected with BmNPV, suggesting that siRNA may play a role in BmNPV (Li et al., 2020). In the 2022 study, BmNPV p35 was shown to promote the replication of the virus by inhibiting host siRNA interference pathways (Zhao S. et al., 2022). By comparing BmNPV infection in RNAi-deficient cells, it was found that BmNPV did not show a difference in these cells, demonstrating the importance of siRNA for viral infection (Zhao S. et al., 2022). Further results showed that pinp35 also inhibited the production of vsiRNA and could disrupt the binding of siRNA to Ago2 protein, resulting in the failure of the host RNAi pathway and promotion of the BmNPV replication (Zhao S. et al., 2022).

Tong et al. (2006) aligned silkworm miRNAs by using mature miRNAs identified from model animals as query sequences, showing that there were 24 miRNAs in the silkworm genome. Wu P. et al. (2016) constructed a miRNA library from the silkworm P50 and applied Solexa sequencing technology to identify 38 differentially expressed miRNAs between infected and uninfected groups (Wu P. et al., 2016). Among the selected differentially expressed miRNAs, miR-282, miR-10-5p, and miR-277-5p were included in the series of genes expressed in cells after BmNPV infection (Wu P. et al., 2016). For the individual miRNA experiment, upregulation of miRNA bmo-miR-3378 could inhibit the replication of the virus, while the inhibited investigation showed the opposite trend (Jiang, 2019). But the validation of downstream targets is still lacking and needs to be improved because of their non-coding features (Jiang, 2019), resulting in unclear mechanisms. As a non-coding RNA, it remains to be clarified how it causes difficulties in viral replication (Table 2).

Using high-throughput sequencing technology, two miRNAs named bmo-miRNA-390 and bmo-miR-217 were predicted from different Bombyx mori. The prediction proved that the above two miRNAs could interact with cg30 and lef-1 proteins, respectively (Wang M. L. et al., 2017; Kang et al., 2018). Both of these BmNPV proteins were indispensable in the replication of the BmNPV (Wang M. L. et al., 2017; Kang et al., 2018). By constructing luciferase reporter plasmids containing genes encoding cg30 and lef-1 proteins, co-transfection with miRNA mimics into cells to investigate the binding of miRNAs and viral proteins. The fluorescence of the experimental group was weaker than the control group, indicating the interaction between miRNA and viral proteins (Wang M. L. et al., 2017; Kang et al., 2018). Meanwhile, in BmNPV-infected BmN cells, transfection of mimics or inhibitors of both miRNAs revealed that bmo-miRNA-390 and bmo-miR-217 could effectively down-regulate the expression of BmNPV-cg30 and BmNPV-lef-1 in infected BmN cells, respectively (Wang M. L. et al., 2017; Kang et al., 2018). It affected the infection and replication of BmNPV in Bombyx mori (Wang M. L. et al., 2017; Kang et al., 2018). Later, Wu et al. (2019a) identified that Bombyx mori miRNA, bmo-miR-2819 was the target of BmNPV-ie-1. BmNPV could reduce the expression of bmo-miR-2819 to enhance the expression of viral ie-1, thus promoting the viral replication (Wu et al., 2019a). Dual luciferase experiments confirmed their interactions between bmo-miR-2819 and ie-1, additional experiments involving overexpression and inhibition of the bmo-miR-2819 demonstrated that the virus could reduce the expression of bmo-miR-2819 to enhance viral ie-1 expression, thus promoting viral replication (Wu et al., 2019a).

In addition to the host miRNAs, viral miRNAs are also significant. Viral miRNAs could regulate host miRNAs and indirectly regulate the expression of the host’s genes or proteins, such as host Ran and Bmo-miR-5738 (Singh et al., 2012; Cao et al., 2017). Thereby realizing the effect of enhancing the self-replication (Tang et al., 2019). In 2004, miRNAs were first discovered in the Epstein–Barr virus (a herpes virus) (Pfeffer et al., 2004). Subsequently, hundreds of viral miRNAs were identified, including baculovirus and adenovirus. In 2009, Shirina successfully predicted the existence of miRNAs in BmNPV, such as bmoNPV-miR-1ph and bmoNPV-miR-2ph, using RNAfold (Shirina et al., 2011). Singh also identified the viral miRNAs, like bmnpv-miR-1 and bmnpv-miR-3 (Singh et al., 2010), and predicted the targets of these miRNAs, which were all considered to be significantly involved in the insect-pathogen interactions by regulating several replicating genes (Singh et al., 2010). It was observed that bmnpv-miR-1 could regulate the miRNAs in host cells, thereby regulating the host-virus defense mechanism affected by the host miRNAs (Singh et al., 2012). Though there is no direct evidence that viral miRNAs can directly regulate host miRNAs, in 2012, the mechanism of bmnpv-miR-1 was proved to affect Ran expression (an essential component of exportin-5 protein) and, in turn, downregulated the expression of bmnpv-miR-1 which resulted in the reduction of BmNPV replication (Singh et al., 2012). Bmnpv-miR-3 is also a BmNPV-produced miRNA that could regulate the P6.9 protein and various viral proteins to escape the immune attack of host cells (Singh et al., 2014). In addition to the direct regulation of multiple proteins, viral miRNAs can also indirectly regulate host miRNAs that can, in turn, regulate the protein expression (Cao et al., 2017). BmNPV-miR-415 was also found that it could down-regulate the expression of Bmo-miR-5738 in BmN cells and indirectly affect the expression of TOR2 protein, which plays an important role in the growth and development of the host (Cao et al., 2017). Although the reasons for how viruses downregulate host miRNAs have not been elaborated, viral miRNAs could regulate host miRNA transport proteins (Singh et al., 2012). BmNPV-miR-415 may likely have a similar function, regulating host proteins to affect the expression of host miRNAs and miRNA-regulated proteins (Cao et al., 2017). miRNAs have enriched the content and depth of BmNPV-resistant research. Still, the specific role of viral miRNAs for BmNPV resistance in silkworms requires further experiments to clarify their requirement for virus replication or their role in the regulation of host immune responses and pathways.

LncRNA is widely found in animals, plants, viruses, and fungi, where they play important roles in epigenetic, transcriptional, and post-transcriptional regulation (Rinn et al., 2007). As a non-coding RNA, its function is as simple as host immunity and viral replication (Fatica and Bozzoni, 2014; Li et al., 2014; Ishihara et al., 2015; Sun and Kraus, 2015; Adelman and Egan, 2017; Zhang et al., 2021). In 2014, Ishihara identified many antisense lncRNAs in the BmNPV genome. There were 15 lncRNAs identified from the sequences of the early and late promoters in BmNPV, and the mutants of these genes were constructed and infected the silkworm (Ishihara et al., 2015). Compared to the ODVs produced by wild types KO122AS, the ODVs produced by the mutant KO122AS, encoded by the Bm122 genome, were significantly reduced, showing that the BmNPV lncRNA might have essential roles in promoting viral replication and spread by regulating its gene expression (Ishihara et al., 2015).

Except for the viral lncRNA in the silkworm, Zhang et al. used whole transcriptomic sequencing technology to analyze 4,450 differentially expressed lncRNAs, 66 DEmiRNAs, and 7,448 DEmRNAs in the BmNPV infected and uninfected BmN cells to construct a lncRNA-miRNA-mRNA network. It was then found that these networks focused on the protein hydrolysis and lysosomal pathway, all of which were significantly associated with the immune pathways in silkworms (Zhang S. et al., 2020). Furthermore, the results indicated that after infection with BmNPV, the silkworm could regulate its self-immunity through the lncRNA-miRNA-mRNA pathway (Zhang S. et al., 2020). In the fat body, Lnc_209997 was significantly down-regulated in BmNPV-infected silkworms, and Lnc_209997 could also inhibit BmNPV replication by regulating miR-275-5p expression (Lin et al., 2022). Therefore, to promote its replication in vivo, BmNPV reduced the expression of host Lnc_209997 and indirectly regulated miR-275-5p (Lin et al., 2022). Since these studies, the function of lncRNA is becoming better understood. As ncRNAs, the specific lncRNA-miRNA-mRNA regulatory network suggests that lncRNA may act indirectly through miRNAs to regulate mRNA expression regardless of the non-coding feature of itself (Ahadi, 2021). The subsequent studies of their functions should focus on discovering their target genes, as this can help elucidate the specific functions of lncRNA, thereby elaborating the capability of lncRNA and their regulation of vital silkworm activities.

Yin et al. (2020) constructed a DElncRNA-bmo-miR-278-3p-BmHSC70 regulatory network in GA-treated BmN cells and found that the lncRNA and miR-278-3p were responsible for the regulation of HSP70 expression (Yin et al., 2020). HSP70, on the other hand, was found to have possible protein interactions with host HSP90 and may be influenced and regulated by BmNPV (Wu et al., 2019b). Zhao used the exogenous inhibitors of HSP90, geldanamycin (GA), which could inhibit the expression of HSP90, and finally inhibit the replication and proliferation of the virus (Zhao Z.-M. et al., 2022), which was similar to the previous studies carried out by Wu et al. (2019b). However, Zhao detected 29 DEmiRNAs in the cells after GA treatment and constructed a DElncRNA-DEmiRNA-DEmRNA regulatory network (Zhao Z.-M. et al., 2022). This network included 169 DElncRNAs, 12 DEmiRNAs, and 69 DEmRNAs. GO and KEGG analysis of these DEmRNAs showed that these DEmRNAs are mainly involved in ubiquitin-mediated proteolytic, phagosomal, proteasomal, and endocytic pathways (Zhao Z.-M. et al., 2022). These functions are all related to membrane transport during virus replication in the silkworm, suggesting that in addition to proteins, ncRNAs may play an important role in the silkworm immune system following the inhibition of HSP90 (Zhao Z.-M. et al., 2022).

CircRNA was first identified by using electron microscopy in 1976 but did not attract much attention then. Only in recent years, with the development of the high-throughput sequencing technique and the discovery of a series of circRNAs, researchers have gradually realized the specificity of circRNA and its functional importance in vivo (Sanger et al., 1976; Jeck et al., 2013; Memczak et al., 2013; Abbas et al., 2023). With the development of technology, the functions of circRNA have been found not only to act as a miRNA sponge that regulates the transcriptional level but also to encode proteins and regulate protein binding at different levels (Qu et al., 2015; Kristensen et al., 2019).

The circRNAs first found in silkworms were in 2017. Researchers discovered specifically expressed circRNAs in the silk glands, and a circRNA-miRNA network was constructed (Gan et al., 2017). Thirteen miRNAs interacted with 193 circRNAs were then identified in BmCPV-infected silkworms by Hu et al. In contrast, in 2018, Hu et al. also identified 353 DEcircRNAs in BmNPV-infected silkworms. GO and KEGG enrichment analysis for the target sites of miRNAs regulated by circRNAs showed ubiquitination, apoptosis, and endocytosis enriched in the silkworms (Gan et al., 2017; Zhang et al., 2022). The signaling pathways of the BmNPV-infected silkworms were determined, and these pathways could regulate the Bombyx mori immune pathways (Gan et al., 2017; Zhang et al., 2022). It was the same as the lncRNA, circRNA-miRNA-mRNA network also found in the silkworm fat body by Yin et al. (2021) and after performing GO and KEGG analysis, it was shown that the fat body is significantly involved in the immune system. These results demonstrated the important role of circRNA in viral infection for the first time.

Although current circRNA research focuses mainly on circRNA as a miRNA sponge to regulate mRNA expression, some studies also demonstrated that circRNA regulated histone modification. Wang et al. found that the cyclized BmEgg gene in B.mori (gene cyclized by the eggless B. mori histone-lysine N-methyltransferase (BmEgg) gene), circEgg positively regulates histone deacetylase (HDAC) expression in B. mori via the sponge with bmo-miR-3391 to inhibit the methylation of lysine 9 of histone H3 and promote acetylation of lysine 9 of histone H3 (Wang et al., 2020). Some sites were essential for the replication of the virus, like the lysine acetylation (Kac) sites, which were found in the late expression factor −3/−4/−6/−11 (Lef-3/−4/−6/−11). In previous studies, 4 Kac sites were required for late viral replication, suggesting the protein acetylation might have unsure functions in the viral replication (Rapp et al., 1998; Mikhailov et al., 2006; Yu and Carstens, 2012; Hu D. et al., 2018). Then in one Kac site, Lef-3, two acetylation modification sites, including K18 and K27, were located in the conserved Lef-3 region. And the K18 acetylation site would cause the binding of Lef-3 and P143 (a DNA helicase protein) by the experiments, making it difficult to replicate the virus in BmN cells (Ito et al., 2004; Gao et al., 2021). Indirectly controlling the protein acetylation modification in infected silkworms through the regulation of miRNA-protein modifying enzymes caused the complex replication of BmNPV in non-resistant silkworms (Wang et al., 2020). These findings suggested that silkworms could be regulated by overexpressing or repressing circRNA to realize the resistance. We believe that by modifying the protein, noncoding RNAs may have a lower impact on the silkworm phenotype than other methods, such as directly knocking out some key genes, because their low expression, non-coding, and target genes are predictable, giving them the unique advantage of having little effect on other none-resistant genes or proteins. All related genes or substances are classified and shown in Figure 3.