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.

Brucella-caused brucellosis is one of the most widespread worldwide zoonoses. Lipopolysaccharide (LPS) of Brucella, which functions as pathogen-associated molecular patterns (PAMPs), is an important virulence factor that elicits protective antibodies. Per of B. melitensis is involved in the biosynthesis of the O-side chain of LPS. Autophagy is a crucial element of the innate immune response against intracellular pathogens including Brucella. In this study, we observed that autophagy was inhibited in RAW264.7 cells infected with Brucella melitensis ∆per. And, a high-throughput array-based screen and qRT-PCR validation were performed to identify the differentially expressed miRNAs in RAW264.7 cells infected with B. melitensis M5-90 ∆per. The results suggested that mmu-miR-146a-5p, mmu-miR-155-5p, mmu-miR-146b-5p, and mmu-miR-3473a were upregulated and mmu-miR-30c-5p was downregulated. During B. melitensis M5-90 ∆per infection, the increased expression of miR-146b-5p inhibited the autophagy activation in RAW264.7 cells. Using a bioinformatics approach, Tbc1d14 was predicted to be a potential target of miR-146b-5p. The results of a luciferase reporter assay indicated that miR-146b-5p directly targeted the 3'-UTR of Tbc1d14, and the interaction between miR-146b-5p and the 3'-UTR of Tbc1d14 was sequence-specific. High-throughput RNA-Seq-based screening was performed to identify differentially expressed genes in Tbc1d14-expressing RAW264.7 cells, and these were validated by qRT-PCR. Among the differentially expressed genes, four autophagy associated genes, IFNγ-inducible p47 GTPase 1 (IIGP1), nuclear receptor binding protein 2 (Nrbp2), transformation related protein 53 inducible nuclear protein 1 (Trp53inp1), and immunity-related GTPase family M member 1 (Irgm1), were obtained. Our findings provide important insights into the functional mechanism of LPS of B. melitensis.

The immune-excluded tumors (IETs) show limited response to current immunotherapy due to intrinsic and adaptive immune resistance. In this study, it is identified that inhibition of transforming growth factor-β (TGF-β) receptor 1 can relieve tumor fibrosis, thus facilitating the recruitment of tumor-infiltrating T lymphocytes. Subsequently, a nanovesicle is constructed for tumor-specific co-delivery of a TGF-β inhibitor (LY2157299, LY) and the photosensitizer pyropheophorbide a (PPa). The LY-loaded nanovesicles suppress tumor fibrosis to promote intratumoral infiltration of T lymphocytes. Furthermore, PPa chelated with gadolinium ion is capable of fluorescence, photoacoustic and magnetic resonance triple-modal imaging-guided photodynamic therapy, to induce immunogenic death of tumor cells and elicit antitumor immunity in preclinical cancer models in female mice. These nanovesicles are further armored with a lipophilic prodrug of the bromodomain-containing protein 4 inhibitor (i.e., JQ1) to abolish programmed death ligand 1 expression of tumor cells and overcome adaptive immune resistance. This study may pave the way for nanomedicine-based immunotherapy of the IETs.

Melioidosis is a severe and fatal tropical zoonosis, which is triggered by Burkholderia pseudomallei. To better understand the host's response to infection of B. pseudomallei, an RNA-Seq technology was used to confirm differentially expressed genes (DEGs) in RAW264.7 cells infected with B. pseudomallei. In total, 4668 DEGs were identified across three time points (4, 8, and 11 hours after infection). Short Time-Series Expression Miner (STEM) analysis revealed the temporal gene expression profiles and identified seven significant patterns in a total of 26 profiles. Kyoto Encyclopedia of Genes and Genomes (KEGG) was utilized to confirm significantly enriched immune process-associated pathways, and 10 DEGs, including Ccl9, Ifnb1, Tnfα, Ptgs2, Tnfaip3, Zbp1, Ccl5, Ifi202b, Nfkbia, and Nfkbie, were mapped to eight immune process-associated pathways. Subsequent quantitative real-time PCR assays confirmed that the 10 DEGs were all upregulated during infection. Overall, the results showed that B. pseudomallei infection can initiate a time-series upregulation of immune process-associated DEGs in RAW264.7 macrophage cells. The discovery of this article helps us better understand the biological function of the immune process-associated genes during B. pseudomallei infection and may aid in the development of prophylaxis and treatment protocols for melioidosis.

Brucellosis is a worldwide zoonosis caused by Brucella species and represents a serious threat to both human and animal health. Omp25 is an important immunogenic and protective antigen in Brucella species; however, the functional mechanism of Omp25 in macrophages has not yet been elucidated. Here, we constructed a Brucella melitensis omp25 deletion mutant (M5-90-Δ omp25) and performed microRNA (miRNA) profiling of infected RAW264.7 cells. Eight differentially expressed miRNAs ( mmu-miR-146a-5p, mmu-miR-155-5p, mmu-miR-3473a, mmu-miR-149-3p, mmu-miR-671-5p, mmu-miR-1224-5p, mmu-miR-1895, and mmu-miR-5126) were identified, with quantitative real-time PCR (qRT-PCR) analysis confirming the up-regulation of mmu-miR-146-a-5p and mmu-miR-155-5p and down-regulation of mmu-miR-149-3p and mmu-miR-5126. mRNA profiling of B. melitensis M5-90-Δo mp25-infected RAW264.7 cells identified 967 differentially expressed genes (DEGs) (fold change ≥ 2). Among these, we focused on genes that were predicted by TargetScan, miRanda, and PicTar to be the potential targets of the differentially expressed miRNAs. The results suggested that 17 separate genes are potentially targeted by mmu-miR-149-3p, with one of these genes, Tbr1, also targeted by mmu-miR-5126. qRT-PCR analysis confirmed the up-regulation of nine of the predicted target genes. Our findings provide important information about the functional molecules in host cells, including miRNA and their target genes, affected by Omp25 from Brucella. This information is particularly valuable for the prophylaxis and treatment of brucellosis.

Pasteurella multocida (P. multocida) is a zoonotic bacterium that predominantly colonizes the respiratory tract and lungs of a variety of farmed and wild animals, and causes severe respiratory disease. To investigate the characteristics of the host immune response induced by P. multocida strains of serotype A and D, high-throughput mRNA-Seq and miRNA-Seq were performed to analyze the changes in goat bronchial epithelial cells stimulated by these two serotypes of P. multocida for 4 h. Quantitative RT-PCR was used to validate the randomly selected genes and miRNAs. The results revealed 204 and 117 differentially expressed mRNAs (|log2(Fold-change)| ≥ 1, p-value < 0.05) in the P. multocida serotype A and D stimulated groups, respectively. Meanwhile, the number of differentially expressed miRNAs (|log2(Fold-change)| > 0.1, p-value < 0.05) were 269 and 290, respectively. GO and KEGG enrichment analyses revealed 13 GO terms (p-value < 0.05) and four KEGG pathways (p-value < 0.05) associated with immunity. In the serotype A-stimulated group, the immune-related pathways were the GABAergic synapse and Toll-like receptor signaling pathways, while in the serotype D-stimulated group, the immune-related pathways were the phagosome and B cell receptor signaling pathways. Based on the predicted results of TargetScan and miRanda, the differentially expressed mRNA-miRNA network of immune-related GO terms and KEGG pathways was constructed. According to the cell morphological changes and the significant immune-related KEGG pathways, it was speculated that the P. multocida serotype D strain-stimulated goat bronchial epithelial cells may induce a cellular immune response earlier than serotype A-stimulated cells. Our study provides valuable insight into the host immune response mechanism induced by P. multocida strains of serotype A and D.

Swine hepatitis E (SHE) is a new type of zoonotic infectious disease caused by swine hepatitis E virus (SHEV). Open reading frame 3 (ORF3) is a key regulatory and virulent protein of SHEV. Circular RNAs (circRNAs) are a special kind of non-coding RNA molecule, which has a closed ring structure. In this study, to identify the circRNA profile in host cells affected by SHEV ORF3, adenovirus ADV4-ORF3 mediated the overexpression of ORF3 in HepG2 cells, whole genome sequencing was used to investigate the differentially expressed circRNAs, GO and KEGG were performed to enrichment analyze of differentially expressed circRNA-hosting gene, and Targetscan and miRanda softwares were used to analyze the interaction between circRNA and miRNA. The results showed adenovirus successfully mediated the overexpression of ORF3 in HepG2 cells, 1,105 up-regulation circRNAs and 1,556 down-regulation circRNAs were identified in ADV4-ORF3 infection group compared with the control. GO function enrichment analysis of differentially expressed circRNAs-hosting genes classified three main categories (cellular component, biological process and molecular function). KEGG pathway enrichment analysis scatter plot showed the pathway term of top20. The circRNAs with top10 number of BS sites for qRT-PCR validation were selected to confirmed, the results indicated that the up-regulated hsa_circ_0001423 and hsa_circ_0006404, and down-regulated of hsa_circ_0004833 and hsa_circ_0007444 were consistent with the sequencing data. Our findings first preliminarily found that ORF3 protein may affect triglyceride activation (GO:0006642) and riboflavin metabolism (ko00740) in HepG2 cells, which provides a scientific basis for further elucidating the effect of ORF3 on host lipid metabolism and the mechanism of SHEV infection.

Pasteurella multocida (P. multocida) is an important zoonotic pathogen. In addition to lung lesions, necropsies have revealed macroscopic lesions in the heart in clinical cases. However, most previous studies focused on lung lesions while ignoring heart lesions. Therefore, to investigate the immune response of the P. multocida-infected heart, two murine infection models were established by using P. multocida serotype A (Pm HN02) and D (Pm HN01) strains. Histopathological examination revealed heterogeneous inflammatory responses, including immune cell infiltration in the epicardial and myocardial areas of the heart. Transcriptome sequencing was performed on infected cardiac tissues. To explore the traits of immune responses, we performed the functional enrichment analysis of differentially expressed genes, gene set enrichment analysis and gene set variation analysis. The results showed that the innate immune pathways were significantly regulated in both groups, including the NOD-like receptor signaling pathway, the complement and coagulation cascade and cytokine-cytokine receptor interaction. The Toll-like receptor signaling pathway was only significantly activated in the Pm HN02 group. For the Pm HN02 group, immunohistochemistry analysis further verified the significant upregulation of the hub component MyD88 at the protein level. In conclusion, this study reveals critical pathways for host heart recognition and defense against P. multocida serotypes A and D. Moreover, MyD88 was upregulated by P. multocida serotype A in the heart, providing a theoretical basis for future prevention, diagnosis and treatment research.

Lipopolysaccharide (LPS) stimulates macrophages to release proinflammatory cytokines. MicroRNAs (miRNAs) are short noncoding RNAs that are involved in inflammatory reaction. Our previously study identified the downregulated expression of mmu-miR-27a-5p in RAW267.4 cells treated with LPS. To dissect the mechanism that mmu-miR-27a-5p regulates target genes and affects proinflammatory cytokine secretion more clearly, based on previous bioinformatics prediction data, one of the potential target genes, MCPIP1 was observed to be upregulated with qRT-PCR and western blot. Luciferase reporter assays were performed to further confirm in silico prediction and determine that MCPIP1 is the target of mmu-miR-27-5p. The results suggested that mmu-miR-27a-5p directly targeted the 3'-UTR of MCPIP1 and the interaction between mmu-miR-27-5p and the 3'-UTR of MCPIP1 is sequence-specific. MCPIP1 overexpression decreased the secretion of IL-6, IL-1β, and IL-10 in macrophage cells stimulated with LPS. Our findings may provide the important information for the precise roles of mmu-miR-27a-5p in the macrophage inflammatory response to LPS stimulation in the future.

Hepatitis E virus- (HEV-) mediated hepatitis has become a global public health problem. An important regulatory protein of HEV, ORF3, influences multiple signal pathways in host cells. In this study, to investigate the function of ORF3 from the swine form of HEV (SHEV), high-throughput RNA-Seq-based screening was performed to identify the differentially expressed genes in ORF3-expressing HepG2 cells. The results were validated with quantitative real-time PCR and gene ontology was employed to assign differentially expressed genes to functional categories. The results indicated that, in the established ORF3-expressing HepG2 cells, the mRNA levels of CLDN6, YLPM1, APOC3, NLRP1, SCARA3, FGA, FGG, FGB, and FREM1 were upregulated, whereas the mRNA levels of SLC2A3, DKK1, BPIFB2, and PTGR1 were downregulated. The deregulated expression of CLDN6 and FREM1 might contribute to changes in integral membrane protein and basement membrane protein expression, expression changes for NLRP1 might affect the apoptosis of HepG2 cells, and the altered expression of APOC3, SCARA3, and DKK1 may affect lipid metabolism in HepG2 cells. In conclusion, ORF3 plays a functional role in virus-cell interactions by affecting the expression of integral membrane protein and basement membrane proteins and by altering the process of apoptosis and lipid metabolism in host cells. These findings provide important insight into the pathogenic mechanism of HEV.

Diabetic nephropathy (DN) is one of the most serious complications in the development of diabetes mellitus, which has become the main cause of end-stage renal disease and one of the main causes of death in diabetic patients. With the prevalence of diabetes, the number of patients at risk for developing DN is increasing, with 20-40 percent of all patients with diabetes at risk for developing DN. Acupuncture and Chinese herbal medicine treatments are often combined to treat DN; however, there has been no meta-analysis on their synergistic effects. Therefore, we aimed to perform a systematic review and meta-analysis to estimate the effectiveness of acupuncture combined with Chinese herbal medicine for DN treatment.

Diabetic peripheral neuropathy (DPN) is one of the most common microvascular complications of diabetes mellitus, with an incidence ranging from 60% to 90%. With the change in modern dietary structure, the incidence of diabetes is increasing year by year, and DPN is also on the rise. Tuina therapy has been widely used in the treatment of DPN, but there is no systematic review on the treatment of DPN. Therefore, this study aimed to conduct a meta-analysis of Tuina in the treatment of DPN to clarify its efficacy.

Mycoplasma capricolum subsp. Capricolum (Mcc) is an important member of the Mycoplasma mycoides cluster (Mm cluster) and causes caprine contagious agalactia. Mcc can infect goats of all age groups, especially pregnant ewes and kids. It can cause the abortion in pregnant ewes and the death of goat kids, leading to enormous losses in the goat breeding industry. To date, the prevalence of epidemic Mcc strains on Hainan Island, China, remains unclear. This study aimed to isolate and identify Mcc strains endemic to Hainan Island, China. Genome sequencing and comparative genomic analysis were performed to reveal the molecular characteristics and evolutionary relationships of the isolated strain. Mcc HN-B was isolated and identified in Hainan Island, China. The Mcc HN-B genome consists of a 1,117,925 bp circular chromosome with a 23.79% G + C content. It contains 912 encoding genes, 3 gene islands, and 14 potential virulence genes. The core genome with the features of the Mm cluster and the specific genes of Mcc HN-B were identified by comparative genomic analysis. These results revealed the evolutionary relationship between Mcc HN-B and other members of the Mm cluster. Our findings provide a reference for further studies on the pathogenic mechanism and local vaccine development of Mcc.

Pasteurella multocida is a highly versatile pathogen that infects a wide range of animals, including goats, causing pneumonia and hemorrhagic septicemia. Circular RNA (circRNA) is a type of non-coding RNA that plays an important role in regulating cellular metabolism. However, whether and how circRNA is involved in regulating immune responses in the goat lung has not been reported. Thus, this study was designed to examine the function of circRNA in goats infected with Pasteurella multocida. Goats were assigned into one of two groups: an uninfected control group (CK) and an infected group challenged with P. multocida. Compared with the CK group, which remained healthy, the infected goats showed clinical signs of infection, including depression, cough, nasal discharge, and dyspnea, along with elevated body temperature and lesions in the lung. Whole-transcriptome sequencing and small RNA sequencing were then performed using lung samples from goats from each group. A total of 138 circRNA, 56 microRNAs (miRNA), and 2,673 messenger RNA (mRNA) molecules were significantly differentially expressed in the P. multocida-infected group compared with the CK group. Randomly selected differentially expressed circRNA, miRNA, and mRNA molecules (n = 5 per group) were then validated by quantitative reverse-transcriptase polymerase chain reaction analysis. Gene ontology (GO) analysis of the source genes indicated that six immune-related terms were enriched among the differentially expressed cirRNA molecules, including inflammatory response, immune effector process, cell activation involved in immune response, cytokine-mediated signaling pathway, response to endogenous stimulus, and immune response. The corresponding circRNA molecules were then selected for construction of a competitive endogenous RNA network to identify networks that may be involved in the immune response to P. multocida infection. The results indicated that P. multocida HN01 may cause pneumonia and stimulate an immune response in goats via regulation of circRNA expression. This study presents the first comprehensive circRNA profile in response to P. multocida infection in goats, thus, providing a basis for understanding the function of circRNA in the host immune response to P. multocida infection.

(1) Background: Mannheimia haemolytica (M. haemolytica) is an opportunistic pathogen and is mainly associated with respiratory diseases in cattle, sheep, and goats. (2) Methods: In this study, a mouse infection model was established using a M. haemolytica strain isolated from goats. Histopathological observations were conducted on various organs of the mice, and bacterial load determination and RNA-seq analysis were specifically performed on the spleens of the mice. (3) Results: The findings of this study suggest that chemokines, potentially present in the spleen of mice following a M. haemolytica challenge, may induce the migration of leukocytes to the spleen and suppress the release of pro-inflammatory factors through a negative feedback regulation mechanism. Additionally, an interesting observation was made regarding the potential of hematopoietic stem/progenitor cells congregating in the spleen to differentiate into immune cells, which could potentially collaborate with leukocytes in their efforts to counteract M. haemolytica invasion. (4) Conclusions: This study revealed the immune regulation mechanism induced by M. haemolytica in the mouse spleen, providing valuable insights into host-pathogen interactions and offering a theoretical basis for the prevention, control, and treatment of mannheimiosis.

Amino acid catabolism is frequently executed in mitochondria; however, it is largely unknown how aberrant amino acid metabolism affects mitochondria. Here we report the requirement for mitochondrial saccharopine degradation in mitochondrial homeostasis and animal development. In Caenorhbditis elegans, mutations in the saccharopine dehydrogenase (SDH) domain of the bi-functional enzyme α-aminoadipic semialdehyde synthase AASS-1 greatly elevate the lysine catabolic intermediate saccharopine, which causes mitochondrial damage by disrupting mitochondrial dynamics, leading to reduced adult animal growth. In mice, failure of mitochondrial saccharopine oxidation causes lethal mitochondrial damage in the liver, leading to postnatal developmental retardation and death. Importantly, genetic inactivation of genes that raise the mitochondrial saccharopine precursors lysine and α-ketoglutarate strongly suppresses SDH mutation-induced saccharopine accumulation and mitochondrial abnormalities in C. elegans Thus, adequate saccharopine catabolism is essential for mitochondrial homeostasis. Our study provides mechanistic and therapeutic insights for understanding and treating hyperlysinemia II (saccharopinuria), an aminoacidopathy with severe developmental defects.

Orf is a zoonotic disease that has caused huge economic losses globally. Systematical analysis of dysregulated cellular micro RNAs (miRNAs) in response to Orf virus (ORFV) infection has not been reported. In the current study, miRNA sequencing and RNA sequencing (RNA-seq) were performed in goat skin fibroblast (GSF) cells at 0, 18, and 30 h post infection (h.p.i). We identified 140 and 221 differentially expressed (DE) miRNAs at 18 and 30 h.p.i, respectively. We also identified 729 and 3961 DE genes (DEGs) at 18 and 30 h.p.i, respectively. GO enrichment analysis indicates enrichment of apoptotic regulation, defense response to virus, immune response, and inflammatory response at both time points. DE miRNAs and DEGs with reverse expression were used to construct miRNA-gene networks. Seven DE miRNAs and seven DEGs related to "negative regulation of viral genome replication" were identified. These were validated by RT-qPCR. Cfa-let-7a, a significantly upregulated miRNA, was found to repress Thrombospondin 1 (THBS1) mRNA and protein expression by directly targeting the THBS1 3' untranslated region. THBS1 has been reported to induce apoptosis; therefore, the cfa-let-7a-THBS1 axis may play an important role in cellular apoptosis during ORFV infection. This study provides new insights into ORFV and host cell interaction mechanisms.

K. pneumoniae is an opportunistic pathogen that leads to widespread infection in humans and animals, seriously threatening human health and animal husbandry development. In our research, we investigated the biological characteristics of the isolate by using a 16S rRNA gene sequencing, biochemical assay, and drug sensitivity test. Moreover, the pathogenicity study, including the bacteria load determination, the histopathology examination, and the RNA sequencing was carried out to explore whether the isolate could cause lung injury in mice through intraperitoneal injection. The results indicated that the isolate was identified as K. pneumoniae and named as KPHN001. The drug susceptibility test showed that KPHN001 was only sensitive to polymyxin B and furazolidone, and was resistant to other 28 antibiotics. In the bacteria load determination, the highest bacterial load of the organs was found in the spleen, and abundant bacterial colonization was also found in the lung. The histopathology showed the mainly acute inflammations in the lung were due to congestion, edema, and exudation. RNA-seq analysis revealed that the differentially expressed genes (DEGs) of inflammatory cytokines and chemokines were expressed massively in mice. In the present research, the biological characteristics and pathogenesis of clinically isolated K. pneumoniae were systematically studied, revealing the pathogenic mechanism of K. pneumoniae to animals, and providing a theoretical basis for the following prevention, control, and diagnosis research.

Delta bovine papillomaviruses (δBPVs) causes fibropapillomas or bladder cancer in cattle. E5 is the major oncogene of δBPVs; however, the influence that E5 oncogene has on host microRNA (miRNA) and mRNA expression profiles remains little elucidated. In the present study, small RNA sequencing and RNA sequencing were used to explore alterations in miRNAs and mRNAs in E5 over-expressing Madin-Darby bovine kidney (MDBK) cells compared with controls. In total, 77 miRNAs (including 30 bovine-derived miRNAs) and 223 genes were differentially expressed (DE) following E5 overexpression. The dysregulated genes were mainly involved in metabolic and biosynthetic processes. We constructed a potential miRNA-gene regulatory network from the differentially expressed genes (DEGs) and DE miRNAs. Finally, 22 DEGs and nine DE miRNAs were selected for RT-qPCR validation. Of these, downregulation of six miRNAs, bta-miR-34c, bta-miR-122, bta-miR-195, bta-miR-449b, bta-miR-2425-5p, and bta-miR-2428-3p were confirmed; In addition, upregulation of 16 genes, ACSS2, DDIT4, INHBE, INSIG1, PNRC1, PSAT1, PSPH, PYCR1, SC4MOL, SLC34A2, SCD, SPARC, IDI1, PCK2, HMGCS1, and SMIM14 and downregulation of two genes, BATF3 and WFDC2 were confirmed. Specially, bta-miR-34c and bta-miR-449b potentially regulated PYCR1 and DDIT4, which were involved in cancer progression and angiogenesis. Our study presented for the first time the comprehensive miRNA and mRNA alterations in MDBK cells expressing the BPV E5 oncogene, providing new insights into the tumorigenesis induced by BPV E5.