Abstract
Enzootic bovine leukosis, a neoplastic disease caused by the bovine leukemia virus (BLV), was the primary cancer affecting cattle in China before 1985. Although its prevalence decreased significantly between 1986 and 2000, enzootic bovine leukosis has been re-emerging since 2000. This re-emergence has been largely overlooked, possibly due to the latent nature of BLV infection or the perceived lack of sufficient evidence. This study investigated the molecular epidemiology of BLV infections in dairy cattle in Henan province, Central China. Blood samples from 668 dairy cattle across nine farms were tested using nested polymerase chain reaction assays targeting the partial envelope (env) gene (gp51 fragment). Twenty-three samples tested positive (animal-level prevalence of 3.4%; 95% confidence interval: 2.2, 5.1). The full-length env gene sequences from these positive samples were obtained and phylogenetically analyzed, along with previously reported sequences from the GenBank database. The sequences from positive samples were clustered into four genotypes (1, 4, 6, and 7). The geographical annotation of the maximum clade credibility trees suggested that the two genotype 1 strains in Henan might have originated from Japan, while the genotype 7 strain is likely to have originated from Moldova. Subsequent Bayesian stochastic search variable selection analysis further indicated a strong geographical association between the Henan strains and Japan, as well as Moldova. The estimated substitution rate for the env gene ranged from 4.39 × 10(-4) to 2.38 × 10(-3) substitutions per site per year. Additionally, codons 291, 326, 385, and 480 were identified as positively selected sites, potentially associated with membrane fusion, epitope peptide vaccine design, and transmembrane signal transduction. These findings contribute to the broader understanding of BLV epidemiology in Chinese dairy cattle and highlight the need for measures to mitigate further BLV transmission within and between cattle herds in China.