Skip to main content

Advertisement

You are viewing the new article page. Let us know what you think. Return to old version

RNA-seq comparative analysis of Peking ducks spleen gene expression 24 h post-infected with duck plague virulent or attenuated virus

Abstract

Duck plague virus (DPV), a member of alphaherpesvirus sub-family, can cause significant economic losses on duck farms in China. DPV Chinese virulent strain (CHv) is highly pathogenic and could induce massive ducks death. Attenuated DPV vaccines (CHa) have been put into service against duck plague with billions of doses in China each year. Researches on DPV have been development for many years, however, a comprehensive understanding of molecular mechanisms underlying pathogenicity of CHv strain and protection of CHa strain to ducks is still blank. In present study, we performed RNA-seq technology to analyze transcriptome profiling of duck spleens for the first time to identify differentially expressed genes (DEGs) associated with the infection of CHv and CHa at 24 h. Comparison of gene expression with mock ducks revealed 748 DEGs and 484 DEGs after CHv and CHa infection, respectively. Gene pathway analysis of DEGs highlighted valuable biological processes involved in host immune response, cell apoptosis and viral invasion. Genes expressed in those pathways were different in CHv infected duck spleens and CHa vaccinated duck spleens. The results may provide valuable information for us to explore the reasons of pathogenicity caused by CHv strain and protection activated by CHa strain.

Introduction

After firstly reported in Netherlands at 1923, duck plague (DP) was rapidly spread around the world [1]. In China, the DP was not suspected until 1958 and the outbreak of DP results in a devastated hit to duck industries in China, one of the largest duck breeding countries in the world, due to DP’s high mortality [2, 3]. The disease is caused by duck plague virus (DPV), a member of alphaherpesvirus subfamily, which is a double-stranded DNA virus composed with capsid, tegument and envelope [4]. DPV Chinese virulent strain (CHv), which was isolated from dead infected ducks in China, is highly pathogenic and could induce massive petechial hemorrhages in parenchymal organ, lymphoid and digestive tract and causes massive ducks death [5,6,7,8]. To control this disease on duck farms, attenuated DPV vaccines, such as CHa, have been put into service against duck plague with billions of doses in China each year [9].

Compared to other Herpesviruses, the development of DPV researches, particularly in viral gene functions and virus-host interplay, undergoes a slow process. After the first reports of DPV highly virulent strain and attenuated strain genomes, the studies of gene structure and functions have been spring up. DPV highly virulent strain genes, such as capsid protein genes UL19, UL35 tegument protein genes UL16, UL51 and envelope glycoprotein genes like UL44, US7, US8 were well studied [10,11,12,13,14,15,16]. Furthermore, it has been reported that the five ORFs, including UL2, UL12, US10, UL47 and UL41, might represent DPV virulence factors [17]. However, although exploring basic gene mechanisms of highly virulent DPV strains and attenuated DPV strains on genomic level is indeed helpful, it still lacks systemic and comprehensive studies on nature host duck biological processes when infect with DP virus. Performed works reported that DPV virulent strain conducted a latent infection following a primary infection in trigeminal ganglion and lymphoid tissues including PBL and spleen [18], whereas how the virus processes a latent infection and what the mechanisms of latency infection have not been clear. Meanwhile, although Yuan observed programmed cell death on lymphocytes after infection of highly virulent DPV strain, and the result indicated membrane proteins of DPV could alter characteristics of apoptotic cells in duck embryo fibroblast culture [19]. However, no further studies about DPV-induced apoptosis were observed. Moreover, as we all know, the processes for virus to entry into a cell is critical to initiate their whole life cycles, however, we have little information about DPV entry [20,21,22]. Although genes involved in other herpesvirus entry such as glycoprotein gB and gD were well studied, it is still a puzzle part that whether these genes are necessary in DPV entry or not [23, 24].

In this study, with the help of RNA-seq platform we explored a bioinformatics analysis on duck spleen gene expression for the first time after highly virulent DPV CHv strain and attenuated DPV CHa strain infection at post-infected 24 h. Previous study indicated that attenuated DPVs were rapid distribution and proliferation in duck’s tissues 12 h post-inoculation, however, copies of attenuated DPV were dramatic decline in duck’s tissues 24 h post-infection (hpi), which means that attenuated DPV could rapidly invoke host immune response and inhibit DPV proliferation in the early stage of inoculation [6]. Thus, 24 h was chosen as a samples collection time point to explore reasons of protection stimulated by attenuated DPV in short time after infection. In order to detect accurate biological differences from sequencing data, the setting of multiple biological replicates and the meticulous analysis of statistic were performed. Furthermore, high quality reads mapping was presented by using a paired-end mRNA sequencing approach. Eventually, 748 differentially expressed genes (DEGs) in CHv infected group and 484 DEGs in CHa vaccinated group were selected. Gene pathway analysis of DEGs highlighted several valuable biological processes involved in host immune response, cell apoptosis and viral invasion. Genes expressed in those pathways were different in CHv infected duck spleens and CHa vaccinated duck spleens. Overall, our findings generated in this work firstly provided a detailed view of global changes in duck spleen gene expression when infected with either highly virulent DPV CHv strain or attenuated DPV CHa strain at 24 hpi and the results may provide valuable information for us to explore the reasons of pathogenicity caused by CHv strain and protection activated by CHa strain.

Materials and methods

Virus and experiment design

Highly virulence DPV CHv strain and attenuated modified vaccine DPV CHa strain used in this study were obtained from the Key Laboratory of Animal Disease and Human Health of Sichuan Province. Twenty-seven Peking ducklings were conducted from a DPV-free farm where vaccination against DPV was not implementation. Three experimental groups, including CHv infected group, CHa vaccinated group and normal saline intramuscular injected control group, were explored in this study. Each experimental group contained three biology repeats and each repeat included three ducks. CHv-infected and CHa-vaccinated groups were inoculated using intramuscular administration (4.1 × 108 copies CHv and CHa, respectively). Ducks were euthanized to collect spleen samples at 24 hpi. The samples were rapidly kept in liquid nitrogen for further experiments.

Library preparation for RNA-seq

Total RNA was isolated from spleen samples by TRIzol extraction (Invitrogen, CA, USA) following the manufacturer’s instructions. RNA purity was checked by using the NanoPhotometer® spectrophotometer (IMPLEN, CA, USA) and its concentration was measured using Qubit® RNA Assay Kit in Qubit® 2.0 Flurometer (Life Technologies, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA). A total amount of 3 μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, USA) following manufacturer’s recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5×). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase (RNase H-). Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNaseH were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3′ ends of DNA fragments, NEBNext Adaptors with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150–200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μL USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37 °C for 15 min followed by 5 min at 95 °C before PCR. Then PCR was performed with phusion high-fidelity DNA polymerase, universal PCR primers and index (X) primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kit v3-cBot-HS (Illumia) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on Hiseq-PE150 platform and 125 bp/150 bp paired-end reads were generated.

Bioinformatics analysis of RNA-seq data

Clean reads were obtained to remove reads containing adapters, poly-N and low quality sequences in raw reads. Clean reads were calculated through Q20, Q30 and GC content to make sure data with high quality. All the downstream analyses were based on clean data with high quality. Paired-end clean reads were mapped to duck reference genome using TopHat v2.0.12. Reads numbers of each gene were counted through the HTSeq v0.6.1. In order to normalize the expression of each gene from different samples, FPKM (expected number of Fragments per kilobase of transcript sequence per Millions base pairs sequenced) was introduced. Genes with an expression level more than 1 FPKM were used for downstream analysis. To test the differential expression genes (DEGs) of two groups, CHv infected ducks and CHa vaccinated ducks (three biological replicates per group), the DESeq R package (1.18.0) was performed and genes with a p-value < 0.05 found by DESeq were assigned as differentially expressed.

Gene ontology enrichment analysis

GOseq R package was used to perform gene ontology (GO) enrichment analysis of differential expression genes. GO terms with p-value < 0.05 were selected as significant enrichment.

KEGG pathway analysis

KEGG is a database resource for understanding functions and utilities of the biological system, such as the cell, the organism and the ecosystem, from molecular-level information. We used KOBAS software to test the statistical enrichment of differential expression genes in KEGG pathways. KEGG pathways with p-value < 0.05 were selected as significant enrichment.

Results

Collection of duck spleen samples and construction of cDNA libraries after CHv and CHa infection

Transcriptome sequencing by RNA-seq was used to detect global changes of Peking duck spleen gene expression after CHv and CHa infection. Three experimental groups, including CHv infected group, CHa vaccinated group and normal saline intramuscular injected control group, were explored in this study. Each experimental group contained three biology repeats and each repeat included three ducks. Nine Peking duck spleen samples were collected and stored at liquid nitrogen until total RNA extraction and quality testing was done. Finally, nine cDNA libraries were performed.

Differentially expressed genes after CHv and CHa infection

By using Hiseq-PE150 platform, a total of 571 million (571 134 310) raw data were produced across nine samples. To guarantee ideal results for following genomic mapping and differential gene change analysis, raw reads were filtered to remove low quality data with a total of 545 million (545 498 360) clean reads acquired, a mean of 61.23% of which mapped to duck reference genome (Table 1 ).

Table 1 Number of reads of all bases detected using RNA-seq in DPV-infected and control ducks

Heat map visually compares the whole level of different expression genes (DEGs) and classifies gene expression patterns in different treatment experimental groups under different experimental conditions. As showed in Figure 1A, notwithstanding genes harbored the same or similar expression pattern in CHv, CHa and control treatment groups, the regulation of gene expression was in an entirely different level. For exploring a global view on the change of duck gene expression between different experimental treatment groups, two paired comparisons (CHv vs control, CHa vs control) were performed. In all, RNA-seq analysis detected 748 and 484 genes, which were expressed at significantly different levels in CHv and CHa groups compared to control at a p value < 0.05, respectively (Figure 1B). The purpose of this study is to explore differences of duck spleen gene expression at 24 h between high-virulence DPV CHv strain infected group and attenuated DPV CHa strain vaccinated group through high throughput RNA-seq analysis. In order to learn clearly and precisely about differentially expressed genes, venn analysis was utilized (Figure 1B) and there were three main parts of it. From the green part, there were 561 genes represented the numbers of genes only expressed in CHv infected ducks when compared to control and as for the red part, there were 297 DEGs that only expressed in CHa vaccinated ducks when compared to control. The intersection part with 187 DE genes were co-expressed genes in both CHv infected and CHa vaccinated groups when compared to control.

Figure 1
figure1

Differentially expressed genes in different experimental conditions. Heatmap is used to classify gene expression patterns under different experimental conditions and Venn diagram displays a global view on the numbers of differentially expressed genes. A Genes with similar expressed patterns were clustered and showed above in the heatmap. Intensity of color indicates gene expression levels that were normalized according to log10 (FPKM + 1) values. Red color represent high expression level genes and blue color represent low expression level genes. The three major clusters represent CHv infected duck group, CHa vaccinated ducks group and control group. B The overlap of differentially expressed genes of CHv vs Control and CHa vs Control. The number in the diagram indicated gene number refers to each comparison.

Annotation of duck spleen DEGs based on GO analysis after CHv infection or CHa vaccination

To ensure differentially expressed genes function, gene ontology (GO) analysis was performed to categorize and annotate DEGs of all of three parts in Venn diagram into three groups including biological process (BP), cellular component (CC) and molecular function (MF), respectively. As shown in Table 2, a total of 2654 GO terms were assigned to DE genes only expressed in CHv infected group including 1806 BP terms, 511 MF terms and 341 CC terms. As for DE genes only expressed in CHa vaccinated group, a total of 1479 GO terms were annotated including 976 BP terms, 277 MF terms and 226 CC terms. Furthermore, DE genes expressed in both CHv infected group and CHa vaccinated group were categorized into GO terms with a total numbers of 1098 including 697 BP terms, 215 MF terms and 186 CC terms.

Table 2 A summary of GO term numbers after CHv and CHa infection

And then, in order to pick out the helpful and useful genes for further exploration, 30 significant GO terms were listed and each terms contained more than 2 genes (Figure 2). Top 30 GO terms were selected according p-value < 0.05. Top three significant GO terms of DEGs expressed in both CHv infected group and CHa vaccinated group were protein localization to centrosome (GO:0071539), smooth muscle contractile fiber (GO:0030485) and extracellular matrix structural constituent (GO:0005201). Top three significant GO terms of DEGs only expressed in CHv infected group were paraxial mesoderm formation (GO:0048341), MHC class I protein complex (GO:0042612) and transmembrane receptor protein tyrosine kinase adaptor activity (GO:0005068). Top three significant GO terms of DEGs expressed only in CHa vaccinated group were positive regulation of cellular extravasation (GO:0002693), clathrin coat of trans-Golgi network vesicle (GO:0030130) and dipeptidase activity (GO:0016805).

Figure 2
figure2

Top 30 Gene ontology (GO) terms of DEGs expressed in CHv infected or CHa vaccinated group. GO-terms were processed under three categories including cellular component (CC), molecular function (MF) and biological process (BP). Top 30 GO terms were selected according p-value < 0.05 A GO annotation of DEGs expressed only in CHv infected group B GO annotation of DEGs expressed only in CHa vaccinated group C GO annotation of DEGs expressed in both CHv infected group and CHa vaccinated group.

Pathway analysis of DEGs based on KEGG after CHv and CHa infection

KEGG pathway analysis was used as a combination way to explore the function of DEGs. Top 20 enrichment KEGG pathways were listed in Figure 3 according to p value < 0.05. Four functional categories were identified to potentially play important roles associated with CHv infection including phagosome, antigen processing and presentation, leukocyte transendothelial migration and Fc gamma R-mediated phagocytosis. Three functional categories were identified to potentially play important roles associated with CHa vaccination including cytokine–cytokine receptor interaction, antigen processing and presentation as well as phagosome. Three functional categories were identified to potentially play important roles associated with both CHv infection and CHa vaccination including ECM-receptor interaction, focal adhesion as well as protein digestion and absorption. Putative functional roles and interactions of these genes involved in mediating host response of ducks are discussed in detail.

Figure 3
figure3

Top 20 KEGG pathways in CHv infected or CHa vaccinated group. A KEGG pathways of DEGs expressed in CHv infected group. B KEGG pathways of DEGs expressed in CHa vaccinated group. C KEGG pathways of DEGs expressed after CHv and CHa infection.

Discussion

Duck plague (DP) is a disease that has been reported to cause significant economic losses on duck farms per year in China [3]. Although attenuated vaccine inoculation is an effective way for preventing this disease on duck industries, the potential risk of live attenuated vaccines of virulence return has raised more and more concern. Furthermore, the lack of effective therapeutic tools is another worry to control this disease spreading, indicating that it is urgent to not only give a deeply insight into pathogenesis mechanism of DPV highly pathogenic strain, but also give a profound understanding of protective mechanism of DPV attenuated vaccine strain. From our results, we performed RNA-seq technology for the first time and explored important information about host-viral interaction at post 24 h infection in which the ducks were infected by DPV highly pathogenic strain CHv and DPV attenuated vaccine strain CHa.

Antigen processing and presentation

Cytotoxic T lymphocytes (CTLs) recognized virus-derived peptides presenting by MHC class I and MHC class II at cell surface is a critical immune monitor system and this high sensitive antigen processing and presentation system also becomes a prime target to escape host immune surveillance deployed by many viruses, particularly large DNA virus like herpesvirus [25]. Thus, antigen processing and presentation pathway annotated by KEGG analysis evoked our strong interests.

From results, we got differentially expressed genes related to antigen processing and presentation in KEGG pathway (Table 3). It was no surprise that at 24 hpi, duck spleen genes expression related to antigen processing and presentation had been changed when ducks infected no matter by CHv strain or by CHa strain. However, the results showed that there was huge diversity on the types of differentially expressed genes. CHv infection resulted in genes differential expression located to both MHC class I and MHC class II pathway, however, as for CHa infection, specific gene expression only changed in MHC class II pathway. Furthermore, from the results we could find that the MHC class II genes expressed were up-regulated in CHv infection group whereas down-regulated in CHa vaccinated group. It means that under the experiment condition in this study, the antigen processing and presentation process presented by highly virulence DPV strain CHv and attenuated DPV strain CHa is totally different in vivo.

Table 3 Genes listed involved in multiple cell biological processes after CHv and CHa infection

Virus antigens are degraded by proteasome in cytosolic origin in MHC class I presentation pathway and then are translocated into the endoplasmic reticulum (ER) by Transporter associated with Antigen Processing (TAP) to load onto newly synthesis MHC class I molecules [26]. In the ER, the MHC class I molecule is assembled from a heavy chain and a light chain called β2-microglobulin (β2m) [27]. The peptide inserts itself deep into the MHC class I peptide-binding groove and then ensure the stability of the complex [27]. The fully assembled molecule loaded with virus protein segments then release from the ER and travels via the Golgi apparatus to the plasma membrane for antigen processing and preparation to CD8+ T cells [28]. From our results, in CHv infected group, genes correlated to MHC class I molecule composition, such as MHC1 and β2m, and antigen processing transporters, such as TAP1 and TAP2, are up-regulated (Table 3). However, in CHa vaccinated group, only the expression of β2m gene was changed. Previous studies reported that herpesviruses could target each steps of MHC class I antigen presenting pathway to evade host immune response. A large number of proteins are encoded by herpesviruses to disturb MHC I antigen presentation and then prevent elimination of the target cell by CTLs. The mRNA of MHC class I high chain is degraded by host shut-off protein of herpes simplex virus (HSV) and bovine herpes virus (BHV-1) [29, 30]. Furthermore, immature MHC class I molecules are directed into cytosol from ER with the use of US10 protein of human cytomegalovirus (HCMV) and finally degraded by proteinase [31]. Moreover, with the help of ICP47, HSV could inhibit the transport of peptides through TAP [32]. Thus, we conclude that there may be some critical mechanisms used by DPV CHa strain to evade MHC class I molecules identify and presentation and then finally evade immune response.

Although MHC class I and MHC class II are similar in presenting virus antigen peptides to CD8+ and CD4+, the basic antigen presentation pathways of them are extremely diverse [28]. Firstly, these molecules have different cell distributions. MHC class I molecules are expressed at kinds of nucleated cells, however, MHC class II molecules are specifically found at antigen presenting cells (APCs), such as dendritic cells macrophages and B cells [33]. Unlike the source of viruses degraded in MHC class I molecules presentation pathway, virus protein presented in MHC II pathway is degraded in endosome. As for MHC class II molecules, after assembling with the transmembrane chains and the invariant chain (Ii, also named CD74) in the ER, MHC class II is transported to the MHC class II compartment (MIIC). In this part, Ii is digested by endosomal proteases such as asparaginyl endopeptidase (APE) and then leave a residual class II-associated Ii peptide (CLIP) occupying the peptide-binding groove of the MHC class II molecules [34]. In the MIIC, the CLIP fragment is displaced by viral protein segment through participation of the H2-DM chaperone [35]. The MHC II molecules then traffic to the plasma membrane and present virus peptides to CD4+ cells [35]. From our results, extreme diversity of different genes expression of MHC-II pathway has been present in CHv infection group and CHa vaccinated group (Table 3). The expression of genes related to basic composition of MHC-II like CD74 was up-regulated in CHv infected group which means the infection of CHv strain processes the potential to be presented to CD4+ cells and then evoke duck adaptive immune response in 24 hpi. Yet,genes like MHC2 and GILT related to CHa vaccinated group were down-regulated. GILT, a gamma-interferon inducible lysosomal thiol reductase, could catalyze antigenic disulfide bond reduction in endocytosis pathway and this function is essential in MHC-II antigen processing [36]. From our results, we conclude that in order to evade immune response, CHa depressed the expression of GILT and escape the degradation of disulfide bond. It may be an effective way for CHa to avoid identification by CD4+ cells. Meanwhile, we also found that, no matter in CHv infected group or CHa vaccinated group, the expression of Class II Transcriptional Activator (CIITA) gene was down-regulated. Previous works reported that, CIITA is a co-activator and tightly controls expression of all components of MHC-II at the transcriptional level, such as MHC-II and Ii [37]. It means that CIITA plays an important role in MHC-II antigen processing in both CHv and CHa infection.

From this part, we conclude that CHv could be present by MHC-I and MHC-I as well as identified by both CD8+ T cells and CD4+ T cells, and as for CHa, there may be some critical ways to escape the immunosurveillance at 24 h post DPV infection.

Cytokine–cytokine receptor interaction

Cytokines are soluble extracellular proteins or glycoproteins including chemokines, interferons (IFN), interleukins (IL), lymphokines, and tumor necrosis factors (TNF) [38]. Cytokines play crucial roles in many cell biological processes. They are engaged in humoral and cell-based immune responses, cell growth, differentiation, maturation, death, angiogenesis and homeostasis [38]. Through binding through corresponding cytokines receptors, cytokines transmit signals from outside into inside of cells. In cells, through phosphorylation and dephosphorylation events, signal transduction pathways are activated [38]. In our results, from cytokine–cytokine receptor KEGG pathway, we found that the expression of many important cytokines and cytokine receptors genes were changed in CHv infected or CHa vaccinated group. For example, the expression of IFN-γ receptor IFNGR1, the common receptor CSF2RB of IL3, IL5 as well as tumor necrosis factors receptor SF18 were depressed after CHv infection. The expression of PDGFRB genes involved in tyrosine kinase receptors and SF21 related to TNF receptors were increased in CHv infected group. Meanwhile, CX3CL1 and CCL19 chemotactic factor as well as transforming growth factor TGFB3 were up-regulated in CHa vaccinated group. Cytokines receptors including IL-6 receptor CNTFR as well as tyrosine kinase receptors FLT1 and FLT4 were also increased after CHa vaccination. Yet, the expression of IL22 receptor, IL22RA2 was down regulation in CHa vaccination. Furthermore, the expression of IL15 and TNFSF13B were down regulation and epidermal growth factor (EGF) receptor EGFR was up regulation in both CHv and CHa infected group. Some significant genes related to viral immune evasion were picked out and performed detailed illustrations.

Viral invasion of a host cell triggers immune responses and the major protection against viral infection implemented by host is activation of the interferon (IFN)-mediated antiviral pathway. IFNs, one of a family of cytokines, contain two subtypes including type I IFNs (IFN-α and IFN-β) and type II IFN (IFN-γ). IFN-α/β is responsible for the recognition and clearance of virus as the first innate immune defense. IFN-γ is produced by activate T cells and nature kill (NK) cells. It is involved in many biological processes including immunosurveillance and cell-mediated immune response [39]. From our results, we found that the expression of IFN-γ receptor IFNGR1 was down-regulated after CHv infection. It has been reported that Kaposi’s sarcoma-associated herpesvirus (KSHV) has been shown to inhibit the function of IFN-γ through depressing the expression of IFNGR1 with the help of K3 and K5 protein [40]. Similarly, another herpesvirus, Epstein–Barr virus (EBV), also block the process of IFN response through down regulation of the IFN-γ receptor gene expression by BZLF1 [41]. It means that just like another herpesvirus, highly virulence DPV CHv strain suppress both cytokines-induced and cell-mediated immunity response through depressing the expression of IFN-γ receptor. It may be an efficient way for CHv to avoid host immune control.

Chemokines are pivotal regulatory factors on the aspects of the accumulation, activation and movement of leukocytes into inflamed cells [38]. CX3CL1 is a cell surface-expressed chemokine and plays a distinct role in several aspects of leukocyte activation and motion during inflammation and viral infection [42]. KSHV encodes CX3CL1 receptor homologs protein US28 to prevent the function of CX3CL1 and then serve the immune evasion purpose [42]. In CHa vaccinated group, CX3CL1 gene was increased in expression. We conclude that CXC3L1 may perform an important role no matter in host immune defense or in viral immune evasion during CHa infection. Further studies need to be done to uncover the exact functions of CXC3L1 in CHa infection.

Nature killing (NK) cells are Innate Lymphoid Cells (ILC) and are capable of killing virus-infected cells in the early stage of infections. IL15, type I cytokines, possesses a pivotal biological activity in stimulating NK cells and controlling various aspects of NK cell immunologic surveillance process [43]. For instance, previous studied reported that IL15 stimulated NK cells cytotoxicity and then severely depress human herpesvirus-6 (HHV-6) expression [44]. This study showed that IL-5 gene expression was reduced after CHv and CHa infection. It means that depression of the expression of IL15 may be an effective way for DPV to evade NK cell immunologic monitor at 24 hpi.

Apoptosis

Cells infected with viruses would occur apoptosis to protect other cells from infection and reduce production of new virus from these infected cells [45]. Therefore, in order to survive long enough in cells, viruses have evolved apoptotic inhibitory mechanisms. Successful inhibitory mechanisms alter host cell signaling pathways to govern cell survival and the viral life cycle.

Viruses selectively inhibit apoptotic pathways through blocking the synthetic processes such as depressing the synthesis of pro-apoptotic protein like caspases and p53 or the expression of anti-apoptotic protein like Bcl-2 [46, 47]. Furthermore, PI3K-Akt, an important anti-apoptosis signaling pathway, is commonly exploited by viruses to regulate cell survival [48]. Our results showed that the genes expression in PI3K-Akt signaling pathway was changed in CHv infected group (Table 3). It means that PI3K-Akt signaling pathway may be a key mechanism for highly virulence DPV CHv strain to regulate anti-apoptosis.

Phosphatidylinositol 3-kinase (PI3K) phosphorylates the lipid substrate PtdIns(4,5)P2 to PtdIns(3,4,5)P3 on the plasma membrane [49]. The phosphorylated lipid acts as a second messenger to activate Akt protein on plasma membrane. Activated Akt phosphorylates a number of target proteins, including transcription factor like NF-κB and pro-apoptotic protein caspase-9 and then lead to suppression of apoptosis [50]. Previous studied reported that PIK3/Akt pathway is also utilized by herpesvirus to maintain cellular survival and viral replication. For example, herpes simplex virus I (HSV-1) induced the phosphorylation of Akt and glycogen synthase kinase 3 (GSK-3) to regulate apoptosis and viral gene expression [48]. From our results, the PIK3CB gene expression was down-regulated after the infection of CHv. PIK3CB, one of the four catalytic subunits of PI3K, plays an important role in cell growth [51]. It means that regulating the expression of PIK3CB may be useful for CHv to anti-apoptosis. Furthermore, some other genes were also found related to apoptosis. IKKA, one of subunits of IkB kinase complex (IKK), from our results, was down regulated in both CHv infected group and CHa vaccinated group. The IKKA is necessary for IkB phosphorylation and NF-kB activation [52]. Moreover, the expression of CASP6, apoptosis-related cysteine protease, was also depressed in CHv infected group and CHa vaccinated group [53]. It means that IKKA and CASP6 genes play a critical role in modulating cell apoptosis in highly virulence CHv infected group and attenuated CHa vaccinated group.

Viral invasion

With no autonomous metabolic and motile capacity, viruses must rely on host cells to achieve genome replication, protein synthesis and structural assembling [54]. In order to finish above processes, viruses must suffer from four steps of cell entry processes that is adsorption to the cell surface, penetration, intracellular transport and egress from nucleus [55].

Adsorption to the cell surface is regarded as a prelude to viral infection and cell surface specific structures serve as viral binding locus. Heparan sulfate proteoglycans (HSPGs) are glycoproteins and participate many biological systems, and one of the important roles for HSPGs is adsorption receptors and facilitates viruses concentrated on the cell surface. HSPGs are the first receptor for many herpes viruses entry process, such as HSV, varicella zoster virus (VZV) and pseudorabies virus (PRV), to attach around host cell surface through their envelope glycoproteins like gB and/or gC [56, 57]. However, the types of adsorption receptors of DPV on duck cells surface and the specific interactive envelope proteins of DPV have not been reported in recent years. In our results, the different gene expression of HSPG2 gene was detected in CHv infected ducks indicating that HSPG perhaps plays an important role in the cell surface adsorption process of DPV (Table 3). Interestingly, although Yong reported that DPV glycoproteins gC played an important role in the initial adsorption on chicken embryo fibroblasts (CEF) cell surface, however, they detected that HSPGs was not interactive sites to DPV gC [58]. Therefore, it is necessary to perform further experiments to explore the functions of HSPGs on the process of DPV entry and the types of DPV glycoprotein interacted with HSPGs.

The next necessary step for enveloped viral entry is penetration reaction and this membrane fusion process is promoted by cues of specific envelope-membrane receptor bindings [54]. Specific cellular receptors participated in penetration process of herpes virus, such as HSV, fall into three classes, and one of important receptors are specific sites in heparan sulphate catalyzed by certain heparan sulfate glucosamine 3-O-sulfotransferase (HS3ST) [59]. Nevertheless, so far there are still no relevant reports paid attention on the types of envelope-membrane binding receptors when ducks infected by DPV. In our study, we detected that the expression of HS3ST genes HS3ST5 was changed in CHv infected ducks (Table 3). It suggests that for highly virulence strain of DPV, HS3ST protein may exercise as the same functions as other herpes virus.

In order to finish viral replication, assembly and maturation, intracellular complexes of virus utilize host intracellular cytoskeletal components and move from cell periphery to synthesis sites [54]. After assembly and maturation in cell, viral particles move from cell center to the cell membrane [54, 60]. Microtubules and actins, two important parts of cell cytoskeleton, play significant roles in herpesvirus transport within cells [61]. For example, microtubules provide long-distance transport for HSV with the assistance of their motor protein, kinesin. Our results showed that the expression of two types of alpha-tubulin (TUBA1C and TUBA1B) gene, whose products are basic components of microtubules, was changed in CHv infected ducks. Notably, the expression of kinesin-related genes were also found changed in CHa vaccinated ducks. N-type kinesin, one of a microtubules-binding plus-end-directed motors, could transport viral particles from their replication and transcription sites to cell membrane [62]. KIF1A, one of a member of N-types of kinesins, plays a key role in HSV transport, and its co-localization with herpes simplex virus II (HSV-2) UL56 protein were also reported [61]. In our study, a N-types kinesin gene KIF14 was detected 
different expression in CHa infected ducks (Table 3). It suggests that CHa strain could rely on important motors of microtubules to finish its cellular long-distance movement just like other herpes virus. In addition to microtubules, another medium, actin filament, is essential for herpesvirus to finish short-range transport [61]. It also contains fast-growing plus-ends and slow-growing minus-ends to link cell center and cell surface [61]. Myosin is a crucial motor to mediate the movement of HSV along actin filaments [61]. Our results also found the expression of some actin-related genes and myosin-related genes to give expression changes in both CHv and CHa infected duck (Table 3), indicating that short-distance transport is an important way to complete intercellular movement process of DPV and some relevant proteins are also used to finish this process.

Although the mechanisms by which DPV interacts with cell surface receptor to achieve viral adsorption, membrane fusion and motors have not been clearly studied, several significant genes related to DPV invasion were detected in this results and it may provide an initial view and supply a foundational data support for future experimental researches in DPV invasion aspect.

Abbreviations

DPV:

duck plague virus

CHv:

DPV Chinese virulent strain

CHa:

attenuated DPV vaccines

HSV:

herpes simplex virus

BHV-1:

bovine herpes virus

HCMV:

human cytomegalovirus

KSHV:

Kaposi’s sarcoma-associated herpesvirus

EBV:

Epstein–Barr virus

HHV-6:

human herpesvirus-6

VZV:

varicella zoster virus

PRV:

pseudorabies virus

DEGs:

differentially expressed genes

GO:

gene ontology

BP:

biological proces

CC:

cellular component

MF:

molecular function

CTLs:

cytotoxic T lymphocytes

ER:

endoplasmic reticulum

β2m:

β2-microglobulin

TAP:

transporter associated with Antigen Processing

APCs:

antigen presenting cells

MIIC:

MHC class II compartment

APE:

asparaginyl endopeptidase

CLIP:

class II-associated Ii peptide

CIITA:

class II transcriptional activator

IFN:

interferon

IL:

interleukin

TNF:

tumor necrosis factors

EGF:

epidermal growth factor

ILC:

innate lymphoid cell

PI3K:

phosphatidylinositol 3-kinase

GSK-3:

glycogen synthase kinase 3

HS3ST:

heparan sulfate glucosamine 3-O-sulfotransferase

References

  1. 1.

    Converse KA, Kidd GA (2001) Duck plague epizootics in the United States, 1967-1995. J Wildl Dis 37:347–357

  2. 2.

    Hansen WR (2008) Duck virus enteritis. In: Gough RE (ed) Duck plague, vol 4. Blackwell Publishing Professional, Hoboken

  3. 3.

    Wang G, Qu Y, Wang F, Hu D, Liu L, Li N, Yue R, Li C, Liu S (2013) The comprehensive diagnosis and prevention of duck plague in northwest Shandong province of China. Poult Sci 92:2892–2898

  4. 4.

    Roizmann B, Desrosiers R, Fleckenstein B, Lopez C, Minson A, Studdert M (1992) The family Herpesviridae: an update. Arch Virol 123:425–449

  5. 5.

    Dhama K, Kumar N, Saminathan M, Tiwari R, Karthik K, Kumar MA, Palanivelu M, Shabbir MZ, Malik YS, Singh RK (2017) Duck virus enteritis (duck plague)—a comprehensive update. Vet Q 37:57–80

  6. 6.

    Qi X, Yang X, Cheng A, Wang M, Guo Y, Jia R (2009) Replication kinetics of duck virus enteritis vaccine virus in ducklings immunized by the mucosal or systemic route using real-time quantitative PCR. Res Vet Sci 86:63–67

  7. 7.

    Qi X, Yang X, Cheng A, Wang M, Zhu D, Jia R, Luo Q, Chen X (2009) Intestinal mucosal immune response against virulent duck enteritis virus infection in ducklings. Res Vet Sci 87:218–225

  8. 8.

    Guo Y, Shen C, Cheng A, Wang M, Zhang N, Chen S, Zhou Y (2009) Anatid herpesvirus 1 CH virulent strain induces syncytium and apoptosis in duck embryo fibroblast cultures. Vet Microbiol 138:258

  9. 9.

    Yang X, Qi X, Cheng A, Wang M, Zhu D, Jia R, Chen X (2010) Intestinal mucosal immune response in ducklings following oral immunisation with an attenuated duck enteritis virus vaccine. Vet J 185:199–203

  10. 10.

    Cai MS, Cheng AC, Wang MS, Chen WP, Zhang X, Zheng SX, Pu Y, Lou KP, Zhang Y, Sun L et al (2010) Characterization of the duck plague virus UL35 gene. Intervirology 53:408–416

  11. 11.

    Chang H, Cheng A, Wang M, Zhu D, Jia R, Liu F, Chen Z, Luo Q, Chen X, Zhou Y (2010) Cloning, expression and characterization of gE protein of duck plague virus. Virol J 7:120

  12. 12.

    Li L, Cheng A, Wang M, Xiang J, Yang X, Zhang S, Zhu D, Jia R, Luo Q, Zhou Y (2011) Expression and characterization of duck enteritis virus gI gene. Virol J 8:1–9

  13. 13.

    Lian B, Xu C, Cheng A, Wang M, Zhu D, Luo Q, Jia R, Bi F, Chen Z, Zhou Y, Yang Z, Chen X (2010) Identification and characterization of duck plague virus glycoprotein C gene and gene product. Virol J 7:1–11

  14. 14.

    Dai B, Cheng A, Wang M (2012) Bioinformatics Analysis and characteristics of UL19 protein of duck plague virus. Appl Mech Mater 236–237:55–60

  15. 15.

    He Q, Yang Q, Cheng A, Wang M, Xiang J, Zhu D, Jia R, Luo Q, Chen Z, Zhou Y, Chen X (2011) Expression and characterization of UL16 gene from duck enteritis virus. Virol J 8:413

  16. 16.

    Shen CJ, Cheng AC, Wang MS, Guo YF, Zhao LC, Wen M, Xie W, Xin HY, Zhu DK (2009) Identification and characterization of the duck enteritis virus UL51 gene. Arch Virol 154:1061–1069

  17. 17.

    Wu Y, Cheng A, Wang M, Yang Q, Zhu D, Jia R, Chen S, Zhou Y, Wang X, Chen X (2012) Complete genomic sequence of Chinese virulent duck enteritis virus. J Virol 86:5965

  18. 18.

    Shawky S, Schat KA (2002) Latency sites and reactivation of duck enteritis virus. Avian Dis 46:308–313

  19. 19.

    Guiping Y, Anchun C, Mingshu W, Xiaoying H, Yi Z, Fei L (2007) Preliminary study on duck enteritis virus-induced lymphocyte apoptosis in vivo. Avian Dis 51:546–549

  20. 20.

    Tawar RG, Schuster C, Baumert TF (2016) HCV receptors and virus entry. In: Miyamura T, Lemon S, Walker C, Wakita T (eds) Hepatitis C virus, I. Springer, Tokyo, pp 81–103

  21. 21.

    Jardetzky TS, Lamb RA (2014) Activation of paramyxovirus membrane fusion and virus entry. Curr Opin Virol 5:24–33

  22. 22.

    Vanderlinden E, Naesens L (2013) Emerging antiviral strategies to interfere with influenza virus entry. Med Res Rev 34:301–339

  23. 23.

    Perelygina L, Patrusheva I, Vasireddi M, Brock N, Hilliard J (2015) B virus (Macacine herpesvirus 1) glycoprotein D is functional but dispensable for virus entry into macaque and Human skin cells. J Virol 89:5515–5524

  24. 24.

    Chesnokova LS, Ahuja MK, Hutt-Fletcher LM (2014) Epstein–Barr virus glycoprotein gB and gH/gL can mediate fusion and entry in trans, and heat can act as a partial surrogate for gH/gL and trigger a conformational change in gB. J Virol 88:12193–12201

  25. 25.

    Horst D, Verweij MC, Davison AJ, Ressing ME, Wiertz EJ (2011) Viral evasion of T cell immunity: ancient mechanisms offering new applications. Curr Opin Immunol 23:96–103

  26. 26.

    Reits E, Griekspoor A, Neijssen J, Groothuis T, Jalink K, Veelen PV, Janssen H, Calafat J, Drijfhout JW, Neefjes J (2003) Peptide diffusion, protection, and degradation in nuclear and cytoplasmic compartments before antigen presentation by MHC Class I. Immunity 18:97–108

  27. 27.

    Ardeniz Ö, Unger S, Onay H, Ammann S, Keck C, Cianga C, Gerçeker B, Martin B, Fuchs I, Salzer U et al (2015) β2-Microglobulin deficiency causes a complex immunodeficiency of the innate and adaptive immune system. J Allergy Clin Immunol 136:392–401

  28. 28.

    Neefjes J, Jongsma ML, Paul P, Bakke O (2011) Towards a systems understanding of MHC class I and MHC class II antigen presentation. Nat Rev Immunol 11:823–836

  29. 29.

    Powers CJ, Früh K (2008) Signal peptide-dependent inhibition of MHC class I heavy chain translation by rhesus cytomegalovirus. PLoS Pathog 4:e1000150

  30. 30.

    Koppers-Lalic D, Rijsewijk FA, Verschuren SB, van Gaans-Van den Brink JA, Neisig A, Ressing ME, Neefjes J, Wiertz EJ (2001) The UL41-encoded virion host shutoff (vhs) protein and vhs-independent mechanisms are responsible for down-regulation of MHC class I molecules by bovine herpesvirus 1. J Gen Virol 82:2071–2081

  31. 31.

    Park B, Spooner E, Houser BL, Strominger JL, Ploegh HL (2010) The HCMV membrane glycoprotein US10 selectively targets HLA-G for degradation. J Exp Med 207:2033–2041

  32. 32.

    Cresswell P, Ackerman AL, Giodini A, Peaper DR, Wearsch PA (2005) Mechanisms of MHC class I-restricted antigen processing and cross-presentation. Immunol Rev 207:145–157

  33. 33.

    Paul P, van den Hoorn T, Jongsma ML, Bakker MJ, Hengeveld R, Janssen L, Cresswell P, Egan DA, van Ham M, Ten Brinke A, Ovaa H, Beijersbergen RL, Kuijl C, Neefjes J (2011) A genome-wide multidimensional RNAi screen reveals pathways controlling MHC class II antigen presentation. Cell 145:268–283

  34. 34.

    Cresswell P, Roche PA (2014) Invariant chain-MHC class II complexes: always odd and never invariant. Immunol Cell Biol 92:471–472

  35. 35.

    Kenneth M, Casey W (2017) Janeway’s immunobiology. Garland Science. http://www.garlandscience.com/product/isbn/9780815344452?fromSearchResults=fromSearchResults, https://scholar.google.com/scholar?hl=zh-CN&q=Janeway%27s+Immunobiology+&btnG=&lr=. Accessed 28 Mar 2017

  36. 36.

    Huang WS, Duan LP, Huang B, Zhou LH, Liang Y, Tu CL, Zhang FF, Nie P, Wang T (2015) Identification of three IFN-γ inducible lysosomal thiol reductase (GILT)-like genes in mud crab Scylla paramamosain with distinct gene organizations and patterns of expression. Gene 570:78–88

  37. 37.

    Reith W, Leibundgutlandmann S, Waldburger JM (2005) Regulation of MHC class II gene expression by the class II transactivator. Nat Rev Immunol 5:793–806

  38. 38.

    Turner MD, Nedjai B, Hurst T, Pennington DJ (2014) Cytokines and chemokines: at the crossroads of cell signalling and inflammatory disease. Biochim Biophys Acta 1843:2563–2582

  39. 39.

    Pestka S, Krause CD, Walter MR (2004) Interferons, interferon-like cytokines, and their receptors. Immunol Rev 202:8–32

  40. 40.

    Li Q, Means R, Lang S, Jung JU (2007) Downregulation of gamma interferon receptor 1 by Kaposi’s sarcoma-associated herpesvirus K3 and K5. J Virol 81:2117–2127

  41. 41.

    Morrison TE, Mauser A, Wong A, Ting JP, Kenney SC (2001) Inhibition of IFN-gamma signaling by an Epstein–Barr virus immediate-early protein. Immunity 15:787–799

  42. 42.

    Mclean KA, Holst PJ, Martini L, Schwartz TW, Rosenkilde MM (2004) Similar activation of signal transduction pathways by the herpesvirus-encoded chemokine receptors US28 and ORF74. Virology 325:241–251

  43. 43.

    Marçais A, Cherfils-Vicini J, Viant C, Degouve S, Viel S, Fenis A, Rabilloud J, Mayol K, Tavares A, Bienvenu J, Gangloff YG, Gilson E, Vivier E, Walzer T (2014) The metabolic checkpoint kinase mTOR is essential for IL-15 signaling during the development and activation of NK cells. Nat Immunol 15:749–757

  44. 44.

    Arena A, Merendino RA, Bonina L, Iannello D, Stassi G, Mastroeni P (2004) Role of IL-15 on monocytic resistance to human herpesvirus 6 infection. New Microbiol 23:105–112

  45. 45.

    Goodkin ML, Morton ER, Blaho JA (2015) Herpes simplex virus infection and apoptosis. Int Rev Immunol 23:141–172

  46. 46.

    Maruzuru Y, Koyanagi N, Takemura N, Uematsu S, Matsubara D, Suzuki Y, Arii J, Kato A, Kawaguchi Y (2016) p53 is a host cell regulator during herpes simplex encephalitis. J Virol 90:6738–6745

  47. 47.

    Kvansakul M, Hinds MG (2015) The Bcl-2 family: structures, interactions and targets for drug discovery. Apoptosis 20:136–150

  48. 48.

    Hsu MJ, Wu CY, Chiang HH, Lai YL, Hung SL (2010) PI3K/Akt signaling mediated apoptosis blockage and viral gene expression in oral epithelial cells during herpes simplex virus infection. Virus Res 153:36–43

  49. 49.

    Eramo MJ, Mitchell CA (2016) Regulation of PtdIns(3,4,5)P3/Akt signalling by inositol polyphosphate 5-phosphatases. Biochem Soc Trans 44:240–252

  50. 50.

    Cooray S (2004) The pivotal role of phosphatidylinositol 3-kinase-Akt signal transduction in virus survival. J Gen Virol 85:1065–1076

  51. 51.

    Lin Z, Zhou P, von Gise A, Gu F, Ma Q, Chen J, Guo H, van Gorp PR, Wang DZ, Pu WT (2014) Pi3kcb links Hippo-YAP and PI3K-AKT signaling pathways to promote cardiomyocyte proliferation and survival. Circ Res 116:35–45

  52. 52.

    Hinz M, Scheidereit C (2014) The IκB kinase complex in NF-κB regulation and beyond. EMBO Rep 15:46–61

  53. 53.

    Wu H, Che X, Zheng Q, Wu A, Pan K, Shao A, Wu Q, Zhang J, Hong Y (2014) Caspases: a molecular switch node in the crosstalk between autophagy and apoptosis. Int J Biol Sci 10:1072–1083

  54. 54.

    Marsh M, Helenius A (2006) Virus entry: open sesame. Cell 124:729–740

  55. 55.

    Greber UF, Way M (2006) A superhighway to virus infection. Cell 124:741–754

  56. 56.

    Sarrazin S, Lamanna WC, Esko JD (2011) Heparan sulfate proteoglycans. Cold Spring Harb Perspect Biol 3:pii:a004952

  57. 57.

    Rue CA, Ryan P (2002) Characterization of pseudorabies virus glycoprotein C attachment to heparan sulfate proteoglycans. J Gen Virol 83:301–309

  58. 58.

    Hu Y, Liu X, Zou Z, Jin M (2013) Glycoprotein C plays a role in the adsorption of duck enteritis virus to chicken embryo fibroblasts cells and in infectivity. Virus Res 174:1–7

  59. 59.

    Spear PG (2004) Herpes simplex virus: receptors and ligands for cell entry. Cell Microbiol 6:401–410

  60. 60.

    Taylor MP, Koyuncu OO, Enquist LW (2011) Subversion of the actin cytoskeleton during viral infection. Nat Rev Microbiol 9:427–439

  61. 61.

    Lyman MG, Enquist LW (2009) Herpesvirus interactions with the host cytoskeleton. J Virol 83:2058–2066

  62. 62.

    Dodding MP, Way M (2011) Coupling viruses to dynein and kinesin-1. EMBO J 30:3527–3539

Download references

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AC and MW conceived and designed the experiments; TL performed the experiments; TL, RJ, QY, YW, KS, DZ, SC, ML, XZ and XC analyzed the data and contributed analysis tools; TL wrote the paper. All authors read and approved the final manuscript.

Acknowledgements

This research was supported by the National Key Research and Development Program of China (2016YFD0500800), China Agricultural Research System (CARS-42-17) and Special Fund for Key Laboratory of Animal Disease and Human Health of Sichuan Province (2016JPT0004).

Animal ethics

The animal studies were approved by Institutional Animal Care and Use Committee of Sichuan Agricultural University, Sichuan, China and followed the National Institutes of Health guidelines for the performance of animal experiments.

Consent for publication

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Correspondence to Anchun Cheng or Mingshu Wang.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark