A new (old) bovine viral diarrhea virus 2 subtype: BVDV-2e

Prevalence and genetic characterization of bovine viral diarrhea virus in dairy cattle in northern China

BACKGROUND

Bovine viral diarrhea virus (BVDV) is one of the major viral pathogens responsible for respiratory disease complexes in cattle and other ruminants; it has spread worldwide and poses a significant threat to the cattle industry. To understand the prevalence and genetic diversity of BVDV in northern China, this study conducted an epidemiological survey of BVDV in dairy cows across 13 provinces in northern China from June 2022 to June 2024. A total of 2,199 nasal swab samples were analyzed by RT-PCR.

RESULTS

The results revealed an overall positive rate of 6.05% for BVDV, with values of 6.47% from June 2022 to June 2023 and 5.59% from July 2023 to June 2024. Notably, the positive rate varied by region, with the highest prevalence in Shandong (9.62%) and the lowest in Hebei (1.61%). Phylogenetic analysis of 53 positive samples revealed that all belonged to BVDV-1, with the predominant sub-genotypes being 1a (47.17%), 1 m (28.30%), and 1c (9.43%). No BVDV-2 or BVDV-3 was detected, indicating that BVDV-1a is the most prevalent strain in northern China. This study also highlighted the genetic diversity of BVDV, with nucleotide homology among the sub-genotypes ranging from 70.2 to 92.1%.

CONCLUSIONS

An epidemiological survey of BVDV conducted in 13 provinces (regions) in northern China between 2022 and 2024 revealed a positive rate of 6.05%. The prevalent genotype identified was BVDV-1, with the predominant sub-genotypes being BVDV-1a, BVDV-1 m, and BVDV-1c. The findings of this study provide new evidence for the molecular epidemiology and genetic evolution of BVDV transmission in northern China, laying a foundational basis for the development of vaccination and control strategies against BVDV.

Respiratory illness in young and adult cattle caused by bovine viral diarrhea virus subgenotype 2b in singular and mixed bacterial infection in a BVDV-vaccinated dairy herd

Bovine respiratory disease (BRD) is a common global health problem in dairy cattle. The definitive diagnosis of BRD is complex because its etiology involves several predisposing and determining factors. This report describes the etiology of a BRD outbreak in a dairy herd in the mesoregion of Central Eastern Paraná, which simultaneously affected young (calves and heifers) and adult (cows) Holstein-Friesian cattle. Nine biological samples, consisting of five lung samples from two cows and three suckling calves, and four nasal swab samples from heifers, were used for etiological diagnosis. The nucleic acids extracted from lung fragments and nasal swabs were subjected to PCR and RT-PCR assays for partial amplification of the genes of five viruses [bovine viral diarrhea virus (BVDV), bovine alphaherpesvirus 1 (BoAHV1), bovine respiratory syncytial virus (BRSV), bovine parainfluenza virus 3 (BPIV-3), and bovine coronavirus (BCoV)] and four bacteria (Mycoplasma bovis, Mannheimia haemolytica, Pasteurella multocida, and Histophilus somni) involved in the etiology of BRD. All nine biological samples from the animals with BRD tested negative for BoAHV1, BRSV, BPIV-3, BCoV, and H. somni. Therefore, the involvement of these microorganisms in the etiology of BRD outbreak can be ruled out. It was possible to identify the presence of BVDV and M. bovis in singular and mixed infections of the lower respiratory tract in cattle. BVDV was also identified in two nasal swabs: one as a single etiological agent and the other in association with two bacteria (P. multocida and M. haemolytica). The phylogenetic analysis conducted in the nucleotide sequence of the 5'UTR region and Npro gene of the BVDV amplicons demonstrated that the BVDV field strains of this BRD outbreak belong to subgenotype 2b. To the best of our knowledge, this is the first report of BVDV-2b involvement in the etiology of BRD in Brazil. Finally, it is necessary to highlight that the cattle were obtained from an open dairy herd with biannual vaccinations for BVDV-1a and - 2a.

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.

Outbreak of persistently infected heifer calves with bovine viral diarrhea virus subgenotypes 1b and 1d in a BVDV-vaccinated open dairy herd

Bovine viral diarrhea virus (BVDV) infection has a significant economic impact on beef and dairy industries worldwide. Fetal infection with a non-cytopathic strain may lead to the birth of persistently infected (PI) offspring, which is the main event in the epidemiological chain of BVDV infection. This report describes the birth of 99 BVDV-PI heifer calves within 52 days of birth in a regular BVDV-vaccinated Brazilian dairy cattle herd and the subgenotypes of the infecting field strains. This study was conducted in a high-yielding open dairy cattle herd that frequently acquired heifers from neighboring areas for replacement. The farm monitors the birth of PI calves by screening all calves born using an ELISA (IDEXX) for BVDV antigen detection. All calves aged 1-7 days were evaluated. For positive and suspected results, the ELISA was repeated when the calves were close to one month old. A total of 294 heifer calves were evaluated between February and March 2021. Of these, 99 (33.7 %) had positive ELISA results and were considered PI calves. To evaluate the predominant BVDV species and subgenotypes in this outbreak, whole blood samples were collected from 31 calves born during the study period. All samples were submitted to the RT-PCR assay for the partial amplification of the BVDV 5'-UTR region, and these amplicons were subjected to nucleotide sequencing. Phylogenetic analysis identified BVDV-1b and BVDV-1d in 16 and 13 heifer calves, respectively. In two calves, it was not possible to determine the BVDV-1 subgenotype. Detection of PI animals and monitoring of circulating BVDV subgenotype strains are central to disease control. This study shows that regular BVDV vaccination alone may be insufficient to prevent BVDV infection in high-yielding open dairy cattle herds. Other biosecurity measures must be adopted to avoid the purchase of cattle with acute infections by BVDV or BVDV-PI, which can cause a break in the health profile of the herd and economic losses.

Phylogenetic analysis of papillomaviruses in dogs from southern Brazil: molecular epidemiology and investigation of mixed infections and spillover events

Papillomaviruses (PVs) have been identified in several animal species, including dogs (canine papillomaviruses, CPVs) and cattle (bovine papillomaviruses, BPVs). Although some BPVs may occasionally infect species other than cattle, to the best of our knowledge, BPVs have not been reported in dogs to date. Herein, we carried out a retrospective phylogenetic study of PVs circulating in dogs from southern Brazil between 2017 and 2022, also investigating possible mixed infections and spillover events. For this, we screened 32 canine papilloma samples by PCR using the degenerate primers FAP59/64 and/or MY09/11, which amplify different regions of the L1 gene; the genomic target often used for PV classification/typing. Out these, 23 PV DNA samples were successfully amplified and sequenced. All PVs amplified by FAP59/64 (n = 22) were classified as CPV-1. On the other hand, PVs amplified by MY09/11 (n = 4) were classified as putative BPV-1. Among these, three samples showed mixed infection by CPV-1 and putative BPV-1. One of the putative BPV-1 detected in co-infected samples had the L1 gene full-sequenced, confirming the gene identity. Furthermore, the phylogenetic classifications from the FAP59/64 and/or MY09/11 amplicons were supported by a careful in silico analysis, which demonstrated that the analysis based on them matches to the classification from the complete L1 gene. Overall, we described CPV-1 circulation in southern Brazil over the years and the potencial BPV infection in dogs (potential spillover event), as well as possible CPV/1/BPV-1 co-infections. Finally, we suggest the analysis of the complete genome of the putative BPVs detected in dogs in order to deepen the knowledge about the PV-host interactions.

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.

Genomic evolution of bovine viral diarrhea virus based on complete genome and individual gene analyses

Bovine viral diarrhea virus (BVDV) genome consists of a single-stranded, positive-sense RNA with high genetic diversity. In the last years, significant progress has been achieved in BVDV knowledge evolution through phylodynamic analysis based on the partial 5'UTR sequences, whereas a few studies have used other genes or the complete coding sequence (CDS). However, no research has evaluated and compared BVDV evolutionary history based on the complete genome (CG), CDS, and individual genes. In this study, phylodynamic analyses were carried out with BVDV-1 (Pestivirus A) and BVDV-2 (Pestivirus B) CG sequences available on the GenBank database and each genomic region: CDS, UTRs, and individual genes. In comparison to the CG, the estimations for both BVDV species varied according to the dataset used, pointing out the importance of considering the analyzed genomic region when concluding. This study may provide new insight into BVDV evolution history while highlighting the need to increase the available BVDV CG sequences to perform more comprehensive phylodynamic studies in the future.

HoBi-like Pestivirus Is Highly Prevalent in Cattle Herds in the Amazon Region (Northern Brazil)

Pestiviruses are globally distributed and cause substantial economic losses to the cattle industry. In Brazil, the country with the world's largest cattle population, pestivirus infections are well described in some regions, such as in the south, where a high frequency of BVDV-2 is described and contrasts with the high prevalence of HoBi-like pestivirus (HoBiPeV) in the northeast. However, there is a lack of information about pestiviruses in the Amazon Region, in northern Brazil, with a cattle population estimated at 55.7 million head, which has a significant impact on the international livestock market. Therefore, this study investigated the seroprevalence and genetic variability of ruminant pestiviruses in 944 bovine serum samples from four states in northern Brazil: Pará (PA), Amapá (AP), Roraima (RR), and Amazonas (AM). Our results showed that 45.4% of the samples were seropositive (19.8% for BVDV-1, 14.1% for BVDV-2, and 20.9% for HoBiPeV). All samples were tested by RT-qPCR, and three were positive and classified as HoBiPeV in a phylogenetic analysis. These serological and molecular results contrast with those from other regions of the world, suggesting that the northern Brazilian states have a high prevalence of all bovine pestiviruses including HoBiPeV.

An appraisal of gene targets for phylogenetic classification of canine distemper virus: Is the hemagglutinin the best candidate?

Sequence analysis of the canine distemper virus (CDV) hemagglutinin (H) gene may provide important insights on virus-host interactions and has also been frequently used for CDV phylogenetic classification. Herein, we performed an in silico analysis of CDV complete genomes (CGs) available in GenBank in order to investigate the suitability of H for CDV classification into lineages/genotypes. In addition, we analyzed the other viral genes for their potential use in CDV classification. Initially, we collected 116 CDV CGs from GenBank and compared their phylogenetic classification with that of their respective H nucleotide (nt) and amino acid (aa) sequences. Subsequently, we calculated the geodesic distance between the CG and H phylogenetic trees. These analyses were later performed with other CDV genes. All CDV CGs were also evaluated for possible recombination events. Nucleotide and aa analyses of H misclassified some Vaccine/America 1/Asia 3 lineage sequences compared to CG analysis, finding supported by both Maximum Likelihood (ML) and Bayesian Markov Chain Monte Carlo (B-MCMC) methods. Moreover, aa-based H analysis showed additional disagreements with the classification obtained by CG. The geodesic distance between the H and CG trees was 0.0680. Strong recombination signals were identified in the H gene, including Vaccine/America 1/Asia 3 lineage sequences. In contrast, C and P were the only genes that fully reproduced the CG classification (by ML and/or B-MCMC) and that did not show strong recombination signals. Furthermore, the P phylogenetic tree showed the lowest geodesic distance from the CG tree (0.0369). These findings suggest C and P as potential targets for CDV phylogenetic classification, especially when full genome sequencing is not possible. Finally, since our results were obtained considering the CDV CGs available to date, future analyses performed as more CDV sequences become available will be useful to assess probable issues of H-based phylogeny and to consolidate the suitability of the C and P genes for CDV classification.

Phylogeny and amino acid analysis in single and mixed bovine papillomavirus infections in Southern Brazil, 2016-2020

Bovine papillomaviruses (BPVs) exhibit a high degree of genetic variability, and several viral types have been identified based on analysis of the L1 gene. The L1 is the main capsid protein and the main target for neutralizing antibodies. We performed a retrospective study on BPVs circulating in Rio Grande do Sul state, Southern Brazil, in 2016-2020. DNA from 43 bovine papilloma samples were amplified using two degenerate primer sets - FAP59/64 and MY09/11 - targeting the L1 region, and analyzed for phylogeny, mixed BPV infections (coinfections) and amino acid (aa) sequences. We also performed an in silico analysis with 114 BPV L1 sequences from the GenBank database to assess the agreement between the phylogeny obtained based on complete L1 sequences versus that based on the region amplified using the FAP59/64 and MY09/11 primer sets. Considering single and coinfections, we identified 31 BPV-1 (31/43; 72.1%), 27 BPV-2 (27/43; 62.8%) and 4 BPV-6 (4/43; 9.3%). Coinfections with BPV-1 and BPV-2 were observed in 61.3% of the samples. Our results are supported by in silico analyses that demonstrate that the classification using FAP59/64 or MY09/11 matches the complete L1 results, except for BPV-17 and -18, which may be mistakenly classified depending on the primers used. Furthermore, we found unique or rare amino acids in at least one L1 sequence of each BPV type identified in our study, some of which have been identified previously in papillomavirus epitopes, suggesting immune-mediated selection. Finally, our study provides an overview of BPVs circulating in Southern Brazil over the last five years and point to the combined use of primers FAP59/64 and MY09/11 for analysis of BPV coinfections and putative epitopes.