Expression profile of microRNAs in bovine lymphocytes infected with Theileria annulata and treated with buparvaquone.

ABSTRACT

Several miRNA-based studies on Theileria-transformed bovine cells have been conducted; however, the mechanism by which transformed cells exhibit uncontrolled proliferation is not yet fully understood. Therefore, it is necessary to screen more microRNAs that may play a role in the transformation process of host cells infected with Theileria annulata to better understand the transformation mechanisms of Theileria-infected cells. RNA sequencing was used to analyze miRNAs expression in the host bovine lymphocytes infected with T. annulata at different time points after buparvaquone (BW720) treatment and DMSO treatment (control groups). Differential miRNAs related to cell proliferation and apoptosis were identified through comparison with gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, and a regulatory network of miRNA-mRNA was constructed. In total, 272 differentially expressed miRNAs were found at 36, 60 and 72 h. The miRNAs change of bta-miR-2285t, novel-miR-622, bta-miR-2478, and novel-miR-584 were significant. Analysis of 27 of these co-differential expressed miRNAs revealed that 15 miRNAs were down-regulated and 12 miRNAs were up-regulated. A further analysis of the changes in the expression of each of these 27 miRNAs in the three datasets suggested that bta-miR-2285t, bta-miR-345-5p, bta-miR-34a, bta-miR-150, and the novel-miR-1372 had significantly changed. Predicted target genes for these 27 miRNAs were analyzed by KEGG and the results demonstrated that EZR, RASSF, SOCS1 were mainly enriched in the signaling pathway microRNAs in cancer. MAPKAPK2, RELB, FLT3LG, and GADD45B were mainly enriched in the MAPK signaling pathway, and some genes were enriched in Axon guidance. This study has provided valuable information to further the understanding of the regulatory function of miRNAs in the host microenvironment and host-parasite interaction mechanisms.