Comparative Transcriptome Analysis of Bombyx mori (Lepidoptera) Larval Midgut Response to BmNPV in Susceptible and Near-Isogenic Resistant Strains

Abstract

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the primary pathogens causing severe economic losses in sericulture. However, the molecular mechanism of silkworm resistance to BmNPV remains largely unknown. Here, the recurrent parent P50 (susceptible strain) and the near-isogenic line BC9 (resistance strain) were used in a comparative transcriptome study examining the response to infection with BmNPV. A total of 14,300 unigenes were obtained from two different resistant strains; of these, 869 differentially expressed genes (DEGs) were identified after comparing the four transcriptomes. Many DEGs associated with protein metabolism, cytoskeleton, and apoptosis may be involved in the host response to BmNPV infection. Moreover, some immunity related genes were also altered following BmNPV infection. Specifically, after removing genetic background and individual immune stress response genes, 22 genes were found to be potentially involved in repressing BmNPV infection. These genes were related to transport, virus replication, intracellular innate immune, and apoptosis. Our study provided an overview of the molecular mechanism of silkworm resistance to BmNPV infection and laid a foundation for controlling BmNPV in the future.

About the Authors:

Xue-Yang Wang

Contributed equally to this work with: Xue-Yang Wang, Hai-Zhong Yu

Affiliation: School of Life Sciences, Anhui Agricultural University, Hefei, People’s Republic of China

Hai-Zhong Yu

Contributed equally to this work with: Xue-Yang Wang, Hai-Zhong Yu

Affiliation: School of Life Sciences, Anhui Agricultural University, Hefei, People’s Republic of China

Lei Geng

Affiliation: School of Life Sciences, Anhui Agricultural University, Hefei, People’s Republic of China

Jia-Ping Xu

* E-mail: [email protected]

Affiliation: School of Life Sciences, Anhui Agricultural University, Hefei, People’s Republic of China

Dong Yu

Affiliation: School of Life Sciences, Anhui Agricultural University, Hefei, People’s Republic of China

Shang-Zhi Zhang

Affiliation: School of Life Sciences, Anhui Agricultural University, Hefei, People’s Republic of China

Yan Ma

Affiliation: School of Life Sciences, Anhui Agricultural University, Hefei, People’s Republic of China

Dong-Qiong Fei

Affiliation: School of Life Sciences, Anhui Agricultural University, Hefei, People’s Republic of China

Introduction

The silkworm, Bombyx mori L. (Lepidoptera: Bombycidae) has been domesticated for production of cocoons for more than 5000 years. Silkworm rearing and the silk industry still play an important role in China, India and many other developing countries. B. mori is also a good model for the study of insect genetics and immunology [1–4]. Bombyx mori nucleopolyhedrovirus (BmNPV) is the principal silkworm pathogen and causes serious economic losses in sericulture every year. Among numerous silkworm strains, most are susceptible to BmNPV infection, although a few strains exhibit resistance [5]. The heredity of silkworm resistance against BmNPV infection is a relatively complicated process because resistance is controlled both by major dominant genes and multiple genes of micro-effect [6].

A series of studies have made significant progress in understanding silkworm resistance against BmNPV infection. Xu et al. reported that Bms3a was potentially involved in resistance to BmNPV infection [7, 8]. B. mori lipase-1, serine protease-2 and alkaline trypsin protein extracted from the digestive juice of larvae midguts showed strong antiviral activity in vitro [9–11]. Using comparative proteomics, arginine kinase was found to be involved in the antiviral process of different resistant strains of silkworm [12]. In our laboratory, a total of 12 proteins that are potentially involved in viral infection were identified using one- and two-dimensional electrophoresis followed by virus overlay assays. These proteins could be categorized into the following groups: endocytosis, intracellular transportation, and host responses [13]. Immune responses were found to be synergistically regulated by the Toll, Janus kinase/signal transducers and activators of the transcription (JAK/STAT) and immune deficiency (IMD) pathways, which could act as an important defense against exogenous pathogenic infection in conjunction with subsistent pathogen recognition receptors and response proteins [14–18]. However, the molecular mechanisms of silkworm resistance to BmNPV infection are still not fully elucidated.

In recent years, the high-throughput nature of next generation sequencing (NGS), using platforms such as Illumina HiSeq^TM 2500 have provided fascinating opportunities in the life sciences and dramatically improved the efficiency and speed of gene discovery, especially in the research of host cell responses to exogenous pathogenic infection [19]. For example, Hu et al. obtained numerous differentially expressed genes (DEGs) involved in metabolism, immunity, and inflammatory responses in Microtus fortis following infection with Schistosoma japonicum based on comparative transcriptome analysis [20]. Diege et al. examined different fish tissues infected with salmon anemia virus (ISAV) using high-throughput transcriptomics and found a strong correlation between functional modules and viral-segment transcription [16]. NGS technology was also used to explore the molecular mechanism of silkworm resistance against exogenous pathogens. Kolliopoulou et al. reported that several genes related to physical barriers, immune response, proteolytic/metabolic enzymes, heat-shock proteins, and hormonal signaling were possibly involved in silkworm resistance against Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) infection; although these genes might be induced by the virus in order to increase infectivity [21]. Additionally, several candidate genes, such as BmEts, BmToll10-3 and Hsp20-1, have been identified in the initial stage of BmNPV infection by analyzing the global transcriptional profile of silkworm cell lines and heads following BmNPV infection [22, 23].

In order to gain a global view of the molecular changes in silkworms during BmNPV infection, we selected near-isogenic line BC9 and recurrent parent P50 for transcriptome sequencing. Through comparative analysis of the transcriptomes from these two strains, a total of 869 DEGs were obtained, which included many genes potentially related to BmNPV-resistance. Our results may provide some reliable evidence to clarify the BmNPV-resistance molecular mechanism in silkworm.

Materials and Methods

Virus and Silkworm

BmNPV (T3 strain) was maintained in the Key Laboratory of Sericulture, Anhui Agricultural University, Hefei, China. The virus was obtained from the hemolymph of infected larvae and purified by repeated and differential centrifugation according to the protocol developed by Rahman et al. [24]. The concentration of the virus (OB/mL) was determined by hemocytometer.

The recurrent parent P50 (susceptible strain), the donor parent A35 (resistant strain) and the near-isogenic line BC9 were maintained in our laboratory. The near-isogenic line was constructed according to the protocol used by Yao et al. [6]. In brief, the recurrent parent P50 was crossed to the donor parent A35; progeny were repeatedly backcrossed with the recurrent parent for nine generations and each progeny was screened by BmNPV.

The first three instar larvae were reared on a fresh artificial diet at 26±1°C, 75±5% relative humidity, and a 12 hours day/night cycle. The rearing temperature for the last two instars was reduced to 24±1°C; other conditions were unchanged. On the first day of fifth instar, all larvae were starved for 24 hours and then fed with 5 μL BmNPV suspended in sterile water (1.0×10⁵ OB/mL) per larva orally; the control group was treated with sterile water. BmNPV occlusion bodies (OB) began fast proliferation at approximately 24 hours post inoculation (hpi) [25]; thus, this time was considered optimal for sample collection. Silkworm larvae were dissected and the midgut tissues were removed and then washed in PBS (137 mM NaCl, 2.7 mM KCl, 4.3 mM Na₂HPO₄, and 1.4 mM KH₂PO₄, pH 7.4) prepared with diethy pyrocarbonate (DEPC) (Sangon, China) treated H₂O. Thirty larvae midguts were mixed together to minimize individual genetic differences. Samples were flash-frozen in liquid nitrogen and pulverized, and 100 mg of sample were added directly into a RNAase free microcentrifuge tube containing 1.0 mL TRIzol Reagent (Invitrogen, USA) and stored at -80°C for later use.

Silkworm strain resistance level bioassays

The level of silkworm resistance to BmNPV was tested following the protocol developed by Cheng et al. [26]. The fourth instar larvae were inoculated with BmNPV at different concentrations; inoculations were conducted in triplicate. The level of silkworm resistance was calculated using IBM SPSS Statistics 20 (IBM, USA).

RNA extraction

The silkworm midguts dissolved in TRIzol Reagent were homogenized. Total RNA of midguts were extracted according to the manufacturer’s protocol. Concentrations were quantified using a NanoDrop 2000 spectrophotometer (Thermo Scientific, USA). The purity of all RNA samples were assessed at an absorbance ratio of A_260/280 and A_260/230, and the integrity of the RNA was confirmed by 1% agarose gel electrophoresis.

Library construction, Illumina sequencing and read assembly

Fragment interruption, cDNA synthesis, addition of adapters, PCR amplification and RNA-Seq were performed by Beijing BioMarker Technologies (Beijing, China). The standard Illumina methods and protocols were adopted to prepare and sequence the cDNA libraries. NEBNext Poly(A) mRNA Magnetic Isolation Module (NEB, USA) was used to enrich mRNA, and then the cDNA library was constructed using the NEBNext mRNA Library Prep Master Mix Set for Illumina (NEB, USA) and NEBNext Multiplex Oligos for Illumina (NEB, USA). The size of the library insert fragments was determined by 1.8% agarose gel electrophoresis, and the fragments were quantified using a Library Quantification Kit/Illumina GA Universal (Kapa, USA). Suitable fragments were selected as templates and sequenced on an Illumina HiSeq^TM 2500 using paired-end technology. Three biological replicates were used to minimize sample differences.

In order to obtain clean and high-quality reads for sequence assembly, the raw reads were filtered by removing adaptor sequences, low-quality sequences (reads with ambiguous bases ‘N’) and reads with 10% > Q < 20% bases [27, 28]. The Trinity assemble program was used to assemble the clean reads into contigs, which covered more full-length transcripts through a broad range of expression levels [29]. The resultant contigs were added to transcripts based on paired-end information. The longest transcript from alternative splicing transcripts was selected as the unigene. These unigenes were combined to produce the final assembly and used for annotation.

Functional annotation

To annotate unigenes, different sequences were searched by BLASTx against the NCBI non-redundant protein (nr) database and other databases, including Swiss-Prot protein database, the Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Cluster of Orthologous Groups (COG) database. Gene Ontology (GO) annotations, including molecular functions, biological processes, and cellular components, were obtained using the Blast2GO program (https://www.blast2go.com/) [30, 31]. All searches were performed with an E-value < 10⁻⁵. Fragments per kilobase of transcript per million fragments mapped (FPKM) was calculated to represent the expression abundance of the unigenes [32]. FPKM may reflect the molar concentration of a transcript by normalizing for RNA length and for the total read number.

Differentially expressed genes (DEGs) analysis

After normalizing genes expression levels, DEGs were obtained by pair-wise comparison of the four transcriptome libraries using IDEG6 software [33]. An FPKM fold change of ≥ 1.2 or ≤ 0.83 between two libraries was defined as the reference standard, with the Benjamini–Hochberg false discovery rate (FDR < 0.01) used to adjust the p-values.

Real-time quantitative PCR analysis

In order to validate the results from our transcriptome sequencing analysis, the relative expression levels of 15 randomly selected genes were confirmed by reverse transcription quantitative PCR (RT-qPCR). Additionally, 9 genes with well reported previously were selected to further validate the genes of interest that might be involved in BmNPV resistance. All the Primers are listed in Table 1. RT-qPCR reactions were prepared with the SYBR Premix Ex Taq^TM Kit (TaKaRa), following the manufacturer’s instruction. Reactions were carried out in Bio-Rad CFX96TM Real-Time System (Bio-Rad, USA). The thermal cycling profile consisted of an initial denaturation at 95°C for 30 s, 40 cycles at 95°C for 5 s, and 60°C for 30 s. All samples were performed in triplicate. Relative expression levels were calculated using the 2^-△△Ct method following the protocol of Livak et al. [34]. In this study, B. mori ribosomal protein s3 (BmRPS3) gene was used as a reference gene. Statistical analysis was conducted using the SPSS software (IBM, www.ibm.com).

[Figure omitted. See PDF.]

Table 1. Primers used in RT-qPCR for validation of DEGs.

https://doi.org/10.1371/journal.pone.0155341.t001

Results

BmNPV infectivity in different silkworm strains

The LC₅₀ value was used to evaluate the resistant level of silkworm to BmNPV infection. The LC₅₀ value of A35 was approximately 26-fold greater than that of BC9 and over 500-fold greater than that of P50. The value of BC9 was 23-fold greater than that of P50 (Table 2).

[Figure omitted. See PDF.]

Table 2. The LC₅₀ value of B. mori larvae infected with BmNPV.

https://doi.org/10.1371/journal.pone.0155341.t002

Overview of the silkworm transcriptome

Transcriptome sequencing is an efficient technology for comparing gene expression levels in different samples, and in our study, it was used to search and analyze DEGs among P50, BC9 control and treatment groups. A total of four cDNA libraries were sequenced: P50- (P50 treated with sterile water), P50+ (P50 infected with BmNPV), BC9- (BC9 treated with sterile water), BC9+ (BC9 infected with BmNPV), with each group created in triplicate. After removing the adaptors and low quality sequences, 144,439,382 sequence reads were obtained (Table 3). The GC content of each of the four libraries was approximately 50%, and CycleQ30% was greater than 89.91% for each library. Thus, the quality and accuracy of the sequencing data were sufficient for further analysis. Most of the reads matched silkworm genomic locations. All the unigenes matched previously described sequences with approximately 70% coverage. The length distribution of unigenes had similar patterns among the four libraries, suggesting that there was little bias in the construction of the four cDNA libraries (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Length distribution of unigenes in the assembled transcriptomes.

The x axis represents different groups and treatments and the y axis shows the number of unigenes.

https://doi.org/10.1371/journal.pone.0155341.g001

[Figure omitted. See PDF.]

Table 3. Summary statistics for silkworm genes based on the RNA-seq data.

https://doi.org/10.1371/journal.pone.0155341.t003

Unigenes annotation and classification

In order to annotate the unigenes, reference sequences were searched using BLASTX against the NCBI nr database (E-value < 10⁻⁵). A total of 12,591 of 13,342 unigenes provided a BLAST result (S1 Table). S2 Table shows the species with the closest match for each unigene. Most of the annotated sequences had the greatest homology with other B. mori sequences (87%), followed by Danaus plexippus (9%).

RT-qPCR validation of differentially expressed transcripts

In order to determine the reliability of the transcriptome sequencing, the relative expression levels of 15 randomly selected genes were analyzed by RT-qPCR (Fig 2A). The results were consistent with the transcriptome data. For example, the gene encoding peptidoglycan-recognition protein was down-regulated in both RNA-seq and RT-qPCR analyses, with a similar fold change. The lipase 1 gene was significantly up-regulated in the resistant strain BC9 after infection with BmNPV, which was also consistent with previously reported results [10]. Linear regression analysis of the correlation between RT-qPCR and RAN-seq (Fig 2B) showed an R² (goodness of fit) value of 0.9169 and a corresponding slope of 1.5281, suggesting a strong positive correlation between RT-qPCR and transcriptome data. Therefore, the transcriptome data were satisfied for further analysis.

[Figure omitted. See PDF.]

Fig 2. Correlation between gene expression ratios obtained from transcriptome data and RT-qPCR.

(A) Expression ratios (FPKM fold change) obtained from transcriptome data (red) and RT-qPCR (blue). (B) Lineage analysis between the transcriptome and RT-qPCR. The ratios obtained by RT-qPCR (x-axis) were plotted against the ratios obtained by RNA-Seq (y-axis).

https://doi.org/10.1371/journal.pone.0155341.g002

Differentially expressed genes (DEGs) and their possible roles in host response to BmNPV

To further elucidate which DEGs had a potential role in antiviral response, a Venn diagram was constructed. A total of 285 DEGs were found to be differentially regulated when comparing P50+ and P50-, of which 122 were up-regulated and 163 were down-regulated. Similarly, 193 DEGs were found to be differentially regulated in the comparison of BC9+ vs. BC9-, with 56 genes up-regulated and 137 down-regulated. In addition, there were 154 DEGs differentially regulated in the comparison of BC9- and P50-, among which 78 genes were up-regulated and 76 genes were down-regulated (Fig 3, Table 4, S3 Table).

[Figure omitted. See PDF.]

Fig 3. Venn diagram showing the DEGs related to BmNPV infection in different resistant strains.

https://doi.org/10.1371/journal.pone.0155341.g003

[Figure omitted. See PDF.]

Table 4. The number of up-regulated, down-regulated and unique DEGs of the pairwise comparisons of the experimental groups.

Shared numbers of DEGs are also shown in Fig 3.

https://doi.org/10.1371/journal.pone.0155341.t004

There were 197, 119, 82 unique DEGs in P50+ vs. P50-, BC9+ vs. BC9- and BC9- vs. P50-, respectively. GO assignments were used to assign a functional classification to these DEGs. For cellular components, the number of unique DEGs fell into the macromolecular complex classification was distinct in BC9 and P50 following BmNPV infection. For molecular functions, the number of transporter activity related unique DEGs was distinct in BC9 and P50 following BmNPV infection. For biological progresses, the number of unique DEGs involved in metabolism processes and localization was distinct in BC9 and P50 following BmNPV infection (Fig 4).

[Figure omitted. See PDF.]

Fig 4. Gene ontology (GO) analysis of DEGs in different comparable groups.

These genes were divided into groups based on cellular component, molecular function, and biological process. The percent of DEGs that could be assigned to the different categories are indicated.

https://doi.org/10.1371/journal.pone.0155341.g004

Analysis of DEGs associated with protein metabolism, cytoskeleton, and apoptosis

The comparisons of BC9+ vs. BC9- and P50+ vs. P50- identified many DEGs that might either be involved in silkworm defense against BmNPV or facilitate BmNPV infection. These genes could be divided into three categories: protein metabolism, cytoskeleton, and apoptosis.

Most of the DEGs (70%) associated with protein metabolism were down-regulated in BC9 following BmNPV infection. In contrast, 80% of the genes in P50 were up-regulated after BmNPV infection. Most of the DEGs (87.5%) associated with the cytoskeleton were up-regulated in P50 after BmNPV infection, but the number of up-regulated genes (50%) decreased in BC9. A total of 19 DEGs associated with apoptosis were identified following BmNPV infection. Ten of these DEGs were altered in BC9 after BmNPV infection. The other DEGs were all over-expressed in P50 following BmNPV infection (Table 5).

[Figure omitted. See PDF.]

Table 5. DEGs involved in protein metabolism, cytoskeleton, and apoptosis after BmNPV infection in different resistant strains.

https://doi.org/10.1371/journal.pone.0155341.t005

Alteration in immunity related gene expression after BmNPV infection in different resistant strains

Pathogen infection stimulated both cellular and humoral responses of insects [33–36]. Genes participating in innate immunity pathways were identified and analyzed in regards to their potential role in BmNPV infection (Table 6). Thirty DEGs were identified, which could be classified into Toll pathway, Imd pathway, polyphenol oxidase (PPO) pathway, pattern recognition receptor, and antimicrobial peptide. As shown in Table 6, 14 (47%) of these genes were down-regulated and 9 (30%) were up-regulated in BC9 after infection with BmNPV, while 17 (57%) were down-regulated and 6 (20%) were up-regulated in P50 after BmNPV infection.

[Figure omitted. See PDF.]

Table 6. Relative FPKMs of innate immune-related genes in silkworm after BmNPV infection in different resistant strains.

https://doi.org/10.1371/journal.pone.0155341.t006

Role of selected DEGs in BmNPV resistance

After removing the genetic background and strain specific immune stress response genes, a total of 22 DEGs were identified as possibly being involved in silkworm resistance to BmNPV (Fig 3, Table 7). For BC9- vs. P50-, 13 genes were up-regulated and the rest were down-regulated. However, the 22 genes were all down-regulated in BC9 following BmNPV infection (Fig 5). Some genes, including prostatic acid phosphatase, protease inhibitor 6, actin cytoskeleton-regulatory complex protein PAN1, and EF-hand domain-containing protein, were further analyzed by RT-qPCR (Fig 6). After BmNPV infection, the expression levels of 4 genes were down-regulation in BC9 and A35 (resistant strain) (Fig 6), which was consistent with the transcriptome data. Thus, we deduced that the 22 genes were possibly involved in resistance to BmNPV infection. Gene functions fell into the following categories: transport, virus replication, intracellular innate immunity, and apoptosis.

[Figure omitted. See PDF.]

Fig 5. Expression patterns of selected genes related to BmNPV resistance in different resistant strains.

Each row represents a different gene, with blue, yellow and red indicating low, medium and high levels of gene expression, respectively.

https://doi.org/10.1371/journal.pone.0155341.g005

[Figure omitted. See PDF.]

Fig 6. Real-time PCR analysis of expression profiles of several genes in B. mori midgut.

After removing genetic background and individual immune stress response genes, 22 differentially expressed genes of interest potentially involved in resistance to NPV infection were obtained. Additionally, 119 genes unique differentially expressed in the isogenic-line BC9 (resistant strain) following BmNPV infection were observed. In order to further conformed the function of these genes, 4 genes from 22 DEGs and 5 genes from 119 DEGs with well reported previously were selected to conduct RT-qPCR.

https://doi.org/10.1371/journal.pone.0155341.g006

[Figure omitted. See PDF.]

Table 7. Relative FPKMs of 22 DEGs of interest in BC9 following BmNPV infection.

https://doi.org/10.1371/journal.pone.0155341.t007

Discussion

Despite the confirmation of an association between many genes and proteins and resistance to BmNPV, the molecular mechanism of antiviral activities was still unclear. Here, transcriptome sequencing was carried out to identify genes related to BmNPV-resistance in silkworm across the genome. By using the near-isogenic line BC9 (resistant strain) and the recurrent parent P50 (susceptible strain) to study silkworm antiviral mechanisms, some DEGs responding to BmNPV infection were successfully identified after comparing infected groups and controls in the two strains.

Overview of specific DEGs in two strains following BmNPV infection

Based on the GO analysis, more DEGs were found to be involved in metabolic processes in BC9 (17.8%) than in P50 (9.2%) following BmNPV infection (Fig 4), which was consistent with previous reports for Sogatella furcifera (Hemiptera: Delphacidae) and Campoletis sonorensis (Hymenoptera: Ichneumonidae) [35, 36]. Approximately twice as many DEGS related to transporter activity were identified in BC9 (12.2%) than in P50 (6.3%) (Fig 4). Down-regulation of transporter related genes, such as lactase-phlorizin hydrolase (LPH), B(0,+)-type amino acid transporter 1 (BAT1), actin cytoskeleton-regulatory complex protein PAN1 (PAN1), MFS-type transporter (MFS), and otoferlin, could repress virus infection in host cells [37–42]. Therefore, the increased number of these genes with altered expression levels in BC9 might be related to BmNPV infection. Moreover, the number of macromolecular complex genes in BC9 (14%) following BmNPV infection that showed an increase in expression compared with P50 (5.1%) were similar to previous reports [43, 44].

Protein metabolism, cytoskeleton, and apoptosis may play an important role in host response to BmNPV infection

Cellular and metabolic processes will be dramatically changed after viral infection [35]. In our study, several DEGs that participate in protein metabolism were found to be of interest. For example, solute carrier family 12 is involved in transporting extraordinarily diverse solutes [45], and cystathionine gamma-lyase participates in hydrolysis of cystathionine [46]. Viruses may have to rely on cell proteins to accomplish replication in intercellular regions [47], therefore, the down-regulation of the genes involved in protein metabolism could inhibit the replication of BmNPV in the host cell. We speculated that the down-regulation of these genes affected virus replication.

Cytoskeleton-dependent intracellular transport is an important strategy for transport of viral particles to different destinations [48]. In this study, some cytoskeleton related genes were found to be of interest, including actin cytoskeleton-regulatory complex protein PAN1 and actin-binding protein. These genes related to actin-coupled endocytosis could promote viral transport [49, 50]. We speculated that the down-regulation of the cytoskeleton genes might affect BmNPV transport.

Apoptosis plays a vital role in regulating cell response in Lepidopteran insects during viral infections, where larvae use selective apoptosis and subsequent sloughing of the infected cells in the midgut epithelium to resist virus infection [51, 52]. In this study, some genes related to activation of apoptosis were found to be of interest, including cytochrome c, inhibitor of caspase-activated DNase, amyloid precursor protein and B-cell lymphoma protein 2. The activity of caspase-activated DNase was blocked by Hepatitis C virus core at physiological levels, resulting in the inhibition of apoptotic cell death [53]. Amyloid precursor protein is a member of several signaling pathways that are involved in abnormal cell cycles, subsequently leading to apoptosis [54]. B-cell lymphoma protein 2 could bind to BH3 domains of various pro-apoptotic regulators to activate apoptosis [55]. The overexpression of cytochrome c in rabies virus showed a decreased pathogenicity in vitro and in vivo [56]. Based on their role in apoptosis activation, hosts need to increase the expression level of these genes to promote apoptosis when exposed to a virus; this supposition explains the up-regulation of genes involved in apoptosis in the transcriptome following BmNPV infection. Additionally, a significantly higher relative transcriptional level of cytochrome c was detected by RT-qPCR in BC9 and A35 (Fig 6), which indicated the up-regulation of this gene could activate apoptosis to repress further BmNPV infection. However, some genes were down-regulated, including cysteine aspartic acid specific protease 9L, protein kinase A, apoptosis-inducing factor, apoptosis signal-regulating kinase 1, TNF-receptor-associated factor 6, TGF-beta-activated kinase 1, and p90 ribosomal S6 kinase. The inhibition of these genes may strongly impair viral infectivity and virus-induced apoptosis [57–62].

Changes to immune gene expression following infection with BmNPV

Although some differential expression of immune genes was observed, this might not be considered biologically important due to low expression levels; low expression levels were also found after BmCPV infection [21]. In this study, most immune related genes were down-regulated, with a few exceptions, such as the 2-fold up-regulation of toll-like protein and lysozyme. The expression of several other genes, including 18 wheeler, macrophage mannose receptor 1, and slit homolog 3 protein, were also up-regulated during the BmNPV infection. Unfortunately, the relationship between immune genes and BmNPV remains unclear and requires further study. We presumed that the low expression levels of immunity related genes may be associated with the disruption of the immune system by BmNPV, similar to the pathogenicity of human immunodeficiency virus (HIV).

Multiple genes have potential roles in repressing BmNPV infection

Based on the Venn diagram, 22 genes of interest were identified, all of which were down-regulated in BC9 following BmNPV infection. These genes were grouped based on their functions, as reported in the literature. These groups, including transport, virus replication, intracellular innate immune, and apoptosis, may play an important role in the process of silkworm resistance to BmNPV.

Several genes related to virus transport, including lactase-phlorizin hydrolase (LPH) [37], B(0,+)-type amino acid transporter 1 (BAT1) [38, 46], actin cytoskeleton-regulatory complex protein PAN1 (PAN1) [63], and otoferlin [41, 42] were identified in this study. The expression level of these genes were all down-regulated in BC9 following BmNPV infection (Fig 5). The resistant strain A35, a donor parent, was used to validate our results. The expression level of PAN1 was chosen for further testing by RT-qPCR. Expression levels of PAN1 were down-regulated in both BC9 and A35 following BmNPV infection (Fig 6). The down-regulation of virus transport-related genes could inhibit the transmembrane and intracellular transport of BmNPV thereby preventing further infection. MFS-type transporter (MFS) [39, 40], a transmembrane facilitator, was induced to increase intracellular monovalent ion concentrations, which led to lysis and cell death after HIV-1 infection. The down-regulation of this gene in silkworm could potentially block BmNPV infection.

Several genes related to virus replication were found to be of interest, including engrailed nuclear homeoprotein-regulated protein [64], 1-acyl-sn-glycerol-3-phosphate acyltransferase alpha (ASGPA) [65, 66], alanine aminotransferase 2 (ALT2) [67, 68], U2AF domain protein (U2AF) [69, 70], adenylate cyclase type 5 (ACT5) [71], EF-hand domain-containing protein (EFHP) [72, 73], prostatic acid phosphatase (PAP) [74], and zinc ribbon domain protein (ZRDP) [75]. In this study, the expression level of these genes were all down-regulated in the transcriptome of BC9 following BmNPV infection (Fig 5). In order to validate the results, RT-qPCR was conducted as described above. The expression levels of EFHP and PAP were lower following BmNPV infection in both BC9 and A35 (Fig 6). The down-regulation of these genes could directly or indirectly participate in repressing BmNPV replication in host cells.

Of the genes related to apoptosis, the expression level of pyruvate dehydrogenase kinase (PDK) was obviously inhibited after treatment with apoptosis-inducing agents [76]; such down-regulation might activate apoptosis in response to BmNPV infection.

Protease inhibitor 6 is a member of the protein superfamily that contains TIL functional domain. Zhao et al. used genome sequences to demonstrate that the expression level of the TIL superfamily were down-regulated following BmNPV infection [77], which was consistent with our results (Fig 6). Furthermore, peptidoglycan-recognition protein 2, hemolin, facilitated trehalose transporter Tret1, integument esterase 2, defense protein precursor, and antimicrobial protein 6Tox precursor also showed differential expression following BmNPV infection. The relationship of these genes to BmNPV is an important area for further study.

Four DEGs were down-regulated in BC9 following BmNPV infection, including ras-responsive element-binding protein 1, GTPase-activating protein, trypsin alkaline A and peroxidase (Fig 5). Previous studies revealed that these proteins play a role in virus infections. Mdv1-miR-M4 encoded by Marek's disease virus efficiently targeted the 3' untranslated regions of ras-responsive element-binding protein 1 (RREB1) [78]. TBC domain proteins belonging to a GTPase-activating protein were knocked out by double stranded RNA interference (RNAi), which led to a decrease in the level of transcripts of white spot syndrome virus genes [79]. Trypsin in the myocardium was able to trigger acute myocarditis following influenza A virus infection [80]. Overexpression and RNA silencing studies revealed that peroxidase was involved in the production of hepatitis C virus particles [81]. Moreover, RT-qPCR results indicated that the expression levels of all genes were down-regulated in BC9 and A35 following BmNPV infection (Fig 6). Therefore, we speculated that the down-regulation of these genes might be involved resistance to BmNPV infection.

Hypothesized modal analysis of the roles of the screened DEGs in silkworm resistance to BmNPV infection pathway

We hypothesized that the 22 DEGs discussed above played a role in the process of host response to BmNPV infection. The endocytosis process is triggered when the BmNPV nucleocapsid containing envelope binds to the cytomembrane. Vacuolar ATP synthase is activated by LPH to promote the fusion of the envelope and endosome thereby releasing the nucleocapsid into the cytoplasm. This process can be promoted by PAN1 and otoferlin. However, the transmembrane transport channel is an alternative pathway for virus to enter the cytoplasm, a process which can be facilitated by BAT1. The released nucleocapsid is transported into the nucleus with the help of the cytoskeleton (EFP). Once viral DNA is released into the nucleus, it will utilize host nucleotides to complete replication. In the final step of replication, viral DNA has to rely on host cell amino acids for assembly on the cytoskeleton (Table 5) [82]. In the cytoplasm, EFHP, ASGPA, ALT2, U2AF, ACT5 and ZRDP play an important role in facilitating virus replication, although the exact mechanism is still unclear. Moreover, transmembrane protein, MFS, is induced by BmNPV to increase cell volume, leading to lysis and cell death. We speculated that the down-regulation of these genes may affect the entry of virus into host cells and virus replication. The apoptosis process could also be triggered by PDK to inhibit BmNPV from further infecting other cells once BmNPV entered a host cell (Fig 7).

[Figure omitted. See PDF.]

Fig 7. Hypothesized modal analysis of the roles of the screened DEGs in BmNPV infection pathway.

V-ATPase is activated by LPH to promote viral entry into the cytoplasm, a process which was also effected by PAN1 and otoferlin. BAT1 related channel could serve as an alternative pathway for virus transmembrane transport. The released nucleocapsid is transported into the nucleus with the help of EFP. During replication, EFHP, ASGPA, ALT2, U2AF, ACT5 and ZRDP play an important role in facilitating virus replication. MFS is induced by the virus to increase cell volume leading to cell death. At the same time, the apoptosis process could be triggered by PDK to inhibit BmNPV infection.

https://doi.org/10.1371/journal.pone.0155341.g007

Taken together, our results provide useful information on silkworm resistance to BmNPV infection.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. The annotation results for all unigenes.

https://doi.org/10.1371/journal.pone.0155341.s001

(XLSX)

S2 Table. The species distribution for the BLAST matches.

https://doi.org/10.1371/journal.pone.0155341.s002

(XLSX)

S3 Table. Genes that were differentially expressed in P50+ vs. P50-, BC9+ vs. BC9-, and BC9- vs. P50-.

https://doi.org/10.1371/journal.pone.0155341.s003

(XLSX)

Author Contributions

Conceived and designed the experiments: JPX. Performed the experiments: XYW LG DY SZZ YM DQF. Analyzed the data: XYW HZY JPX. Contributed reagents/materials/analysis tools: JPX. Wrote the paper: XYW HZY JPX.

Citation: Wang X-Y, Yu H-Z, Geng L, Xu J-P, Yu D, Zhang S-Z, et al. (2016) Comparative Transcriptome Analysis of Bombyx mori (Lepidoptera) Larval Midgut Response to BmNPV in Susceptible and Near-Isogenic Resistant Strains. PLoS ONE 11(5): e0155341. https://doi.org/10.1371/journal.pone.0155341

References

1. Goldsmith MR, Shimada T, Abe H. The genetics and genomics of the silkworm, Bombyx mori. Annu Rev Entomol. 2005; 50:71–100. pmid:15355234

2. Nagaraju J, Goldsmith MR. Silkworm genomics—progress and prospects. Curr Sci India. 2002; 83(4):415–25.

3. Li MW, Shen L, Xu AY, Miao XX, Hou CX, Sun PJ, et al. Genetic diversity among silkworm (Bombyx mori L., Lep., Bombycidae) germplasms revealed by microsatellites. Genome. 2005; 48(5):802–10. pmid:16391686

4. Tanaka H, Sagisaka A, Fujita K, Kaneko Y, Imanishi S, Yamakawa M. Lipopolysaccharide elicits expression of immune-related genes in the silkworm, Bombyx mori. Insect Mol Biol. 2009; 18(1):71–5. pmid:19196348

5. Bao YY, Tang XD, Lv ZY, Wang XY, Tian CH, Xu YP, et al. Gene expression profiling of resistant and susceptible Bombyx mori strains reveals nucleopolyhedrovirus-associated variations in host gene transcript levels. Genomics. 2009; 94(2):138–45. pmid:19389468

6. Yao Q, Li MW, Wang Y, Wang WB, Lu J, Dong Y, et al. Screening of molecular markers for NPV resistance in Bombyx mori L. (Lep., Bombycidae). J. Appl. Entomol. 2003; 127(3):134–6.

7. Xu JP, Chen KP, Yao Q, Liu XY. Fluorescent differential display analysis of gene expression for NPV resistance in Bombyx mori L. J. Appl. Entomol. 2005; 129(1):27–31.

8. Xu JP, Chen KP, Yao Q, Liu MH, Gao GT, Zhao Y. Identification and characterization of an NPV infection-related gene Bmsop2 in Bombyx mori L. J. Appl. Entomol. 2005; 129(8):425–31.

9. Nakazawa H, Tsuneishi E, Ponnuvel KM, Furukawa S, Asaoka A, Tanaka H, et al. Antiviral activity of a serine protease from the digestive juice of Bombyx mori larvae against nucleopolyhedrovirus. Virology. 2004; 321(1):154–62. pmid:15033574

10. Ponnuvel KM, Nakazawa H, Furukawa S, Asaoka A, Ishibashi J, Tanaka H, et al. A lipase isolated from the silkworm Bombyx mori shows antiviral activity against nucleopolyhedrovirus. J Virol. 2003; 77(19):10725–9. pmid:12970462

11. Ponnuvel KM, Nithya K, Sirigineedi S, Awasthi AK, Yamakawa M. In Vitro Antiviral Activity of an Alkaline Trypsin from the Digestive Juice of Bombyx Mori Larvae against Nucleopolyhedrovirus. Arch Insect Biochem. 2012; 81(2):90–104. pmid:22898997

12. Kang LQ, Shi HF, Liu XY, Zhang CY, Yao Q, Wang Y, et al. Arginine kinase is highly expressed in a resistant strain of silkworm (Bombyx mori, Lepidoptera): Implication of its role in resistance to Bombyx mori nucleopolyhedrovirus. Comp Biochem Phys B. 2011; 158(3):230–4. pmid:21146627

13. Cheng Y, Wang XY, Hu H, Killiny N, Xu JP. A Hypothetical Model of Crossing Bombyx mori Nucleopolyhedrovirus through Its Host Midgut Physical Barrier. PloS One. 2014; 9(12):e115032. pmid:25502928

14. Liu W, Liu JB, Lu YH, Gong YC, Zhu M, Chen F, et al. Immune signaling pathways activated in response to different pathogenic micro-organisms in Bombyx mori. Mol Immunol. 2015; 65(2):391–7. pmid:25745806

15. Wu S, Zhang XF, He YQ, Shuai JB, Chen XM, Ling EJ. Expression of antimicrobial peptide genes in Bombyx mori gut modulated by oral bacterial infection and development. Dev Comp Immunol. 2010; 34(11):1191–8. pmid:20600274

16. Valenzuela-Miranda D, Boltana S, Cabrejos ME, Yanez JM, Gallardo-Escarate C. High-throughput transcriptome analysis of ISAV-infected Atlantic salmon Salmo salar unravels divergent immune responses associated to head-kidney, liver and gills tissues. Fish Shellfish Immunol. 2015; 45(2):367–77. pmid:25910847

17. Wu K, Yang B, Huang W, Dobens L, Song H, Ling E. Gut immunity in Lepidopteran insects. Dev Comp Immunol. 2016. pmid:26872544.

18. Liu F, Xu Q, Zhang Q, Lu A, Beerntsen BT, Ling E. Hemocytes and hematopoiesis in the silkworm, Bombyx mori. Isj-Invert Surviv J. 2013; 10(1):102–9.

19. Ansorge WJ. Next-generation DNA sequencing techniques. New Biotechnol. 2009; 25(4):195–203.

20. Hu Y, Xu YX, Lu WY, Yuan ZY, Quan H, Shen YJ, et al. De novo assembly and transcriptome characterization: novel insights into the natural resistance mechanisms of Microtus fortis against Schistosoma japonicum. BMC genomics. 2014, 15(1):1–13. pmid:24886088

21. Kolliopoulou A, Van Nieuwerburgh F, Stravopodis DJ, Deforce D, Swevers L, Smagghe G. Transcriptome analysis of Bombyx mori larval midgut during persistent and pathogenic cytoplasmic polyhedrosis virus infection. PLoS One. 2015; 10(3):e0121447. pmid:25816294

22. Sagisaka A, Fujita K, Nakamura Y, Ishibashi J, Noda H, Imanishi S, et al. Genome-wide analysis of host gene expression in the silkworm cells infected with Bombyx mori nucleopolyhedrovirus. Virus Res. 2010; 147(2):166–75. pmid:19883703.

23. Wang GB, Zhang JJ, Shen YW, Zheng Q, Feng M, Xiang XW, et al. Transcriptome analysis of the brain of the silkworm Bombyx mori infected with Bombyx mori nucleopolyhedrovirus: A new insight into the molecular mechanism of enhanced locomotor activity induced by viral infection. J Invertebr Pathol. 2015; 128:37–43. pmid:25912089

24. Rahman MM, Gopinathan KP. Systemic and in vitro infection process of Bombyx mori nucleopolyhedrovirus. Virus Res. 2004; 101(2):109–18. pmid:15041178

25. Yao Q, Gao L, Cheng KP, Hu ZG. Detection of proliferation of Bombyx mori nucleopolyhedrovirus in its host by fluorescence quantitative PCR. Acta Entomologica Sinica. 2005; 48(6):871–5.

26. Cheng Y, Wang XY, Du C, Gao J, Xu JP. Expression analysis of several antiviral related genes to BmNPV in different resistant strains of silkworm, Bombyx mori. J Insect Sci. 2014; 14(2):191–7. pmid:25373223

27. Xie FL, Burklew CE, Yang YF, Liu M, Xiao P, Zhang BH, et al. De novo sequencing and a comprehensive analysis of purple sweet potato (Impomoea batatas L.) transcriptome. Planta. 2012; 236(1):101–13. pmid:22270559

28. Wang ZY, Fang BP, Chen JY, Zhang XJ, Luo ZX, Huang LF, et al. De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweetpotato (Ipomoea batatas). BMC genomics. 2010; 11(53):726. pmid:21182800

29. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011; 29(7):644–52. pmid:21572440

30. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005; 21(18):3674–6. pmid:16081474

31. Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008; 36(10):3420–35. pmid:18445632

32. Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008; 5(7):621–8. pmid:18516045

33. Romualdi C, Bortoluzzi S, D'Alessi F, Danieli GA. IDEG6: a web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol Genomics. 2003; 12(2):159–62. pmid:12429865

34. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods. 2001; 25(4):402–8. pmid:11846609

35. Xu Y, Zhou WW, Zhou YJ, Wu JX, Zhou XP. Transcriptome and Comparative Gene Expression Analysis of Sogatella furcifera (Horvath) in Response to Southern Rice Black-Streaked Dwarf Virus. PLoS One. 2012; 7(4):e36238. pmid:22558400

36. Luan JB, Li JM, Varela N, Wang YL, Li FF, Bao YY, et al. Global Analysis of the Transcriptional Response of Whitefly to Tomato Yellow Leaf Curl China Virus Reveals the Relationship of Coevolved Adaptations. J Virol. 2011; 85(7):3330–40. pmid:21270146

37. Zhou J, Chen Y, Wang Y, Gao X, Qu D, Liu CY. A Comparative Study on the Metabolism of Epimedium koreanum Nakai-Prenylated Flavonoids in Rats by an Intestinal Enzyme (Lactase Phlorizin Hydrolase) and Intestinal Flora. Molecules. 2014; 19(1):177–203. pmid:24368601

38. Rajan DP, Huang W, Kekuda R, George RL, Wang J, Conway SJ, et al. Differential influence of the 4F2 heavy chain and the protein related to b(0,+) amino acid transport on substrate affinity of the heteromeric b(0,+) amino acid transporter. J Biol Chem. 2000; 275(19):14331–5. pmid:10799513

39. Marger MD, Saier MH Jr. A major superfamily of transmembrane facilitators that catalyse uniport, symport and antiport. Trends Biochem Sci. 1993; 18(1):13–20. pmid:8438231.

40. Voss TG, Gatti PJ, Fermin CD, Garry RF. Reduction of human immunodeficiency virus production and cytopathic effects by inhibitors of the Na+/K+/2Cl- cotransporter. Virology. 1996; 219(1):291–4. pmid:8623543.

41. Johnson CP, Chapman ER. Otoferlin is a calcium sensor that directly regulates SNARE-mediated membrane fusion. Mol Biol Cell. 2010, 191(1):187–97. pmid:20921140

42. Lek A, Evesson FJ, Sutton RB, North KN, Cooper ST. Ferlins: Regulators of Vesicle Fusion for Auditory Neurotransmission, Receptor Trafficking and Membrane Repair. Traffic. 2012; 13(2):185–94. pmid:21838746

43. Sun JF, Han ZX, Shao YH, Cao ZZ, Kong XG, Liu SW. Comparative proteome analysis of tracheal tissues in response to infectious bronchitis coronavirus, Newcastle disease virus, and avian influenza virus H9 subtype virus infection. Proteomics. 2014; 14(11):1403–23. pmid:24610701

44. Zhang Y, Liu JF, Liu HL, He YL, Yi RT, Niu YH, et al. Comparative study of the different activities of hepatitis B virus whole-X protein and HBx in hepatocarcinogenesis by proteomics and bioinformatics analysis. Arch Virol. 2015; 160(7):1645–56. pmid:25913689

45. Hediger MA, Romero MF, Peng JB, Rolfs A, Takanaga H, Bruford EA. The ABCs of solute carriers: physiological, pathological and therapeutic implications of human membrane transport proteins—Introduction. Pflug Arch Eur J Phy. 2004; 447(5):465–8. pmid:14624363

46. Barry C, Key T, Haddad R, Duncan R. Features of a Spatially Constrained Cystine Loop in the p10 FAST Protein Ectodomain Define a New Class of Viral Fusion Peptides. J Biol Chem. 2010; 285(22):16424–33. pmid:20363742

47. Emmett SR, Dove B, Mahoney L, Wurm T, Hiscox JA. The cell cycle and virus infection. Methods Mol Biol. 2005; 296:197–218. pmid:15576934

48. Harries PA, Schoelz JE, Nelson RS. Intracellular Transport of Viruses and Their Components: Utilizing the Cytoskeleton and Membrane Highways. Mol Plant Microbe In. 2010; 23(11):1381–93. pmid:20653412

49. Martin R, Hellwig D, Schaub Y, Bauer J, Walther A, Wendland J. Functional analysis of Candida albicans genes whose Saccharomyces cerevisiae homologues are involved in endocytosis. Yeast. 2007; 24(6):511–22. pmid:17431925

50. Barrientos T, Frank D, Kuwahara K, Bezprozvannaya S, Pipes GCT, Bassel-Duby R, et al. Two novel members of the ABLIM protein family, ABLIM-2 and -3, associate with STARS and directly bind F-actin. J Biol Chem. 2007; 282(11):8393–403. pmid:17194709

51. Wu Y, Hui T, Wu H, Wang W. Reaper homologue IBM1 in silkworm Bombyx mori induces apoptosis upon baculovirus infection. FEBS letters. 2013; 587(6):600–6. pmid:23416298

52. Narayanan K. Insect defence: its impact on microbial control of insect pests. Curr Sci India. 2004; 86(6):800–14.

53. Sacco R, Tsutsumi T, Suzuki R, Otsuka M, Aizaki H, Sakamoto S, et al. Antiapoptotic regulation by hepatitis C virus core protein through up-regulation of inhibitor of caspase-activated DNase. Virology. 2003; 317(1):24–35. pmid:14675622

54. Neve RL, McPhie DL. The cell cycle as a therapeutic target for Alzheimer's disease. Pharmacol Therapeut. 2006; 111(1):99–113. pmid:16274748

55. Bhat V, Olenick MB, Schuchardt BJ, Mikles DC, McDonald CB, Farooq A. Biophysical basis of the promiscuous binding of B-cell lymphoma protein 2 apoptotic repressor to BH3 ligands. J Mol Recognit. 2013; 26(10):501–13. pmid:23996493

56. Pulmanausarakul R, Faber M, Morimoto K, Spitsin S, Weihe E, Hooper DC, et al. Overexpression of cytochrome c by a recombinant rabies virus attenuates pathogenicity and enhances antiviral immunity. J Virol. 2001; 75(22):10800–7. pmid:11602721

57. Liu L, Cao ZY, Chen J, Li R, Cao YH, Zhu CL, et al. Influenza A Virus Induces Interleukin-27 through Cyclooxygenase-2 and Protein Kinase A Signaling. J Biol Chem. 2012; 287(15):11899–910. pmid:22343630

58. Kim Y, Lee C. Porcine epidemic diarrhea virus induces caspase-independent apoptosis through activation of mitochondrial apoptosis-inducing factor. Virology. 2014; 460:180–93. pmid:25010284

59. Follstaedt SC, Barber SA, Zink MC. Mechanisms of minocycline-induced suppression of simian immunodeficiency virus encephalitis: inhibition of apoptosis signal-regulating kinase 1. J Neurovirol. 2008; 14(5):376–88. pmid:19003592

60. Lafferty MK, Fantry L, Bryant J, Jones O, Hammoud D, Weitzmann MN, et al. Elevated suppressor of cytokine signaling-1 (SOCS-1): a mechanism for dysregulated osteoclastogenesis in HIV transgenic rats. Pathog Dis. 2014; 71(1):81–9. pmid:24376119

61. Mogensen TH, Melchjorsen J, Hollsberg P, Paludan SR. Activation of NF-kappa B in virus-infected macrophages is dependent on mitochondrial oxidative stress and intracellular calcium: Downstream involvement of the kinases TGF-beta-activated kinase 1, mitogen-activated kinase/extracellular signal-regulated kinase kinase 1, and I kappa B kinase. J Immunol. 2003; 170(12):6224–33. pmid:12794154

62. Avey D, Tepper S, Li WW, Turpin Z, Zhu FX. Phosphoproteomic Analysis of KSHV-Infected Cells Reveals Roles of ORF45-Activated RSK during Lytic Replication. Plos Pathog. 2015; 11(7). pmid:26133373

63. Blissard GW, Wenz JR. Baculovirus gp64 envelope glycoprotein is sufficient to mediate pH-dependent membrane fusion. J Virol. 1992; 66(11):6829–35. pmid:1404622

64. Wagner LM, Bayer A, DeLuca NA. Requirement of the N-Terminal Activation Domain of Herpes Simplex Virus ICP4 for Viral Gene Expression. J Virol. 2013; 87(2):1010–8. pmid:23135715

65. Prasad SS, Garg A, Agarwal AK. Enzymatic activities of the human AGPAT isoform 3 and isoform 5: localization of AGPAT5 to mitochondria. J Lipid Res. 2011; 52(3):451–62. pmid:21173190

66. Hyodo K, Taniguchi T, Manabe Y, Kaido M, Mise K, Sugawara T, et al. Phosphatidic Acid Produced by Phospholipase D Promotes RNA Replication of a Plant RNA Virus. Plos Pathog. 2015; 11(5). pmid:26020241

67. Yang RZ, Blaileanu G, Hansen BC, Shuldiner AR, Gong DW. cDNA cloning, genomic structure, chromosomal mapping, and functional expression of a novel human alanine aminotransferase. Genomics. 2002; 79(3):445–50. pmid:11863375

68. Salgado MC, Meton I, Anemaet IG, Baanante IV. Activating transcription factor 4 mediates up-regulation of alanine aminotransferase 2 gene expression under metabolic stress. Bba-Gene Regul Mech. 2014; 1839(4):288–96. pmid:24418603

69. Wu TB, Fu XD. Genomic functions of U2AF in constitutive and regulated splicing. Rna Biol. 2015; 12(5):479–85. pmid:ISI:000354042400001.

70. Zhang L, Stoltzfus CM. A suboptimal src 3' splice site is necessary for efficient replication of Rous sarcoma virus. Virology. 1995; 206(2):1099–107. pmid:7856084

71. Danaher RJ, Savells-Arb AD, Black SA Jr., Jacob RJ, Miller CS. Herpesvirus quiescence in neuronal cells IV: virus activation induced by pituitary adenylate cyclase-activating polypeptide (PACAP) involves the protein kinase A pathway. J Neurovirol. 2001; 7(2):163–8. pmid:11517389

72. de Nijs L, Leon C, Nguyen L, LoTurco JJ, Delgado-Escueta AV, Grisar T, et al. EFHC1 interacts with microtubules to regulate cell division and cortical development. Nat Neurosci. 2009; 12(10):1266–74. pmid:19734894

73. Taterka J, Sutcliffe M, Rubin DH. Selective reovirus infection of murine hepatocarcinoma cells during cell division. A model of viral liver infection. J Clin Invest. 1994; 94(1):353–60. pmid:8040276

74. Arnold F, Schnell J, Zirafi O, Sturzel C, Meier C, Weil T, et al. Naturally occurring fragments from two distinct regions of the prostatic acid phosphatase form amyloidogenic enhancers of HIV infection. J Virol. 2012; 86(2):1244–9. pmid:22090109

75. Ballana E, Senserrich J, Pauls E, Faner R, Mercader JM, Uyttebroeck F, et al. ZNRD1 (Zinc Ribbon Domain-Containing 1) Is a Host Cellular Factor That Influences HIV-1 Replication and Disease Progression. Clin Infect Dis. 2010; 50(7):1022–32. pmid:20192730

76. Dlamini Z, Ntlabati P, Mbita Z, Shoba-Zikhali L. Pyruvate dehydrogenase kinase 4 (PDK4) could be involved in a regulatory role in apoptosis and a link between apoptosis and insulin resistance. Exp Mol Pathol. 2015; 98(3):574–84. pmid:25794976

77. Zhao P, Dong ZM, Duan J, Wang GH, Wang LY, Li YS, et al. Genome-Wide Identification and Immune Response Analysis of Serine Protease Inhibitor Genes in the Silkworm, Bombyx mori. PLoS One. 2012; 7(2):e31168. pmid:22348050

78. Muylkens B, Coupeau D, Dambrine G, Trapp S, Rasschaert D. Marek's disease virus microRNA designated Mdv1-pre-miR-M4 targets both cellular and viral genes. Arch Virol. 2010; 155(11):1823–37. pmid:20680360

79. Yingvilasprasert W, Supungul P, Tassanakajon A. PmTBC1D20, a Rab GTPase-activating protein from the black tiger shrimp, Penaeus monodon, is involved in white spot syndrome virus infection. Dev Comp Immunol. 2014; 42(2):302–10. pmid:24076066

80. Pan HY, Sun HM, Xue LJ, Pan M, Wang YP, Kido H, et al. Ectopic trypsin in the myocardium promotes dilated cardiomyopathy after influenza A virus infection. Am J Physiol-Heart C. 2014; 307(6):H922–32. pmid:25038143

81. Morikawa K, Gouttenoire J, Hernandez C, Thi VLD, Tran HTL, Lange CM, et al. Quantitative Proteomics Identifies the Membrane-Associated Peroxidase GPx8 as a Cellular Substrate of the Hepatitis C Virus NS3-4A Protease. Hepatology. 2014; 59(2):423–33. pmid:23929719

82. Volkman LE. Baculovirus infectivity and the actin cytoskeleton. Curr Drug Targets. 2007; 8(10):1075–83. pmid:17979667.