Metadata beyond the sequence enables the phylodynamic inference of bovine viral diarrhea virus type 1a isolates from Western Canada

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.

Temporal and geographic dynamics of bovine viral diarrhea virus in American countries

Bovine viral diarrhea virus (BVDV) is a worldwide distributed pathogen of livestock classified into three species, BVDV-1 (Pestivirus A), BVDV-2 (Pestivirus B), and HoBi-like pestivirus (HoBiPeV; Pestivirus H). Despite being considered endemic in several regions of the Americas, the spatiotemporal distribution of BVDV is scarcely known. This study aimed to reconstruct the population dynamics of BVDV in American countries. The analyses performed with the partial 5´UTR gene showed that BVDV-1 and -2 would have started their diversification in the 1670s and 1790s in the United States, whereas HoBiPeV probably emerged in the 1980s in Brazil. No evident geographic clustering was observed in the Bayesian trees, which may indicate that multiple introductions events would have occurred following the first introduction. This study provides new insights into BVDV dynamics, although further analyses including sequences from other American countries and continents will help to expand the knowledge of BVDV evolution and transmission.

Isolation of BVDV-1a, 1m, and 1v strains from diarrheal calf in china and identification of its genome sequence and cattle virulence

Bovine viral diarrhea virus (BVDV) is an important livestock viral pathogen responsible for causing significant economic losses. The emerging and novel BVDV isolates are clinically and biologically important, as there are highly antigenic diverse and pathogenic differences among BVDV genotypes. However, no study has yet compared the virulence of predominant genotype isolates (BVDV-1a, 1b, and 1m) in China and the emerging genotype isolate BVDV-1v. The serological relationship among these genotypes has not yet been described. In this study, we isolated three BVDV isolates from calves with severe diarrhea, characterized as BVDV-1a, 1m, and novel 1v, based on multiple genomic regions [including 5-untranslated region (5'-UTR), Npro, and E2] and the phylogenetic analysis of nearly complete genomes. For the novel genotype, genetic variation analysis of the E2 protein of the BVDV-1v HB-03 strain indicates multiple amino acid mutation sites, including potential host cell-binding sites and neutralizing epitopes. Recombination analysis of the BVDV-1v HB-03 strain hinted at the possible occurrence of cross-genotypes (among 1m, 1o, and 1q) and cross-geographical region transmission events. To compare the pathogenic characters and virulence among these BVDV-1 genotypes, newborn calves uninfected with common pathogens were infected intranasally with BVDV isolates. The calves infected with the three genotype isolates show different symptom severities (diarrhea, fever, slowing weight gain, virus shedding, leukopenia, viremia, and immune-related tissue damage). In addition, these infected calves also showed bovine respiratory disease complexes (BRDCs), such as nasal discharge, coughing, abnormal breathing, and lung damage. Based on assessing different parameters, BVDV-1m HB-01 is identified as a highly virulent strain, and BVDV-1a HN-03 and BVDV-1v HB-03 are both identified as moderately virulent strains. Furthermore, the cross-neutralization test demonstrated the antigenic diversity among these Chinese genotypes (1a, 1m, and 1v). Our findings illustrated the genetic evolution characteristics of the emerging genotype and the pathogenic mechanism and antigenic diversity of different genotype strains, These findings also provided an excellent vaccine candidate strain and a suitable BVDV challenge strain for the comprehensive prevention and control of BVDV.

Bovine viral diarrhea virus in China: A comparative genomic and phylogenetic analysis with complete genome sequences

Bovine viral diarrhea virus (BVDV), causing bovine viral diarrhea (BVD) in cattle, is one of the highly contagious and devastating diseases of cattle. Since 1980, BVDV has been identified all-over China in a variety of animal species including cattle, camels, yaks, sheep, water buffalo, goats, Sika deer and pigs. In this study, 31 BVDV complete genomes reported in China (from 2004 to 2020) with other 112 genomes reported around the world were comparatively analyzed. Phylogenetic analysis shows that BVDV genomes reported worldwide clustered in three major clades i.e., BVDV-1, BVDV-2, and BVDV-3. The BVDV-1 is genetically the most diverged genotype and phylogenetically classified into 7 sub-clades in our study based on full-length genomes. The China BVDV genomes fall into all three major clades, e.g., BVDV-1, BVDV-2 and BVDV-3. China BVDV-1 clustered into five sub-clades, e.g., 1, 2, 3, 6 and 7, where sub-clade 7 clustered as a separate sub-clade. Full-length genome recombination analysis reveals that the BVDV-1 reported in China appears to be mainly involved in recombination events. In addition, comparative analysis of E2 proteins between BVDV-1, BVDV-2, and BVDV-3 reveals that the amino acid variations could affect 12 potential linear B cell epitopes, demonstrating a dramatic antigen drift in the E2 protein. These results provide a thorough view of the information about the genetic and antigenic diversity of BVDVs circulating in China and therefore could benefit the development of suitable strategies for disease control.

Temporal analysis of bovine pestivirus diversity in Brazil

In this study, phylogenetic and evolutionary analyses of cattle pestiviruses (BVDV-1, 2 and HoBiPeV) originating in Brazil were used to investigate the temporal diversification of subgenotypes in the country. Inferred dated phylogeny and time of the most recent common ancestor (tMRCA) demonstrated that some BVDV subgenotypes (1a, 1b, 1d, 1e, and 2b) and HoBi-like sequences clustered according to the region in which they were collected and that the diversification of subgenotypes appears to have occurred around the introduction of first Bos taurus and then Bos indicus, followed by expansion to form the adapted Brazilian breeds. The present results help to elucidate the temporal facts that led to diversification of ruminant pestiviruses in cattle in Brazil.

Genomic diversity and phylodynamic of bovine viral diarrhea virus in Argentina

Bovine viral diarrhea virus (BVDV) is an important pathogen of ruminants worldwide and is characterized by high genetic diversity and a wide range of clinical presentations. In Argentina, several studies have evaluated the genetic diversity of BVDV but no phylodynamic study has been published yet. In this study, a comprehensive compilation and update of Argentinean BVDV sequences were performed, and the evolutionary history of BVDV was characterized by phylodynamic analyses based on the 5´UTR. Although BVDV-1b and BVDV-1a were the most frequent subtypes, novel subtypes for Argentina, 1e and 1i, were identified. The phylodynamic analysis suggested that BVDV started its diversification in the mid-1650s with an exponential increase in viral diversity since the late 1990s, possibly related to the livestock expansion and intensification in the country. Evolutionary rate in the 5´UTR was faster for BVDV-1a than for BVDV-1b, and both subtypes presented an endemic nature according to the demographic reconstructions. The current study contributes to clarify the evolutionary history of BVDV in the main cattle region of the country and provides useful information about the epidemiology and future development of diagnostic and control tools in Argentina.

Viral Enteritis in Cattle: To Well Known Viruses and Beyond

Livestock products supply about 13 percent of energy and 28 percent of protein in diets consumed worldwide. Diarrhea is a leading cause of sickness and death of beef and dairy calves in their first month of life and also affecting adult cattle, resulting in large economic losses and a negative impact on animal welfare. Despite the usual multifactorial origin, viruses are generally involved, being among the most important causes of diarrhea. There are several viruses that have been confirmed as etiological agents (i.e., rotavirus and coronavirus), and some viruses that are not yet confirmed as etiological agents. This review summarizes the viruses that have been detected in the enteric tract of cattle and tries to deepen and gather knowledge about them.

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.

Insights into the origin and diversification of bovine viral diarrhea virus 1 subtypes

In this study, we performed phylogenetic and evolutionary analysis on bovine viral diarrhea virus 1 (BVDV-1) sequences to investigate the origin and temporal diversification of different BVDV-1 subtypes. Dated phylogenies using the complete polyprotein sequence were reconstructed, and the time of the most recent common ancestor (tMRCA) was estimated. The results demonstrated that BVDV-1 subtypes clustered into two phylogenetic clades, where the predominant subtypes worldwide grouped together. In the temporal analysis, the tMRCA of BVDV-1 was 1336, and the diversification into different subtypes appears to have occurred around 363 years ago. The present results help to elucidate the origins of BVDV-1 subtypes and the dynamics of ruminant pestiviruses.

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.

An extensive field study reveals the circulation of new genetic variants of subtype 1a of bovine viral diarrhea virus in Uruguay

Bovine viral diarrhea virus (BVDV) is a major pathogen worldwide, causing significant economic losses to the livestock sector. In Uruguay, BVDV seroprevalence at the farm level is >80%. In this work, 2546 serum, blood or tissue samples collected from animals suspected of being affected by BVD between 2015 and 2017 were analyzed by reverse transcription PCR and sequencing. Analysis of the BVDV genomic regions 5'UTR/N^pro, N^pro and E2 revealed that BVDV-1a, 1i and 2b circulate in the country, with BVDV-1a being the most prevalent subtype. Population dynamics studies revealed that BVDV-1a has been circulating in our herds since ~1990. This subtype began to spread and evolve, accumulating point mutations at a rate of 3.48 × 10^-3 substitutions/site/year, acquiring specific genetic characteristics that gave rise to two local genetic lineages of BVDV-1a. These lineages are divergent from those circulating worldwide, as well as the vaccine strain currently used in Uruguay. The most notable differences between field and vaccine strains were found in the E2 glycoprotein, suggesting that the amino acid substitutions could result in failure of cross-protection/neutralization after vaccination. This is the first study that compares Uruguayan BVDV field and vaccine strains with other BVDV strains from throughout the world. The results obtained in this study will be very useful for developing a suitable immunization program for BVDV in Uruguay by identifying local field strains as candidates for vaccine development.

Assessing biases in phylodynamic inferences in the presence of super-spreaders

Phylodynamic analyses using pathogen genetic data have become popular for making epidemiological inferences. However, many methods assume that the underlying host population follows homogenous mixing patterns. Nevertheless, in real disease outbreaks, a small number of individuals infect a disproportionately large number of others (super-spreaders). Our objective was to quantify the degree of bias in estimating the epidemic starting date in the presence of super-spreaders using different sample selection strategies. We simulated 100 epidemics of a hypothetical pathogen (fast evolving foot and mouth disease virus-like) over a real livestock movement network allowing the genetic mutations in pathogen sequence. Genetic sequences were sampled serially over the epidemic, which were then used to estimate the epidemic starting date using Extended Bayesian Coalescent Skyline plot (EBSP) and Birth–death skyline plot (BDSKY) models. Our results showed that the degree of bias varies over different epidemic situations, with substantial overestimations on the epidemic duration occurring in some occasions. While the accuracy and precision of BDSKY were deteriorated when a super-spreader generated a larger proportion of secondary cases, those of EBSP were deteriorated when epidemics were shorter. The accuracies of the inference were similar irrespective of whether the analysis used all sampled sequences or only a subset of them, although the former required substantially longer computational times. When phylodynamic analyses need to be performed under a time constraint to inform policy makers, we suggest multiple phylodynamics models to be used simultaneously for a subset of data to ascertain the robustness of inferences.

Highlighting priority areas for bovine viral diarrhea control in Italy: A phylogeographic approach

The prevalence and genetic diversity of bovine viral diarrhea virus (BVDV) in a geographic area are largely influenced by live animal trade and management practices. Despite control and eradication programs currently underway in several European countries, the risk of BVDV spread within and among countries is still present. BVDV-1 is the predominant type circulating in European cattle population. In this study, a phylogeographic analysis was applied to the BVDV-1 highest prevalent subtypes in Italy to reconstruct the origin and spatial-temporal distribution and to trace main viral flows between different locations to highlight priority areas for BVDV control. A comprehensive dataset of BVDV-1b (n = 173) and 1e (n = 172) 5' UTR sequences was analysed, including both novel and published sequences from Italy and from European countries bordering and/or with commercial cattle flows with Italy. A common phylogeographic pattern was observed for BVDV-1b and 1e subtypes: interspersion from multiple Italian areas and European countries was widespread until the end of the last century, whereas significant local clusters were observed starting from 2000. These findings support a continuous viral flow among different areas over long time scales with no evidence of significant geographical structure, while local transmission networks are limited to more recent years. Northern Italy has been confirmed as the area of origin of the main clades of both BVDV subtypes at national level, acting both as a crucial area for introduction and a maintenance source for other areas. Piedmont, Central and Southern Italian regions contributed to limited geographical distribution and local BVDV-1b and 1e persistence. On the whole, priority control measures for BVDV-1b and 1e in Italy should be focused on: i) implementation of BVDV systematic control in all Northern Italian regions to break the viral flow from larger to smaller animal populations; and ii) breaking the dynamics of infections in regions with self-maintenance of BVDV by voluntary control programs.

Carbon monoxide and biliverdin suppress bovine viral diarrhoea virus replication

Bovine viral diarrhoea virus (BVDV) causes significant economic losses to the cattle industry worldwide. Previously, we demonstrated that heme oxygenase-1 (HO-1) can inhibit BVDV replication via an unknown molecular mechanism. To elucidate the mechanism involved, we assess whether the HO-1 downstream metabolites carbon monoxide (CO), biliverdin (BV) and iron affect BVDV replication. We treated Madin-Darby bovine kidney (MDBK) cells with an exogenous CO donor, CORM-2. We found that CORM-2 but not its inactive form (iCORM-2) inhibited BVDV replication in a dose-dependent and time duration-dependent manner, suggesting a CO-specific mediation of the CORM-2 antiviral effect. Direct incubation of BVDV with high-dose CORM-2 reduced virus titres, suggesting that CORM-2 attenuates BVDV growth by both physically inactivating virus particles in the extracellular environment and affecting intracellular BVDV replication, but mainly via an intracellular mechanism. Exogenous BV treatment, both post-infection and co-incubation with BVDV, inhibited BVDV replication in a dose-dependent manner, indicating that BV has potent antiviral activity against BVDV. Direct incubation of BVDV with BV had no significant effect on virus titres, indicating that BV is not virucidal and attenuates BVDV growth by affecting intracellular BVDV replication. Furthermore, BV was found to affect BVDV penetration but not attachment. However, increased iron via addition of FeCl3 did not interfere with BVDV replication. Collectively, the results of the present study demonstrate that the HO-1 metabolites BV and CO, but not iron, inhibit BVDV replication. These findings not only provide new insights into the molecular mechanism of HO-1 inhibition of BVDV replication but also suggest potential new control measures for future BVDV infection.

Evolution of Bovine viral diarrhea virus in Canada from 1997 to 2013

Bovine viral diarrhea virus (BVDV) is a rapidly evolving, single-stranded RNA virus and a production limiting pathogen of cattle worldwide. 79 viral isolates collected between 1997 and 2013 in Canada were subjected to next-generation sequencing. Bayesian phylogenetics was used to assess the evolution of this virus. A mean substitution rate of 1.4×10^-3 substitutions/site/year was found across both BVDV1 and BVDV2. Evolutionary rates in the E2 gene were slightly faster than other regions. We also identified population structures below the sub-genotype level that likely have phenotypic implications. Two distinct clusters within BVDV2a are present and can be differentiated, in part, by a tyrosine to isoleucine mutation at position 963 in the E2 protein, a position implicated in the antigenicity of BVDV1 isolates. Distinct clustering within all sub-genotypes, particularly BVDV2a, is apparent and could lead to new levels of genotypic classification. Continuous monitoring of emerging variants is therefore necessary.

Variability and Global Distribution of Subgenotypes of Bovine Viral Diarrhea Virus

Bovine viral diarrhea virus (BVDV) is a globally-distributed agent responsible for numerous clinical syndromes that lead to major economic losses. Two species, BVDV-1 and BVDV-2, discriminated on the basis of genetic and antigenic differences, are classified in the genus Pestivirus within the Flaviviridae family and distributed on all of the continents. BVDV-1 can be segregated into at least twenty-one subgenotypes (1a–1u), while four subgenotypes have been described for BVDV-2 (2a–2d). With respect to published sequences, the number of virus isolates described for BVDV-1 (88.2%) is considerably higher than for BVDV-2 (11.8%). The most frequently-reported BVDV-1 subgenotype are 1b, followed by 1a and 1c. The highest number of various BVDV subgenotypes has been documented in European countries, indicating greater genetic diversity of the virus on this continent. Current segregation of BVDV field isolates and the designation of subgenotypes are not harmonized. While the species BVDV-1 and BVDV-2 can be clearly differentiated independently from the portion of the genome being compared, analysis of different genomic regions can result in inconsistent assignment of some BVDV isolates to defined subgenotypes. To avoid non-conformities the authors recommend the development of a harmonized system for subdivision of BVDV isolates into defined subgenotypes.

Phylogeography, phylodynamics and transmission chains of bovine viral diarrhea virus subtype 1f in Northern Italy

Bovine viral diarrhea virus (BVDV) type 1 in Italy is characterized by high genetic diversity, with at least 20 subtypes. Subtype 1f is endemic in a restricted geographic area, meaning that it has local distribution. We investigated the population dynamics of BVDV-1f in Northern Italy and characterized the transmission chains of a subset of samples from Piedmont and Aosta Valley regions. A total of 51 samples from 1966 to 2013 were considered and 5' UTR sequences were used for phylogeography. A subset of 12 samples was selected for Npro gene sequencing and further characterization of the transmission chains using both molecular and epidemiological data. Phylogeography estimated the root of BVDV-1f tree in Veneto in 1965. Four significant subclades included sequences clustering by region: Lombardy (n=3), Lombardy and Emilia-Romagna (n=7), Piedmont (n=17), Piedmont and Aosta Valley (n=21). The Piedmont-only subclade has a ladder-like branching structure, while the Piedmont and Aosta Valley subclade has a nearly complete binary structure. In the subset, the outbreak reconstruction identified one sample from Piedmont as the most probable source of infection for the Aosta Valley cases. An ad hoc questionnaire submitted to public veterinarians revealed connections between sampled and non-sampled farms by means of trades, exhibitions and markets. According to the phylogeography, BVDV-1f moved westward, entering from Veneto, and spreading to Lombardy and Emilia-Romagna in the early 1990s, and finally to Piedmont and Aosta Valley in the first decade of 2000s. Both phylogeographic analyses on the whole dataset and on the selection of Npro dataset pointed out that subtype 1f entered Aosta Valley from Piedmont. The integration of molecular and epidemiological data revealed connections between farms, and such approach should be considered in any control plan. In Aosta Valley, the study showed that BVDV1f can be controlled only monitoring the introduction of cattle from Piedmont region.

Heme Oxygenase-1 Suppresses Bovine Viral Diarrhoea Virus Replication in vitro

Viral cycle progression depends upon host-cell processes in infected cells, and this is true for bovine viral diarrhoea virus (BVDV), the causative agent of BVD that is a worldwide threat to the bovine industry. Heme oxygenase-1 (HO-1) is a ubiquitously expressed inducible isoform of the first and rate-limiting enzyme for heme degradation. Recent studies have demonstrated that HO-1 has significant antiviral properties, inhibiting the replication of viruses such as ebola virus, human immunodeficiency virus, hepatitis C virus, and porcine reproductive and respiratory syndrome virus. However, the function of HO-1 in BVDV infection is unclear. In the present study, the relationship between HO-1 and BVDV was investigated. In vitro analysis of HO-1 expression in BVDV-infected MDBK cells demonstrated that a decrease in HO-1 as BVDV replication increased. Increasing HO-1 expression through adenoviral-mediated overexpression or induction with cobalt protoporphyrin (CoPP, a potent HO-1 inducer), pre- and postinfection, effectively inhibited BVDV replication. In contrast, HO-1 siRNA knockdown in BVDV-infected cells increased BVDV replication. Therefore, the data were consistent with HO-1 acting as an anti-viral factor and these findings suggested that induction of HO-1 may be a useful prevention and treatment strategy against BVDV infection.

Genetic Variability of Bovine Viral Diarrhea Virus and Evidence for a Possible Genetic Bottleneck during Vertical Transmission in Persistently Infected Cattle

Bovine viral diarrhea virus (BVDV), a Pestivirus in the family Flaviviridae, is an economically important pathogen of cattle worldwide. The primary propagators of the virus are immunotolerant persistently infected (PI) cattle, which shed large quantities of virus throughout life. Despite the absence of an acquired immunity against BVDV in these PI cattle there are strong indications of viral variability that are of clinical and epidemiological importance. In this study the variability of E2 and NS5B sequences in multiple body compartments of PI cattle were characterized using clonal sequencing. Phylogenetic analyses revealed that BVDV exists as a quasispecies within PI cattle. Viral variants were clustered by tissue compartment significantly more often than expected by chance alone with the central nervous system appearing to be a particularly important viral reservoir. We also found strong indications for a genetic bottleneck during vertical transmission from PI animals to their offspring. These quasispecies analyses within PI cattle exemplify the role of the PI host in viral propagation and highlight the complex dynamics of BVDV pathogenesis, transmission and evolution.