Virulent duck enteritis virus infected DEF cells generate a unique pattern of viral microRNAs and a novel set of host microRNAs
BMC Veterinary Research volume 14, Article number: 144 (2018)
Duck enteritis virus (DEV) belongs to the family Herpesviridae and is an important epornitic agent that causes economic losses in the waterfowl industry. The Chinese virulent (CHv) and attenuate vaccines (VAC) are two different pathogenic DEV strains. MicroRNAs (miRNAs) are a class of non-coding RNAs that regulate gene expression in viral infection. Nonetheless, there is little information on virulent duck enteritis virus (DEV)-encoded miRNAs.
Using high-throughput sequencing, we identified 39 mature viral miRNAs from CHv-infected duck embryo fibroblasts cells. Compared with the reported 33 VAC-encoded miRNAs, only 13 miRNA sequences and 22 “seed sequences” of miRNA were identical, and 8 novel viral miRNAs were detected and confirmed by stem-loop RT-qPCR in this study. Using RNAhybrid and PITA software, 38 CHv-encoded miRNAs were predicted to target 41 viral genes and formed a complex regulatory network. Dual luciferase reporter assay (DLRA) confirmed that viral dev-miR-D8-3p can directly target the 3’-UTR of CHv US1 gene (p < 0.05). Gene Ontology analysis on host target genes of viral miRNAs were mainly involved in biological regulation, cellular and metabolic processes. In addition, 598 novel duck-encoded miRNAs were detected in this study. Thirty-eight host miRNAs showed significant differential expression after CHv infection: 13 miRNAs were up-regulated, and 25 miRNAs were down-regulated, which may affect viral replication in the host cell.
These data suggested that CHv encoded a different set of microRNAs and formed a unique regulatory network compared with VAC. This is the first report of DEF miRNAs expression profile and an analysis of these miRNAs regulatory mechanisms during DEV infection. These data provide a basis for further exploring miRNA regulatory roles in the pathogenesis of DEV infection and contribute to the understanding of the CHv-host interaction at the miRNA level.
Duck viral enteritis, also called as duck plague, is an acute, contagious and fatal disease of duck and geese, resulting in considerable economic losses in the waterfowl breeding industry [1,2,3,4]. The causative agent of this disease is duck enteritis virus (DEV) which belongs to the species Anatid herpesvirus I, genus Mardivirus, subfamily Alphaherpesvirinae, family Herpesviridae . Many countries, such as China, Britain, the United States, Germany, and Netherlands have reported the prevalence of this virus [6,7,8]. The genome of DEV is a linear double-stranded DNA molecule composed of a unique long region (UL) and a unique short region (US) flanked by an internal repeat sequence (IRS) and a terminal repeat sequence (TRS). Its genomic arrangement pattern (UL-IRS-US-TRS) is consistent with the members of Marek’s disease virus 1 and 2 (MDV-1 and MDV-2), herpes simplex virus types 1 and 2 (HSV-1 and HSV-2) and Pseudorabies virus (PRV) [7, 8].
MicroRNAs (miRNAs) are small (18–24 nt), endogenous non-coding RNAs that widely found in plant, animal and viral genomes and are now increasingly recognized as important regulators of gene expression through post-transcriptional mechanisms, leading to mRNA degradation or translational inhibition by binding to fully or partially complementary 3′ untranslated regions (3’UTR) . These small miRNAs participate in a variety of biological processes, including cellular proliferation, differentiation, apoptosis, signal transduction and the process of virus-host interactions [10,11,12,13,14].
Over 300 virus-encoded miRNAs have been identified (miRBase 22.0). They were encoded by multiple virus families [15, 16], such as herpesviruses adenoviruses, polyomaviruses and retroviruses [17,18,19]. approximately 95% of viral miRNAs were encoded by herpesvirus families . This phenomenon suggested the importance of miRNA-mediated gene regulation in the biology of herpesvirus infections. Some functions of viral miRNAs were validated by experiments in the pathogenesis of herpesvirus infection [21, 22].
As with many other miRNA-encoding α-herpesviruses [23,24,25,26,27,28], DEV-encoded miRNAs were identified from VAC-infected chicken embryo fibroblast (CEF) by deep sequencing technology . This research identified 24 pre-miRNAs in VAC genome producing 33 mature miRNAs. The VAC strain was attenuated and was widely used against duck viral enteritis , while the CHv strain (Chinese virulent DEV strain) can cause epidemical and fatal disease in waterfowl . CHv and VAC are two different pathogenic DEV strains . The mechanism of the two viruses causing different pathogenesis is not well understood. Our aim was to confirm whether the CHv encoded the same miRNAs as VAC and explore those miRNAs regulatory roles in CHv infection. Moreover, recent studies have demonstrated that host miRNAs play crucial roles in viral infection [20, 21], but DEF-encoded miRNAs have not been reported until now. For the above purposes, we constructed and analysed the miRNA expression profile from CHv-infected and uninfected DEF cells using high-throughput sequencing. The potential targets of viral and host miRNAs were predicted by RNAhybrid and PITA software. These data may contribute to the understanding of CHv pathogenesis and the CHv virus-host interaction at the overall miRNA level.
Virus and cells
CHv (GenBank accession No. JQ647509), a classic Chinese virulent strain, was isolated from an infected duck farm and kept in our laboratory. Primary duck embryo fibroblast (DEF) cells were made using 10-day-old embryonated duck eggs (Chengdu Egg & Poultry Co. China) for virus propagation. The use of duck embryos in this study was approved by the Animal Ethics Committee of Sichuan Agricultural University (approval No. XF2014–18). Cell monolayers were cultured in Dulbecco’s Modified Eagle’s Medium (DMEM, Gibco, Grand Island, NY USA) supplemented with 8% foetal bovine serum (FBS, Gibco, USA) and 1% penicillin-streptomycin (Gibco, USA) at 37 °C in a 5% CO2 humidified incubator.
Isolation and sequence of RNA
Duck embryo fibroblasts (DEF) cells (80% confluency per dish) were infected with CHv at a multiplicity of infection (MOI) of 1.0, with mock-infected DEF as a control. Cells were harvested at 2, 4, 6, 8, 12, 18, 24 and 30 h post-infection (hpi) and resuspended in TRIzol (TIANGEN, Beijing, China). Total RNAs from DEV-infected and uninfected DEF cells at the above time points were extracted according to the manufacturer’s directions (TIANGEN, Beijing, China) and quantified using a NanoDrop 2000 Spectrophotometer (Thermo, Carlsbad, CA, USA). The RNA (0.125 μg) extracted from the eight time points was mixed as a group. Our experiments were performed in triplicate and all the infected and control samples were subjected to Huada (Guangdong, China) for high-throughput sequencing of small RNAs (sRNAs). The same mixed RNA samples were used in the subsequent stem-loop RT-qPCR experiments.
The CHv genome has been sequenced and the total size is 162,175 bp. The annotated VAC-encoded miRNAs were from miRBase 22.0 (http://www.mirbase.org/). Duck genomic sequences and the 3’UTR of duck genes were downloaded from the Ensembl database (http://www.ensembl.org). The annotated chicken and Zebra Finch mature miRNAs were from miRBase 22.0 (http://www.mirbase.org/).
Analysis of viral small RNAs
The total raw small RNA (sRNA) reads were detected by an Illumina Genome Analyser. The cleaned sequence reads were obtained after the filtering procedure as previously described [32, 33]. Using the Bowtie algorithm , the filtered sRNA reads were aligned to the known DEV pre-miRNA sequences in miRBase 22.0 with no mismatch and then aligned to the corresponding mature miRNA with at least 16 nts overlap allowing offsets. The known CHv-encoded miRNAs including the pre-miRNA sequences, length and count of reads would be obtained. The remaining sRNA reads mapped to genome were subjected for novel miRNA prediction. Mireap software (http://sourceforge.net/projects/mireap/) was used to predict novel miRNA by exploring the secondary structure. Dicer cleavage sites and predicted minimum free energies of unannotated sRNA reads.
Analysis of host small RNAs
There are not any Anas platyrhynchos miRNAs annotated in the miRBase 22.0. All host small RNA sequences were aligned with known mature miRNAs of two reference species (Gallus and Taeniopygia guttata) and Anas platyrhynchos genome by the Bowtie algorithm . Different miRNA expression levels were normalized to get the number of transcripts per million (TPM) in two samples (CHv-infected and uninfected). Normalization formula: Normalized expression = Actual miRNA count/Total count of clean reads*1000000. A change of at least 2-fold between libraries was considered significant. Fold-change formula: Fold-change = log2 (treatment/control). P-value was set as the reported formula . P-value < 0.05 indicated significance differentially expressed miRNA.
Target prediction and GO analysis of viral and host miRNAs
Target genes of viral and host miRNAs were predicted using RNAhybrid and PITA software, and the parameters were strictly set as a previously reported program in the seed sequence . The potential host target genes were analysed using the Gene Ontology (GO) program (http://www.geneontology.org). Gene Ontology enrichment analysis of the target genes was performed using Goseq  to detect the significantly enriched GO terms of the host target. The GO terms with p < 0.05 were considered significant. The WEGO software (http://wego.genomics.org.cn) was used to produce histograms of the GO annotations, including three fields: cellular component, biological process and molecular function.
The stem-loop RT-qPCR was conducted as previously described [36, 38]. Briefly, 1000 ng of RNA mixture were reverse-transcribed to cDNA and then 2 μL cDNA was used for Real-time PCR amplification according to the company kit instructions (Thermo, Carlsbad, CA, USA). All primers used are listed in (Additional file 1: Table S1). The reaction conditions were as follows: reverse transcription was incubated at 50 °C for 45 min and kept at 85 °C for 5 min. Next, real-time PCR was 95 °C for 5 min, 39 cycles of denaturing at 95 °C for 15 s, annealing and extending 60 °C for 60 s, and the cellular miRNA U6 was used as an internal control. The relative expression values were calculated using the comparative 2-ΔΔCt method .
Vector constructs and luciferase assay
The dev-miR-D8-3p mimic and negative control mimic (miR-NC) were synthesized by Ribobio (Guangzhou, China). The CHv US1 gene 3’UTR (nt 136,085–136,248) including the predicted dev-miR-D8-3p binding sites were synthesised and cloned into a pmirGLO vector (Promega, Madison, WI, USA) with SacI and XhoI sites and named pmirGLO-WT-US1, Accordingly, the mutant 3’UTR of the US1 vector was constructed and named pmirGLO-MU-US1. For luciferase assay, COS7 cells were seeded in 96-well plates and co-transfected with dev-miR-D8-3p mimic, miR-NC, pmirGLO-WT-US1 and pmirGLO-MU-US1 with Lipofectamine 3000 (Invitrogen, Carlsbad, CA, USA). We performed site-directed dual luciferase reporter assay (DLRA), and luciferase activity was measured at 36 h post-transfection according to the manufacturer’s protocol (Promega, Madison, WI, USA).
Each experiment was performed in triplicate and the data were presented as the means (M) ± standard deviations (SD) by the software GraphPad Prism (version7.0). The significance of the variability between different treatment groups was determined by one-way analysis of variance (ANOVA) tests of variance using the GraphPad Prism software (version 7.0). P-values < 0.05 was considered statistically significant.
Analysis of sRNA libraries by deep sequencing
In this study, we obtained 12,088,641 and 12,263,713 sRNA reads of 18–30 nucleotides from CHv-infected and uninfected DEF cells. After filtering adapter sequences and low-quality sequences. 11,462,557 (94.82%) and 11,836,099 (96.51%) high quality reads from infected and uninfected sample were obtained, respectively. Among each sample, approximately 89.36% and 92.85% sRNAs ranged from 20 to 24 nt respectively, and most of the sRNA reads were 22 nt in length (Fig. 1a). In addition to miRNAs, other noncoding sRNAs were also detected and categorized by following the priority rule: microRNA (miRNA) > repeat > rRNA > tRNA > snoRNA > snRNA (Additional file 2: Table S2). Ultimately, 7,446,931 (64.97%) and 7,995,424 (67.55%) miRNA reads from CHv-infected and uninfected libraries respectively were matched to the annotated miRNAs of VAC and the two reference species (Gallus gallus and Taeniopygia guttata), and remaining 3,158,331 (27.55%) and 3,085,287 (26.07%) unannotated sRNA reads from two libraries were matched to CHv and the duck genome for predicting novel miRNAs (Fig. 1b, c).
Conservation analysis of miRNAs in CHv and VAC
In our study, we obtained 29 pre-miRNAs (Additional file 3) and 39 mature miRNAs from the CHv strain by deep sequencing. The names, sequences, length and location of 39 mature miRNAs are listed in Table 1. Compared with previously reported 33 mature VAC-encoded miRNAs , 31 of 33 reported miRNAs were detected and were shown in Table 2. The remaining two miRNAs, dev-miR-D2–3p and dev-miR-D10-3p were not detected in our study. Among 31 detected miRNAs, only 13 miRNA sequences were identical, and 18 were different in contrast to VAC-encoded miRNAs (Table 2). Twenty-two miRNAs were identical in the “seed sequence” and the other 9 were not identical. The difference of the “seed sequence” mostly occurs in 2–8 nucleotides at the 5′ end of miRNAs. For example, dev-miR-D19-5p and dev-miR-D21-5p had one deleted base, dev-miR-D7-5p, dev-miR-D11-3p, dev-miR-D13-5p, dev-miR-D14-3p and dev-miR-D23-3p had two deleted bases, the dev-miR-D4-3p had four deleted bases, and the dev-miR-D17-3p had three inserted bases. In addition, 8 novel CHv-encoded miRNAs were identified and were named from dev-miR-D25-5p to dev-miR-D31-3p (Table 1). The pre-miRNA hairpin structures and isoform expression profile of these novel miRNAs are shown in Additional file 3. Thirty-nine CHv-encoded miRNAs were distributed mostly the unique long region (UL) and the repeat region (IRS and TRS) of the genome (Fig. 2). This result was consistent with the previous report about distribution of VAC-encoded miRNAs . We found that 7 miRNAs were present in two copies, which were located in two loci in the CHv genome. Including dev-miR-D20 to dev-miR-D24 (Table 1). Those miRNAs mapped in the internal repeat sequence (IRS) were marked as ‘a’ and the homologous miRNAs in terminal repeat sequences (TRS) were marked as ‘b’ (Fig. 2). including dev-miR-D20a/b-5p, dev-miR-D21a/b-3P, dev-miR-D21a/b-5p dev-miR-D22a/b-3p, dev-miR-D22a/b-5P, dev-miR-D23a/b-3p and dev-miR-D24a/b-3p. This “two-copy” phenomenon seems to be a common feature in α-herpesviruses.
Self-regulation analysis of viral miRNAs
Prediction results showed that 41 viral genes were targeted by 38 viral miRNAs. Some novel viral miRNAs (like dev-miR-D27-5p and dev-miR-D28-3p) could target multiple CHv genes, and the some CHv genes (like UL24, UL28 and UL52) could be targeted by multiple viral miRNAs. A complex regulatory network was formed according to the regulation interaction between viral miRNAs and target genes (Additional file 4: Figure S2).
Regulatory analysis of viral miRNAs on host genes
Analysis results showed that the 3’UTRs of 4703 host genes were targeted by 39 viral miRNAs using the intersection of the two software programs (Additional file 5: Table S3). Gene Ontology (GO) annotation was performed to analyze biological function of the host target genes. The results reflected that these host target genes were mainly concentrated in the cellular process, metabolic process, signal-organism process, biological regulation process and others (Additional file 6: Table S4). Among of these host target genes, GO enrichment analysis showed that 236 genes were related to signaling processes (p < 0.05) and 66 genes were related to immune system processes (p < 0.05) (Fig. 3, Additional file 6: Table S4), which implied that viral miRNAs may play important regulatory function during viral infection and immune evasion.
Expression and differential analysis for host miRNAs
Alignment results showed that 598 mature host miRNAs were detected in this study (Fig. 4a). Among these, 386 (64.5%) miRNAs (264 aligned and 122 novel) were co-expressed in both libraries (Additional file 7: Table S5), 108 (18.1%) miRNAs were unique to the DEV-infected group and 104 (17.4%) miRNAs were unique to the uninfected group (Additional file 8: Table S6). Among the co-expressed host miRNAs, 38 miRNAs were differentially expressed between the CHv-infected sample and uninfected sample (Additional file 9: Table S7). Thirteen were significantly up-regulated and 25 were significantly down-regulated after CHv infection (Fig. 4b). Thirty-eight differentially expressed host miRNAs were predicted to target viral genes using the RNAhybrid and PITA software, and the results showed that the 3’UTRs of 40 viral genes were targeted by 36 host miRNAs by the intersection of two software (Additional file 10: Figure S3).
Stem-loop RT-qPCR for miRNAs confirmation
To further validate deep sequencing results, 8 novel viral miRNAs and 10 randomly differentially expressed host miRNAs were confirmed using stem-loop RT-qPCR. The results obtained by RT-qPCR were highly consistent with the deep sequencing data (Fig. 4c, d).
Dev-miR-D8-3p target the 3’UTR of US1 gene
Dual luciferase reporter assay (DLRA) showed that the luciferase level of the pmirGLO-WT-US1 was significantly repressed by dev-miR-D8-3p compared to the negative control miR-NC (p < 0.05) (Fig. 5a, b). To further ascertain that the down-regulation of targets by dev-miR-D8-3p is binding sites dependent, the binding sites of US1 were mutated and constructed as pmirGLO-MU-US1 vector (Fig. 5a). As expected, the dev-miR-D8-3p lost its repression effect on the mutant vector of pmirGLO-MU-US1. These results indicated that the dev-miR-D8-3p can directly target the CHv US1 gene by 7 nucleotide complementary seed sequence.
Previous research has reported that the VAC encoded 33 mature miRNAs in the viral genome . We obtained 39 mature viral miRNAs from CHv-infected DEF cells, 22 of 39 CHv-encoded miRNAs share identical “seed sequence” with VAC-encoded miRNAs. Another 17 miRNAs (9 different “seed sequence” miRNAs and 8 novel miRNAs) were different in the “seed sequence”. As we know, target-gene recognition of viral miRNA is strictly dependent on the full base complementarity of the “seed sequence”, which covered 2 to 8 nucleotides from the 5′ end of the miRNA . Thus, the stability of the “seed sequence” of viral miRNA is crucial for target-gene discrimination. We speculate that the 22 miRNAs of the identical “seed sequence” play same regulatory roles in DEV-infection. Overall, the data analysis showed that CHv encoded a different pattern of miRNA Compared with VAC, which might form a complex regulatory network between viral miRNAs and their target genes. The differences of miRNAs regulatory network might lead to the differential pathogenesis of these two viruses.
Studies have revealed that viral and host miRNAs play important roles in host-virus interactions [20, 22, 40]. CHv is a virulent herpesvirus that can mainly cause contagious lethal disease in ducks [30, 31] and the VAC is an avirulent virus and has been reported to encode 33 mature miRNAs from VAC-infected CEF cells . However, a precise regulatory network analysis of DEV miRNAs is unlikely to be achieved without the discovery of the virulent DEV miRNAome. In this research, we collected sRNA samples from CHv-infected DEF cells at eight time points to detect as many viral miRNAs as possible. Using High-throughput sequencing technology, we obtained 29 pre-miRNA sequences with 39 mature miRNAs from CHv-infected DEF cells. Eight novel viral miRNAs were predicted and were confirmed by stem-loop RT-qPCR (Fig. 4c). In addition, we also first made a repertoire of DEF cells miRNAs transcriptome in CHv-infected and uninfected cells and have performed a preliminary analysis of the functions of these miRNAs. These data provide a foundation for further investigations on host-herpesvirus interactions.
Among the 8 novel viral miRNAs, the dev-miR-D27-5p, dev-miR-D28-3p, dev-miR-D29-5p and dev-miR-D30-3p were located in coding region of UL26, UL24, UL19 and UL9 gene, respectively. The remaining four miRNAs were encoded in the in noncoding regions of CHv genome. Several reports revealed that most α-herpesvirus-encoded miRNAs were found clustered in the repeat or other adjacent regions of the viral genome [20, 23,24,25,26,27, 41]. However, the 39 CHv-encoded miRNAs were distributed mostly in the unique long region (UL) and the repeat region (IRS and TRS) of the genome (Fig. 2). This result was consistent with the previous report about the distribution of VAC-encoded miRNAs . Moreover, of the seven miRNAs detected with two copies, miR-D22b-3p, miR-D22b-5p and miR-D23b-3p were located in the coding region of ICP4 in an antisense orientation, which could theoretically lead to the cleavage of the transcript and negative regulation of the gene like siRNAs [42,43,44,45,46].
Several studies have confirmed that herpesvirus-encoded miRNAs can target viral immediate-early (IE) genes to regulate viral latent and lytic infection [20,21,22]. The hsv1-miR-H2–3p and hsv1-miR-H6 target the ICP0 and ICP4 genes of HSV-1 respectively , the hsv2-miR-H2–3 target the ICP0 gene of HSV-2 [46, 48], the mdv1-miR-M7-5P target the ICP4 and ICP27 genes of MDV-1  and the litv-miR-I5 target the ICP4 gene of LITV . The above target genes acted as viral IE genes which upregulate early and late genes of herpesvirus subfamilies and downregulate latency-associated transcript (LAT), inducing the virus towards lytic infection [20, 21]. The targeting of IE genes by viral miRNAs was thought to inhibit entry into viral replication and maintain the latent infection state . In our study, dev-miR-D4-3p, dev-miR-D11-3p, and dev-miR-D20-5p were predicted to target the 3’UTR region of the CHv ICP4 gene. Dev-miR-D1-5p, dev-miR-D8-3p, dev-miR-D12–5p, dev-miR-D17-3p, dev-miR-D26-3p, dev-miR-D28-3p and dev-miR-D30-3p were predicted to target the 3’UTR region of the CHv US1 gene. Our results confirmed that dev-miR-D8-3p could directly target the 3’-UTR of the US1 gene. Dev-miR-D13-5p and dev-miR-D14-3p are predicted to target the 3’UTR region of the CHv UL54 gene (Additional file 4: Figure S2b). The ICP4, US1 and UL54 of CHv were considered the functional equivalent of the immediate-early (IE) genes ICP4, ICP22 and ICP27 of HSV-1 [7, 30, 51, 52]. Thus, we speculate that these viral miRNAs may play key roles in the regulation of the CHv lytic and latent infection.
Some virus-encoded miRNAs could regulate the cellular signal pathway to evade the immune response. For example, the hcmv-miR-UL112–3p was reported to target toll-like receptors 2 (TLR2), inhibiting IRAK1/NFκB signaling and avoiding the related inflammatory response . The mdv1-miR-M4-5P regulated the endogenous TLR3 gene that repressed IFN-β production expression and facilitates virus replication . Three virus-encoded miRNAs (e.g., hcmv-miR-UL112–1, ebv-miR-BART2-5p and kshv-miR-K12–7) could repress identical target gene MICB and lead to a similar outcome, evading NK cell recognition and immune response . The kshv-miR-K12–9 and kshv-miR-K12–5 could target IRAK1 and MYD88, respectively, which repressed TLR/IL-1R signaling, resulting in reduced inflammation . The kshv-miR-K12–11 could target IκB kinase epsilon (IKKε), inhibiting type I interferon signal pathway . Moreover, viral miRNAs could regulate cell growth and survival to favour viral replication. For example, kshv-miR-K12–10 could inhibit TWEAK-induced apoptosis by targeting the cellular TWEAKR , which contributed to cell survival. In addition, mdv1-miR-M4-5p could target LTBP1, which suppressed the TGF-β signaling . Kshv-miR-K12–11 targeted SMAD5 which interfered with the TGF-β pathway . The suppression of TGF-β signaling ultimately result in increased cell survival and virally induced oncogenesis [58, 59]. In our study, GO analysis on the cellular targets of viral miRNAs showed that these targets were involved in complex cellular processes, including signal-organism processes, the metabolic pathway, biological regulation, immune response and signaling process.
The virus could alter host miRNA expression profiles to favour viral replication. In our study, 38 cellular miRNAs were expressed differentially in both the CHv-infected library and mock library. These dysregulated host miRNAs were identified to play crucial roles in other viral infections. For example, miR-let-7a was downregulated in NPC cells after EBV-infection, which in turn promoted viral replication by targeting the dicer gene . The gga-miR-26a was downregulated in MDV-infected spleens at cytolytic infection, latency and tumour transformation phases. Decreasing the expression of gga-miR-26a had been shown to contribute to MDV-induced lymphomagenesis upon regulation of NEK6 proteins . The previous research showed that the differential expression of gga-miR-181a contributed to MDV-induced lymphomagenesis by targeting IGF2BP3/MYBL1 genes . Cellular microRNA miR-181b inhibited replication of mink enteritis virus (MEV) by repression of non-structural protein 1(NS1) translation . The gga-miR-15b was downregulated in splenic tumours after MDV infection and had a negative effect on the expression of ATF2, facilitating viral replication by increasing the expression of the ATF2 . Expression of miR-146 was upregulated after EBV infection, which could downregulate levels of IRAK1 and TRAF6 proteins, reducing the activity of host immune and inflammatory response . Recently, miR-148 was reported as a novel biomarker in non-small-cell lung cancer screening . In our prediction results, miR-148a-5p could target UL1, UL2 UL3, UL24 and UL25 genes of CHv. MiR-181a-3p could target UL24, UL54, US3, US5 and US8 genes of CHv (Additional file 10: Figure S3a). The ICP4 gene of CHv was targeted by miR-135a-1-3p and miR-135a-2-3p, while the UL54 gene was targeted by miR-124a-3p, miR-135a-1-3p, miR-135a-2-3p, miR-15b-3p, miR-181a-3p and miR-181b-1-3p. A complex regulatory network was formed between 36 differentially expressed host miRNAs and their 40 viral target genes (Additional file 10: Figure S3b). However, the regulatory functions of these dysregulated cellular miRNAs in the process of CHv replication need further analysis.
In this study, we obtained 39 DEV-encoded miRNAs from CHv-infected DEF cells by high-throughput sequencing. Of these, 8 novel viral miRNAs were detected and confirmed through stem-loop RT-qPCR. Conservative analysis showed that CHv encoded a different set of miRNAs and formed a unique regulatory network compared with VAC. In addition, a total of 598 novel duck-encoded miRNAs were detected by aligning with known mature miRNAs of Gallus gallus and Taeniopygia guttata. This is the first report of a DEF miRNA expression profile and an analysis of these miRNAs regulatory mechanisms during DEV infection.
Duck enteritis virus
Herpes simplex virus 1
Herpes simplex virus 2
Inverted repeated sequences
Inverted repeated sequences
Marek’s disease virus 1
Marek’s disease virus 2
Real-time quantitative reverse transcriotion PCR
terminal repeated sequences
Unique long region
Unique long region
Unique short region
Metwally SA SD, Glisson JR, McDougald LR, Nolan, LK, Suarez DL, Nair V,. Duck Virus Enteritis (Duck Plague). In: Dis poultry, 13th edition 2013:431–440.
Converse KA, Kidd GA. Duck plague epizootics in the United States, 1967-1995. J Wildl Dis. 2001;37(2):347–57.
Kaleta EF, Kuczka A, Kühnhold A, Bunzenthal C, Bönner BM, Hanka K, et al. Outbreak of duck plague (duck herpesvirus enteritis) in numerous species of captive ducks and geese in temporal conjunction with enforced biosecurity (in-house keeping) due to the threat of avian influenza a virus of the subtype Asia H5N1. Dtw. Dtsch Tierarztl Wochenschr. 2007;114(1):3–11.
Abend JR, Uldrick T, Ziegelbauer JM. Regulation of tumor necrosis factor-like weak inducer of apoptosis receptor protein (TWEAKR) expression by Kaposi's sarcoma-associated herpesvirus MicroRNA prevents TWEAK-induced apoptosis and inflammatory cytokine expression. J Virol. 2010;84(23):12139–51.
King AMQ. Virus taxonomy: classification and nomenclature of viruses: ninth report of the international committee on taxonomy of viruses: Elsevier/academic press; 2012.
Wang G, Qu Y, Wang F, Hu D, Liu L, Li N, et al. The comprehensive diagnosis and prevention of duck plague in Northwest Shandong province of China. Poultry Sci. 2013;92(11):2892–8.
Li Y, Huang B, Ma X, Wu J, Li F, Ai W, et al. Molecular characterization of the genome of duck enteritis virus. Virology. 2009;391(2):151–61.
Wang J, Höper D, Beer M, Osterrieder N. Complete genome sequence of virulent duck enteritis virus (DEV) strain 2085 and comparison with genome sequences of virulent and attenuated DEV strains. Virus Res. 2011;160(1):316–25.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116(2):281–97.
Boss IW, Plaisance KB, Renne R. Role of virus-encoded microRNAs in herpesvirus biology. Trends in Microbiol. 2009;17(12):544–53.
Alvarezgarcia I, Miska EA. MicroRNA functions in animal development and human disease. Develop. 2005;132(21):4653–62.
Cullen BR. Viral and cellular messenger RNA targets of viral microRNAs. Nature. 2009;457(7228):421–5.
Chen Z, Zeng H, Guo Y, Liu P, Pan H. Deng a, et al. miRNA-145 inhibits non-small cell lung cancer cell proliferation by targeting c-Myc. J Experi & Clin Cancer Res. 2010;29(1):1–10.
Huang S, He X. The role of microRNAs in liver cancer progression. Bri. J Cancer. 2011;104(2):235–40.
Kincaid RP, Sullivan CS. Virus-encoded microRNAs: an overview and a look to the future. PLoS Pathog. 2012;8(12):e1003018.
Chen CJ, Cox JE, Kincaid RP, Martinez A, Sullivan CS. Divergent MicroRNA Targetomes of closely related circulating strains of a polyomavirus. J Virol. 2013;87(20):11135–47.
Yao Y, Smith LP, Nair V, Watson M. An avian retrovirus uses canonical expression and processing mechanisms to generate viral microRNA. J Virol. 2014;88(1):2–9.
Kincaid RP, Burke JM, Cox JC, Villiers EMD, Sullivan CS. A human torque Teno virus encodes a MicroRNA that inhibits interferon signaling. PLoS Pathog. 2013;9(12):e1003818.
Kozomara A, Griffithsjones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42:68–73.
Yao Y, Nair V. Role of virus-encoded microRNAs in avian viral diseases. Viruses. 2014;6(3):1379–94.
Piedade D, Azevedo-pereira JM. The role of microRNAs in the pathogenesis of herpesvirus infection. Viruses. 2016;8(6):156.
Grey F. Role of microRNAs in herpesvirus latency and persistence. J Gen Virol. 2015;96(4):739–51.
Burnside J, Bernberg E, Anderson A, Lu C, Meyers BC, Green PJ, et al. Marek's disease virus encodes MicroRNAs that map to meq and the latency-associated transcript. J Virol. 2006;80(17):8778–86.
Yao Y, Zhao Y, Xu H, Smith LP, Lawrie CH, Watson M, et al. MicroRNA profile of Marek's disease virus-transformed T-cell line MSB-1: predominance of virus-encoded microRNAs. J Virol. 2008;82(8):4007–15.
Yao Y, Zhao Y, Xu H, Smith LP, Lawrie CH, Sewer A, et al. Marek's disease virus type 2 (MDV-2)-encoded microRNAs show no sequence conservation with those encoded by MDV-1. J Virol. 2007;81(13):7164–70.
Waidner LA, Morgan RW, Anderson AS, Bernberg EL, Kamboj S, Garcia M, et al. MicroRNAs of Gallid and Meleagrid herpesviruses show generally conserved genomic locations and are virus-specific. Virology. 2009;388(1):128–36.
Yao YX, Zhao YG, Smith LP, Watson M, Nair V. Novel microRNAs (miRNAs) encoded by herpesvirus of turkeys: evidence of miRNA evolution by duplication. J Virol. 2009;83(13):6969–73.
Rachamadugu R, Lee JY, Wooming A, Kong BW. Identification and expression analysis of infectious laryngotracheitis virus encoding microRNAs. Virus Genes. 2009;39(3):301–8.
Yao Y, Smith LP, Petherbridge L, Watson M, Nair V. Novel microRNAs encoded by duck enteritis virus. J Gen Virol. 2012;93(7):1530–6.
Wu Y, Cheng A, Wang M, Yang Q, Zhu D, Jia R, et al. Complete genomic sequence of Chinese virulent duck enteritis virus. J Virol. 2012;86(10):5965.
Wu Y, Cheng A, Wang M, Zhu D, Jia R, Chen S, et al. Comparative genomic analysis of duck enteritis virus strains. J Virol. 2012;86(24):13841–2.
Liu X, Liao S, Xu Z, Zhu L, Yang F, Guo W. Identification and analysis of the porcine MicroRNA in porcine cytomegalovirus-infected macrophages using deep sequencing. PLoS One. 2016;11(3):e0150971.
Liu X, Zhu L, Liao S, Xu Z, Zhou Y. The porcine microRNA transcriptome response to transmissible gastroenteritis virus infection. PLoS One. 2015;10(3):e0120377.
Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.
Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7(10):986–95.
Liu F, Zheng H, Tong W, Li GX, Tian Q, Liang C, et al. Identification and analysis of novel viral and host dysregulated MicroRNAs in variant pseudorabies virus-infected PK15 cells. PLoS One. 2016;11(3):e0151546.
Young MD, Wakeeld MJ. Smyth GK. Oshlack A. goseq: Gene Ontology testing for RNA-seq datasets; 2012.
Varkonyi-Gasic E, Hellens RP. Quantitative stem-loop RT-PCR for detection of microRNAs. Meth mole. Biol. 2011;744(9):145–57.
Bartel DP. MicroRNA target recognition and regulatory functions. Cell. 2009;136(2):215–33.
Frappier L. Regulation of herpesvirus reactivation by host microRNAs. J Virol. 2015;89(5):2456.
Morgan R, Anderson A, Bernberg E, Kamboj S, Huang E, Lagasse G, et al. Sequence conservation and differential expression of Marek's disease virus microRNAs. J Virol. 2008;82(24):12213–20.
Barth S, Pfuhl T, Mamiani A, Ehses C, Roemer K, Kremmer E, et al. Epstein-Barr virus-encoded microRNA miR-BART2 down-regulates the viral DNA polymerase BALF5. Nucleic Acids Res. 2008;36(2):666–75.
Seo GJ, Fink LH, O'Hara B, Atwood WJ, Sullivan CS. Evolutionarily conserved function of a viral microRNA. J Virol. 2008;82(20):9823–8.
Seo GJ, Chen CJ, Sullivan CS. Merkel cell polyomavirus encodes a microRNA with the ability to autoregulate viral gene expression. Virology. 2009;383(2):183–7.
Tang S, Bertke AS, Patel A, Wang K, Cohen JI, Krause PR. An acutely and latently expressed herpes simplex virus 2 viral microRNA inhibits expression of ICP34.5, a viral neurovirulence factor. Proc Natl Acad Sci U S A. 2008;105(31):10931–6.
Shuang Tang AP, Krause PR. Novel less-abundant viral MicroRNAs encoded by herpes simplex virus 2 latency-associated transcript and their roles in regulating ICP34.5 and ICP0 mRNAs. J Virol. 2009;83(3):1433–42.
Umbach JL, Kramer MF, Jurak I, Karnowski HW, Coen DM, Cullen BR. MicroRNAs expressed by herpes simplex virus 1 during latent infection regulate viral mRNAs. Nature. 2008;454(7205):780–3.
Tang S, Bertke AS, Patel A, Margolis TP, Krause PR. Herpes simplex virus 2 MicroRNA miR-H6 is a novel latency-associated transcript-associated MicroRNA, but reduction of its expression does not influence the establishment of viral latency or the recurrence phenotype. J Virol. 2011;85(9):4501–9.
Strassheim S, Stik G, Rasschaert D, Laurent S. mdv1-miR-M7-5p, located in the newly identified first intron of the latency-associated transcript of Marek’s disease virus, targets the immediate-early genes ICP4 and ICP27. J Gen Virol. 2012;93(8):1731–42.
Waidner LA, Burnside J, Anderson AS, Bernberg EL, German MA, Meyers BC, et al. A microRNA of infectious laryngotracheitis virus can downregulate and direct cleavage of ICP4 mRNA. Virology. 2011;411(1):25–31.
Brehm M, Samaniego LA, Bonneau RH, Deluca NA, Tevethia SS. Immunogenicity of herpes simplex virus type 1 mutants containing deletions in one or more α-genes: ICP4, ICP27, ICP22, and ICP0. Virology. 1999;256(2):258–69.
Liu C, Cheng A, Wang M, Chen S, Jia R, Zhu D, et al. Duck enteritis virus UL54 is an IE protein primarily located in the nucleus. Virol J. 2015;12(1):1–8.
Landais I, Pelton C, Streblow D, Defilippis V, Mcweeney S, Nelson JA. Human cytomegalovirus miR-UL112-3p targets TLR2 and modulates the TLR2/IRAK1/NFκB signaling pathway. PLoS Pathog. 2015;11(5):e1004881.
Hu X, Ye J, Qin A, Zou H, Shao H, Qian K. Both MicroRNA-155 and virus-encoded MiR-155 Ortholog regulate TLR3 expression. PLoS One. 2015;10(5):e0126012.
Stern-Ginossar N, Elefant N, Zimmermann A, Wolf DG, Saleh N, Biton M, et al. Host immune system gene targeting by a viral miRNA. Science. 2007;317(5836):376–81.
Abend JR, Ramalingam D, Kieffer-Kwon P, Uldrick TS, Yarchoan R, Ziegelbauer JM. Kaposi's sarcoma-associated herpesvirus MicroRNAs target IRAK1 and MYD88, two components of the toll-like receptor/interleukin-1R signaling Cascade, to reduce inflammatory-cytokine expression. J Virol. 2012;86(21):11663–74.
Liang D, Gao Y, Lin X, He Z, Zhao Q, Deng Q, et al. A human herpesvirus miRNA attenuates interferon signaling and contributes to maintenance of viral latency by targeting IKKε. Cell Res. 2011;21(5):793–806.
Chi JQ, Teng M, Yu ZH, Xu H, Su JW, Zhao P, et al. Marek′s disease virus-encoded analog of microRNA-155 activates the oncogene c-Myc by targeting LTBP1 and suppressing the TGF-β signaling pathway. Virology. 2015;476:72–84.
Yunhua L, Rui S, Xianzhi L, Deguang L, Qiang D, Ke L. Kaposi's sarcoma-associated herpesvirus-encoded microRNA miR-K12-11 attenuates transforming growth factor beta signaling through suppression of SMAD5. J Virol. 2012;86(3):1372–81.
Mansouri S, Pan Q, Blencowe BJ, Claycomb JM, Frappier L. Epstein-Barr virus EBNA1 protein regulates viral latency through effects on let-7 microRNA and dicer. J Virol. 2014;88(19):11166–77.
Li X, Lian L, Zhang D, Qu L, Yang N. Gga-miR-26a targets NEK6 and suppresses Marek's disease lymphoma cell proliferation. Poultry Sci. 2014;93(5):1097–105.
Lian L. Qu L, Chen Y, Lamont SJ, Yang N, Lian L, et al. a systematic analysis of miRNA transcriptome in Marek’s disease virus-induced lymphoma reveals novel and differentially expressed miRNAs. PLoS One. 2011;7(11):e51003.
Sun JZ, Wang J, Yuan D, Wang S, Li Z, Yi B, et al. Cellular microRNA miR-181b inhibits replication of mink enteritis virus by repression of non-structural protein 1 translation. PLoS One. 2013;8(12):e81515.
Tian F, Luo J, Zhang H, Chang S, Song J. MiRNA expression signatures induced by Marek's disease virus infection in chickens. Genomics. 2012;99(3):152–9.
Taganov KD, Boldin MP, Chang KJ, Baltimore D. NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc Natl Acad Sci U S A. 2006;103(33):12481–6.
Yang JS, Li BJ, Lu HW, Chen Y, Lu C, Zhu RX, et al. Serum miR-152, miR-148a, miR-148b, and miR-21 as novel biomarkers in non-small cell lung cancer screening. Tumour Biology the Journal of the International Society for Oncodevelopmental Biology & Medicine. 2014;36(4):3035–42.
We would like to thank Jue Wang (BGI Genomics Co, Ltd., China) for technical assistance.
This work was supported by National Science and Technology Support Program (2015BAD12B05), Integration and Demonstration of Key Technologies for Duck Industrial in Sichuan Province (2014NZ0030), China Agricultural Research System (CARS-43-8), and Sichuan Province Research Programs (2017JY0014/2014–002).
Availability of data and materials
The datasets used or analysed during this study are included within the article.
Ethics approval and consent to participate
The study was approved by the Animal Ethics Committee of Sichuan Agricultural University (approval No. XF2014-18).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primers used to amplify virus and host miRNAs by stem-loop RT-qPCR. Stem-loop RT-qPCR was conducted using miRNA specific stem-loop RT primers together with corresponding miRNA specific forward (F) primers and universal reverse (UR) primer. (DOCX 15 kb)
Table S2. Distribution of sRNAs in DEV-infected and uninfected samples. (DOCX 15 kb)
The expression profiling of CHv miRNAs and pre-miRNA secondary structures. MiRNA sequences and their coresponding reads mapped on the precursors of CHv miRNA genes. Opening parentheses indicate pairing nucleotides. Inside the closed parentheses indicate the minimum free energy for the secondary structure of the miRNA. The number of reads mapped to the miRNA precursors is indicated in the right side. Mature miRNAs are denoted in red. The hairpin structures of pre-miRNA is shown at the back. (PDF 649 kb)
Figure S2. Regulatory network of CHv-encoded miRNAs and CHv genes. a Gene regulatory network formed by CHv-encoded miRNAs (red ellipses) and their target genes (yellow ectangles). b Gene regulatory network of CHv-encoded miRNAs (red ellipses) and their target immediate-early(IE) genes (yellow rectangles) (PDF 401 kb)
Table S3. Predicted host target genes of DEV miRNAs. Host target genes were predicted using PITA and RNAhybrid software. (XLSX 832 kb)
Table S4. Gene ontology analysis on host targets of the viral miRNAs (XLSX 23 kb)
Table S5. Novel duck miRNA expressed in DEF cells. The 386 mature sequence and sequence read count of novel duck miRNAs predicted with miRDeep in each sequenced sample of DEF cells. (XLSX 26 kb)
Table S6. Novel duck miRNA expressed in DEF cells. The 212 mature sequence and sequence read count of novel duck miRNAs predicted with miRDeep in each sequenced sample of DEF cells. (XLSX 15 kb)
Table S7. Differentially expressed host alignment miRNA in CHv-infected DEF cells. Fold-change = log2 (infected/mock in expression) > 1 or < − 1, and p-value < 0.05 indicated significance differentially expressed miRNA. Fold-change < − 1 indicates down-regulated, Fold-change > 1 indicates up-regulated. (XLSX 12 kb)
Figure S3. Regulatory network of DEF miRNAs and CHv genes. a Gene regulatory network formed by differentially-expressed DEF miRNAs (blue ellipses) and their target genes (yellow rectangles). b Gene regulatory network of differentially-expressed DEF miRNAs (blue circles) and target immediate-early (IE) genes (yellow rectangles). (PDF 353 kb)
About this article
Cite this article
Wu, X., Jia, R., Zhou, J. et al. Virulent duck enteritis virus infected DEF cells generate a unique pattern of viral microRNAs and a novel set of host microRNAs. BMC Vet Res 14, 144 (2018). https://doi.org/10.1186/s12917-018-1468-2