Bovine leukemia virus (BLV) is a retrovirus that causes enzootic bovine leukosis (EBL) and has worldwide distribution. Infections with BLV have been reported in cattle from Kazakhstan but the virus has not yet been thoroughly characterized. In this study, we detect and estimate the level of BLV proviral DNA by qPCR in DNA samples from 119 cattle naturally infected with BLV, from 18 farms located in four different geographical regions of Kazakhstan. Furthermore, we conducted the phylogenetic and molecular analysis of 41 BLV env-gp51 gene sequences from BLV infected cattle. Phylogenetic analysis showed the affiliation of sequences to two already known genotypes G4 and G7 and also to a new genotype, classified as genotype G12. In addition, a multivariate method was employed for analysis of the association between proviral load and different variables such as the geographical location of the herd, cattle breeds, age of animals, and the presence of particular BLV genotypes. In summary, the results of this study provide the first evidence on molecular characterization of BLV circulating in cattle from Kazakhstan.

Bovine viral diarrhea virus (BVDV) belongs to the Flaviviridae family and the Pestivirus genus. Infection with BVDV causes a disease with a wide spectrum of clinical symptoms, most often mild, although infections with this virus constitute a serious economic problem all over the world. The virus is characterized by a high genetic variability, while the accumulation of single mutations leads to the formation of its new variants. The aim of this study was to better understand the complicated pathogenesis of this disease at the molecular level via the analysis of the transcriptome of cells infected with this virus. The bovine kidney cell line (MDBK), the cytopathic (cp) reference strain, and two non-cytopathic (ncp) BVD virus field strains were used in transcriptomic studies. The cell transcriptome was tested 24 and 72 h after infection. The results of the microarray analysis revealed changes in the expression levels of numerous genes. Genes with changed expression as a result of infection with the cp strain caused changes in the expression levels of a large number of genes and enriched a number of pathways. Genes with increased expression levels were enriched among other pathways involved in the cell cycle, while genes with reduced expression levels enriched pathways mostly related to metabolism. Genes with increased expression levels as a result of infection with ncp strains enriched a much smaller number of pathways, among them, pathways related to signaling activity 24 h post-infection and serine biosynthetic pathways both 24 and 72 h post-infection. Pathways enriched by genes with reduced expression levels were related to the innate immune response (72 h post-infection) or metabolism (24 and 72 h post-infection). The results of microarray studies can help us to better understand the host’s response to BVDV infection.

Characterization of the global genetic diversity of the bovine leukemia virus (BLV) is an ongoing international research effort. Up to now BLV sequences have been classified into eleven distinct genotypes. Although BLV genotyping and molecular analysis of field isolates were reported in many countries, there is no report describing BLV genotypes present in cattle from Pakistan. In this study we examined 27 env gene sequences from BLV-infected cattle coming from four farms located in Khyber Pakhtunkwa, Gilgit Baltisan and Punjab provinces. Phylogenetic analyses revealed the classification of Pakistani sequences into genotypes G1 and G6. The alignment with the FLK-BLV sequence revealed the presence of 45 mutations, namely, seven in genotype G1 and 33 in genotype G6. Five mutations were found in both, G1 and G6 genotypes. Twelve amino acid substitutions were found in the analyzed sequences, of which only one P264S was specific for sequences from Pakistan. Furthermore, a certain degree of nucleotide heterogeneity was identified by NGS. These results highlight the need for further study on the importance of genetic variability of BLV, especially in the context of its pathogenicity and potential effect on serological detection.

Limited data are available on the incidence of variations in nucleotide sequences of long terminal repeat (LTR) regions of Bovine Leukemia Virus (BLV). Consequently, the possible impact of SNPs on BLV LTR function are poorly elucidated. Thus, a detailed and representative study of full-length LTR sequences obtained from sixty-four BLV isolates from different geographical regions of Poland, Moldova, Croatia, Ukraine and Russia were analyzed for their genetic variability. Overlap extension PCR, sequencing and Bayesian phylogenetic reconstruction of LTR sequences were performed. These analyses were followed by detailed sequence comparison, estimation of genetic heterogeneity and identification of transcription factor binding site (TFBS) modifications. Phylogenetic analysis of curated LTR sequences and those available in the GenBank database reflected the acknowledged env gene classification of BLV into 10 genotypes, and further clustered analysed sequences into three genotypes - G4, G7 and G8. Additional molecular studies revealed the presence of 97 point mutations distributed at 89 positions throughout all 64 LTR sequences. The highest rate of variability was noted in U3 and U5 subregions. However, the variability in regulatory sequences (VR) was assessed as lower than the variability within non-regulatory sequences (VNR) for both, U3 and U5 subregions. In contrast, VR value for R subregion, as well as for the total LTR, was higher than the VNR suggesting the existence of positive selection. Twelve unique SNPs for these LTR sequences localized in regulatory and non-regulatory elements were identified. The presence of different types of substitutions lead to the abrogation of present or to the creation of additional TFBS. This study represents the largest study of LTR genetic variability of BLV field isolates from Eastern part of Europe. Phylogenetic analysis of LTRs supports the clustering BLV variants based on their geographic origin. The SNP screening showed variations modifying LTR regulatory sequences, as well as altering TFBS. These features warrant further exploration as they could be related to proviral load and distinctive regulation of BLV transcription and replication.

Small ruminant lentiviruses (SRLVs) are a group of highly divergent viruses responsible for global infection in sheep and goats. In a previous study we showed that SRLV strains found in mixed flocks in Poland belonged to subtype A13 and A18, but this study was restricted only to the few flocks from Małopolska region. The present work aimed at extending earlier findings with the analysis of SRLVs in mixed flocks including larger numbers of animals and flocks from different part of Poland. On the basis of gag and env sequences, Polish SRLVs were assigned to the subtypes B2, A5, A12, and A17. Furthermore, the existence of a new subtypes, tentatively designed as A23 and A24, were described for the first time. Subtypes A5 and A17 were only found in goats, subtype A24 has been detected only in sheep while subtypes A12, A23, and B2 have been found in both sheep and goats. Co-infection with strains belonging to different subtypes was evidenced in three sheep and two goats originating from two flocks. Furthermore, three putative recombination events were identified within gag and env SRLVs sequences derived from three sheep. Amino acid (aa) sequences of immunodominant epitopes in CA protein were well conserved while Major Homology Region (MHR) had more alteration showing unique mutations in sequences of subtypes A5 and A17. In contrast, aa sequences of surface glycoprotein exhibited higher variability confirming type-specific variation in the SU5 epitope. The number of potential N-linked glycosylation sites (PNGS) ranged from 3 to 6 in respective sequences and were located in different positions. The analysis of LTR sequences revealed that sequences corresponding to the TATA box, AP-4, AML-vis, and polyadenylation signal (poly A) were quite conserved, while considerable alteration was observed in AP-1 sites. Interestingly, our results revealed that all sequences belonging to subtype A17 had unique substitution T to A in the fifth position of TATA box and did not have a 11 nt deletion in the R region which was noted in other sequences from Poland. These data revealed a complex picture of SRLVs population with ovine and caprine strains belonging to group A and B. We present strong and multiple evidence of dually infected sheep and goats in mixed flocks and present evidence that these viruses can recombine in vivo.