Orf, caused by Orf virus (ORFV), is a globally distributed zoonotic disease responsible for serious economic losses in the agricultural sector. However, the mechanism underlying ORFV infection remains largely unknown. Circular RNAs (circRNAs), a novel type of endogenous non-coding RNAs, play important roles in various pathological processes but their involvement in ORFV infection and host response is unclear. In the current study, whole transcriptome sequencing and small RNA sequencing were performed in ORFV-infected goat skin fibroblast cells and uninfected cells. A total of 151 circRNAs, 341 messenger RNAs (mRNAs), and 56 microRNAs (miRNAs) were differently expressed following ORFV infection. Four circRNAs: circRNA1001, circRNA1684, circRNA3127 and circRNA7880 were validated by qRT-PCR and Sanger sequencing. Gene ontology (GO) analysis indicated that host genes of differently expressed circRNAs were significantly enriched in regulation of inflammatory response, epithelial structure maintenance, positive regulation of cell migration, positive regulation of ubiquitin-protein transferase activity, regulation of ion transmembrane transport, etc. The constructed circRNA-miRNA-mRNA network suggested that circRNAs may function as miRNA sponges indirectly regulating gene expression following ORFV infection. Our study presented the first comprehensive profiles of circRNAs in response to ORFV infection, thus providing new clues for the mechanisms of interactions between ORFV and the host.

Delta bovine papillomaviruses (δBPVs) mainly infect cattle and cause fibropapillomas. δBPVs encode three oncogenes, E5, E6 and E7. The effect of E6 on microRNA (miRNA) and mRNA expression profiles is not well characterized. In this study, RNA sequencing and small RNA sequencing were used to explore alterations in mRNAs and miRNAs in E6 over-expressing NIH/3T3 cells (NH-E6) compared with control cells (NH-GFP). We found that 350 genes (181 upregulated and 169 downregulated) and 54 miRNAs (26 upregulated and 28 downregulated) were differentially expressed (DE) following E6 expression. The top 20 significantly enriched GO terms in "biological process" included inflammatory response, innate immune response, immune response, immune system process, positive regulation of apoptotic process, cell adhesion, and angiogenesis. We constructed a potential miRNA-gene regulatory network from the differentially expressed genes (DEGs) and DE miRNAs. Finally, we selected 19 immune-response related DEGs and 11 DE miRNAs for qPCR validation. Of these, upregulation of 12 genes, Ccl2, Ccl7, Cxcl1, Cxcl5, Tlr2, Nfkbia, Fas, Il1rl1, Ltbp1, Rab32, and Zc3h12a, Dclk1 and downregulation of four genes, Agtr2, Ptx3, Sfrp1, and Thbs1 were confirmed. Ccl2, Ccl7, Cxcl1 and Cxcl5 were upregulated more than ten-fold in NH-E6 compared with NH-GFP. Also, upregulation of three miRNAs, mmu-miR-129-2-3p, mmu-miR-149-5p-R-2 and mmu-miR-222-3p, and downregulation of five miRNAs, mmu-miR-582-3p-R+1, mmu-miR-582-5p, mmu-miR-708-3p, mmu-miR-708-5p and mmu-miR-1197-3p, were confirmed. Our study describes changes in both mRNA and miRNA profiles in response to BPV E6 expression, providing new insights into BPV E6 oncogene functions.