Identification and Analysis of Novel Viral and Host Dysregulated MicroRNAs in Variant Pseudorabies Virus-Infected PK15 Cells

Integrated analysis of microRNA and messenger RNA expression profiles reveals functional microRNA in infectious bovine rhinotracheitis virus-induced mitochondrial damage in Madin-Darby bovine kidney cells

BACKGROUND

Studies have confirmed that Infectious bovine rhinotracheitis virus (IBRV) infection induces mitochondrial damage. MicroRNAs (miRNAs) are a class of noncoding RNA molecules, which are involved in various biological processes and pathological changes associated with mitochondrial damage. It is currently unclear whether miRNAs participate in IBRV-induced mitochondrial damage in Madin-Darby bovine kidney (MDBK) cells.

RESULTS

In the present study, we used high-throughput sequencing technology, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to screen for mitochondria-related miRNAs and messenger RNAs (mRNAs). In total, 279 differentially expressed miRNAs and 832 differentially expressed mRNAs were identified in 6 hours (IBRV1) versus 24 hours (IBRV2) after IBRV infection in MDBK cells. GO and KEGG enrichment analysis revealed that 42 differentially expressed mRNAs and 348 target genes of differentially expressed miRNAs were correlated with mitochondrial damage, and the miRNA-mitochondria-related target genes regulatory network was constructed to elucidate their potential regulatory relationships. Among the 10 differentially expressed miRNAs, 8 showed expression patterns consistent with the high-throughput sequencing results. Functional validation results showed that overexpression of miR-10a and miR-182 aggravated mitochondrial damage, while inhibition of miR-10a and miR-182 alleviated mitochondrial damage.

CONCLUSIONS

This study not only revealed the expression changes of miRNAs and mRNAs in IBRV-infected MDBK cells, but also revealed possible biological regulatory relationship between them. MiR-10a and miR-182 may have the potential to be developed as biomarkers for the diagnosis and treatment of IBRV. Together, Together, these data and analyses provide additional insights into the roles of miRNA and mRNA in IBRV-induced mitochondria damage.

The Role of Latency-Associated Transcripts in the Latent Infection of Pseudorabies Virus

Pseudorabies virus (PRV) can cause neurological, respiratory, and reproductive diseases in pigs and establish lifelong latent infection in the peripheral nervous system (PNS). Latent infection is a typical feature of PRV, which brings great difficulties to the prevention, control, and eradication of pseudorabies. The integral mechanism of latent infection is still unclear. Latency-associated transcripts (LAT) gene is the only transcriptional region during latent infection of PRV which plays the key role in regulating viral latent infection and inhibiting apoptosis. Here, we review the characteristics of PRV latent infection and the transcriptional characteristics of the LAT gene. We also analyzed the function of non-coding RNA (ncRNA) produced by the LAT gene and its importance in latent infection. Furthermore, we provided possible strategies to solve the problem of latent infection of virulent PRV strains in the host. In short, the detailed mechanism of PRV latent infection needs to be further studied and elucidated.

Functional Analysis of a Frontal miRNA Cluster Located in the Large Latency Transcript of Pseudorabies Virus

MicroRNAs (miRNAs) have been identified as a class of crucial regulators of virus-host crosstalk, modulating such processes as viral replication, antiviral immune response, viral latency, and pathogenesis. Pseudorabies virus (PRV), a model for the study of alphaherpesvirus biology, codes for 11 distinct miRNAs mapped to the ~4.6 kb intron of Large Latency Transcript (LLT). Recent studies have revealed the role of clusters consisting of nine and eleven miRNA genes in the replication and virulence of PRV. The function of separate miRNA species in regulating PRV biology has not been thoroughly investigated. To analyze the regulatory potential of three PRV miRNAs located in the frontal cluster of the LLT intron, we generated a research model based on the constitutive expression of viral miRNAs in swine testis cells (ST_LLT [1–3] cell line). Using a cell culture system providing a stable production of individual miRNAs at high levels, we demonstrated that the LLT [1–3] miRNA cluster significantly downregulated IE180, EP0, and gE at the early stages of PRV infection. It was further determined that LLT [1–3] miRNAs could regulate the infection process, leading to a slight distortion in transmission and proliferation ability. Collectively, our findings indicate the potential of LLT [1–3] miRNAs to retard the host responses by reducing viral antigenic load and suppressing the expansion of progeny viruses at the early stages of infection.

Tandem Mass Tag-Based Quantitative Proteomic Analysis of ISG15 Knockout PK15 Cells in Pseudorabies Virus Infection

Pseudorabies virus (PRV) is recognized as one of the most important pathogens of swine and poses a serious threat to the swine industry worldwide. Available commercial vaccines fail to protect against the emergence of new PRV strains. Therefore, the new protein targets against PRV highlight the urgent need for uncovering the molecular determinants of host cellular proteins following PRV infection. Interferon-stimulated gene 15 (ISG15) demonstrates an outstanding antiviral response. However, the molecular mechanism of ISG15 that affects PRV replication is incompletely known. Here, we performed a tandem mass tag (TMT)-based approach to quantitatively identify protein expression changes in PRV-infected ISG15 knockout PK15 (ISG15^−/−-PK15) cells. In total, 4958 proteins were identified by using TMT coupled with LC-MS/MS in this study. In the PRV- and mock-infected groups, 241 differentially expressed proteins (DEPs) were identified, 162 upregulated and 79 downregulated proteins at 24 h post-infection (hpi), among which AFP, Vtn, Hsp40, Herc5, and Mccc1 may play important roles in PRV propagation. To ensure the validity and reliability of the proteomics data, the randomly selected DEPs were verified by RT-qPCR and Western blot analysis, and the results were consistent with the TMT results. Bioinformatics analyses further demonstrated that the DEPs are mainly involved in various biological processes and signaling pathways, such as signal transduction, the digestive system, and the PI3K-AKT pathway. These findings may provide new insight into molecular mechanisms for PRV infection, which is helpful for identifying potential protein targets for antiviral agents.

Long non-coding RNA LNC_000641 regulates pseudorabies virus replication

Long non-coding RNAs (lncRNAs) are a new arm of gene regulatory mechanism as discovered by sequencing techniques and follow-up functional studies. The lncRNAs regulation of pseudorabies virus (PRV) infection has rarely been reported so far. Using RNA sequencing analysis, 225 lncRNAs with significant altered expressions in 3D4/21 cells infected with PRV (ZJ01) were identified. Five lncRNAs upregulated in PRV-infected cells were verified in cells infected with different PRV strains by qRT-PCR. By down- and up-regulation of lnc641, the accelerating effect of lnc641 on PRV replication was confirmed. Furthermore, we found that lnc641 regulated PRV replication by inhibiting the JAK-STAT1 pathway. This study suggests that lnc641 could be a new host factor target for developing antiviral therapies against PRV infection.

Quantitative Proteomic Analysis Uncovers the Mediation of Endoplasmic Reticulum Stress-Induced Autophagy in DHAV-1-Infected DEF Cells

Autophagy is a tightly regulated catabolic process and is activated in cells in response to stress signals. Despite extensive study, the interplay between duck hepatitis A virus type 1 (DHAV-1) and the autophagy of host cells is not clear. In this study, we applied proteomics analysis to investigate the interaction mechanism between DHAV-1 and duck embryo fibroblast (DEF) cells. In total, 507 differentially expressed proteins (DEPs) were identified, with 171 upregulated proteins and 336 downregulated proteins. The protein expression level of heat shock proteins (Hsps) and their response to stimulus proteins and zinc finger proteins (ZFPs) were significantly increased while the same aspects of ribosome proteins declined. Bioinformatics analysis indicated that DEPs were mainly involved in the “response to stimulus”, the “defense response to virus”, and the “phagosome pathway”. Furthermore, Western blot results showed that the conversion of microtubule-associated protein 1 light chain 3-I (LC3-I) to the lipidation form of LC3-II increased, and the conversion rate decreased when DEF cells were processed with 4-phenylbutyrate (4-PBA). These findings indicated that DHAV-1 infection could cause endoplasmic reticulum (ER) stress-induced autophagy in DEF cells, and that ER stress was an important regulatory factor in the activation of autophagy. Our data provide a new clue regarding the host cell response to DHAV-1 and identify proteins involved in the DHAV-1 infection process or the ER stress-induced autophagy process.

Expression of pseudorabies virus-encoded long noncoding RNAs in epithelial cells and neurons

Long noncoding RNAs (lncRNAs) play important roles in regulating eukaryotic genome replication and gene expression in diverse biological systems. Here, we identified lncRNAs transcribed from pseudorabies virus (PRV)-infected PK-15 cells. Based on high-throughput sequencing data, we obtained 87,263,926 and 93,947,628 clean reads from mock-infected and PRV-infected PK-15 cells, respectively. Through a normalized analytic protocol, we identified three novel viral lncRNAs. According to an analysis of differential expression between the mock-infected and PRV-infected cells, 4151 host lncRNAs were significantly upregulated and 2327 host lncRNAs were significantly downregulated in the latter group. Viral lncRNAs and several host lncRNAs were verified by northern blotting and real-time PCR. The findings showed that the viral lncRNA LDI might regulate the expression of IE180, a potent transcriptional activator of viral genes. Furthermore, we characterized the expression of viral lncRNAs in a culture of infected primary chicken dorsal root ganglia (DRG). Collectively, the obtained data suggest that PRV generates lncRNAs in both epithelial cells and chick DRG neurons.

Utilizing host endogenous microRNAs to negatively regulate the replication of porcine reproductive and respiratory syndrome virus in MARC-145 cells

MicroRNAs (miRNAs) contribute to gene regulation at the post-transcriptional level and are capable of mRNA silencing by binding to target sites exhibiting high degrees of complementarity. Therefore, cloning host miRNA-recognition sequences into the genome of RNA viruses represents a rational strategy for manipulating viral replication. Here, we performed deep sequencing to obtain small-RNA (sRNA)-expression profiles from in vitro-cultured MARC-145 cells post infection with porcine reproductive and respiratory syndrome virus (PRRSV) and chose six candidate miRNAs of different abundance (miR-21, miR-140-3p, miR-185, miR-26a, miR-505, and miR-199a) for further study. Based on the full-length cDNA clone p7USC, we constructed a number of PRRSV mutants that provided complementary base-pairing target sites for the miRNAs in 3′ untranslated regions. Our results showed that all low- and moderate- abundant miRNA-target mutants showed similar growth properties, whereas the highest-abundant miRNA-target mutant blocked both viral transcription and replication. Discontinuous mutations in high-abundant miRNA-target sites subsequently recovered viral viability and propagation. These results demonstrated the copy number of endogenous miRNAs and the extent of sRNA complementarity were key factors to silence potential mRNA expression/translation, thereby determining PRRSV viability. Interestingly, the mutant containing miR-140-target sites (v140-t) showed strong suppression of viral replication from P1 to P3 in vitro, as shown by virus titer, plaque morphology, and qRT-PCR assays. To assess genetic stability, sequencing of v140-t (P1, P3, P5 and P10) revealed spontaneous mutations preferentially located among several nucleotides near the 3′ end of the insertion region and corresponding to the “seed region” of miR-140-3p, explaining the induced viral repression and the direction of virus evolution. This approach provided a general silencing strategy for limiting PRRSV replication by endogenous miRNAs in MARC-145 cells.

Virulent duck enteritis virus infected DEF cells generate a unique pattern of viral microRNAs and a novel set of host microRNAs

Background

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.

Results

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.

Conclusions

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.

Electronic supplementary material

The online version of this article (10.1186/s12917-018-1468-2) contains supplementary material, which is available to authorized users.

iTRAQ-based Proteomic Analysis of Porcine Kidney Epithelial PK15 cells Infected with Pseudorabies virus

Pseudorabies virus (PRV) is one of the most important pathogens of swine, resulting in severe economic losses to the pig industry. To improve our understanding of the host responses to PRV infection, we applied isobaric tags for relative and absolute quantification (iTRAQ) labeling coupled with liquid chromatography-tandem mass spectrometry to quantitatively identify the differentially expressed cellular proteins in PRV-infected PK15 cells. In total, relative quantitative data were identified for 4333 proteins in PRV and mock- infected PK15 cells, among which 466 cellular proteins were differentially expressed, including 234 upregulated proteins and 232 downregulated proteins. Bioinformatics analysis disclosed that most of these differentially expressed proteins were involved in metabolic processes, cellular growth and proliferation, endoplasmic reticulum (ER) stress response, cell adhesion and cytoskeleton. Moreover, expression levels of four representative proteins, beta-catenin, STAT1, GRB2 and PCNA, were further confirmed by western blot analysis. This is the first attempt to analyze the protein profile of PRV-infected PK15 cells using iTRAQ technology, and our findings may provide valuable information to help understand the host response to PRV infection.

MicroRNA analysis in mouse neuro-2a cells after pseudorabies virus infection

Pseudorabies virus (PRV), an alpha herpesvirus can enter the mammalian nervous system, causing Aujezsky's disease. Previous studies have reported an alteration of microRNA (miRNA) expression levels during PRV infections. However, knowledge regarding miRNA response in nervous cells to PRV infection is still unknown. To address this issue, small RNA libraries from infected and uninfected mouse neuroblastoma cells were assessed after Illumina deep sequencing. A total of eight viral miRNA were identified, and ten host miRNAs showed significantly different expression upon PRV infection. Among these, five were analyzed by stem-loop RT-qPCR, which confirmed the above data. Interestingly, these viral miRNAs were mainly found in the large latency transcript region of PRV, and predicted to target a variety of genes, forming a complicated regulatory network. Moreover, ten cellular miRNAs were expressed differently upon PRV infection, including nine upregulated and one downregulated miRNAs. Host targets of these miRNAs obtained by bioinformatics analysis belonged to large signaling networks, mainly encompassing calcium signaling pathway, cAMP signaling pathway, MAPK signaling pathway, and other nervous-associated pathways. These findings further highlighted miRNA features in nervous cells after PRV infection and contributed to unveil the underlying mechanisms of neurotropism as well as the neuropathogenesis of PRV.