Bovine viral diarrhea virus (BVDV) significantly impacts the bovine industries, both dairy and beef sectors. BVDV can infect various domestic and wild animals, most notably cattle. The dynamic variations among BVDV serotypes due to the continuous genetic diversity, especially in BVDV1 (BVDV1), reduce the effectiveness of the currently available vaccines and reduce the specificity/sensitivity of the diagnostic assays. The development of novel, safe, and effective vaccines against BVDV requires deep knowledge of the antigenicity and virulence of the virus. Previous studies on the antigenicity and the virulence of BVDV serotypes have been mainly focused on one or a few BVDV proteins. While however, little is known about the orchestration of all BVDV in the context of viral virulence and immunogenicity. The main aim of the current study was to do a comparative computational evaluation of the immunogenicity, and virulence for all the encoded proteins of both BVDV1 and BVDV2 and their sub-genotypes.
To achieve this goal, 11,737 protein sequences were retrieved from Virus Pathogen Resource. The analysis involved a total of 4,583 sequences after the removal of short sequences and those with unknown collection time. We used the MP3 tool to map the pathogenic proteins across different BVDV strains. The potential protective and the epitope motifs were predicted using the VaxiJen and EMBOSS antigen tools, respectively.
The virulence prediction revealed that the NS4B proteins of both BVDV1 and BVDV2 likely have essential roles in BVDV virulence. Similarly, both the capsid (C) and the NS4-A proteins of BVDV1 and the Npro and P7 proteins of BVDV2 are likely important virulent factors. There was a clear trend of increasing predicted virulence with the progression of time in the case of BVDV1 proteins, but that was not the case for the BVDV2 proteins. Most of the proteins of the two BVDV serotypes possess antigens predicted immunogens except Npro, P7, and NS4B. However, the predicted antigenicity of the BVDV1 was significantly higher than that of BVDV2. Meanwhile, the predicted immunogenicity of the immunodominant-E2 protein has been decreasing over time. Based on our predicted antigenicity and pathogenicity studies of the two BVDV serotypes, the sub-genotypes (1a, 1f, 1k, 2a, and 2b) may represent ideal candidates for the development of future vaccines against BVDV infection in cattle. In summary, we identified some common differences between the two BVDV genotypes (BVDV1 and BVDV2) and their sub-genotypes regarding their protein antigenicity and pathogenicity. The data presented here will increase our understanding of the molecular pathogenesis of BVDV infection in cattle. It will also pave the way for developing some novel diagnostic assays and novel vaccines against BVDV in the near future.