Novel genomic targets for proper subtyping of bovine viral diarrhea virus 1 (BVDV-1) and BVDV-2

Whole-genome phylogenetic analysis, the most suitable strategy for subtyping bovine viral diarrhea virus 1 (BVDV-1) and BVDV-2, is not feasible for many laboratories. Consequently, BVDV isolates/strains have been frequently subtyped based on analysis of single genomic regions, mainly the 5’ untranslated region (UTR). This approach, however, may lead to inaccurate and/or poorly statistically supported viral classification. Herein, we describe novel primer sets whose amplicons may be easily sequenced and used for BVDV subtyping. Initially, genomic regions previously described as the most suitable targets for BVDV subtyping were analyzed for design of high-coverage primers. The putative amplicons were analyzed in silico for their suitability to reproduce the phylogenetic classification of 118 BVDV-1 and 88 BVDV-2 complete/near-complete genomes (CNCGs) (GenBank). This analysis was also performed considering the region amplifiable by primers HCV90-368, 324-326 and BP189-389 (5’UTR), which have been used for BVDV diagnosis and/or classification. After confirming the agreement between the analyses of our primers’ amplicon versus the CNCGs, we optimized the RT-PCRs and evaluated their performance for amplification of BVDV isolates/strains (n = 35 for BVDV-1; n = 33 for BVDV-2). Among the potential targets for BVDV subtyping, we designed high-coverage primers for NS3-NS4A (BVDV-1) (526 bp amplicon) and NS5B (BVDV-2) (728 bp). The classification based on these regions fully reproduced the subtyping of all CNCGs. On the other hand, subtyping based on the putative amplicons from primers HCV90-368, 324-326 and BP189-389 showed disagreements in relation the CNCG analysis. The NS3-NS4A and NS5B primers also allowed the amplification of all BVDV isolates/strains tested. Finally, we suggest the use of these primers in future phylogenetic and epidemiological studies of BVDVs.