Resolving Bovine viral diarrhea virus subtypes from persistently infected U.S. beef calves with complete genome sequence

Association of ARRDC3 and NFIA variants with bovine congestive heart failure in feedlot cattle

BACKGROUND

Bovine congestive heart failure (BCHF) has become increasingly prevalent among feedlot cattle in the Western Great Plains of North America with up to 7% mortality in affected herds. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. Genes associated with BCHF in feedlot cattle have not been previously identified. Our aim was to search for genomic regions associated with this disease.

METHODS

A retrospective, matched case-control design with 102 clinical BCHF cases and their unaffected pen mates was used in a genome-wide association study. Paired nominal data from approximately 560,000 filtered single nucleotide polymorphisms (SNPs) were analyzed with McNemar's test.

RESULTS

Two independent genomic regions were identified as having the most significant association with BCHF: the arrestin domain-containing protein 3 gene ( ARRDC3), and the nuclear factor IA gene ( NFIA, mid- p-values, 1x10 -8 and 2x10 -7, respectively). Animals with two copies of risk alleles at either gene were approximately eight-fold more likely to have BCHF than their matched pen mates with either one or zero risk alleles at both genes (CI 95 = 3-17). Further, animals with two copies of risk alleles at both genes were 28-fold more likely to have BCHF than all others ( p-value = 1×10 -7, CI 95 = 4-206). A missense variant in ARRDC3 (C182Y) represents a potential functional variant since the C182 codon is conserved among all other jawed vertebrate species observed. A two-SNP test with markers in both genes showed 29% of 273 BCHF cases had homozygous risk genotypes in both genes, compared to 2.5% in 198 similar unaffected feedlot cattle. This and other DNA tests may be useful for identifying feedlot animals with the highest risk for BCHF in the environments described here.

CONCLUSIONS

Although pathogenic roles for variants in the ARRDC3 and NFIA genes are unknown, their discovery facilitates classifying animals by genetic risk and allows cattle producers to make informed decisions for selective breeding and animal health management.

Association of ARRDC3 and NFIA variants with bovine congestive heart failure in feedlot cattle

Background

Bovine congestive heart failure (BCHF) has become increasingly prevalent among feedlot cattle in the Western Great Plains of North America with up to 7% mortality in affected herds. BCHF is an untreatable complex condition involving pulmonary hypertension that culminates in right ventricular failure and death. Genes associated with BCHF in feedlot cattle have not been previously identified. Our aim was to search for genomic regions associated with this disease.

Methods

A retrospective, matched case-control design with 102 clinical BCHF cases and their unaffected pen mates was used in a genome-wide association study. Paired nominal data from approximately 560,000 filtered single nucleotide polymorphisms (SNPs) were analyzed with McNemar's test.

Results

Two independent genomic regions were identified as having the most significant association with BCHF: the arrestin domain-containing protein 3 gene ( ARRDC3), and the nuclear factor IA gene ( NFIA, mid- p-values, 1x10 -8 and 2x10 -7, respectively). Animals with two copies of risk alleles at either gene were approximately eight-fold more likely to have BCHF than their matched pen mates with either one or zero risk alleles at both genes (CI 95 = 3-17). Further, animals with two copies of risk alleles at both genes were 28-fold more likely to have BCHF than all others ( p-value = 1×10 -7, CI 95 = 4-206). A missense variant in ARRDC3 (C182Y) represents a potential functional variant since the C182 codon is conserved among all other jawed vertebrate species observed. A two-SNP test with markers in both genes showed 29% of 273 BCHF cases had homozygous risk genotypes in both genes, compared to 2.5% in 198 similar unaffected feedlot cattle. This and other DNA tests may be useful for identifying feedlot animals with the highest risk for BCHF in the environments described here.

Conclusions

Although pathogenic roles for variants in the ARRDC3 and NFIA genes are unknown, their discovery facilitates classifying animals by genetic risk and allows cattle producers to make informed decisions for selective breeding and animal health management.

Two bovine hepacivirus genome sequences from U.S. cattle

Bovine hepacivirus (BoHV) is closely related to the hepatitis C virus (HCV) in humans and can cause both acute and chronic liver infections in cattle. BoHV was first identified in Ghana and Germany in 2015 and since then it has been detected and characterized in other countries around the world, but no strains have been sequenced from U.S. cattle. To date, BoHV has been classified into 2 genotypes (1 and 2), with genotype 1 being further divided into 11 subtypes (A-K). However, the true genetic diversity of BoHV is likely underestimated given limited surveillance and a lack of published genome sequences. Here, we sequenced 2 nearly complete BoHV genomes from serum samples collected in 2019 from beef cattle in Missouri. Sequence comparisons and phylogenetic analysis showed that isolate MARC/2019/60 had high sequence homology with genotype 1, subtype E isolates from China. In contrast, isolate MARC/2019/50 represented a novel BoHV subtype within genotype 2. Thus, we report the first genomic characterization of BoHV isolates from U.S. cattle, and the second complete BoHV2 genome worldwide. This work increases our knowledge of the global genetic diversity of BoHV and demonstrates the co-circulation of divergent BoHV strains in U.S. cattle.

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.

Multivariate Analysis as a Method to Evaluate Antigenic Relationships between Bovine Viral Diarrhea Virus 1b Isolates and Vaccine Strains

The antigenicity of bovine viral diarrhea virus (BVDV) has been evaluated using virus-neutralizing titer data analyzed by principal component analysis (PCA) and has demonstrated numerous isolates to be antigenically divergent from US vaccine strains. The lack of BVDV-1b strains in currently licensed vaccines has raised concerns regarding the lack of protection against BVDV-1b field strains. The aim of this study was to evaluate the antigenic diversity of BVDV-1b strains and better understand the breadth of antigenic relatedness using BVDV-1b antisera and antisera from vaccine strains. Results from this analysis demonstrate the antigenic diversity observed among BVDV-1b isolates and genetic assignment into the BVDV-1b subgenotype is not representative of antigenic relatedness. This is demonstrated by BVDV-1b isolates (2280N, SNc, Illc, MSU, and 2337) observed to be as antigenically dissimilar as BVDV-2a isolates when using BVDV-1b antisera. Additionally, when BVDV-1a vaccine antisera was used for comparisons, a greater percentage of BVDV-1b isolates clustered with BVDV-1a vaccine strains as part of PC1, suggesting antigenic relatedness and potentially partial protection. Collectively, data from this study would suggest that while most BVDV-1b isolates are antigenically similar, there are antigenically dissimilar BVDV-1b isolates as determined by the lack of cross-reactivity, which may contribute to the lack of protection.

First report of bovine viral diarrhea virus subgenotypes 1d and 1e in southern Chile

Bovine viral diarrhea virus (BVDV) affects cattle worldwide causing severe productive and economic loss. In this study, we investigated the subgenotypes of BVDV circulating in cattle samples from the Aysén region, an active cattle breeding area located in southern Chile. Partial amplification of the 5' untranslated region (UTR) was performed by polymerase chain reaction (PCR), and twelve samples were analyzed by Sanger sequencing and phylogenetic analysis. Eight samples were identified as belonging to Pestivirus bovis subgenotype 1e, three to 1-b, and one to 1-d. The phylogenetic analyses performed revealed a marked distance between these now-identified strains and those previously reported in the country. These findings support the need to continually expand the analysis of the variability of the viral phylogeny for the currently circulating BVDV strains and to update the vaccines recommended for this livestock area and surrounding areas.

First gene-edited calf with reduced susceptibility to a major viral pathogen

Bovine viral diarrhea virus (BVDV) is one of the most important viruses affecting the health and well-being of bovine species throughout the world. Here, we used CRISPR-mediated homology-directed repair and somatic cell nuclear transfer to produce a live calf with a six amino acid substitution in the BVDV binding domain of bovine CD46. The result was a gene-edited calf with dramatically reduced susceptibility to infection as measured by reduced clinical signs and the lack of viral infection in white blood cells. The edited calf has no off-target edits and appears normal and healthy at 20 months of age without obvious adverse effects from the on-target edit. This precision bred, proof-of-concept animal provides the first evidence that intentional genome alterations in the CD46 gene may reduce the burden of BVDV-associated diseases in cattle and is consistent with our stepwise, in vitro and ex vivo experiments with cell lines and matched fetal clones.

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.

Recent Emergence of Bovine Coronavirus Variants with Mutations in the Hemagglutinin-Esterase Receptor Binding Domain in U.S. Cattle

Bovine coronavirus (BCoV) has spilled over to many species, including humans, where the host range variant coronavirus OC43 is endemic. The balance of the opposing activities of the surface spike (S) and hemagglutinin-esterase (HE) glycoproteins controls BCoV avidity, which is critical for interspecies transmission and host adaptation. Here, 78 genomes were sequenced directly from clinical samples collected between 2013 and 2022 from cattle in 12 states, primarily in the Midwestern U.S. Relatively little genetic diversity was observed, with genomes having >98% nucleotide identity. Eleven isolates collected between 2020 and 2022 from four states (Nebraska, Colorado, California, and Wisconsin) contained a 12 nucleotide insertion in the receptor-binding domain (RBD) of the HE gene similar to one recently reported in China, and a single genome from Nebraska collected in 2020 contained a novel 12 nucleotide deletion in the HE gene RBD. Isogenic HE proteins containing either the insertion or deletion in the HE RBD maintained esterase activity and could bind bovine submaxillary mucin, a substrate enriched in the receptor 9-O-acetylated-sialic acid, despite modeling that predicted structural changes in the HE R3 loop critical for receptor binding. The emergence of BCoV with structural variants in the RBD raises the possibility of further interspecies transmission.

A Serosurvey for Ruminant Pestivirus Exposure Conducted Using Sera From Stray Mexico Origin Cattle Captured Crossing Into Southern Texas

The US Department of Agriculture (USDA), Animal Plant Health Inspection Service (APHIS), Cattle Fever Tick Eradication Program (CFTEP) monitor a quarantine zone along the Texas border to prevent the introduction of stray livestock carrying cattle fever ticks entering the United States from Mexico. Stray cattle collected by CFTEP are checked for ticks and several infectious disease-causing pathogens, but not for bovine viral diarrhea virus (BVDV). BVDV is one of the most economically impactful viruses affecting US cattle producers. BVDV is present in all parts of the world, but it has been demonstrated that another distantly related pestivirus, HoBi-like pestivirus (HoBiPev), can also cause BVD. To date, HoBiPev has not been detected in the United States, but is commonly found in Brazil, and sporadically in Europe and Asia. The objective of the current study was to evaluate the seroprevalence of pestiviruses, with a specific focus on HoBiPev, in stray cattle. Virus neutralization (VN) assay was used to determine seroprevalence (or antibody titers) of BVDV-1, BVDV-2, and HoBiPev. Approximately 50% (67 of 134) of the samples were seropositive for pestiviruses; all 67 positive samples were positive (50%) for BVDV-1, 66 samples of the 67 were positive (49.3%) for BVDV-2, and the same 66 samples of the 67 were also positive (49.3%) for HoBiPev. Due to the antigenic cross-reactivity among Pestiviruses, the comparative antibody against each pestivirus was calculated from all VN-positive samples. Titers were clearly higher against BVDV-1, and only one sample had a titer clearly higher against BVDV-2. No sample had an antibody titer higher for HoBiPev, and while this does not prove the absence of HoBiPev, it does provide evidence that the prevalence of HoBiPev is less predominant than BVDV-1. Additionally, data from these samples provide evidence on the susceptibility of animals that may enter into the United States, with ~50% of the animals seronegative for bovine pestiviruses. This cattle population provides a unique opportunity to evaluate and monitor changes in seroprevalence of economically important cattle diseases affecting the cattle industry.

Use of multivariate analysis to evaluate antigenic relationships between US BVDV vaccine strains and non-US genetically divergent isolates

Bovine viral diarrhea virus (BVDV) comprises two species, BVDV-1 and BVDV-2. But given the genetic diversity among pestiviruses, at least 22 subgenotypes are described for BVDV-1 and 3-4 for BVDV-2. Genetic characterization is generally accomplished through complete or partial sequencing and phylogeny, but it is not a reliable method to define antigenic relationships. The traditional method for evaluating antigenic relationships between pestivirus isolates is the virus neutralization (VN) assay, but interpretation of the data to define antigenic relatedness can be difficult to discern for BVDV isolates within the same BVDV species. Data from this study utilized a multivariate analysis for visualization of VN results to analyze the antigenic relationships between US vaccine strains and field isolates from Switzerland, Italy, Brazil, and the UK. Polyclonal sera were generated against six BVDV strains currently contained in vaccine formulations, and each serum was used in VNs to measure the titers against seven vaccine strains (including the six homologous strains) and 23 BVDV field isolates. Principal component analysis (PCA) was performed using VN titers, and results were interpreted from PCA clustering within the PCA dendrogram and scatter plot. The results demonstrated clustering patterns among various isolates suggesting antigenic relatedness. As expected, the BVDV-1 and BVDV-2 isolates did not cluster together and had the greatest spatial distribution. Notably, a number of clusters representing antigenically related BVDV-1 subgroups contain isolates of different subgenotypes. The multivariate analysis may be a method to better characterize antigenic relationships among BVDV isolates that belong to the same BVDV species and do not have distinct antigenic differences. This might be an invaluable tool to ameliorate the composition of current vaccines, which might well be important for the success of any BVDV control program that includes vaccination in its scheme.

Evaluating Large Spontaneous Deletions in a Bovine Cell Line Selected for Bovine Viral Diarrhea Virus Resistance

Bovine viral diarrhea virus’s (BVDV) entry into bovine cells involves attachment of virions to cellular receptors, internalization, and pH-dependent fusion with endosomal membranes. The primary host receptor for BVDV is CD46; however, the complete set of host factors required for virus entry is unknown. The Madin-Darby bovine kidney (MDBK) cell line is susceptible to BVDV infection, while a derivative cell line (CRIB) is resistant at the level of virus entry. We performed complete genome sequencing of each to identify genomic variation underlying the resistant phenotype with the aim of identifying host factors essential for BVDV entry. Three large compound deletions in the BVDV-resistant CRIB cell line were identified and predicted to disrupt the function or expression of the genes PTPN12, GRID2, and RABGAP1L. However, CRISPR/Cas9 mediated knockout of these genes, individually or in combination, in the parental MDBK cell line did not impact virus entry or replication. Therefore, resistance to BVDV in the CRIB cell line is not due to the apparent spontaneous loss of PTPN12, GRID2, or RABGAP1L gene function. Identifying the functional cause of BVDV resistance in the CRIB cell line may require more detailed comparisons of the genomes and epigenomes.

Cattle connection: molecular epidemiology of BVDV outbreaks via rapid nanopore whole-genome sequencing of clinical samples

BACKGROUND

As a global ruminant pathogen, bovine viral diarrhea virus (BVDV) is responsible for the disease Bovine Viral Diarrhea with a variety of clinical presentations and severe economic losses worldwide. Classified within the Pestivirus genus, the species Pestivirus A and B (syn. BVDV-1, BVDV-2) are genetically differentiated into 21 BVDV-1 and four BVDV-2 subtypes. Commonly, the 5' untranslated region and the N^pro protein are utilized for subtyping. However, the genetic variability of BVDV leads to limitations in former studies analyzing genome fragments in comparison to a full-genome evaluation.

RESULTS

To enable rapid and accessible whole-genome sequencing of both BVDV-1 and BVDV-2 strains, nanopore sequencing of twelve representative BVDV samples was performed on amplicons derived through a tiling PCR procedure. Covering a multitude of subtypes (1b, 1d, 1f, 2a, 2c), sample matrices (plasma, EDTA blood and ear notch), viral loads (Cq-values 19-32) and species (cattle and sheep), ten of the twelve samples produced whole genomes, with two low titre samples presenting 96 % genome coverage.

CONCLUSIONS

Further phylogenetic analysis of the novel sequences emphasizes the necessity of whole-genome sequencing to identify novel strains and supplement lacking sequence information in public repositories. The proposed amplicon-based sequencing protocol allows rapid, inexpensive and accessible obtainment of complete BVDV genomes.

Evaluation of Antigenic Comparisons Among BVDV Isolates as it Relates to Humoral and Cell Mediated Responses

Antigenic differences between bovine viral diarrhea virus (BVDV) vaccine strains and field isolates can lead to reduced vaccine efficacy. Historically, antigenic differences among BVDV strains were evaluated using techniques based on polyclonal and monoclonal antibody activity. The most common method for antigenic comparison among BVDV isolates is determination of virus neutralization titer (VNT). BVDV antigenic comparisons using VNT only account for the humoral component of the adaptive immune response, and not cell mediated immunity (CMI) giving an incomplete picture of protective responses. Currently, little data is available regarding potential antigenic differences between BVDV vaccine strains and field isolates as measured by CMI responses. The goal of the current paper is to evaluate two groups of cattle that differed in the frequency they were vaccinated, to determine if similar trends in CMI responses exist within each respective group when stimulated with antigenically different BVDV strains. Data from the current study demonstrated variability in the CMI response is associated with the viral strain used for stimulation. Variability in IFN-γ mRNA expression was most pronounced in the CD4⁺ population, this was observed between the viruses within each respective BVDV subgenotype in the Group 1 calves. The increase in frequency of CD25⁺ cells and IFN-γ mRNA expression in the CD8⁺ and CD335⁺ populations were not as variable between BVDV strains used for stimulation in the Group 1 calves. Additionally, an inverse relationship between VNT and IFN-γ mRNA expression was observed, as the lowest VNT and highest IFN-γ mRNA expression was observed and vice versa, the highest VNT and lowest IFN-γ mRNA expression was observed. A similar trend regardless of vaccination status was observed between the two groups of calves, as the BVDV-1b strain had lower IFN-γ mRNA expression. Collectively, data from the current study and previous data support, conferring protection against BVDV as a method for control of BVDV in cattle populations is still a complex issue and requires a multifactorial approach to understand factors associated with vaccine efficacy or conversely vaccine failure. Although, there does appear to be an antigenic component associated with CMI responses as well as with humoral responses as determined by VNT.

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.

Persistent infection of American bison (Bison bison) with bovine viral diarrhea virus and bosavirus

Bovine viral diarrhea viruses (BVDV) are significant pathogens of cattle, leading to losses associated with reproductive failure, respiratory disease and immune dysregulation. While cattle are the reservoir for BVDV, a wide range of domestic and wild ruminants are susceptible to infection and disease caused by BVDV. Samples from four American bison (Bison bison) from a captive herd were submitted for diagnostic testing due to their general unthriftiness. Metagenomic sequencing on pooled nasal swabs and serum identified co-infection with a BVDV and a bovine bosavirus. The BVDV genome was more similar to the vaccine strain Oregon C24 V than to other BVDV sequences in GenBank, with 92.7 % nucleotide identity in the open reading frame. The conserved 5'-untranslated region was 96.3 % identical to Oregon C24 V. Bosavirus has been previously identified in pooled fetal bovine serum but its clinical significance is unknown. Sequencing results were confirmed by virus isolation and PCR detection of both viruses in serum and nasal swab samples from two of the four bison. One animal was co-infected with both BVDV and bosavirus while separate individuals were positive solely for BVDV or bosavirus. Serum and nasal swabs from these same animals collected 51 days later remained positive for BVDV and bosavirus. These results suggest that both viruses can persistently infect bison. While the etiological significance of bosavirus infection is unknown, the ability of BVDV to persistently infect bison has implications for BVDV control and eradication programs. Possible synergy between BVDV and bosavirus persistent infection warrants further study.

Multivariate analysis as a method to evaluate antigenic relationships between BVDV vaccine and field strains

Bovine viral diarrhea virus (BVDV) is comprised of two species, BVDV-1 and BVDV-2, but given the genetic diversity among pestiviruses, at least 21 subgenotypes are described for BVDV-1 and 4 for BVDV-2. Genetic characterization can be achieved through complete or partial sequencing and phylogeny, but antigenic characterization can be difficult to determine due to the antigenic diversity and cross-neutralization that exists among isolates. The traditional method for evaluating antigenic relationships between pestivirus isolates is the virus neutralization (VN) assay, but interpretation of the data to determine antigenic difference can be unclear. Data from this study utilized a multivariate analysis for visualization of VN results to analyze the antigenic relationships between vaccine strains and multiple field isolates. Polyclonal sera were generated against 6 BVDV strains currently contained in vaccine formulations, and each serum was used in VN's to measure the neutralizing antibody titers against 15 BVDV field isolates characterized as prevalent and divergent subgenotypes in the USA. Principal component analysis (PCA) were performed on the VN assay datasets, and results were interpreted from PCA clustering within the PCA dendrogram and scatter plot. The results demonstrated clustering patterns among isolates suggestive of antigenic differences. While expected, the BVDV-1 and BVDV-2 isolates did not cluster together and had the greatest spatial distribution. In addition, other BVDV isolates had distinct spatial patterns suggesting antigenically divergent isolates. This analysis provides an alternative and more efficient means to analyze large VN datasets to visualize antigenic relationships between pestivirus isolates. This analysis could be beneficial for vaccine development and evaluation of efficacy, since most vaccines cannot fully protect animals from the broad range diversity of BVDV viruses.

Bovine viral diarrhea virus: An updated American College of Veterinary Internal Medicine consensus statement with focus on virus biology, hosts, immunosuppression, and vaccination

Control of bovine viral diarrhea virus (BVDV) in cattle populations across most of the world has remained elusive in spite of advances in knowledge about this viral pathogen. A central feature of virus perseverance in cattle herds is the unique mechanism of persistent infection. Managing BVDV infection in herds involves controlling persistently infected carrier animals using a multidimensional approach of vaccination, biosecurity, and identification of BVDV reservoirs. A decade has passed since the original American College of Veterinary Internal Medicine consensus statement on BVDV. While much has remained the same with respect to clinical signs of disease, pathogenesis of infection including persistent infection, and diagnosis, scientific articles published since 2010 have led to a greater understanding of difficulties associated with control of BVDV. This consensus statement update on BVDV presents greater focus on topics currently relevant to the biology and control of this viral pathogen of cattle, including changes in virus subpopulations, infection in heterologous hosts, immunosuppression, and vaccination.

Viruses in Bovine Respiratory Disease in North America: Knowledge Advances Using Genomic Testing

Advances in viral detection in bovine respiratory disease (BRD) have resulted from advances in viral sequencing of respiratory tract samples. New viruses detected include influenza D virus, bovine coronavirus, bovine rhinitis A, bovine rhinitis B virus, and others. Serosurveys demonstrate widespread presence of some of these viruses in North American cattle. These viruses sometimes cause disease after animal challenge, and some have been found in BRD cases more frequently than in healthy cattle. Continued work is needed to develop reagents for identification of new viruses, to confirm their pathogenicity, and to determine whether vaccines have a place in their control.

Bovine viral diarrhoea virus loses quasispecies diversity rapidly in culture

Bovine viral diarrhoea (BVD) is an important disease of cattle, with significant impacts on animal health and welfare. The wide host range of the causative pestiviruses may lead to formation of virus reservoirs in other ruminant or wildlife species, presenting a concern for the long-term success of BVD eradication campaigns. It is likely that the quasispecies nature of these RNA viruses contributes to their interspecies transmission by providing genetic plasticity. Understanding the spectrum of sequence variants present in persistently infected (PI) animals is, therefore, essential for studies of virus transmission. To analyse quasispecies diversity without amplification bias, we extracted viral RNA from the serum of a PI cow, and from cell culture fluid after three passages of the same virus in culture, to produce cDNA without amplification. Sequencing of this material using Illumina 250 bp paired-read technology produced full-length virus consensus sequences from both sources and demonstrated the quasispecies diversity of this pestivirus A genotype 1a field strain within serum and after culture. We report the distribution and diversity of over 800 SNPs and provide evidence for a loss of diversity after only three passages in cell culture, implying that cultured viruses cannot be used to understand quasispecies diversity and may not provide reliable molecular markers for source tracing or transmission studies. Additionally, both serum and cultured viruses could be sequenced as a set of 25 overlapping PCR amplicons that demonstrated the same consensus sequences and the presence of many of the same quasispecies variants. The observation that aspects of the quasispecies structure revealed by massively parallel sequencing are also detected after PCR and Sanger sequencing suggests that this approach may be useful for small or difficult to analyse samples.

Upregulation of the type I interferon pathway in feedlot cattle persistently infected with bovine viral diarrhea virus

Bovine viral diarrhea virus (BVDV) has a profound economic impact on the cattle industry. Calves infected in utero and born persistently infected (PI) with BVDV have increased morbidity, mortality, and reduced productivity. Further, they serve as a continual source of viral exposure to herd mates and thereby pose a significant risk to animal wellbeing and production efficiency. Understanding the mechanisms through which PI is established and maintained is therefore important in working toward finding means to prevent or mitigate losses due to infection. Early studies of acute infection suggested BVDV infection alters the host's ability to mount a type I interferon (IFN) response, thereby allowing for the establishment of PI. More recently, however, animals experimentally challenged with the virus demonstrated a chronic yet modest upregulation of the IFN pathway. To identify if the IFN or other pathways are altered due to PI by BVDV in a natural infection, the circulating blood transcriptome was analyzed from PI feedlot cattle (N = 10 BVDV1a, 8 BVDV1b, 8 BVDV2), cattle co-mingling with PI cattle but not themselves infected (N = 9), and a group of unrelated, unexposed controls (N=10). Differential expression analyses included contrasts among BVDV subtypes, and all pair-wise comparisons of PI, co-mingled non-PI, and unexposed animals. Analyses in limma-voom revealed no difference in the transcriptome based upon the BVDV genotype with which the animal was infected. However, gene expression did differ (adj P < 0.05 and |logFC|> 1) at 175 loci between the PI and co-housed, non-PI contemporaries and when compared to the unexposed controls, 489 loci were differentially expressed. Pathway analyses predict that alterations in the transcriptome of the PI cattle indicate significant upregulation of innate immune function including IFN signaling. These data support prior work suggesting IFN signaling is not completely suppressed in cattle naturally PI with BVDV.

In vitro method to evaluate virus competition between BVDV-1 and BVDV-2 strains using the PrimeFlow RNA assay

Bovine viral diarrhea viruses (BVDV), segregated in BVDV-1 and BVDV-2 species, lead to substantial economic losses to the cattle industry worldwide. It has been hypothesized that there could be differences in level of replication, pathogenesis and tissue tropism between BVDV-1 and BVDV-2 strains. Thus, this study developed an in vitro method to evaluate virus competition between BVDV-1 and BVDV-2 strains. To this end the competitive dynamics of BVDV-1a, BVDV-1b, and BVDV-2a strains in cell cultures was evaluated by a PrimeFlow RNA assay. Similar results were observed in this study, as was observed in an earlier in vivo transmission study. Competitive exclusion was observed as the BVDV-2a strains dominated and excluded the BVDV-1a and BVDV-1b strains. The in vitro model developed can be used to identify viral variations that result in differences in frequency of subgenotypes detected in the field, vaccine failure, pathogenesis, and strain dependent variation in immune responses.

High prevalence of persistently infected animals from bovine viral diarrhea in Colombian cattle

Background

Bovine Viral Diarrhea Virus (BVDV) is associated with gastrointestinal, respiratory and reproductive diseases of livestock across the world that causes continuous economic losses in the cattle industry. This virus can establish a persistent infection (PI) in calves after the fetal infection, making BVDV positive catle carriers and primary reservoirs which will constantly transmit the virus to healthy and new-born animals. For this reason, the detection of the PI animals in herds is the first line of prevention of the viral infection.

Results

In this study, PI animals were detected in five different regions of Colombia through RT-PCR techniques and confirmed by sequencing. BVDV genotypes were determined using one fragment of the 5’UTR. It was found a 7% BVDV prevalence in animals and 22% in farms; and genotype 1 was identified as a single genotype for all of the samples. All samples were BVDV 1a.

Conclusion

This is the first report in Colombia with higher prevalence rates compared with other places in the world, turned out to be of great importance for the ranchers, the vaccine producers and animal health control parties.

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.

Selection and characterization of specific nanobody against bovine virus diarrhea virus (BVDV) E2 protein

Bovine viral diarrhea-mucosal disease (BVD-MD) is caused by bovine viral diarrhea virus (BVDV), and results in abortion, stillbirth, and fetal malformation in cows. Here, we constructed the phage display vector pCANTAB 5E-VHH and then transformed it into Escherichia coli TG1-competent cells, to construct an initial anti-BVDV nanobody gene library. We obtained a BVDV-E2 antigen epitope bait protein by prokaryotic expression using the nucleotide sequence of the E2 gene of the BVDV-NADL strain published in GenBank. Phage display was used to screen the anti-BVDV nanobody gene library. We successfully constructed a high quality phage display nanobody library, with an initial library capacity of 4.32×105. After the rescue of helper phage, the titer of the phage display nanobody library was 1.3×1011. The BVDV-E2 protein was then expressed in Escherichia coli (DE3), and a 49.5 kDa band was observed with SDS-PAGE analysis that was consistent with the expected nanobody size. Thus, we were able to isolate one nanobody that exhibits high affinity and specificity against BVDV using phage display techniques. This isolated nanobody was then used in Enzyme Linked Immunosorbent Assay and qRT-PCR, and ELISA analyses of BVDV infection of MDBK cells indicated that the nanobodies exhibited good antiviral effect.