Development of a loop-mediated isothermal amplification (LAMP) assay for rapid and visual detection of Anguillid herpesvirus 1.

China has the largest aquaculture eel production in the world. High-density cultivation pattern often results in an outbreak of epidemic diseases. Since the 1990s, eel "mucus sloughing and hemorrhagic septicemia disease" was often broke out in China, and brought huge economic losses to eel breeders. Anguillid herpesvirus 1 (AngHV) was detected and isolated from the diseased eel, and proved to be the pathogen of the disease. In this study, a loop-mediated isothermal amplification (LAMP) assay was developed for rapid, sensitive, and specific detection of AngHV. A set of six primers targeting the ORF51 gene of AngHV was designed, which could effectively detect purified AngHV virions, AngHV-infected cells, or eel tissue samples. The suitable reaction temperature is 63℃, and the reaction time is 40 min. There was no cross-reaction with eel and other fish viruses, including Infectious pancreatic necrosis virus (IPNV), Marine birnavirus (MABV), Rana grylio virus (RGV), Cyprinid herpesvirus 3 (CyHV-3), and Eel iridovirus (EIV). The lower detection limit of the AngHV LAMP assay is 10 copies of AngHV genome DNA, which is at least 100 times more sensitive than conventional PCR in detecting AngHV. The assay could effectively detect AngHV from collected samples with typical clinical symptoms of AngHV infection. It suggested that the LAMP assay could be used in specific detection of AngHV and has great potential for early diagnosis of AngHV infection in the farm.

Immune responses and protective efficacy of American eel (Anguilla rostrata) immunized with a formalin-inactivated vaccine against Anguillid herpesvirus.

Anguillid herpesvirus 1 (AngHV), the causative agent of "mucus sloughing and hemorrhagic septicemia disease", causes serious infectious diseases in farmed eel. Among the effective prevention and control strategies, vaccination is one of the most effective approaches. However, no vaccine for AngHV is available. Our study developed a formalin-inactivated AngHV vaccine and evaluated its performance in American eels. Initially, AngHV-FJ, a strain of AngHV, was inactivated completely by 0.1 % formaldehyde, mixed with adjuvant Montanide ISA 763 A VG (763A). Then, vaccines containing different amount of antigen (3 × 106 PFU, 3 × 105 PFU, 3 × 104 PFU, 3 × 103 PFU) were immunized in each American eels. The results showed that the 3 × 105 PFU/fish was the proper dose. The inactivated AngHV vaccine was proven safe for American eels by back intramuscular injection. The results of twice immunization showed that antibody production peaked in the 8th week after the first immunization, and the antibody titer was 1:64,000. Furthermore, the immunized fishes challenged with AngHV (105 PFU/ml immersion) showed a significantly lower incidence rate (33.33 %) than the control group (95.65 %). The survival of the fish in the vaccine group (94.44 %) was significantly higher than the control group (60.87 %). The relative survival rate of the vaccinated group was 85.80 %. Also, vaccine group tissue collected at 7th d post-challenge showed reduced tissue damage and a lower virus load than the control group. The expression of cytokines of IL-1ß, IFN-α, IFN-γ, Mx1, RIG-1, and IRF-3, were significantly lower in the vaccine group than the control group at the 7th and 14th d post-challenge. Overall, the formalin-inactivated AngHV vaccine was safe and had immune protective effects against AngHV infection.

Development and application of a recombinase-aided amplification combined with a lateral flow dipstick assay for rapid visual detection of anguillid herpesvirus 1.

Eel (Anguilla sp.) is an important freshwater-cultured species with high economic value in China. Anguillid herpesvirus 1 (AngHV-1) has been proven to be the pathogen of "mucus sloughing and haemorrhagic septicaemia disease" in eels, resulting in significant mortality and substantial losses to the eel industry. Current diagnostic methods for detecting AngHV-1 are limited to laboratory-based tests, for example, conventional end-point PCR and qPCR. Therefore, there is an urgent need to develop an accurate, rapid, and simple detection method for on-site diagnosis of AngHV-1. In this study, we developed a recombinase-aided amplification combined lateral flow dipstick (RAA-LFD) assay for the detection of AngHV-1. The RAA-LFD assay can be performed within a temperature range of 18-45°C, with a reaction time of just 10 min for amplification. Importantly, the established RAA-LFD assay exhibited no reactivity with other common aquatic viral pathogens, indicating its high specificity. The limit of detection for this method is 102 copies of AngHV-1, which is more sensitive than the established conventional end-point PCR method similarly targeting ORF95. Clinical detection of the diseased samples demonstrated that the accuracy of RAA-LFD was significantly higher than that of the conventional end-point PCR. In conclusion, the developed RAA-LFD assay has proven to be a convenient, rapid, sensitive, and reliable tool for on-site diagnosis of AngHV-1. This advancement will be invaluable for the prevention and control of AngHV-1 in the eel farming industry.

Metagenomic analysis fecal microbiota of dysentery-like diarrhoea in a pig farm using next-generation sequencing.

Porcine enteric diseases including swine dysentery involves a wide range of possible aetiologies and seriously damages the intestine of pigs of all ages. Metagenomic next-generation sequencing is commonly used in research for detecting and analyzing pathogens. In this study, the feces of pigs from a commercial swine farm with dysentery-like diarrhea was collected and used for microbiota analysis by next-generation sequencing. While Brachyspira spp. was not detected in diarrheal pig fecal samples, indicating that the disease was not swine dysentery. The quantity of microbial population was extremely lowered, and the bacterial composition was altered with a reduction in the relative abundance of the probiotics organisms, Firmicutes and Bacteroidetes, with an increase in pathogens like Fusobacterium and Proteobacteria, in which the specific bacteria were identified at species-level. Viral pathogens, porcine circovirus type 2, porcine lymphotropic herpesviruses 1, and porcine mastadenovirus A were also detected at pretty low levels. Carbohydrate-active enzymes (CAZy) analysis indicated that the constitute of Firmicutes and Bacteroidete were also changed. Further, the Kyoto Encyclopedia of Genes and Genomes (KEGG) alignment analysis indicated that the microbiota of diarrheal pigs had a lower ability in utilizing energy sources but were enriched in multi-drug resistance pathways. Comprehensive Antibiotic Resistance Database (CARD) and Virulence Factors of Pathogenic Bacteria (VFDB) analysis indicated that genes for elfamycin and sulfonamide resistance and the iron uptake system were enriched in diarrheal pigs. This revealed potential bacterial infection and can guide antibiotic selection for treating dysentery. Overall, our data suggested that alterations in both the population and functional attributes of microbiota in diarrheal pigs with decreased probiotic and increased pathogenic microorganisms. These results will help elucidate the mechanism of dysentery-like diarrhea and the development of approaches to control the disease.

Evaluation of housekeeping genes as references for quantitative real-time PCR analysis of European eel, Anguilla anguilla.

Eels are important aquaculture species for which an increasing number of reference genes are being identified and applied. In this study, five housekeeping genes [RPL7 (ribosomal protein L7), 18 S (18 S ribosomal RNA), EF1A (elongation factor 1α), ACTB (ß-actin) and GAPDH (glyceraldehyde-3-phosphate dehydrogenase)] were chosen to evaluate their reliability as reference genes for quantitative real-time PCR (qPCR) for the study of Anguilla anguilla. The expression of the selected genes in different eel tissues was determined using qPCR at different growth stages or upon challenge by Anguillid herpesvirus (AngHV), and the expression levels of these genes were then compared and evaluated using the geNorm and NormFinder algorithms. Then, RefFinder was used to comprehensively rank the examined housekeeping genes. Interestingly, the expression of the evaluated housekeeping genes exhibited tissue-dependent and treatment-dependent variations. In different growth periods A. anguilla tissues, the most stable genes were the following: ACTB in mucus; 18 S in skin and kidney; RPL7 in muscle, gill, intestine and brain; EF1A in heart and liver; and GAPDH in spleen. In contrast, in AngHV-challenged A. anguilla tissues, the most stable genes were the following: 18 S in mucus; RPL7 in skin, gill, heart, spleen, kidney and intestine; EF1A in muscle and liver; and ACTB in brain. Further comparison analysis indicated that the expression of RPL7 and EF1A was stable in multiple A. anguilla tissues in different growth periods and in eels challenged by AngHV. Nonetheless, the expression level of GAPDH in eel tissues was lower, and it was unstable in several tissues. These results indicated that the selection of reference genes for qPCR analysis in A. anguilla should be made in accordance with experimental parameters, and both RPL7 and EF1A could be used as reference genes for qPCR study of A. anguilla at different growth stages or upon challenge by AngHV. The reference genes identified in this study could improve the accuracy of qPCR data and facilitate further studies aimed at understanding the biology of eels.

Expression Profiles of Zebrafish (Danio rerio) Lysozymes and Preparation of c-Type Lysozyme with High Bacteriolytic Activity against Vibrio vulnificus.

Lysozyme is a crucial component of the host's innate immune system. Due to its natural non-toxic and harmless characteristics, lysozyme is considered to be an ideal antibiotic substitute. In this study, we analyzed the expression profiles of lysozymes from zebrafish (Danio rerio) in uninfected or V. vulnificus-infected tissues using real-time quantitative PCR (qPCR). Furthermore, lysozymes that might be involved in the defense against V. vulnificus were selected for over-expression, and the antibacterial activity of over-expressed lysozyme proteins were evaluated using V. vulnificus. The results showed that three types of zebrafish lysozyme, i.e., c-type lysozyme (DrLysC), g1-type lysozyme (DrLysG1), and g2-type lysozyme (DrLysG2), were identified, and V. vulnificus infection significantly changed the expression levels of DrLysC and DrLysG1. Then, DrLysC and DrLysG1 were over-expressed in E. coli, and the purified recombinant DrLysC (rDrLysC) showed more potent antibacterial activity against V. vulnificus. This finding lays the foundation for further application of rDrLysC to treat V. vulnificus infection.

Proteomic Profiling Skin Mucus of European Eel Anguilla anguilla Infected with Anguillid Herpesvirus.

Anguillid herpesvirus 1 (AngHV) is an important viral pathogen affecting eel. This study was designed to investigate the potential molecular mechanisms and immune response elicited at the protein levels in the skin mucus of AngHV-infected Anguilla anguilla. Tandem mass tag (TMT)-labelling proteomics with the liquid chromatography tandem mass spectrometry (LC-MS/MS) was used for performing quantitative identification of the proteins. In addition, the quantitative protein amount was detected by parallel reaction monitoring (PRM) analysis. A total of 3486 proteins were identified, of which 2935 were quantified. When a protein fold change was greater than 1.3 or less than 0.76, it indicated a differentially expressed protein (DEP). Overall, 187 up-regulated proteins and 126 down-regulated proteins were detected, and most of the DEPs were enriched in the CAMs pathway, intestinal immune pathway, herpes simplex virus 1 infection pathway, phagosome pathway and p53 signaling pathway. The results of the DEPs detected by PRM were highly consistent with the results of the TMT-labelled quantitative proteomic analysis. The findings of this study provide an important research basis for further understanding the pathogenesis of AngHV.

Bombyx mori Nucleopolyhedrovirus p26 Is Associated with Viral Late Stage Replication.

Bombyx mori nucleopolyhedrovirus (BmNPV) p26 is conserved among all Lepidoptera baculoviruses that have been completely sequenced thus far, and some baculoviruses even have two copies of p26, which suggested that p26 may play an important role in the virus infection cycle. This study aimed to characterize BmNPV p26. We found that BmNPV p26 transcripts were detectable as early as 3 h post-infection (hpi), and the transcript levels rapidly increased starting from 12 hpi. Western blot analysis using an anti-p26 polyclonal antibody demonstrated that the corresponding protein was also detectable from 6 hpi in BmNPV-infected cell lysates. Immunofluorescence analysis demonstrated that p26 was mainly dispersed in the infected cell cytoplasm, whereas the over-expressed fusion protein EGFP-p26 also accumulated in the nucleus. These results indicated that p26 is an early BmNPV gene and has functions both in the cytoplasm and the nucleus. RNAi-based knockdown of p26 could produce infectious virus and normal-appearing virions but decreased budded virus (BV) production in BmNPV-infected cells at 72 hpi. Moreover, the results of further quantitative PCR (Q-PCR) analysis indicated that the gp64 and p74 transcripts levels decreased significantly. These results indicated that BmNPV p26 may be associated with BmNPV replication during the late infection stage.

Isolation of a novel adomavirus from cultured American eels, Anguilla rostrata, with haemorrhagic gill necrosis disease.

Recently, the culture of American eels (Anguilla rostrate) in China has been impacted by emergence of a disease with signs of haemorrhagic gill necrosis. The gills of diseased eels are covered with petecchia and they bleed when the operculum is pressed. In this study, a novel American eel adomavirus (AEAdoV) was isolated from the diseased eels using the eel ovary cell line (EO). The virus proliferated in the EO cells with a maximum TCID50 /ml of 106.29 ± 0.23 at 6 days post-infection. The virions were non-enveloped with a diameter of 75-85 nm and shown to be a DNA virus upon 5-iodo-2'-deoxyuridine (IDU) treatment. PCR assays showed that AEAdoV encodes a superfamily 3 helicases (S3H) replicase and shared high similarities with Anguilla marmorata adomavirus (MEAdoV). Although no clinical signs or mortality was observed among the eels injected with AEAdoV, the virus was reisolated from livers, kidneys and gills of injected eels at 35 days post-injection. Our results suggested that AEAdoV exhibited a latent infection in A. rostrata. The pathogenicity of the AEAdoV needs to be confirmed further.

Anguillid herpesvirus 1 (AngHV) ORF95 encodes a late, structural envelope protein.

Anguillid herpesvirus 1 (AngHV) is one of the vital pathogenic agents found in the wild and cultured eel populations, which has brought significant losses to eel culture industry in China. In this study, AngHV ORF95 was characterized. Bioinformatics analysis showed that ORF95 putatively encodes a structural protein that is homologous to hemagglutinin-esterase (HE) protein of infectious salmon anemia virus (ISAV). Temporal transcription and expression analysis indicated that ORF95 is a viral late gene. Subcellular localization analysis revealed that ORF95 was predominantly localized in the cytoplasm. Further, western blot analysis indicated that ORF95 is a structural protein of virion envelope. These results provide a novel basis to make further efforts to clarify the function of ORF95 in the process of AngHV infection and the possibility to use ORF95 as antigen to develop AngHV subunit vaccine.

Development and characterization of a continuous cell line (EL) from the liver of European eel Anguilla anguilla.

In the present study, a new hepatic tissue-origin cell line from European eel Anguilla anguilla has been developed and characterized. This cell line designated EL has been maintained in Leibovitz L-15 supplemented with 10% fetal bovine serum over 72 months, and subcultured more than 90 times. The EL cell line consisted predominantly of fibroblast-like cells, which could survive over 100 days in vitro, and could grow at 15-32°C. The optimum temperature for growth was 27°C. The chromosome analysis revealed a modal diploid karyotype of 2n = 38. The origin of this cell line was confirmed by the 18S recombinant (r)RNA sequencing. The susceptibility test indicated significant cytopathic effects in the EL cells with regard to the Rana grylio virus and the Herpesvirus anguillae. The viral replication was confirmed by transmission electron microscopy and polymerase chain reaction analysis. Following poly (I:C) exposure, the expression levels of the immune-related molecules interferon regulatory factor-7 (irf7) and transforming growth factor-ß (TGF-ß) were downregulated in EL cells, whereas the expression levels of the rf3 and the cytochrome P450 (CYP450) were upregulated. All four genes were significantly upregulated following inflammation by lipopolysaccharide (LPS). These data suggested the application of EL cell line for viral identification, as well as for immunodiagnosis and pharmacological targeting.

An insect nucleoside diphosphate kinase (NDK) functions as an effector protein in wheat - Hessian fly interactions.

Like pathogens, galling insects deliver effectors into plant tissues that induce gall formation. The gall midge Mayetiola destructor, also called Hessian fly, can convert a whole wheat seedling into a gall by inducing the formation of nutritive cells at the feeding site, inhibiting wheat growth, and reprogramming metabolic pathways of the attacked plants. Here we demonstrated the identification of a secreted Hessian fly protein, the nucleoside diphosphate kinase (NDK), in infested wheat plants through liquid chromatography-tandem mass spectrometry (LC-MS/MS) and western blots. In association with the NDK presence, enzymatic activity of NDK increased significantly in wheat tissues at the feeding site. In addition, there was a sudden increase in ATP abundance at the feeding site of infested susceptible wheat seedlings 24 h following Hessian fly larval infestation. Even though a direct link between cause and effect remains to be established, our data points to the direction that Hessian fly larvae injected NDK into wheat tissues as an effector protein, which plays a role in manipulating host plants and converting the plants into galls.

Development of an in vitro model system to study the interactions between Mycobacterium marinum and teleost neutrophils.

The lack of a reliable mammalian neutrophil in vitro culture system has restricted our ability to examine their precise roles in mycobacterial infections. Previously, we developed the procedures for the isolation and culture of primary kidney-derived neutrophil-like cells from goldfish that are functionally and morphologically similar to mammalian neutrophils. The cultured primary goldfish neutrophils exhibited prolonged viability and functional effector responses. In this study, we demonstrate that when exposed to live or heat-killed Mycobacterium marinum, goldfish neutrophils increased their mRNA levels for several pro-inflammatory cytokines (il-1ß1, il-1ß2, tnfα-1, tnfα-2) and the cytokine receptors (ifngr1-1, tnfr1, tnfr2). These neutrophils also exhibited chemotaxis toward live mycobacteria, internalized the bacilli, and produced reactive oxygen intermediates (ROI) in response to pathogen exposure. The survival of intracellular mycobacteria was significantly reduced in activated neutrophils, indicating a robust killing response by these teleost granulocytes. We suggest that this goldfish primary neutrophil in vitro model system will provide important information regarding neutrophil-mediated host defense mechanisms against mycobacteria in teleosts as well as in higher vertebrates.

Sequence heterogeneity of the ORF3 gene of porcine epidemic diarrhea viruses field samples in Fujian, China, 2010-2012.

Twenty-seven field samples that showed positive in PEDV detection were collected from different farms of Fujian province from 2010 to 2012. Their heterogeneity was investigated by analysis of the ORF3 gene because of its potential function as a representation of virulence. According to the results, six Fujian strains in Group 1 showed a different genotype with unique point mutations, which might be used in differentiation between PEDV groups and brought potential antigenic variation. P55 and five reference strains in Group 2 had a long length deletion, showing another genotype and might be involved in the variation of virulence. Phylogenetic analysis revealed that the collected Fujian strains were very distant from the vaccine development strain CV777, which might be the reason why the vaccine was inefficient to control the disease. The results can help to reconsider the strategy of PEDV vaccine management and prevent outbreaks of PEDV-induced diarrhea more efficiently.

Analysis of the immune response in infections of the goldfish (Carassius auratus L.) with Mycobacterium marinum.

The rapid doubling time and genetic relatedness of the fish pathogen Mycobacterium marinum to Mycobacterium tuberculosis has rendered the former an attractive model for investigating mycobacterial host-pathogen interactions. We employed the M. marinum-goldfish infection model to investigate the in vivo immune responses to this pathogen in the context of a natural host. Histological analysis revealed mycobacterial infiltrates in goldfish kidney and spleen tissues, peaking 28 days post infections (dpi). Quantitative gene expression analysis showed significant increases of mRNA levels of pro-inflammatory cytokines (IFNγ, IL-12p40, IL-1ß1) and cytokine receptors (IFNGR1-1, TNFR2) at 7 dpi. Conversely, the gene expression levels of key anti-inflammatory cytokines TGFß and IL-10 were elevated at 14 dpi. Furthermore, M. marinum infections markedly increased the cytokine-primed oxidative burst responses of isolated kidney phagocytes at 7 but not 56 dpi. We believe that the M. marinum-goldfish infection model will be invaluable in furthering the understanding of the mycobacterium host-pathogen interface.

Molecular characterization and phylogenetic analysis of porcine epidemic diarrhea virus (PEDV) samples from field cases in Fujian, China.

The outbreak of porcine epidemic diarrhea virus (PEDV) has been a big problem of swine industry in China in recent years. In this study, we investigated molecular diversity, phylogenetic relationships, and protein characterization of Fujian field samples with other PEDV reference strains. Sequence analysis of the S1 and sM genes showed that each sample had unique characteristics, and the sample P55 may be differentiated from the others by the unique deletions and insertions of sM gene. Phylogenetic analysis based on S1 or sM gene, which have high levels of variations, indicated that each sample was related to the specific reference strain, and this finding was consistent with the protein characterization prediction analysis. The study is useful to better understand the prevalence of PEDV and its prevention and control in Fujian.

Characterization of a late gene, ORF75 from Bombyx mori nucleopolyhedrovirus.

Open reading frame 75 (Bm-p33) of Bombyx mori nucleopolyhedrovirus (BmNPV) is a homologue of Autographa californica multiple nucleopolyhedrovirus ORF92. The gene is conserved among all baculoviruses that have been completely sequenced to date and is considered to be a baculovirus core set gene. No amino acid mutation was found in Bm-p33 sequences among six BmNPV strains differing in geography, phenotype, or host. The Bm-p33 transcript can be detected as early as 12 h post infection (h p.i.) and remains detectable until 96 h p.i. The Bm-p33 protein was detected in cell lysates from 18 h p.i. through 96 h p.i., and no positive band could be detected in budded viruses (BVs) and occlusion-derived viruses (ODVs) by western blot using anti-Bm-p33 serum. Immunofluorescence microscopy indicated that Bm-p33 accumulated in the nuclear membrane and the intranuclear region, especially near the nuclear membrane of the virus-infected cells. Bm75 RNAi significantly decreased the mRNA level. However, no obvious effects on ODV formation and BV production in BmNPV-infected cells could be detected. Bm-p33 is a BmNPV late gene encoding a nonstructural protein which may function mainly in the nucleus of the infected cells.

Characterization of an early gene orf122 from Bombyx mori nucleopolyhedrovirus.

The open reading frame 122 of Bombyx mori nucleopolyhedrovirus (BmNPV) (Bm122) has been observed to be a conserved gene in the lepidopteran baculoviruses that have been completely sequenced so far. Its transcript was detected at 3 h post infection (h p.i.) and remained detectable at up to 96 h p.i. Temporal transcription analysis indicated that Bm122 is transcribed by host RNA polymerase. The size of the translational product of the Bm122 gene in Tn5B-1-4 cells was approximately 23 kDa, which is in agreement with the predicted value of 22.9 kDa, suggesting that no major posttranslational modification occurred in the primary protein product. The subcellular localization of Bm122 was studied using EGFP-Bm122, which revealed that Bm122 protein was accumulated within the nuclear region of virus-infected BmN cells. All these results suggest that Bm122 is an early gene encoding a protein that functions in the nucleus.

Bombyx mori nucleopolyhedrovirus ORF9 is a gene involved in the budded virus production and infectivity.

The ORF9 of Bombyx mori nucleopolyhedrovirus (BmNPV) (Bm9) is conserved in all completely sequenced lepidopteran nucleopolyhedroviruses. RT-PCR analysis demonstrated that Bm9 is an early and late transcribed gene that is initiated at 3 h post-infection, and immunofluorescence microscopy showed that Bm9 is localized mainly in the cytoplasm of infected cells. To determine the role of Bm9 during virus infection, Bm9 was knocked out by recombination in a BmNPV genome propagated as a bacmid in Escherichia coli. The budded virus (BV) production of Bm9-deleted bacmids was reduced more than 10-fold compared with wild-type (wt) bacmid; however, the kinetics of viral DNA replication were unaffected. The defect in BV production was recovered by the Bm9 rescue bacmid. In addition, electron microscope observations revealed that polyhedra formation was not affected by the deletion of Bm9. Bioassays showed that the Bm9-deleted bacmid took approximately 14-22 h longer to kill fifth instar B. mori larvae than wt bacmid, and the LD(50) was about 15 times higher than that of the wt bacmid. In conclusion, Bm9 is an important but not essential factor in virus production and infectivity in vivo and in vitro.

Differentially expressed genes in resistant and susceptible Bombyx mori strains infected with a densonucleosis virus.

We investigated variations in the gene expression of Bombyx mori following infection with a densonucleosis virus (BmDNV-Z). Two B. mori near-isogenic lines, Jingsong and Jingsong.nsd-Z.NIL, which are highly susceptible and completely resistant to BmDNV-Z, respectively, were used in this study. The infection profiles of BmDNV-Z in the midguts of the B. mori Jingsong and Jingsong.nsd-Z.NIL larvae revealed that the virus invaded the midguts of both of these strains. However, its proliferation was notably inhibited in the midgut of the resistant strain. By using the suppression subtractive hybridization method, three cDNA libraries were constructed to compare BmDNV-Z responsive gene expression between the two silkworm lines. In total, 151 differentially expressed genes were obtained. Real-time qPCR analysis confirmed that 11 genes were significantly up-regulated in the midgut of the Jingsong.nsd-Z.NIL strain following BmDNV-Z infection. Our results imply that these up-regulated genes might be involved in B. mori immune responses against BmDNV infection.