Differential CircRNA Expression Profiles in PK-15 Cells Infected with Pseudorabies Virus Type II

Circular RNAs (circRNAs) belong to a class of non-coding RNAs with diverse biological functions. However, little is known about their roles in case of pseudorabies virus (PrV) infection. Here, we analyzed the expression profile of host circRNAs from a virulent PrV type II strain DX (PrV-DX) infected and an attenuated gE/TK deficient (gE^-TK^-PrV) strain of PrV infected PK-15 cells. CircRNAs were identified by find_circ and analyzed with DESeq 2. Compared with the mock cells, 449 differentially expressed (DE) circRNAs (233 down-regulated and 216 up-regulated) from PrV-DX infected and 578 DE circRNAs (331 down-regulated and 247 up-regulated) from gE^-TK^- PrV infected PK-15 cells were identified. In addition, 459 DE circRNAs (164 down-regulated and 295 up-regulated) between the PrV-DX and gE^-TK^-PrV infected cells were identified. The expression patterns of 13 circRNAs were validated by reverse transcription quantitative real-time PCR (RT-qPCR) and results were similar as of RNA-seq. The putative target miRNA binding sites of DE circRNAs were predicted by using miRanda and psRobot. The circRNA-miRNA-mRNA network was constructed and certain miRNAs that have possible roles in antiviral immune response, such as miR-210 and miR-340, were predicted. GO and KEGG pathway analysis demonstrated that DE circRNAs were enriched in the processes such as cellular metabolism, protein binding, RNA degradation and regulation of actin cytoskeleton. Collectively, these findings might provide the useful information for a better understanding of mechanisms underlying the interaction between PrV-II and host cells.