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

• Published in: Veterinary World 2024-07, Vol.17 (7), p.1562-1574
• Main Authors: Khan, S U, Wuryastuty, Hastari, Wibowo, M H, Sarmin, Sarmin, Irianingsih, S H
• Format: Article
• Language: English
• Subjects: bovine viral diarrhea virus; cysteine mutation; e2 protein; serine; v11 bovine viral diarrhea virus1; v16 bovine viral diarrhea virus1
• ISSN: 0972-8988; 2231-0916
• DOI: 10.14202/vetworld.2024.1562-1574

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. 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. 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