Ping-Zhen Xu
Mei-Rong Zhang1,2
Xue-Yang Wang
Yang-Chun Wu- 1Jiangsu Key Laboratory of Sericultural Biology and Biotechnology, School of Biotechnology, Jiangsu University of Science and Technology, Zhenjiang, China
- 2Key Laboratory of Silkworm and Mulberry Genetic Improvement, Ministry of Agriculture, Sericultural Research Institute, Chinese Academy of Agricultural Sciences, Zhenjiang, China
The mulberry silkworm (Bombyx mori) is a model organism, and BmNPV is a typical baculovirus. Together, these organisms form a useful model to investigate host–baculovirus interactions. Prothoracic glands (PGs) are also model organs, used to investigate the regulatory effect of synthetic ecdysone on insect growth and development. In this study, day-4 fifth instar silkworm larvae were infected with BmNPV. Wandering silkworms appeared in the infected groups 12 h earlier than in the control groups, and the ecdysone titer in infected larvae was significantly higher than that of the control larvae. We then used RNA sequencing (RNA-seq) to analyze silkworm PGs 48 h after BmNPV infection. We identified 15 differentially expressed genes (DEGs) that were classified as mainly being involved in metabolic processes and pathways. All 15 DEGs were expressed in the PGs, of which Novel01674, BmJing, and BmAryl were specifically expressed in the PGs. The transcripts of BmNGDN, BmTrypsin-1, BmACSS3, and BmJing were significantly increased, and BmPyd3, BmTitin, BmIGc2, Novel01674, and BmAryl were significantly decreased from 24 to 72 h in the PGs after BmNPV infection. The changes in the transcription of these nine genes were generally consistent with the transcriptome data. The upregulation of BmTrypsin-1 and BmACSS3 indicate that these DEGs may be involved in the maturation process in the latter half of the fifth instar of silkworm larvae. These findings further our understanding of silkworm larval development, the interaction between BmNPV infection and the host developmental response, and host–baculovirus interactions in general.
Introduction
The mulberry silkworm (Bombyx mori) has been reared for the past 5,000 years in China due to its importance for silk production. In addition to this economic importance, B. mori has recently played an essential role as a model organism in scientific research, molecular biology, and genetics studies (Mita et al., 2004; Xia et al., 2004). B. mori undergoes complete (egg–larva–pupa–adult) metamorphosis within each generation; however, only the larval stage feeds. In general, silkworm larvae are tetramolters that proceed through four instars, molting between each instar. The durations of the larval instar stages are as follows: 3–4 days in the first instar, 2–3 days in the second instar, 3–4 days in the third instar, 5–6 days in the fourth instar, and 6–8 days in the fifth instar. However, the duration of the larval stage depends on the silkworm strain and rearing temperature. Silkworm larvae grow rapidly, and the weight of a terminal fifth instar larva is ∼10,000 times that of a newly hatched larva (Xu et al., 2019). In particular, day 3 of the fifth instar larvae is the boundary for the larval stage, where the larvae feed and grow quickly from day 1 to day 3 of the fifth instar. The larvae then (in the gluttonous stage) greatly synthesize silk proteins in the silk gland (Xu et al., 2019), which indicates maturation and leads to spinning in the terminal fifth instar stage. After the completion of silk spinning, the silkworms proceed with larval–pupal metamorphosis.
Sericulture is one of the main sources of income for farmers in many developing countries, such as China, India, Brazil, Vietnam, and Thailand (Jiang and Xia, 2014). China produces almost 80% of cocoons worldwide. Sericulture faces biological challenges from pathogenic fungi, bacteria, and viruses, which can cause annual cocoon production losses of 20–30% (Jiang et al., 2013a). Although antibiotics are administered to silkworms to prevent and treat bacterial diseases, and fresh lime and chlorine-containing preparations are used to disinfect the rearing seat to prevent fungal diseases, there are no effective prevention and treatment methods for viral diseases. Viral diseases are responsible for almost 80% of the annual cocoon production losses, and B. mori nucleopolyhedrovirus (BmNPV) is one of the major pathogens and the most prevalent threat to sericulture in almost all countries worldwide (Jiang et al., 2013a; Jiang and Xia, 2014). BmNPV is an enveloped double-stranded DNA virus that presents a biphasic infection process throughout its viral life cycle, generating progeny with two different phenotypes, namely, occlusion-derived virus (ODV) and budded virus (BV) (Gomi et al., 1999). ODVs are packaged in occlusion bodies (OBs). Both forms play different roles during pathogenesis. The alkalinity of the silkworm midgut triggers the dissolution of OBs and the release of ODV in the midgut lumen. The ODV is responsible for the primary infection through oral transmission of the virus among silkworm larvae, while the BV is responsible for the secondary infection, causing systemic spreading all over the host within the infected silkworm larvae (Jiang, 2021).
Bombyx mori is a model organism, and BmNPV is a typical baculovirus (Gomi et al., 1999; Mita et al., 2004; Xia et al., 2004), and together, they present an important model to assess host–baculovirus interactions (Jiang and Xia, 2014). Insights from previous host studies revealed that innate antiviral immunity in lepidopteran insects plays important roles in host–baculovirus interactions (Jiang, 2021). Antiviral proteins, including red fluorescent proteins (RFPs) (Sunagar et al., 2011; Manjunatha et al., 2018), Bmlipase (Bmlipase-1 and Bmlipase member H-A) (Ponnuvel et al., 2003; Zhang S. Z. et al., 2020), serine proteases (SPs), and serine protease homologs (SPHs) (Nakazawa et al., 2004; Ponnuvel et al., 2012), show strong antiviral activity in the digestive juice of the silkworm. Moreover, heat shock protein 19.9 (Bmhsp19.9) is involved in antiviral immunity against BmNPV function (Jiang et al., 2021b). BmNPV has also evolved diverse mechanisms to counter host responses and ensure its replication. For example, BmNPV activates the expression of BmPGRP2-2 to inhibit phosphatase and tensin homolog (PTEN), which relieves its suppression of the PI3K-Akt pathway and triggers an increase in Akt phosphorylation (p-Akt) to inhibit cell apoptosis; the resulting increased cell survival is beneficial for viral replication (Jiang et al., 2019). BmSpry is upstream of ERK and JNK and is downregulated by BmNPV to elevate p-ERK and ensure viral reproduction in the silkworm (Guo et al., 2019). BmNPV activates the host ERK and JNK signal pathways for efficient replication (Katsuma et al., 2007). The baculovirus ecdysteroid UDP-glucosyltransferase gene (egt) encoding the enzyme ecdysteroid UDP-glucosyltransferase catalyzes the transfer of glucose from UDP-glucose to ecdysteroid molting hormones, and the expression of this enzyme blocks the molting of infected larval insects (O’Reilly and Miller, 1989). The BmNPV egt gene prolongs the survival time of infected silkworms to increase virus reproduction (Katsuma and Shimada, 2015).
The expression of NPV genes occurs in four phases: immediate early phase (0–4 h post-infection, hpi), delayed early stage (5–7 hpi), late stage (8–18 hpi), and very late stage (>18 hpi). Viral DNA replication starts at 8 hpi and represents the transition from the early stage to the late stage (Huh and Weaver, 1990; Jiang et al., 2013b, 2021a). Global shutoff of host gene expression and protein synthesis in insect cells begins at the early stage at around 12–18 h after NPV infection (Du and Thiem, 1997; Shirata et al., 2010; Ikeda et al., 2013). However, previous studies that investigated the interactions between BmNPV and its hosts have mainly focused on newly exuviated fifth instar silkworm larvae infected by BmNPV and the systemic process of infection by BmNPV within 48 hpi (i.e., silkworm larvae in the first half of the fifth instar). Until now, no studies have investigated the interactions between BmNPV and silkworm larvae in the latter half of the fifth instar.
The use of next-generation sequencing technologies in genome-wide studies of silkworms and BmNPV interactions is a recent development and is rapidly advancing. Recently, several studies have reported on the transcriptional response of silkworm larvae against BmNPV infection in the major innate immune tissues of the fat body and midgut (Chen et al., 2019; Huang et al., 2019; Jiang et al., 2019; Toufeeq et al., 2019; Zhang X. et al., 2020). However, the gene expression of prothoracic glands (PGs) infected by BmNPV has not yet been analyzed.
In the present study, we first investigated the precocious molting and metamorphosis of silkworm larvae under BmNPV infection, and the ecdysone titer in infected larvae was significantly higher than that of the control larvae. We then used RNA sequencing (RNA-seq) to analyze silkworm PGs 48 h after BmNPV infection. The classifications of the 15 differentially expressed genes (DEGs) were mainly involved in the metabolic processes and pathways. The reverse transcription quantitative PCR (RT-qPCR) results of the DEGs in the PGs of BmNPV-infected larvae at 24, 48, and 72 h were generally consistent with the transcriptome data. The transcripts of BmTrypsin-1 and BmACSS3 were significantly increased from 24 to 72 h after BmNPV infection, indicating that they may be involved in the maturation process in the latter half of the fifth instar of silkworm larvae. This study was conducted to further our understanding of the complex biological processes in the interactions between BmNPV and its precocious metamorphic insect hosts.
Materials and Methods
Study Animals and Virus
Bombyx mori F50 strain larvae were reared on fresh mulberry leaves under a 12:12 h day/night cycle at 25 ± 1°C and 60% relative humidity. The majority of the fifth instar larvae started wandering on day 8, depending on the batch of the silkworm. The larvae underwent oral inoculation with a wild BmNPV T3 strain, and the OBs were obtained from the larvae hemolymph before the larvae died. The OBs were purified by repeated and differential centrifugation, as previously described (Rahman and Gopinathan, 2004).
Sample Collection
In total, 500 day-4 fifth instar larvae were orally infected with BmNPV using 2.0 × 106 OB/larva. Control larvae (n = 500) were fed the same volume of sterile distilled water. The larvae of the infected and control groups were maintained in isolation and reared under the same conditions. The PGs were entwined in pairs in the tracheal bush of the first spiracle (Supplementary Figure 1). The PGs were carefully removed from the larvae of the infected and control groups after 24, 48, and 72 h (Supplementary Figure 2). Hemolymph was collected from day 6 (48 hpi), day 6.5 (60 hpi), day 7 (72 hpi), and day 7.5 (84 hpi) fifth instar larvae of the infected and control groups for use in assays of the ecdysteroid titers among the different developmental stages.
Statistics of Precocious Maturation of Silkworms After BmNPV Infection
Day-4 fifth instar larvae were divided into 6 groups with 200 in each group. All 200 larvae used in each of the three independent experiments were orally infected with BmNPV using 2.0 × 106 OB/larva. The 200 control larvae used in each of the three independent experiments were fed with the same volume of sterile distilled water. The larvae of the infected and control groups were maintained in isolation and reared under the same conditions. Diseased and dead larvae were removed and counted during rearing. When the proportion of mature silkworms was > 5% (first gate), the statistics was started.
Cholesterol and 7-Dehydrocholesterol Feeding Experiments
As previously described (Wu et al., 2016), silkworm larvae were fed mulberry leaves supplemented with 8,000 mg/L of cholesterol and 7-dehydrocholesterol (7dC). Mulberry leaves supplemented with the same volume of sterile distilled water were used as the control. Day-5 fifth instar larvae were initially fed (first feed session) with cholesterol and 7dC supplemented leaves and then again 24 h later (second feed session). Replacement mulberry leaves were added 6 h after each feeding session. The proportion of mature vs. immature larvae was counted, and the second gate was determined to be the point at which the majority of the larvae had started maturing. The point at which all larvae (100%) had reached maturity defined the third gate. All experiments were repeated three times per group.
Assay of Ecdysteroid Titers in Hemolymph and Examination of Viral DNA in PGs
The hemolymph samples were homogenized in 50% MeOH (800 μl). The resultant homogenates were centrifuged, and the supernatant was used to assay the ecdysteroid titers using an Insect Ecdysone ELISA Kit (Shanghai MEILIAN Biotechnology Co., Ltd.) according to the manufacturer’s instructions. RT-PCR was used to analyze the BmNPV virus replication level. The total DNA was extracted from the PGs of the BmNPV-infected larvae at 24, 48, and 72 h, as well as from the control larvae at 48 h. The DNA templates (10 ng) were PCR amplified using primers for the BmNPV GP41 gene. The silkworm glyceraldehyde-3-phosphate dehydrogenase (BmGAPDH) was used as the internal control. The specific primers for each gene used in the RT-PCR are shown in Supplementary Table 1. The RT-PCR product of each gene was defined as previously described (Zhang et al., 2019).
Transcriptome Analysis
Total RNA was isolated from the PGs of the silkworm larvae using the TRIzol reagent (Invitrogen, New York, NY, United States) according to the manufacturer’s instructions. RNA purity was quantified using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, New York, NY, United States). Poly(A)-tailed RNA prepared using magnetic oligo (dT) beads was broken into short fragments using a fragment buffer and was then reverse transcribed to synthesize first-strand complementary DNA (cDNA) with a random primer. DNA polymerase I was then mixed with RNase H, deoxyribonucleotide triphosphate (dNTP), and the buffer solution to synthesize the complementary strand. The libraries were constructed using the Illumina methods and protocols, following the manufacturer’s instructions. The insert size and concentration of the cDNA library were both checked and quantified by an Agilent Bioanalyzer 2100 (Agilent Technologies, Inc., Santa Clara, CA) and Qubit® RNA Assay Kit (Life Technologies, CA, United States), respectively. RNA-seq was carried out using an Illumina HiSeq 2500 instrument (Illumina, San Diego, CA, United States). To obtain clean reads and ensure the quality of information analysis, the raw reads were filtered by removing the adapter sequences, empty reads, unknown nucleotides (ratio ≥ 10%), and low-quality reads with a basic mass value of Q ≤ 20, which accounted for more than 50% of the whole read length. The clean read assembly was performed according to a previous report (Grabherr et al., 2011). The paired-end clean reads were mapped to the B. mori genome using the software package TopHat2 (version 2.0.12) (Kim et al., 2013). The genome sequences and annotation file were downloaded from SilkDB. The RNA-seq reads were aligned and then used to construct transcripts with Cufflinks (version 2.1.1) (Trapnell et al., 2012). HTSeq (version 0.6.1) was used to count the reads mapped to each gene to quantify the gene expression levels (Anders et al., 2015). The fragments per kilobase of transcript per million mapped reads (FPKM) of each gene were then calculated based on the length of the gene and read count mapped to a given gene. Genes with a FPKM ≥ 1.0 were identified as “expressed.” A ratio (log2 fold change) between the infected and control groups of ≥ 1.5 was identified as the determinant of the DEGs. The raw data have been submitted to the Gene Expression Omnibus (GEO) database with the accession number GSE167875. The functional annotation of DEGs was performed using the Gene Ontology (GO) assignments (Gotz et al., 2008) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichments (Kanehisa et al., 2012).
Tissue Expression Patterns of DEGs
The PGs are important endocrine organs that are significantly different from other tissues in both their morphology and function. In silkworms, the day-3 fifth instar is the boundary for the whole larval development stage (Xu et al., 2019). To analyze the tissue expression patterns of the identified DEGs in the PGs, the PGs, head, integument, midgut, fat body, hemocyte, ovary, testis, Malpighian tubule, trachea, anterior silk gland (ASG), median silk gland (MSG), and posterior silk gland (PSG) of day-3 fifth instar larvae were collected. We detected the expression patterns in the multiple tissues of day-3 fifth instar larvae. Total RNA was extracted using the TRIzol reagent (Invitrogen, Carlsbad, CA, United States). Total RNA concentrations were quantified, and single-stranded cDNAs were synthesized. The BmGAPDH gene was used as an intrinsic control.
RT-qPCR Analysis
The genes selected according to the RNA-seq results were compared by RT-qPCR. Total RNA was extracted from the PGs samples of the infected and control groups at 24, 48, and 72 h. The first-strand cDNA was synthesized using the PrimeScript Reverse Transcriptase kit (TaKaRa, Dalian, China) according to the manufacturer’s instructions. RT-qPCR was performed as previously described (Xu et al., 2019). The BmGAPDH gene was used as an intrinsic control (Guo et al., 2016).
Results
Statistics of Precocious Maturation of Silkworm After BmNPV Infection
The point when the proportion of mature silkworms was >5% was considered as the first gate, the point when the majority of larvae started maturing was considered as the second gate, and the point when all larvae had reached maturity (100%) was considered as the third gate. The duration of the fifth instar larval stage of the B. mori F50 strain was almost 8.5 days (Figure 1). The day-4 fifth instar larvae infected with BmNPV matured early. The times that both the first- and second-gate mature silkworms appeared in the infected groups were 12 h earlier than in the control groups (Figure 1). The weights of mature silkworms in the infected groups were significantly decreased when compared with the control groups (Supplementary Table 2). The spinning process was normal, and there was no difference between the infected groups and the control groups. Approximately half of the larvae in the infected groups died during the larval–pupal stage. Compared with the control groups, the cocoon sizes and the weights of the pupae (female and male) were observably reduced in the BmNPV-infected groups (Supplementary Table 2). The fifth instar larvae underwent precocious maturation after infection with BmNPV. Moreover, the day-5 fifth instar larvae were fed with 8,000 mg/L cholesterol and 7dC (via supplemented leaves). The results of feeding with cholesterol and 7dC were also shown to induce precocious maturation when compared to the control, where a certain number of larvae exhibited an anal prolapse in each group fed with cholesterol and 7dC (Supplementary Table 3).
Figure 1. Statistics of the duration of the fifth instar and precocious maturation for day-4 fifth instar silkworm larvae infected with BmNPV. Values represent means ± SDs of three independent investigations. Significant differences are indicated by (*p < 0.05) or (**p < 0.01).
Assay of Ecdysteroid Titers and RT-PCR Analysis of Viral DNA After BmNPV Infection
Based on the findings of precocious maturation of silkworms after BmNPV infection, the molting and metamorphosis of silkworm requires the presence of 20-hydroxyecdysone (20E). We speculated that the BmNPV infection would have some influences on the ecdysone titer. Thus, the titers of ecdysone in infected and control larvae were determined by ELISA. The results indicated that the ecdysone titers were significantly higher in infected larvae than in the control larvae (Figure 2A). Meanwhile, RT-PCR was used to analyze the virus genomic DNA copies in the PGs at 24, 48, and 72 h after BmNPV infection. The expression of the BmNPV GP41 gene was detected in the PGs from 24 to 72 hpi (Figure 2B). The expression of the GP41 gene was not detected in the PGs from the uninfected larvae at 48 h (Figure 2B). These results were useful for selecting the time point for the RNA-seq experiments.
Figure 2. Assay of ecdysteroid titers in hemolymph and examination of virus genomic DNA copies in the prothoracic glands (PGs). (A) ELISA analysis was repeated three times for each set of protein samples. Values represent the means ± SDs of three independent determinations. Significant difference is indicated by (*p < 0.05). (B) The expression of the BmNPV GP41 gene was detected in the PGs of the BmNPV-infected silkworm larvae at 24, 48, and 72 h, as well as the control larvae at 48 h. M: DL2000 DNA marker; numbers 1–4 indicate the control group at 48 h and the BmNPV-infected groups at 24, 48, and 72 h, respectively.
General Information of RNA-Seq and DEGs
The square of the Pearson correlation coefficient (R2) between the four samples was > 0.938 (Supplementary Figure 3). The Q20 values for the clean reads (for each group) were above 95% (Supplementary Table 4). The percentage of clean sequences located on the genome was > 80%. These results indicated that the transcriptome data were assembled with high quality and can be used for further research. The number of expressed genes was 10,152 in the control groups and 10,404 in the BmNPV-infected groups (Supplementary Table 5). In total, seven upregulated and eight downregulated DEGs were screened out (Table 1 and Supplementary Table 6). The functions of the 15 DEGs were primarily located in the binding proteins of nucleic acids, ions, and proteins, and BmTrypsin-1 had a serine-type endopeptidase activity (Table 1). The DEGs were then annotated by GO analysis to determine their involvement in biological processes, molecular functions, and cellular components (Supplementary Figure 4). The upregulated expression genes were related to biological processes that were mainly focused on metabolic and biological processes (Supplementary Figure 4 and Supplementary Table 7). The downregulated expression genes were focused on biological and metabolic processes and respond to stress factors and stimuli in biological processes (Supplementary Figure 4 and Supplementary Table 7). Regarding the cellular components, only the downregulated expression genes were involved in the membrane and integral components of the membrane, and the upregulated expression genes were not enriched (Supplementary Figure 4 and Supplementary Table 7). Within the molecular function, the upregulated expression genes were primarily located in the catalytic activity, and the down regulated expression genes were involved in protein binding, hydrolase activity, and catalytic activity (Supplementary Figure 4 and Supplementary Table 7). There were some differences in the GO functional annotations between the upregulated and downregulated genes (Supplementary Figure 4). The KEGG pathway enrichment analysis of the identified DEGs showed that the enriched genes were mainly involved in pathways, including metabolic pathways, propanoate metabolism, beta-alanine metabolism, drug metabolism with other enzymes, pyrimidine metabolism, and pantothenate and coenzyme A (CoA) biosynthesis (Table 2 and Supplementary Table 8). The KEGG pathway enrichment analysis of BmACSS3 revealed that acyl-coenzyme A (AcCoA) synthase activity and BmACSS3 could be involved in the acyl-CoA to cholesterogenesis pathways.
Table 1. List of the differentially expressed genes in silkworm prothoracic glands with a 1.5-fold change after BmNPV infection.
Table 2. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis.
Spatial Expression Patterns of the Identified DEGs
We investigated the spatial expression patterns of the identified DEGs in multiple tissues of the PGs, head, integument, midgut, fat body, hemocyte, ovary, testis, Malpighian tubule, trachea, ASG, MSG, and PSG of day-3 fifth instar larvae. Such findings can further our understanding of the PGs and elucidate the expression characteristics of DEGs. Expression signals of all of the 15 DEGs were detected in the PGs (Figure 3). The genes of Novel01674, BmJing, and BmAryl were expressed only in PGs (Figure 3). BmACSS3 was expressed only in the PGs, head, and integument (Figure 3). The other 11 genes were expressed in multiple tissues (Figure 3).
Figure 3. Expression patterns of the differentially expressed genes (DEGs) in multiple tissues of day-3 fifth instar silkworm larvae. RT-PCR was performed, and the BmGAPDH gene was used as the internal control.
Expression Analysis of the Identified DEGs
We used RT-qPCR to investigate the relative expression levels of nine randomly selected genes of the DEGs in the PGs of BmNPV-infected larvae at 24, 48, and 72 h. The expression levels of BmNGDN, BmTrypsin-1, BmACSS3, and BmJing were upregulated in the transcriptome data, i.e., their transcripts were significantly increased from 24 to 72 h after BmNPV infection (Figures 4A–D). Meanwhile, the expression levels of BmPyd3, BmTitin, BmIGc2, Novel01674, and BmAryl were downregulated in the transcriptome data, i.e., their transcripts were significantly decreased from 24 to 72 h after BmNPV infection (Figures 4E–I). The changes in the transcription of the nine genes were generally consistent with the transcriptome data.
Figure 4. Analysis of the nine randomly selected genes of the differentially expressed genes (DEGs) in the prothoracic glands (PGs) of BmNPV-infected silkworm larvae. The expression levels of BmNGDN (A), BmTrypsin-1 (B), BmACSS3 (C), BmJing (D) were upregulated and the expression levels of BmPyd3 (E), BmTitin (F), BmIGc2 (G), Novel01674 (H), BmAryl (I) were downregulated from 24 to 72 h after BmNPV infection, respectively. Here, BmGAPDH was used as the internal control. The experiments were repeated three times. Values represent the means ± SDs of three independent experiments. Significant differences are indicated by (*p < 0.05) or (**p < 0.01).
Discussion
The mulberry silkworm is one of the best models to study insect physiology and biochemistry, especially to better understand the relationship between induction factors (external and internal) and development. In general, the larvae of tetramolter silkworms proceed through five instars and undergo molting between each instar. The last instar larva completes the larval–pupal transition. The developmental speed of silkworm larvae has been shown to be regular and constant within the same silkworm strain and when maintained under the same rearing conditions. The day-3 fifth instar larval stage is considered to represent the boundary for the larval stage.
In this study, the day-4 fifth instar larvae infected with BmNPV (using 2.0 × 106 OB/larva) matured early, and the ecdysone titer in infected larvae was significantly higher than that of the control larvae. In addition, BmNPV infection (using 2.0 × 107 OB/larva) also caused larvae precocious maturation (only at the first gate), followed by illness and death (data not shown). Meanwhile, day-5 fifth instar larvae fed with cholesterol and 7dC also exhibited precocious maturation. Cholesterol and 7dC supplementation in the latter half of the fifth instar shortened the fifth instar period. In contrast, cholesterol and 7dC supplementation in the first half of the fifth instar (days 1–3) prolonged the fifth instar period, but no results were observed for BmNPV infection in the first half of the fifth instar because the larvae died. Briefly, the prothoracicotropic hormone (PTTH) secreted by the brain stimulates the PGs to release ecdysteroid, which in turn induces larval or metamorphic ecdysis depending on the presence of juvenile hormone (JH) secreted by the corpora allata. Throughout the latter half of the fifth instar, the first step of the ecdysteroid titer increased in a stepwise manner, the second step of the ecdysteroid titer showed a small increase (which led to the silkworms wandering) followed by a plateau, and the third step of the ecdysteroid titer showed an initially gradual but then steep increase to reach a peak 1 day later (Mizoguchi et al., 2001).
Insights from host studies reveal that baculoviruses manipulate host behavior to enhance transmission to new victims. For example, baculoviruses enable infected larvae to continue to seek foliage and prolong insect feeding after infection, thus resulting in an increased OBs production (O’Reilly and Miller, 1989). The egt gene of NPV, expressed immediately, encodes an enzyme that inactivates the molting hormone 20E by transferring a sugar moiety from a nucleotide sugar donor to a hydroxyl group on 20E (Hoover et al., 2011). The ecdysone blood level is reduced by up to 90% in silkworms as a result of the transgenic expression of the egt gene of BmNPV and because egt expression in egt-transgenic silkworms prolongs the duration of the larval and pupal stages resulting in the arresting of the pupal-to-adult metamorphosis (Zhang et al., 2012). Interestingly, in the present study, silkworm larvae infected with BmNPV in the latter half of the fifth instar showed precocious molting and metamorphosis and a higher level of hemolymph ecdysone titer. The egt gene of BmNPV is dispensable for normal virus production (Katsuma et al., 2008). The fast-killing phenotype is observed in the three egt-mutated BmNPVs only when the infection process progresses through silkworm larval–larval transition, but under infection in the middle stages of the fifth instar, the slow-killing phenotype is observed than that of the wild-type virus-infected larvae (Katsuma and Shimada, 2015). In particular, in the gluttonous stage, silkworms synthesize an enormous amount of silk proteins in the silk gland (Xu et al., 2019), and silk proteins are dispensable for normal silkworm development. A certain amount of silk proteins can remain in the body (incomplete spinning) that can lead to an incomplete larval–pupal transition. Silkworm and BmNPV interactions are largely dependent on the developmental stage of the host larvae infected by the virus. In addition, the overproduction of silkworm PTTH induces higher than normal levels of hemolymph ecdysteroids, which have been found to inhibit the pathogenicity of the virus, but did not have any observable effects on the development of infected Spodoptera frugiperda larvae (O’Reilly et al., 1995). Moreover, insect innate immunity can be activated by 20E and 20E, which induce antimicrobial peptide (AMP) gene expression and thus act as immune activators (Dimarcq et al., 1997; Roxstrom-Lindquist et al., 2005; Flatt et al., 2008). The ecdysteroid titer showed a small increase followed by a plateau that occurred 1 day before the silkworms started wandering. Thereafter, the titer increased gradually and then steeply to reach a peak (where the majority of silkworms had started wandering) the following day (Mizoguchi et al., 2001). Therefore, we speculated that BmNPV infection in the latter half of the fifth instar of silkworm larvae induced precocious molting and metamorphosis and a higher level of hemolymph ecdysone titer, which would enable infected larvae to complete their larval–pupal transition.
Briefly, PGs are an important endocrine organ with characteristically cholesterol-rich tissue as the main site of synthetic ecdysteroids (Igarashi et al., 2018). In this study, the RT-PCR results confirmed that the silkworm PGs were infected by BmNPV through oral inoculation. Seven upregulated and eight downregulated DEGs were identified from silkworm PGs sequenced by RNA-seq 48 h after BmNPV infection. The RT-qPCR results of the DEGs in the PGs of BmNPV-infected larvae at 24, 48, and 72 h were generally consistent with the transcriptome data. The classifications of the 15 DEGs were primarily located in binding activity of nucleic acids, ions, and proteins that were mainly involved in the metabolic processes and pathways. The spatial expression profiles of Novel01674, BmAryl, and BmJing indicated that they were specifically expressed in the silkworm PGs. The KEGG pathway enrichment analysis of BmACSS3 (BGIBMGA010070) revealed that acyl-coenzyme A (AcCoA) synthase activity and BmACSS3 could be involved in the acyl-CoA to cholesterogenesis pathways. The acyl-coA synthetase catalyzes fatty acids to form a thioester with CoA, which is a common initial step of all fatty acid metabolic processes (Watkins, 1997). Fatty acids are the building blocks of many lipids, including triacylglycerol and cholesteryl esters. RNAi-mediated knockdown of the acyl-coenzyme A synthetase gene prolongs and extends the maximum lifespan (Eisenberg et al., 2014). Trypsin-1 serine protease (BmTrypsin-1) had serine-type endopeptidase activity. SPs play crucial roles in insect development and innate immunity. RNAi-mediated silencing of SPs results in severe molting defects, specifically by reducing the expression of genes in the 20E synthesis and signaling pathway, and increases larval sensitivity to bacteria (Broehan et al., 2010; Yang et al., 2019). In silkworm PGs, the transcripts of BmTrypsin-1 and BmACSS3 were significantly increased from 24 to 72 h after BmNPV infection. BmTrypsin-1 and BmACSS3 may be involved in the maturation process in the latter half of the fifth instar of silkworm larvae.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
P-ZX performed the experiment and wrote the manuscript. M-RZ performed the literature review and analyzed the data. X-YW prepared the illustrations and collected the data. Y-CW suggested important research points. All authors have read and approved the final version of the manuscript.
Funding
This research was financially supported by the Jiangsu Provincial Natural Science Foundation of China (Grant No. BK2012273) and the National Natural Science Foundation of China (Grant Nos. 31302035 and 2735011802).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2021.650972/full#supplementary-material
Supplementary Figure 1 | View of the prothoracic glands (PGs) entwined in the tracheal bush of the first spiracle.
Supplementary Figure 2 | View of the prothoracic glands (PGs) from a fifth instar silkworm larva.
Supplementary Figure 3 | Analysis of the correlation of RNA-seq data. (A) Diagram of the correlation coefficients between samples. (B) Correlation between C1 and C2. (C) Correlation between C1 and N1. (D) Correlation between C1 and N2. (E) Correlation between C2 and N1. (F) Correlation between C2 and N2. (G) Correlation between N1 and N2. C1 and C2 indicate the two independent biological experiments of transcriptome sequencing of prothoracic glands (PGs) in the control groups, respectively. N1 and N2 indicate the two independent biological experiments of transcriptome sequencing of PGs in the control groups, respectively.
Supplementary Figure 4 | GO enrichment analysis of the differentially expressed genes (DEGs). Genes were annotated by the biological process, cellular component, and molecular function.
Supplementary Table 1 | Primers of genes for RT-PCR and RT-qPCR.
Supplementary Table 2 | The weight of mature silkworms and pupae in BmNPV-infected and control groups.
Supplementary Table 3 | The results of feeding experiment with cholesterol, and 7-dehydrocholesterol (7dC), using water as control.
Supplementary Table 4 | Summary statistics of prothoracic gland of Bombyx mori genes based on transcriptome data.
Supplementary Table 5 | A list of genes with FPKM ≥ 1.0 between infected and control groups.
Supplementary Table 6 | A list of differential expression genes between infected and control groups.
Supplementary Table 7 | Infected vs. control DEG_GO_enrichment result.
Supplementary Table 8 | Infected vs. control all DEG_KEGG_pathway_enrichment result.
References
Anders, S., Pyl, P. T., and Huber, W. (2015). HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics 31, 166–169. doi: 10.1093/bioinformatics/btu638
Broehan, G., Arakane, Y., Beeman, R. W., Kramer, K. J., Muthukrishnan, S., and Merzendorfer, H. (2010). Chymotrypsin-like peptidases from Tribolium castaneum: a role in molting revealed by RNA interference. Insect Biochem. Mol. Biol. 40, 274–283. doi: 10.1016/j.ibmb.2009.10.009
Chen, T. T., Tan, L. R., Hu, N., Dong, Z. Q., Hu, Z. G., Qin, Q., et al. (2019). Specific genes related to nucleopolyhedrovirus in Bombyx mori susceptible and near-isogenic resistant strains through comparative transcriptome analysis. Insect Mol. Biol. 28, 473–484. doi: 10.1111/imb.12566
Dimarcq, J. L., Imler, J. L., Lanot, R., Ezekowitz, R. A., Hoffmann, J. A., Janeway, C. A., et al. (1997). Treatment of l(2)mbn Drosophila tumorous blood cells with the steroid hormone ecdysone amplifies the inducibility of antimicrobial peptide gene expression. Insect Biochem. Mol. Biol. 27, 877–886. doi: 10.1016/s0965-1748(97)00072-6
Du, X., and Thiem, S. M. (1997). Responses of insect cells to baculovirus infection: protein synthesis shutdown and apoptosis. J. Virol. 71, 7866–7872. doi: 10.1128/jvi.71.10.7866-7872.1997
Eisenberg, T., Schroeder, S., Andryushkova, A., Pendl, T., Kuttner, V., Bhukel, A., et al. (2014). Nucleocytosolic depletion of the energy metabolite acetyl-coenzyme a stimulates autophagy and prolongs lifespan. Cell Metab. 19, 431–444. doi: 10.1016/j.cmet.2014.02.010
Flatt, T., Heyland, A., Rus, F., Porpiglia, E., Sherlock, C., Yamamoto, R., et al. (2008). Hormonal regulation of the humoral innate immune response in Drosophila melanogaster. J. Exp. Biol. 211, 2712–2724. doi: 10.1242/jeb.014878
Gomi, S., Majima, K., and Maeda, S. (1999). Sequence analysis of the genome of Bombyx mori nucleopolyhedrovirus. J. Gen. Virol. 80(Pt 5), 1323–1337. doi: 10.1099/0022-1317-80-5-1323
Gotz, S., Garcia-Gomez, J. M., Terol, J., Williams, T. D., Nagaraj, S. H., Nueda, M. J., et al. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36, 3420–3435. doi: 10.1093/nar/gkn176
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–U130.
Guo, H., Jiang, L., and Xia, Q. (2016). Selection of reference genes for analysis of stress-responsive genes after challenge with viruses and temperature changes in the silkworm Bombyx mori. Mol. Genet. Genomics 291, 999–1004. doi: 10.1007/s00438-015-1125-4
Guo, H. Z., Sun, Q., Wang, B. B., Wang, Y. M., Xie, E. Y., Xia, Q. Y., et al. (2019). Spry is downregulated by multiple viruses to elevate ERK signaling and ensure viral reproduction in silkworm. Dev. Comp. Immunol. 98, 1–5. doi: 10.1016/j.dci.2019.04.001
Hoover, K., Grove, M., Gardner, M., Hughes, D. P., McNeil, J., Slavicek, J., et al. (2011). A gene for an extended phenotype. Science 333, 1401–1401.
Huang, H. L., Wu, P., Zhang, S. L., Shang, Q., Yin, H. T., Hou, Q. R., et al. (2019). DNA methylomes and transcriptomes analysis reveal implication of host DNA methylation machinery in BmNPV proliferation in Bombyx mori. BMC Genomics 20:736. doi: 10.1186/s12864-019-6146-7
Huh, N. E., and Weaver, R. F. (1990). Categorizing some early and late transcripts directed by the Autographa californica nuclear polyhedrosis virus. J. Gen. Virol. 71(Pt 9), 2195–2200. doi: 10.1099/0022-1317-71-9-2195
Igarashi, F., Ogihara, M. H., Iga, M., and Kataoka, H. (2018). Cholesterol internalization and metabolism in insect prothoracic gland, a steroidogenic organ, via lipoproteins. Steroids 134, 110–116. doi: 10.1016/j.steroids.2018.01.012
Ikeda, M., Yamada, H., Hamajima, R., and Kobayashi, M. (2013). Baculovirus genes modulating intracellular innate antiviral immunity of lepidopteran insect cells. Virology 435, 1–13. doi: 10.1016/j.virol.2012.10.016
Jiang, L. (2021). Insights into the antiviral pathways of the silkworm Bombyx mori. Front. Immunol. 12:639092. doi: 10.3389/fimmu.2021.639092
Jiang, L., Goldsmith, M. R., and Xia, Q. Y. (2021a). Advances in the arms race between silkworm and Baculovirus. Front. Immunol. 12:628151. doi: 10.3389/fimmu.2021.628151
Jiang, L., Liu, W. Q., Guo, H. Z., Dang, Y. H., Cheng, T. C., Yang, W. Y., et al. (2019). Distinct functions of Bombyx mori peptidoglycan recognition protein 2 in immune responses to bacteria and viruses. Front. Immunol. 10:776. doi: 10.3389/fimmu.2019.00776
Jiang, L., and Xia, Q. (2014). The progress and future of enhancing antiviral capacity by transgenic technology in the silkworm Bombyx mori. Insect. Biochem. Mol. Biol. 48, 1–7. doi: 10.1016/j.ibmb.2014.02.003
Jiang, L., Xie, E. Y., Guo, H. Z., Sun, Q., Liuli, H. Y., Wang, Y. M., et al. (2021b). Heat shock protein 19.9 (Hsp19.9) from Bombyx mori is involved in host protection against viral infection. Dev. Comp. Immunol. 114:103790. doi: 10.1016/j.dci.2020.103790
Jiang, L., Zhao, P., Cheng, T. C., Sun, Q., Peng, Z. W., Dang, Y. H., et al. (2013a). A transgenic animal with antiviral properties that might inhibit multiple stages of infection. Antivir. Res. 98, 171–173. doi: 10.1016/j.antiviral.2013.02.015
Jiang, L., Zhao, P., Wang, G. H., Cheng, T. C., Yang, Q., Jin, S. K., et al. (2013b). Comparison of factors that may affect the inhibitory efficacy of transgenic RNAi targeting of baculoviral genes in silkworm, Bombyx mori. Antivir. Res. 97, 255–263. doi: 10.1016/j.antiviral.2012.12.020
Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40, D109–D114.
Katsuma, S., Kawaoka, S., Mita, K., and Shimada, T. (2008). Genome-wide survey for baculoviral host homologs using the Bombyx genome sequence. Insect. Biochem. Mol. 38, 1080–1086. doi: 10.1016/j.ibmb.2008.05.008
Katsuma, S., Mita, K., and Shimada, T. (2007). ERK- and JNK-Dependent signaling pathways contribute to Bombyx mori nucleopolyhedrovirus infection. J. Virol. 81, 13700–13709. doi: 10.1128/jvi.01683-07
Katsuma, S., and Shimada, T. (2015). The killing speed of egt-inactivated Bombyx mori nucleopolyhedrovirus depends on the developmental stage of B-mori larvae. J. Invertebr. Pathol. 126, 64–70. doi: 10.1016/j.jip.2015.01.012
Kim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., and Salzberg, S. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14:R36.
Manjunatha, G. K. S., Peter, A., Naika, M. B. N., Niranjana, P., and Shamprasad, P. (2018). Identification of in-vitro red fluorescent protein with antipathogenic activity from the midgut of the silkworm (Bombyx Mori L.). Protein Pept. Lett. 25, 302–313. doi: 10.2174/0929866525666180115121853
Mita, K., Kasahara, M., Sasaki, S., Nagayasu, Y., Yamada, T., Kanamori, H., et al. (2004). The genome sequence of silkworm, Bombyx mori. DNA Res. 11, 27–35. doi: 10.1093/dnares/11.1.27
Mizoguchi, A., Ohashi, Y., Hosoda, K., Ishibashi, J., and Kataoka, H. (2001). Developmental profile of the changes in the prothoracicotropic hormone titer in hemolymph of the silkworm Bombyx mori: correlation with ecdysteroid secretion. Insect. Biochem. Mol. Biol. 31, 349–358. doi: 10.1016/s0965-1748(00)00127-2
Nakazawa, H., Tsuneishi, E., Ponnuvel, K. M., Furukawa, S., Asaoka, A., Tanaka, H., et al. (2004). Antiviral activity of a serine protease from the digestive juice of Bombyx mori larvae against nucleopolyhedrovirus. Virology 321, 154–162. doi: 10.1016/j.virol.2003.12.011
O’Reilly, D. R., Kelly, T. J., Masler, E. P., Thyagaraja, B. S., Robson, R. M., Shaw, T. C., et al. (1995). Overexpression of Bombyx mori prothoracicotropic hormone using baculovirus vectors. Insect. Biochem. Mol. Biol. 25, 475–485. doi: 10.1016/0965-1748(94)00087-f
O’Reilly, D. R., and Miller, L. K. (1989). A baculovirus blocks insect molting by producing ecdysteroid UDP-glucosyl transferase. Science 245, 1110–1112. doi: 10.1126/science.2505387
Ponnuvel, K. M., Nakazawa, H., Furukawa, S., Asaoka, A., Ishibashi, J., Tanaka, H., et al. (2003). A lipase isolated from the silkworm Bombyx mori shows antiviral activity against nucleopolyhedrovirus. J. Virol. 77, 10725–10729. doi: 10.1128/jvi.77.19.10725-10729.2003
Ponnuvel, K. M., Nithya, K., Sirigineedi, S., Awasthi, A. K., and Yamakawa, M. (2012). In vitro antiviral activity of an alkaline trypsin from the digestive juice of Bombyx mori larvae against nucleopolyhedrovirus. Arch. Insect Biochem. Physiol. 81, 90–104. doi: 10.1002/arch.21046
Rahman, M. M., and Gopinathan, K. P. (2004). Systemic and in vitro infection process of Bombyx mori nucleopolyhedrovirus. Virus Res. 101, 109–118. doi: 10.1016/j.virusres.2003.12.027
Roxstrom-Lindquist, K., Assefaw-Redda, Y., Rosinska, K., and Faye, I. (2005). 20-Hydroxyecdysone indirectly regulates Hemolin gene expression in Hyalophora cecropia. Insect Mol. Biol. 14, 645–652. doi: 10.1111/j.1365-2583.2005.00593.x
Shirata, N., Ikeda, M., and Kobayashi, M. (2010). Identification of a Hyphantria cunea nucleopolyhedrovirus (NPV) gene that is involved in global protein synthesis shutdown and restricted Bombyx mori NPV multiplication in a B. mori cell line. Virology 398, 149–157. doi: 10.1016/j.virol.2009.11.049
Sunagar, S. G., Savanurmath, C. J., and Hinchigeri, S. B. (2011). The profiles of red fluorescent proteins with antinucleopolyhedrovirus activity in races of the silkworm Bombyx mori. J. Insect Physiol. 57, 1707–1714. doi: 10.1016/j.jinsphys.2011.09.009
Toufeeq, S., Wang, J., Zhang, S. Z., Li, B., Hu, P., Zhu, L. B., et al. (2019). Bmserpin2 is involved in BmNPV infection by suppressing melanization in Bombyx mori. Insects 10:399. doi: 10.3390/insects10110399
Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016
Wu, F., Wang, P. Y., Zhao, Q. L., Kang, L. Q., Xia, D. G., Qiu, Z. Y., et al. (2016). Mutation of a cuticle protein gene, BmCPG10, is responsible for silkworm non-moulting in the 2nd instar mutant. PLoS One 11:e0153549. doi: 10.1371/journal.pone.0153549
Xia, Q., Zhou, Z., Lu, C., Cheng, D., Dai, F., Li, B., et al. (2004). sequence for the genome of the domesticated silkworm (Bombyx mori). Science 306, 1937–1940. doi: 10.1126/science.1102210
Xu, P., Zhang, M., Qian, P., Li, J., Wang, X., and Wu, Y. (2019). iTRAQ-based quantitative proteomic analysis of digestive juice across the first 48 hours of the fifth instar in silkworm larvae. Int. J. Mol. Sci. 20:6113. doi: 10.3390/ijms20246113
Yang, W. J., Chen, C. X., Yan, Y., Xu, K. K., and Li, C. (2019). Clip-domain serine protease gene (LsCLIP3) is essential for larval-pupal molting and immunity in Lasioderma serricorne. Front. Physiol. 10:1631. doi: 10.3389/fphys.2019.01631
Zhang, M. R., Xu, P. Z., Pang, H. L., Chen, T., and Zhang, G. Z. (2019). Expression analysis of mRNA decay of maternal genes during Bombyx mori maternal-to-zygotic transition. Int. J. Mol. Sci. 20:5651. doi: 10.3390/ijms20225651
Zhang, S. Z., Zhu, L. B., You, L. L., Wang, J., Cao, H. H., Liu, Y. X., et al. (2020). A novel digestive proteinase lipase member H-A in Bombyx mori contributes to digestive juice antiviral activity against B. mori nucleopolyhedrovirus. Insects 11:154. doi: 10.3390/insects11030154
Zhang, X., Xue, R. Y., Cao, G. L., Hu, X. L., Wang, X. J., Pan, Z. H., et al. (2012). Effects of egt gene transfer on the development of Bombyx mori. Gene 491, 272–277. doi: 10.1016/j.gene.2011.09.026
Keywords: Bombyx mori, Bombyx mori nucleopolyhedrovirus, prothoracic gland, transcriptome, 20-hydroxyecdysone
Citation: Xu P-Z, Zhang M-R, Wang X-Y and Wu Y-C (2021) Precocious Metamorphosis of Silkworm Larvae Infected by BmNPV in the Latter Half of the Fifth Instar. Front. Physiol. 12:650972. doi: 10.3389/fphys.2021.650972
Received: 08 January 2021; Accepted: 06 April 2021;
Published: 10 May 2021.
Edited by:
Liang Jiang, Southwest University, ChinaReviewed by:
Jingchen Sun, South China Agricultural University, ChinaWei Yu, Zhejiang Sci-Tech University, China
Copyright © 2021 Xu, Zhang, Wang and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ping-Zhen Xu, eHB6MTk4MjQ5QGp1c3QuZWR1LmNu; Yang-Chun Wu, MTk5NDAwMDAxOTM5QGp1c3QuZWR1LmNu