Genetic analyses of the structural protein E2 bovine viral diarrhea virus isolated from dairy cattle in Yogyakarta, Indonesia

Background and Aim: Bovine viral diarrhea (BVD), a highly pathogenic ribonucleic acid (RNA) virus, causes devastating financial losses and reproductive deaths among dairy cattle in Yogyakarta and globally. This study aimed to identify point mutations within the E2 structural protein of the acquired BVD virus (BVDV) isolates using genetic analysis. Materials and Methods: The study period shows that we performed the research in 2023. We collected 118 serum samples from 2019 to 2023, among which only 10 BVDV positive were used and 108 were negative lacking the BVDV antigen. An anti-Erns monoclonal antibody-coated protein was used in indirect antigen capture enzyme-linked immunosorbent assay (I-ACE) to detect the BVD antigen present in positive BVDV serum specimens. In the initial step of the two-step reverse transcription polymerase chain reaction, the enzyme (superscript III reverse transcriptase) and the primer (random hexamer) were used to convert the RNA of the BVDV into complementary deoxyribonucleic acid (cDNA) during the process of reverse transcription. The final step involved the amplification of the E2 gene of the resultant BVDV cDNA through gene-specific primers (E2_fwd: 5′-TGGTGGCCTTATGAGAC-3′ and P7_rev: 5′-CCCATCATCACTATTTCACC-3′) and enzyme (platinum taq DNA polymerase high fidelity). For conducting Sanger sequencing, those 3 BVDV-1-positive isolates (about 2.6% of all isolates) were selected as a typical specimen for each site and year between 2019 and 2023 using a proportional computation. Therefore, only two BVDV isolates with complete genomes were chosen to perform their homological and genetic analysis based on the E2 gene by means of Blast and MEGA Version 11 in addition to the Bioedit 7.2.5 program. Results: By applying phylogenetic analysis relying on the E2 gene, a sum of 1011 nucleotides of the BVDV-1 isolates derived from each of the two BVDV-1 Indonesian isolates (n = 2) and its 23 reference BVDV strains were acquired from the National Center for Biotechnology Information (NCBI) database. The findings of the genetic analysis inside the phylogenetic tree revealed that the two BVDV Indonesian isolates were clustered into BVDV-1a subgenotype, while the reference BVDV strains were clustered into the five BVDV subgenotype, BVDV-1a (n = 6), BVDV-1b (n = 3), BVDV-1c (n = 11), BVDV-1m (n = 1), and BVDV-1n (n = 2). The branch exists in phylogenetic tree located before the division of our two BVDV isolates was divided into two branches with the same maximum bootstrap values of 99%, indicating a high degree of confidence, was seen. Next, we observed the branch near our study samples, which displayed the bootstrap value of 100, indicating that our 02 isolates were identical. In both isolates, V11 BVDV1/Indonesia/Yogyakarta/2023 and V16 BVDV1/Indonesia/ Yogyakarta/2023 with GenBank accession numbers PP836388 and PP836389, respectively, conserved D7E residues were mutated as well as cysteine changed/altered into serine (S) was identified at amino acid position 201. Conclusion: We identified two isolates of BVDV belonging to the BVDV-1a subgenotype. Our findings indicate that the conserved D7E residues of isolates V11 BVDV1/Indonesia/Yogyakarta/2023 and V16 BVDV1/Indonesia/Yogyakarta/2023 were altered. The Indonesian BVDV isolates exhibited a cysteine to serine mutation at amino acid position 201, leads to vaccination failure, range of animal’s host will increase, and diagnostic kit will not be effective.