Genoepidemiological evaluation of Bovine viral diarrhea virus 2 species based on secondary structures in the 5' untranslated region

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

Bovine pestiviruses are members of the species Pestivirus A (bovine viral diarrhea virus 1, BVDV-1), Pestivirus B (BVDV-2) or Pestivirus H (HoBiPeV). To date, BVDV-2 isolates/strains have been classified into three subtypes (a-c) by phylogenetic analysis, and an additional subtype (d) has been proposed based on 5' untranslated region (UTR) secondary structures. In a previous study, we identified some BVDV-2 sequences in the GenBank database that could not be classified as subtype a, b or c by phylogenetic analysis of their genomes, UTRs or individual genes. Here, we performed a detailed study of these sequences and assessed whether they might represent a distinct BVDV-2 subtype. Initially, we collected 85 BVDV-2 complete/near-complete genomes (CNCGs) from GenBank and performed a "proof of equivalence" between phylogenetic analyses based on CNCGs and open reading frames (ORFs), which showed that ORFs may be reliably used as a reference target for BVDV-2 phylogeny, allowing us to increase our dataset to 139 sequences. Among these, we found seven sequences that could not be classified as BVDV-2a-c. The same was observed in the phylogenetic analysis of CNCGs and viral genes. In addition, the seven non-BVDV-2a-c sequences formed a distinct cluster in all phylogenetic trees, which we propose to term BVDV-2e. BVDV-2e also showed 44 amino acid changes compared to BVDV-2a-c, 20 of which are in well-defined positions. Importantly, an additional phylogenetic analysis including BVDV-2d and a pairwise comparison of BVDV-2e and BVDV-2d sequences also supported the difference between these subtypes. Finally, we propose the recognition of BVDV-2e as a distinct BVDV-2 subtype and encourage its inclusion in future phylogenetic analyses to understand its distribution and evolution.

Genetic diversity of Bovine Viral Diarrhea Virus in cattle in France between 2018 and 2020

Bovine Viral Diarrhea Virus (BVDV) is one of the main pathogens that affects ruminants worldwide, generating significant economic losses. Like other RNA viruses, BVDV is characterized by a high genetic variability, generating the emergence of new variants, and increasing the risk of new outbreaks. The last report on BVDV genotypes in France was in 2008, since which there have been no new information. The goal of this study is to determine the genetic diversity of BVDV strains currently circulating in France. To this aim, samples of cattle were taken from different departments that are part of the main areas of livestock production during the years 2018 to 2020. Using the partial sequence of the 5'UTR region of the viral genome, we identified and classified 145 samples corresponding to Pestivirus A and one sample corresponding to Pestivirus D. For the Pestivirus A samples, the 1e, 1b, 1d, and 1l genotypes, previously described in France, were identified. Next, the 1r and 1s genotypes, not previously described in the country, were detected. In addition, a new genotype was identified and was tentatively assigned as 1x genotype. These results indicate an increase in the genetic diversity of BVDV in France.

Subtyping bovine viral diarrhea virus (BVDV): Which viral gene to choose?

Bovine viral diarrhea virus-1 (BVDV-1, Pestivirus A) and BVDV-2 (Pestivirus B) have been clustered into 21 and 4 subtypes, respectively. This genetic diversity, in addition to the lack of consensus on which genomic region to use for BVDV subtyping, has resulted in conflicting classifications depending on the target analyzed. Here, we investigated which genes or UTRs would reproduce the phylogeny obtained by complete genome (CG) analyses. The study was carried out with 91 (BVDV-1) and 85 (BVDV-2) CG available on GenBank database. The viruses were subtyped by analyzing their CG, as well as their individual genes and UTRs (complete 3' and 5'UTRs, and partial 5'UTR); and the phylogeny results were compared to each other. The sequences were aligned using the ClustalW multiple method (BioEdit Alignment Editor software, v.7.0.5.3) and the phylogenetic analyses were performed by the Maximum Likelihood method (MEGA-X software, v.10.2.4), with 1000 bootstrap replicates. The best analysis model for each gene/UTR was defined using the jModelTest software. The geodesic distance between the CG (reference) and individual genes/UTRs trees was also calculated (TreeCmp software, v.2.0). In general, 3'UTR-based analyses, followed by 5'UTR, presented the least reliable subtyping results. Regarding BVDV-1, phylogeny based on C, E^rns, E1, E2, p7, NS2, NS3, NS4B, NS5A and NS5B was consistent with that of CG. In contrast, analyses performed with individual BVDV-2 genes showed at least one different clustering from the phylogeny based on the CG. After analyzing the geodesic distance between the CG and genes/UTRs trees, we observed that NS4B (for BVDV-1) and NS5A (BVDV-2) presented the closest topology and edge length to the CG analyses. Finally, comparing the phylogeny performed with the CG and the genes/UTRs, as well as the geodesic distance between them, we understand that NS4B and NS5A represent the most suitable targets for BVDV-1 and -2 subtyping, respectively, and may be considered in future phylogenetic studies.

A genetic profile of bovine pestiviruses circulating in Brazil (1998-2018)

The pestiviruses bovine viral diarrhea virus 1 (BVDV-1), 2 (BVDV-2), and HoBi-like (HoBiPeV) are endemic among Brazilian cattle, the world's largest commercial bovine herd. In the last two decades (1998-2018) over 300 bovine pestiviruses have been partially or fully sequenced in Brazil, including viruses from different regions, different epidemiological backgrounds, and associated with diverse clinical presentations. Phylogenetic analysis of these viruses demonstrated a predominance of BVDV-1 (54.4%), with subgenotypes -1a (33.9% of total) and -1b (16.3%) being more frequent and subgenotypes -1d, -1e, and -1i at very low frequencies. The overall BVDV-2 frequency was 25.7% but it varied largely by region, reaching up to 48% in Southern states. BVDV-2b was the predominant subgenotype (84.8% of BVDV-2), followed by BVDV-2a (8.86%). HoBiPeV accounted for 19.9% (61/307) of the genotyped viruses and were detected at high frequency in cattle from Northeastern states. These findings demonstrate a unique mix of pestivirus species and subgenotypes, unlike that seen in Europe or North America. The design of effective diagnostic tools, vaccines, and control programs for limiting bovine pestivirus infections in Brazil must take into consideration this unique mix of viruses. This article provides a critical review of two decades of genetic identification of pestiviruses in Brazil.

Increased genetic variation of bovine viral diarrhea virus in dairy cattle in Poland

Background

Bovine viral diarrhea virus (BVDV) causes severe economic losses and is one of the most important viral pathogens of ruminants worldwide. The infection manifests itself in a variety of clinical symptoms. Phylogenetic studies based mainly on 5’UTR of its genome, identified many different subtypes of BVDV. Previous study indicated the predominance of BVDV-1b and BVDV-1d in Poland. The aim of this study was to genotype BVDV isolates currently circulating in Polish dairy herds.

Results

BVDV was detected in 30 herds. Viral subtypes were identified using sequences of the 5’UTR fragment and they were confirmed within a fragment of the N^pro region. Seven subtypes of BVDV-1 species have been identified: 1b, 1 g, 1f, 1d, 1r, 1 s and 1e.

Conclusion

The number of subtypes of BVDV in Poland evolves and 2 new subtypes have been identified for the first time. Such studies may have a positive impact on successful eradication of the virus using effective vaccines and diagnostic tests.

Electronic supplementary material

The online version of this article (10.1186/s12917-019-2029-z) contains supplementary material, which is available to authorized users.

Identification of BVDV2b and 2c subgenotypes in the United States: Genetic and antigenic characterization

Bovine viral diarrhea virus (BVDV), a ubiquitous pathogen of cattle, causes subclinical to severe acute disease. Two species of BVDV are recognized, BVDV1 and BVDV2 with BVDV1 divided into at least 21 subgenotypes and BVDV2 into 3-4 subgenotypes, most commonly using sequences from the 5' untranslated region (5' UTR). We report genomic sequencing of 8 BVDV2 isolates that did not segregate into the 2a subgenotype; but represented two additional BVDV2 subgenotypes. One BVDV2 subgenotype was previously recognized only in Asia. The other seven viruses fell into a second subgenotype that was first reported in Brazil and the U.S. in 2002. Neutralization assays using antiserum raised against vaccine strain BVDV2a 296c revealed varying degrees of neutralization of genetically diverse BVDV2 isolates. Neutralization titers decreased from 1.8 to more than a four log(2) decrease. This study illustrated the considerable genetic and antigenic diversity in BVDV2 circulating in the U.S.

Pestivirus infection in cattle dairy farms: E2 glycoprotein ELISA reveals the presence of bovine viral diarrhea virus type 2 in northwestern Italy

Background

Bovine viral diarrhea virus (BVDV) types 1 and 2 are members of the Pestivirus genus of the Flaviviridae family. This genus also includes the HoBi-like virus, tentatively classified as BVDV type 3. BVDV-1 is widely distributed in Italy despite the extensive use of BVDV-1-based vaccines, while BVDV-2 and HoBi-like Pestivirus have been detected occasionally. Monitoring the occurrence of sporadic or atypical pestiviruses is a useful approach to evaluate the need for additional vaccine strains that can be used in BVDV control programs.

Results

In this study we developed a multiwell antibody ELISA based on the recombinant E2 protein of the three bovine pestiviruses. We evaluated the assay’s applicability for surveillance purposes using pooled milk samples, each prepared from a maximum of 35 lactating cows and collected from 176 dairy herds. As expected, the majority of the pooled samples reacted to a greater extent against the BVDV-1 E2 antigen. All three milk pools from a single farm reacted to the BVDV-2 antigen, however. Further analysis using spot tests, antigen detection, and sequence analysis of the 5′-UTR region confirmed the presence of five persistently infected calves carrying a BVDV-2a strain.

Conclusions

This study highlights for the first time that sporadic circulation of BVDV-2 can be predicted by immunoenzymatic methods in the absence of specific vaccination.

Electronic supplementary material

The online version of this article (10.1186/s12917-017-1305-z) contains supplementary material, which is available to authorized users.

Genetic diversity of bovine viral diarrhea virus in cattle from Mexico

Bovine viral diarrhea virus (BVDV) infects cattle populations worldwide, causing significant economic losses though its impact on animal health. Previous studies have reported the prevalence of BVDV species and subgenotypes in cattle from the United States and Canada. We investigated the genetic diversity of BVDV strains detected in bovine serum samples from 6 different Mexican regions. Sixty-two BVDV isolates from Mexico were genetically typed based on comparison of sequences from the 5' untranslated region (5'-UTR) of the viral genome. Phylogenetic reconstruction indicated that 60 of the samples belonged to the BVDV-1 genotype and 2 to the BVDV-2 genotype. Comparison of partial 5'-UTR sequences clustered 49 samples within BVDV-1c, 8 samples within BVDV-1a, 3 samples within BVDV-1b, and 2 samples clustered with the BVDV-2a subgenotypes. Our study, combined with information previously published on BVDV field strain diversity in the United States and Canada, benefits the development of effective detection assays, vaccines, and control programs for North America.

Priming Cross-Protective Bovine Viral Diarrhea Virus-Specific Immunity Using Live-Vectored Mosaic Antigens

Bovine viral diarrhea virus (BVDV) plays a key role in bovine respiratory disease complex, which can lead to pneumonia, diarrhea and death of calves. Current vaccines are not very effective due, in part, to immunosuppressive traits and failure to induce broad protection. There are diverse BVDV strains and thus, current vaccines contain representative genotype 1 and 2 viruses (BVDV-1 & 2) to broaden coverage. BVDV modified live virus (MLV) vaccines are superior to killed virus vaccines, but they are susceptible to neutralization and complement-mediated destruction triggered by passively acquired antibodies, thus limiting their efficacy. We generated three novel mosaic polypeptide chimeras, designated N^proE2¹²³; NS²³¹; and NS²³², which incorporate protective determinants that are highly conserved among BVDV-1a, 1b, and BVDV-2 genotypes. In addition, strain-specific protective antigens from disparate BVDV strains were included to broaden coverage. We confirmed that adenovirus constructs expressing these antigens were strongly recognized by monoclonal antibodies, polyclonal sera, and IFN-γ-secreting T cells generated against diverse BVDV strains. In a proof-of-concept efficacy study, the multi-antigen proto-type vaccine induced higher, but not significantly different, IFN-γ spot forming cells and T-cell proliferation compared to a commercial MLV vaccine. In regards to the humoral response, the prototype vaccine induced higher BVDV-1 specific neutralizing antibody titers, whereas the MLV vaccine induced higher BVDV-2 specific neutralizing antibody titers. Following BVDV type 2a (1373) challenge, calves immunized with the proto-type or the MLV vaccine had lower clinical scores compared to naïve controls. These results support the hypothesis that a broadly protective subunit vaccine can be generated using mosaic polypeptides that incorporate rationally selected and validated protective determinants from diverse BVDV strains. Furthermore, regarding biosafety of using a live vector in cattle, we showed that recombinant human adenovirus-5 was cleared within one week following intradermal inoculation.

Genetic typing of bovine viral diarrhoea virus in cattle on Irish farms

The aim was to carry out a phylogenetic study of bovine viral diarrhoea viruses (BVDV) circulating in Irish cattle herds from 2011 to 2014. Three hundred and twenty five viruses from 267 herds were subtyped by nucleotide sequence analysis of the 5'UTR and/or the N^pro regions. All viruses investigated in this study belonged to species BVDV-1 with BVDV-1a as the prominent subtype (97%). Subtypes BVDV-1b, BVDV-1d and BVDV-1e were also identified for the first time in Ireland. Pairwise alignments of 225 viruses with complete sequences for the 5'UTR and the N^pro regions were performed to determine a low conflict threshold for virus strain demarcation. One hundred and seventy seven unique virus strains were identified. The study revealed significant levels of herd specific clustering of strains but no geographical or temporal clustering. Similar virus strains were identified in different counties, provinces and years indicating the potential to investigate the epidemiology of the disease by combining sequence analysis with animal movement data.

Molecular Characterization of a Novel Bovine Viral Diarrhea Virus Isolate SD-15

As one of the major pathogens, bovine viral diarrhea virus caused a significant economic loss to the livestock industry worldwide. Although BVDV infections have increasingly been reported in China in recent years, the molecular aspects of those BVDV strains were barely characterized. In this study, we reported the identification and characterization of a novel BVDV isolate designated as SD-15 from cattle, which is associated with an outbreak characterized by severe hemorrhagic and mucous diarrhea with high morbidity and mortality in Shandong, China. SD-15 was revealed to be a noncytopathic BVDV, and has a complete genomic sequence of 12,285 nucleotides that contains a large open reading frame encoding 3900 amino acids. Alignment analysis showed that SD-15 has 93.8% nucleotide sequence identity with BVDV ZM-95 isolate, a previous BVDV strain isolated from pigs manifesting clinical signs and lesions resembling to classical swine fever. Phylogenetic analysis clustered SD-15 to a BVDV-1m subgenotype. Analysis of the deduced amino acid sequence of glycoproteins revealed that E2 has several highly conserved and variable regions within BVDV-1 genotypes. An additional N-glycosylation site (240NTT) was revealed exclusively in SD-15-encoded E2 in addition to four potential glycosylation sites (Asn-X-Ser/Thr) shared by all BVDV-1 genotypes. Furthermore, unique amino acid and linear epitope mutations were revealed in SD-15-encoded E^rns glycoprotein compared with known BVDV-1 genotype. In conclusion, we have isolated a noncytopathic BVDV-1m strain that is associated with a disease characterized by high morbidity and mortality, revealed the complete genome sequence of the first BVDV-1m virus originated from cattle, and found a unique glycosylation site in E2 and a linear epitope mutation in E^rns encoded by SD-15 strain. Those results will broaden the current understanding of BVDV infection and lay a basis for future investigation on SD-15-related pathogenesis.

Genetic diversity of ruminant Pestivirus strains collected in Northern Ireland between 1999 and 2011 and the role of live ruminant imports

BACKGROUND

The genus pestivirus within the family Flaviviridae includes bovine viral diarrhoea virus (BVDV) types 1 and 2, border disease virus (BDV) and classical swine fever virus. The two recognised genotypes of BVDV are divided into subtypes based on phylogenetic analysis, namely a-p for BVDV-1 and a-c for BVDV-2.

METHODS

Three studies were conducted to investigate the phylogenetic diversity of pestiviruses present in Northern Ireland. Firstly, pestiviruses in 152 serum samples that had previously tested positive for BVDV between 1999 and 2008 were genotyped with a RT-PCR assay. Secondly, the genetic heterogeneity of pestiviruses from 91 serum samples collected between 2008 and 2011 was investigated by phylogenetic analysis of a 288 base pair portion of the 5' untranslated region (UTR). Finally, blood samples from 839 bovine and 4,437 ovine animals imported in 2010 and 2011 were tested for pestiviral RNA. Analysis of animal movement data alongside the phylogenetic analysis of the strains was carried out to identify any links between isolates and animal movement.

RESULTS

No BVDV-2 strains were detected. All of the 152 samples in the first study were genotyped as BVDV-1. Phylogenetic analysis indicated that the predominant subtype circulating was BVDV-1a (86 samples out of 91). The remaining five samples clustered close to reference strains in subtype BVDV-1b. Out of the imported animals, 18 bovine samples tested positive and 8 inconclusive (Ct ≥36), while all ovine samples were negative. Eight sequences were obtained and were defined as BVDV-1b. Analysis of movement data between herds failed to find links between herds where BVDV-1b was detected.

CONCLUSION

Given that only BVDV-1a was detected in samples collected between 1968 and 1999, this study suggests that at least one new subtype has been introduced to Northern Ireland between 1999 and 2011 and highlights the potential for importation of cattle to introduce new strains.

Complete Genome Sequence of Noncytopathic Bovine Viral Diarrhea Virus 1 Contaminating a High-Passage RK-13 Cell Line

A high-passage rabbit kidney RK-13 cell line (HP-RK-13[KY], originally derived from the ATCC CCL-37 cell line) used in certain laboratories worldwide is contaminated with noncytopathic bovine viral diarrhea virus (ncpBVDV). On complete genome sequence analysis, the virus strain was found to belong to BVDV group 1b.

Genomic characterization of a bovine viral diarrhea virus 1 isolate from swine

The SD0803 strain of the bovine viral diarrhea virus (BVDV) was isolated from a piglet in China in 2008 and has been classified as a novel subgenotype of BVDV-1. To describe the molecular features of this novel subgenotype, we sequenced and characterized the complete genome of the SD0803 virus. The genome is 12,271 bp in length and contains 5' and 3' untranslated regions (UTRs) that flank an open reading frame (ORF) encoding a 3,898-amino-acid polypeptide. The full-length genome of the SD0803 strain shares 78.8% to 83.3% identity with those of other BVDV-1 strains, 70.0% to 70.7% identity with those of BVDV-2 strains, and less than 67.6% identity with those of other pestiviruses. The highest level of shared identity was 83.3% between the complete SD0803 genome and that of the ZM-95 strain of BVDV-1. Phylogenetic analysis of the 5' UTR and the coding sequence for the N-terminal protease fragment of the SD0803 polyprotein indicated that the SD0803 virus is a member of the novel subgenotype BVDV-1q, isolates of which have been identified recently in dairy cattle and camels in China.

Numerical taxonomy of the genus Pestivirus: new software for genotyping based on the palindromic nucleotide substitutions method

The genus Pestivirus from the family Flaviviridae is represented by four established species; Bovine viral diarrhea virus 1 (BVDV-1); Bovine viral diarrhea virus 2 (BVDV-2); Border disease virus (BDV); and Classical swine fever virus (CSFV); as well a tentative species from a Giraffe. The palindromic nucleotide substitutions (PNS) in the 5' untranslated region (UTR) of Pestivirus RNA has been described as a new, simple and practical method for genotyping. New software is described, also named PNS, that was prepared specifically for this PNS genotyping procedure. Pestivirus identification using PNS was evaluated on five hundred and forty-three sequences at genus, species and genotype level using this software. The software is freely available at www.pns-software.com.