A bovine viral diarrhea virus type 1a strain in China: isolation, identification, and experimental infection in calves
Virology Journalvolume 11, Article number: 8 (2014)
Bovine viral diarrhea virus (BVDV) is one of the most important pathogens in cattle. Previously, BVDV sub-genotypes of 1b, 1c, 1d, and 1 m were detected in China. However, isolation of BVDV type 1a from cattle has not been reported in China. In 2010, twenty nasal swabs and blood samples were collected from the cattle suspected BVDV infection in Henan province, China. A BVDV isolate was isolated using cell culture, and the pathogenesis of the virus isolate was studied.
Virus isolation was performed on MDBK cells. The virus identification was conducted by RT-PCR, neutralization test and immunofluorescence assay. In order to determine the genotype of the newly isolated virus, the 5′ un-translated region (5′UTR) of the virus isolate was cloned, sequenced and phylogenetically analyzed. To evaluate the virulence of the virus isolate, four BVDV sero-negative calves were intranasally inoculated with the virus suspension. Rectal temperatures and clinical signs were recorded daily. Blood samples were analyzed for changes in white blood cell counts, and tissue samples were taken for histopathology analysis.
A new isolate of bovine viral diarrhea virus (BVDV), named HN01, was isolated from the nasal swabs using MDBK cell culture. The HN01 strain caused cytopathic effect (CPE) in MDBK cell cultures after two passages. The virus specifically reacted to BVDV1-specific monoclonal antibody in an immunofluorescence assay. A fragment of 288 bp of genome from this isolate was amplified by the RT-PCR. Phylogenetic analysis of 5′UTR indicated that the virus was BVDV 1a. In the pathogenesis study, four calves experimentally infected with the BVDV strain developed depression, cough and other clinical signs. Calves showed high temperature over 40°C, and white blood cell counts dropped more than 40%.
A new subgenotype 1a strain of BVDV was firstly isolated from dairy cattle in China. The experimental infection showed that the virus was moderate pathogenic to cattle and can be used as a BVDV challenge virus to evaluate the efficacy of BVDV vaccines in the target animals.
Bovine viral diarrhea virus (BVDV) is a great economically pathogen in cattle and other ruminants in the world [1–10]. The virus is associated with several clinical symptoms, including diarrhea, respiratory disease, congenital malformations, reproductive disorders and mucosal disease [9, 11–14].
BVDV belongs to the genus pestivirus with classical swine fever virus and border disease virus in the family Flaviviridae. The genome of the BVDV consists of a single positive-stranded RNA, which usually have a length of 12.3 kb . Two biotypes of BVDV categorized as cytopathogenic or noncytopathogenic based on their activity in cell culture have been recognized in the past years [16, 17]. Based on the basis of the nucleotide sequence of 5′-untranslated region (5′UTR), Npro or E2 gene, BVDV strains can be divided into two different genotypes, BVDV1 and BVDV2 . Each genotype can be further divided into different subgroups, and currently at least 11 genetic subgroups of BVDV1 and three genetic subgroups of BVDV2 are identified [3, 6, 18–24]. Recently, a new virus referred to as HoBi-like BVDV3 was identified in Europe, the virus can be divided into two sub-groups, Thai origin and Brazilian origin . BVDV1 spreads worldwide in cattle population [22, 23, 26, 27]. In the case of BVDV2 species, the highest occurrence is reported in the USA and Canada [25, 28], partially in Japan [29–31], Indian , South America , and occasionally in some European countries [6, 15, 32–35].
Virulence is both important to understanding the mechanisms of pathology and selecting the challenge strains for evaluation of a vaccine. Variation in virulence among BVDV2 strains has been extensively reported [36–38], but much less information is available on variation in virulence among BVDV1 strains. To date, the strain NY-1 has been used as a challenge strain for evaluating efficacy of vaccine protection against BVDV1. While it is well characterized, the clinical presentation infected with NY-1 indicates it is more likely a low virulence strain . So investigating an efficacious challenge virus to access the vaccine efficacy is very important.
To the present, many subgenotypes of BVDV1 have been isolated and detected in China [27, 40, 41]. Based on the phylogenetic tree, the clustering of BVDV1b and BVDV 1 m were the major prevalent subgenotypes in China [27, 41, 42]. However, BVDV subgenotype 1a was not isolated from cattle in China. Moreover, pathogenesis of above all strains was seldom reported.
In this study, one virus was isolated from nasal swabs of cattle using MDBK cell cultures, and identified as a BVDV isolate by the virus neutralization test, reverse transcriptase-polymerase chain reaction (RT-PCR) method and immunofluorescence assay. To investigate the genetic subgroup of the strain, the 5′UTR gene of the virus was sequenced and compared with other 13 reference BVDV strains by phylogenetic analysis. The pathogenesis of the virus was evaluated by intranasally inoculating to four susceptible calves to assess the potential endemic risk to the cattle herd in China.
A virus isolate was isolated from the nasal swab samples after two blind passages on MDBK cells. The MDBK cells inoculated with nasal swab samples developed obvious cytopathic effects (CPE) as cells appeared aggregated and formed net-like in the monolayers within 48 hours of incubation. The non-inoculated control cells did not show the CPE. Specific fluorescent staining showed the appearance in form of granules distribution all over the cytoplasm in the infected MDBK cells, when an indirect immunofluorescence test was conducted using BVDV1 type-specific monoclonal antibodies (E2/gp53 IgG2a Isotype, VMRD, USA). The control cells did not show any fluorescence (Figure 1).
In virus neutralization assay using BVDV anti-sera, results revealed that the inoculated MDBK cells did not show CPE after mix of the virus fluids with BVDV1 positive serum, but cells inoculated with virus fluid only developed CPE. These results indicate that the virus isolate was a BVDV 1 strain.
Amplification, sequencing, and analysis of 5′UTR
The fragment of 288 bp of 5′UTR gene was amplified from the isolated virus by RT-PCR. The amplified product was purified, cloned and sequenced. The sequence of the isolated BVDV1 strain was compared to the reference sequences using BLAST program. Result revealed that the isolated BVDV 1 strain was similar with the BVDV strains Singer, with sequence identity about 95%. The sequence of 5′UTR has been deposited in GenBank with the accession No.JX878887.
Phylogenetic analysis based on 5′UTR revealed the strain HN01 was closed to BVDV1a Singer (L32875). The homology for nucleotide of 5′UTR is 95.3%, and with the strain SH1060 (JN248741) was 89.2%. This result indicates that the newly isolated BVDV strain, HN01, belongs to subgroup BVDV 1a of genotype BVDV1 (Figure 2).
Clinical signs were scored based on clinical severity, included depression, loss of appetite, nasal discharge, distress, excessive salivation, and elevated rectal temperatures. Two of the four BVDV infected calves developed significant clinical signs with nasal discharge and depression beginning at 2 days post-challenge (dpi), and subsequently 3 of four infected calves showed asthma and excessive salivation. Calves in the control group had no clinical signs during the experiment course. The rectal temperature of cattle #655 was over 40°C from 6 dpi to 8 dpi with the highest temperature of 41.4°C. The rectal temperature of cattle #642 was over 40°C at 2, 5 and 6 dpi with the highest recorded temperature of 40.1°C. The rectal temperature of animal #675 was over 40°C from 5 to 7 dpi with the highest temperature of 41.1°C. The rectal temperature of cattle #666 was elevated to 40°C from 2 to 4 dpi with the higest temperature of 40.3°C (Figure 3A).
WBC counts of challenged calves started decreasing from day 2 post-challenge. The WBC counts of 655# dropped significantly from 8.6 (1000 per mm3) at day 0 to 3.6 (1000 per mm3) at 4 dpi, a 58% reduction to the initial number. 642# decreased from 7.2 (1000 per mm3) (0 day) to 3.0 (1000 per mm3) (6 dpi), a 58% reduction to the initial number. The cattle of 675# decreased from 10.3 (1000 per mm3) (0 day) to 4.8 (1000 per mm3) (2 dpi), a 53% reduction to the initial number. The cattle of 666# decreased from 10.1 (1000 per mm3) (0 day) to 4.6 (1000 per mm3) (6 dpi), a 54.5% reduction to the initial number. Average greatest decline of WBC counts among the control animal group was less than 20%, which were significantly less (P < 0.05) than those of the challenge group (Figure 3B).
Viral detection was positive from 2 dpi to 8dpi in nasal samples. Viral shedding was detected by virus isolating in MDBK cells from nasal swabs. Two infected cattle had viral shedding as early as 2 dpi. The long shedding period reached to 8 days after inoculation, and the highest shedding peak was reached at 4-8 dpi.
Two randomly chosen calves (#675 and #642) and one mock calf (#799) were euthanized at 8 dpi for histopathology analysis. Gross pathological findings included hemorrhages in the spleen and mesenteric lymph nodes. Histopathologic changes included thymic and spleen atrophy, the lymph node was lymphadenopathy and hemorrhage, and lymphoid depletion of the Peyer′s patches in all the infected calves (Figure 4). All of the samples from the control cattle did not show pathological changes.
This study was the first report on BVDV 1a isolated from cattle in China, and the pathogenesis of the isolate was tested by experimental infection. Similar virus at the genetic subtype was detected in pigs in China . Identification of BVDV 1a in cattle and swine seemed that the virus have evolved well to replicate in different species. To our knowledge, the diversity of the genotype is a dominant character of the BVDV. Phylogenetic analysis is useful for molecular epidemiological studies and vaccine research.as well as tracing the origin of newly introduced viruses [19, 44–46]. The 5′UTR analysis is commonly used to classify the genotypes of BVDV strains. To investigate the molecular feature of the strain, 13 reference sequences (Table 1) in China, USA and Korea field isolates were used to phylogenetic analysis. A high genetically identity with the BVDV1a strain was detected among all the reference sequences. Recently, a BVDV 1a strain named SH1060 was isolate from pigs in Shanghai, China , but the virus only has the 89.2% homology of 5′UTR with the HN01 reported in this study. This may suggest that the virus have the great variation in evolution between the different kind species.
The subgroup of BVDV might be correlated with the geographical origin . Previous study indicated that the BVDV 1a subtype was predominant and widespread in neighboring countries, such as Korea and Japan [47, 48]. In China and India, BVDV 1b was the predominant subgenotypes [27, 49]. Hence, the isolated BVDV1a strain in this study may be originated from the import of infected livestock. The findings of the current study also suggest a possible role of BVDV in the etiology of diarrhea and bovine respiratory disease. In addition, the presence of BVDV could represent a clinical threat to diary and cattle industry in China.
In this study, calves infected strain HN01 of BVDV1a showed the mild clinical signs, high rectal temperature and lymphopenia, which suggested that the strain had the moderate virulence to cattle. Previous studies on the pathogenesis of BVDV1 and BVDV2 demonstrated that the low virulent strains were characterized with a subclinical course, a mild short-term pyrexia, rectal temperatures between 39.2°C and 40.0°C for 1-2 days, and a drop in circulating lymphocytes ranging between 20% and 40% . In contrast, clinical signs infected with a high virulent strain were severe with a longer term and more severe pyrexia, temperatures above 40°C, ranging upwards to 41.7°C for 3 or more days, decreases of greater than 40% in circulating lymphocytes and decreases of greater than 40% in platelets. Death losses ranging between 20% and 50% have been observed in field outbreaks associated with highly virulent BVDV2 strains [36, 39, 51]. Although the experimental animals showed over 40°C for 2 days, over 40% WBC drop, the virus was categorized as moderate virulence due to the lack of severe clinical signs. This study demonstrated that strain HN01 was an alternate challenge strain in addition to the currently available strain NY-1 for vaccine efficacy test from the USDA . Furthermore, much of the research on BVDV virulence demonstrated that the clinical signs, such as depression, febrile response and excessive lacrimation were more reproducible than hemorrhage or diarrhea. It may be suggested that low virulence strains were more dominant in nature [52–54]. Despite the lack of severe clinical signs, infection with BVDV may transiently compromise the immune system and thus possibly predisposes cattle to secondary infections . Measures such as vaccination should be taken to control BVDV1 infection.
In this study, a BVDV strain genotype 1a was firstly isolated from cattle herd in China. The study carried out a molecular epidemiology and pathogenesis on the strain. Calves inoculated with the isolate developed clinical signs, including high temperature, depression, cough and nasal discharge. The results confirmed the existence of the BVDV type 1a in China, and it may be helpful in preventing the BVD in China and a challenge virus strain for efficacy evaluations of vaccines.
Materials and methods
Samples and virus isolation
Twenty nasal swabs and blood samples were collected from the cattle that showed mild respiratory clinical signs, such as nasal discharge and cough in Henan province, China in 2010. The nasal swab samples were collected and put into a tube containing 2 ml DMEM (HyClone, USA) supplemented with 10% horse serum (Hyclone, USA), 150 μg/ml gentamicin sulfate (Sigma, USA), 7.5 μg/ml fungizone (Sigma, USA), and Streptomycin at 100 μg/ml. Three milliliter (3 ml) blood samples were collected using an EDTA vacuum blood tube from the jugular vein. All the samples were kept at 2-8°C and quickly transferred to the laboratory. Nasal swabs were inoculated into MDBK cell monolayers in 24-well tissue culture plates for virus isolation. Briefly, following centrifugation at 1500 rpm for 10 min, the samples were filtered through 0.45 μm membrane (Sigma, USA) and then inoculated onto the MDBK cell monolayer cultured in 0.5 ml DMEM (HyClone, USA) supplemented with 6% horse serum (Hyclone, USA) in 24-well cell culture plates, and incubated at 37°C, with 5% CO2 for 2 hours. Then, the supernatants were discarded, and plates were rinsed twice with PBS (pH7.2, 0.01 mol/L), and 1 ml DMEM (HyClone, USA), with 3.5% horse serum was added. The infected-MDBK cell plates were checked daily and appearance of cytopathic effects (CPE) was observed and recorded. If the CPE was not found, the cultures were frozen and thawed twice and the clarified supernatant was passaged three times in MDBK cells. Un-infected MDBK cultures were included as negative controls and MDBK cells inoculated with BVDV NM01 strain, which was previously isolated and identified by our laboratory, was used as positive control. After 3-4 days of incubation at 37°C, with 5% CO2 supply, the virus was confirmed by immunofluorescence on cell monolayers as described below.
Virus neutralization test
To identify the isolated virus, virus neutralization (VN) test was conducted using antisera raised in cattle against BVDV1a NM01, a vaccine strain in our laboratory. Briefly, virus fluids were mixed with the equal volume hyperimmune antisera (antibody VN titer ≥ 1∶64) and incubated at 37°C for one hour, and then the mixture was transferred into MDBK cell cultures in the 96-well tissue culture plates. After 4 days of incubation at 37°C, with 5% CO2 supply, the plates were observed for CPE to determine if the virus isolate was neutralized by the BVDV antisera. The antisera positive control and negative control were also included in this test.
To detect the BVDV in the infected-MDBK culture, an immunofluorescent assay (IFA) was conducted as follows. Briefly, cell lysate was added to each of four wells of a 96-well MDBK tissue culture plate. Positive virus (NM01 strain) and DMEM media (negative control) were also added to 4 wells respectively in the test. Then plates were incubated for 4 days at 37°C, with 5% CO2 supply. Then, the supernatants were dumped and plates were fixed with the mixture of 80% acetone and 20% methyl alcohol; the fixed plates were incubated with monoclonal antibody specific to anti-BVDV1 (E2/gp53 IgG2a Isotype, VMRD, USA). Finally, the rabbit anti-mouse fluorescein isothiocyanate (FITC)-conjugated immunoglobulin G (IgG) (Sigma, USA) was added to the plates, then incubated in a 37°C humid box for an hour. The plates were examined for fluorescence under a fluorescent microscope (Zeiss Axioskop-40, Germany).
PCR amplification and sequencing
Primers (Forward: 5′ATGCCCTTAGTAGGACTAGCA3′; Reverse: 5′TCAACTCCATGTGCCATGTAC3′) for 5′UTR gene sequence of BVDV were designed according to the BVDV strain Singer (GenBank accession number L32875) and strain NADL (GenBank accession number M31182). Total BVDV RNA was extracted from infected MDBK cell fluids using Trizol reagent (Invitrogen, China) according to the manufacturer′s instructions. The reverse transcription polymerase chain reaction (RT-PCR) was carried out using 5 μl of RNA extraction samples and 1 μl of BVDV reverse primer as reverse transcription primer in a 0.2 ml tube for 5 min at 70°C; and the tube was added with 1 μl (200U) M-MLV reverse transcriptase RNase inhibitor and 2 μl RNAse-free water, and then was incubated for one hour at 42°C. The amplification of cDNA by PCR was carried out in a total volume of 50 μl solution containing 1 μl cDNA, 5 μl 10 × Buffer, 3 μl 2.5 mmol/L dNTP, 0.5 μl each primer, 0.5 μl Platinum Pfx DNA polymerase and 39.5 μl sterile water. Amplification was performed as follows: 95°C for 5 min, and then submitted to 30 cycles of amplification. The conditions for the amplification were 1 min at 94°C, 1 min at 55°C, and 1 min at 72°C, then with a final extension at 72°C for 10 min. The PCR products were detected by electrophoresis through a 1% agarose gel stained with ethydium bromide and visualized under UV light. The target fragment of PCR products were purified and cloned into pMD18-T vector (TaKaRa, Japan), and then transfected into E.coli DH5α cells. The positive clones were selected for sequencing by Shanghai Invitrogen Biotechnology Co., Ltd.
The sequence analysis was performed by computer software DNASTAR (Madison, Wis USA) and MEGA5.1. The 5′UTR gene sequence of the virus isolate was evaluated for further sub-genotype classification and compared to 13 BVDV reference sequences (Table 1) from NCBI database information to develop a phylogenetic tree. Percent identities and divergences were calculated using Neighbor-Joining method of Mega5.1 software.
Calves and housing
Six 4-month old calves were obtained from a calf farm in Inner Mongolian, China. All animals were negative to BVDV and apparently healthy, presenting no signs of depression, cough or other health disorders. The selected animals were transported to an animal facility in Inner Mongolian. All animal experiments were approved by the Institutional Animal Care and Use Committee of Jilin University.
Virus infection and temperature recording
Three days prior to challenge, all animals were transferred to a bio-level 3 safety facility. Four calves were intranasally inoculated with 6 ml of cell culture medium containing 106.8TCID50 per milliliter of the HN01 strain, and other two calves were inoculated with the cell culture medium to serve as the controls. All the culture medium were demonstrated free of adventitious pathogens. The challenge procedure was performed by spraying 3 ml of virus samples into each nostril, using an atomizer, and then the calves were monitored for 8 days. The temperature was taken two times at the same time every day by investigators who were unaware of the treatment codes in each study. Three randomly chosen calves were euthanized at 8 dpi for histopathology analysis.
Calves were clinically observed daily from 1 to 8 days post-challenge. Clinical signs including depression, cough, asthma, and other respiratory disease were recorded and scored using a scale of 0 – 3 (0 = absence; 1 = mild; 2 = moderate; 3 = severe).
White blood cell (WBC) counts were conducted from 2 days pre-challenge through 8 days post-challenge and white blood cell (WBC) counts were conducted by a Vetscan HM5 veterinary hematology system (Abaxis, USA). Deep nasal swab specimens were collected at 1 day prior to challenge through 8 days post-challenge. After collection, swabs were placed into a tube containing 3 ml of transport medium, consisting of DMEM (HyClone, USA) supplemented with 10% horse serum (Hyclone, USA), 150 μg/ml gentamicin sulfate (Sigma, USA), 7.5 μg/ml fungizone and streptomycin at 100 μg/ml. The samples were kept at 2-8°C and quickly transferred to the laboratory. Upon arriving in laboratory, all swab specimens were stored at -70°C or below until they were cultured for virus isolation.
On days 8 dpi, two calves in the infected group and one in the control group were euthanized; tissue samples of liver, spleen, lung, heart, kidney, intestine, mandibular lymph node, and mesenteric lymph node were collected and fixed in 10% buffered formalin for histopathologic analysis in Inner Mongolia agricultural university.
Borsberry S: Challenges of BVDV eradication in ‘closed’ cattle herds. Vet Rec 2012, 170: 652.
Czopowicz M, Kaba J, Schirrmeier H, Bagnicka E, Szalus-Jordanow O, Nowicki M, Witkowski L, Frymus T: Serological evidence for BVDV-1 infection in goats in Poland - short communication. Acta Vet Hung 2011, 59: 399-404. 10.1556/AVet.2011.022
Pogranichniy RM, Schnur ME, Raizman EA, Murphy DA, Negron M, Thacker HL: Isolation and genetic analysis of bovine viral diarrhea virus from infected cattle in indiana. Vet Med Int 2011, 2011: 925910.
Uttenthal A, Stadejek T, Nylin B: Genetic diversity of bovine viral diarrhoea viruses (BVDV) in Denmark during a 10-year eradication period. APMIS 2005, 113: 536-541. 10.1111/j.1600-0463.2005.apm_227.x
Canal CW, Strasser M, Hertig C, Masuda A, Peterhans E: Detection of antibodies to bovine viral diarrhoea virus (BVDV) and characterization of genomes of BVDV from Brazil. Vet Microbiol 1998, 63: 85-97. 10.1016/S0378-1135(98)00232-6
Wolfmeyer A, Wolf G, Beer M, Strube W, Hehnen HR, Schmeer N, Kaaden OR: Genomic (5′UTR) and serological differences among German BVDV field isolates. Arch Virol 1997, 142: 2049-2057. 10.1007/s007050050222
Siegwart N, Hilbe M, Hassig M, Braun U: Increased risk of BVDV infection of calves from pregnant dams on communal Alpine pastures in Switzerland. Vet J 2006, 172: 386-388. 10.1016/j.tvjl.2005.07.018
Otter A, Welchman Dde B, Sandvik T, Cranwell MP, Holliman A, Millar MF, Scholes SF: Congenital tremor and hypomyelination associated with bovine viral diarrhoea virus in 23 British cattle herds. Vet Rec 2009, 164: 771-778. 10.1136/vr.164.25.771
Hessman BE, Sjeklocha DB, Fulton RW, Ridpath JF, Johnson BJ, McElroy DR: Acute bovine viral diarrhea associated with extensive mucosal lesions, high morbidity, and mortality in a commercial feedlot. J Vet Diagn Invest 2012, 24: 397-404. 10.1177/1040638711436244
Yesilbag K, Gungor B: Seroprevalence of bovine respiratory viruses in North-Western Turkey. Trop Anim Health Prod 2008, 40: 55-60. 10.1007/s11250-007-9053-x
Bielanski A, Loewen KS, Del Campo MR, Sirard MA, Willadsen S: Isolation of bovine herpesvirus-1 (BHV-1) and bovine viral diarrhea virus (BVDV) in association with the in vitro production of bovine embryos. Theriogenology 1993, 40: 531-538. 10.1016/0093-691X(93)90406-U
Flores EF, Gil LH, Botton SA, Weiblen R, Ridpath JF, Kreutz LC, Pilati C, Driemeyer D, Moojen V, Wendelstein AC: Clinical, pathological and antigenic aspects of bovine viral diarrhea virus (BVDV) type 2 isolates identified in Brazil. Vet Microbiol 2000, 77: 175-183. 10.1016/S0378-1135(00)00274-1
Fulton RW, Ridpath JF, Saliki JT, Briggs RE, Confer AW, Burge LJ, Purdy CW, Loan RW, Duff GC, Payton ME: Bovine viral diarrhea virus (BVDV) 1b: predominant BVDV subtype in calves with respiratory disease. Can J Vet Res 2002, 66: 181-190.
Vilcek S, Durkovic B, Kolesarova M, Greiser-Wilke I, Paton D: Genetic diversity of international bovine viral diarrhoea virus (BVDV) isolates: identification of a new BVDV-1 genetic group. Vet Res 2004, 35: 609-615. 10.1051/vetres:2004036
Marques Antunes de Oliveira A, Stalder H, Peterhans E, Sauter KS, Schweizer M: Complete genome sequences of both biotypes of a virus pair of bovine viral diarrhea virus subgenotype 1k. Genome Announc 2013., 1:
Corapi WV, Donis RO, Dubovi EJ: Monoclonal antibody analyses of cytopathic and noncytopathic viruses from fatal bovine viral diarrhea virus infections. J Virol 1988, 62: 2823-2827.
Harding MJ, Cao X, Shams H, Johnson AF, Vassilev VB, Gil LH, Wheeler DW, Haines D, Sibert GJ, Nelson LD, et al.: Role of bovine viral diarrhea virus biotype in the establishment of fetal infections. Am J Vet Res 2002, 63: 1455-1463. 10.2460/ajvr.2002.63.1455
Ridpath JF, Bolin SR: Differentiation of types 1a, 1b and 2 bovine viral diarrhoea virus (BVDV) by PCR. Mol Cell Probes 1998, 12: 101-106. 10.1006/mcpr.1998.0158
Flores EF, Ridpath JF, Weiblen R, Vogel FS, Gil LH: Phylogenetic analysis of Brazilian bovine viral diarrhea virus type 2 (BVDV-2) isolates: evidence for a subgenotype within BVDV-2. Virus Res 2002, 87: 51-60. 10.1016/S0168-1702(02)00080-1
Nagai M, Sato M, Nagano H, Pang H, Kong X, Murakami T, Ozawa T, Akashi H: Nucleotide sequence homology to bovine viral diarrhea virus 2 (BVDV 2) in the 5′ untranslated region of BVDVs from cattle with mucosal disease or persistent infection in Japan. Vet Microbiol 1998, 60: 271-276. 10.1016/S0378-1135(98)00158-8
Peletto S, Zuccon F, Pitti M, Gobbi E, Marco LD, Caramelli M, Masoero L, Acutis PL: Detection and phylogenetic analysis of an atypical pestivirus, strain IZSPLV_To. Res Vet Sci 2012, 92: 147-150. 10.1016/j.rvsc.2010.10.015
Gao Y, Wang S, Du R, Wang Q, Sun C, Wang N, Zhang P, Zhang L: Isolation and identification of a bovine viral diarrhea virus from sika deer in china. Virol J 2011, 8: 83. 10.1186/1743-422X-8-83
Li QC, Miao LG, Li HT, Liu YH, Zhang GL, Xiao JM: Genome sequencing and analysis of the bovine viral diarrhea virus-2 strain JZ05-1 isolated in China. Bing Du Xue Bao 2010, 26: 238-243.
Vilcek S, Durkovic B, Kolesarova M, Paton DJ: Genetic diversity of BVDV: consequences for classification and molecular epidemiology. Prev Vet Med 2005, 72: 31-35. discussion 215-219 10.1016/j.prevetmed.2005.08.004
Bauermann FV, Flores EF, Ridpath JF: Antigenic relationships between Bovine viral diarrhea virus 1 and 2 and HoBi virus: possible impacts on diagnosis and control. J Vet Diagn Invest 2012, 24: 253-261. 10.1177/1040638711435144
Zhu LQ, Ren M, Lin YQ, Ding XY, Zhang GP, Zhao X, Zhu GQ: Identification of a bovine viral diarrhea virus 2 isolated from cattle in China. Acta Virol 2009, 53: 131-134. 10.4149/av_2009_02_131
Xue F, Zhu YM, Li J, Zhu LC, Ren XG, Feng JK, Shi HF, Gao YR: Genotyping of bovine viral diarrhea viruses from cattle in China between 2005 and 2008. Vet Microbiol 2010, 143: 379-383. 10.1016/j.vetmic.2009.11.010
Pellerin C, van den Hurk J, Lecomte J, Tussen P: Identification of a new group of bovine viral diarrhea virus strains associated with severe outbreaks and high mortalities. Virology 1994, 203: 260-268. 10.1006/viro.1994.1483
Nagai M, Ito T, Sugita S, Genno A, Takeuchi K, Ozawa T, Sakoda Y, Nishimori T, Takamura K, Akashi H: Genomic and serological diversity of bovine viral diarrhea virus in Japan. Arch Virol 2001, 146: 685-696. 10.1007/s007050170139
Minami F, Nagai M, Ito M, Matsuda T, Takai H, Jinkawa Y, Shimano T, Hayashi M, Seki Y, Sakoda Y, et al.: Reactivity and prevalence of neutralizing antibodies against Japanese strains of bovine viral diarrhea virus subgenotypes. Comp Immunol Microbiol Infect Dis 2011, 34: 35-39. 10.1016/j.cimid.2009.10.007
Shimazaki T, Sekiguchi H, Nakamura S, Taguchi K, Inoue Y, Satoh M: Segregation of bovine viral diarrhea virus isolated in Japan into genotypes. J Vet Med Sci 1998, 60: 579-583. 10.1292/jvms.60.579
Vilcek S, Mojzisova J, Bajova V, Paulik S, Strojny L, Durkovic B, Hipikova V: A survey for BVDV antibodies in cattle farms in Slovakia and genetic typing of BVDV isolates from imported animals. Acta Vet Hung 2003, 51: 229-236. 10.1556/AVet.51.2003.2.11
Pizarro-Lucero J, Celedon MO, Aguilera M, de Calisto A: Molecular characterization of pestiviruses isolated from bovines in Chile. Vet Microbiol 2006, 115: 208-217. 10.1016/j.vetmic.2006.02.009
Yilmaz H, Altan E, Ridpath J, Turan N: Genetic diversity and frequency of bovine viral diarrhea virus (BVDV) detected in cattle in Turkey. Comp Immunol Microbiol Infect Dis 2012, 35: 411-416. 10.1016/j.cimid.2012.03.006
Giammarioli M, Pellegrini C, Casciari C, Rossi E, De Mia GM: Genetic diversity of bovine viral diarrhea virus 1: Italian isolates clustered in at least seven subgenotypes. J Vet Diagn Invest 2008, 20: 783-788. 10.1177/104063870802000611
Liebler-Tenorio EM, Ridpath JE, Neill JD: Distribution of viral antigen and development of lesions after experimental infection with highly virulent bovine viral diarrhea virus type 2 in calves. Am J Vet Res 2002, 63: 1575-1584. 10.2460/ajvr.2002.63.1575
Carman S, van Dreumel T, Ridpath J, Hazlett M, Alves D, Dubovi E, Tremblay R, Bolin S, Godkin A, Anderson N: Severe acute bovine viral diarrhea in Ontario, 1993-1995. J Vet Diagn Invest 1998, 10: 27-35. 10.1177/104063879801000106
Odeon AC, Kelling CL, Marshall DJ, Estela ES, Dubovi EJ, Donis RO: Experimental infection of calves with bovine viral diarrhea virus genotype II (NY-93). J Vet Diagn Invest 1999, 11: 221-228. 10.1177/104063879901100303
Ridpath JF, Neill JD, Peterhans E: Impact of variation in acute virulence of BVDV1 strains on design of better vaccine efficacy challenge models. Vaccine 2007, 25: 8058-8066. 10.1016/j.vaccine.2007.09.014
Zhu LQ, Lin YQ, Ding XY, Ren M, Tao J, Wang JY, Zhang GP, Zhu GQ: Genomic sequencing and characterization of a Chinese isolate of Bovine viral diarrhea virus 2. Acta Virol 2009, 53: 197-202. 10.4149/av_2009_03_197
Zhong F, Li N, Huang X, Guo Y, Chen H, Wang X, Shi C, Zhang X: Genetic typing and epidemiologic observation of bovine viral diarrhea virus in Western China. Virus Genes 2011, 42: 204-207. 10.1007/s11262-010-0558-4
Gao S, Du J, Shao J, Lang Y, Lin T, Cong G, Zhao F, Belak S, Liu L, Chang H, Yin H: Genome analysis reveals a novel genetically divergent subgenotype of bovine viral diarrhea virus in China. Infect Genet Evol 2013.
Deng Y, Sun CQ, Cao SJ, Lin T, Yuan SS, Zhang HB, Zhai SL, Huang L, Shan TL, Zheng H, et al.: High prevalence of bovine viral diarrhea virus 1 in Chinese swine herds. Vet Microbiol 2012, 159: 490-493. 10.1016/j.vetmic.2012.04.023
Nagai M, Hayashi M, Itou M, Fukutomi T, Akashi H, Kida H, Sakoda Y: Identification of new genetic subtypes of bovine viral diarrhea virus genotype 1 isolated in Japan. Virus Genes 2008, 36: 135-139. 10.1007/s11262-007-0190-0
Behera SP, Mishra N, Vilcek S, Rajukumar K, Nema RK, Prakash A, Kalaiyarasu S, Dubey SC: Genetic and antigenic characterization of bovine viral diarrhoea virus type 2 isolated from cattle in India. Comp Immunol Microbiol Infect Dis 2011, 34: 189-196. 10.1016/j.cimid.2010.11.002
Mishra N, Dubey R, Rajukumar K, Tosh C, Tiwari A, Pitale SS, Pradhan HK: Genetic and antigenic characterization of bovine viral diarrhea virus type 2 isolated from Indian goats (Capra hircus). Vet Microbiol 2007, 124: 340-347. 10.1016/j.vetmic.2007.04.023
Yamamoto T, Kozasa T, Aoki H, Sekiguchi H, Morino S, Nakamura S: Genomic analyses of bovine viral diarrhea viruses isolated from cattle imported into Japan between 1991 and 2005. Vet Microbiol 2008, 127: 386-391. 10.1016/j.vetmic.2007.08.020
Oem JK, Hyun BH, Cha SH, Lee KK, Kim SH, Kim HR, Park CK, Joo YS: Phylogenetic analysis and characterization of Korean bovine viral diarrhea viruses. Vet Microbiol 2009, 139: 356-360. 10.1016/j.vetmic.2009.06.017
Mishra N, Vilcek S, Jain P, Pitale SS, Pradhan HK: Genetic analysis of Indian bovine viral diarrhea virus 1 isolates in N(pro) and entire gene region coding structural proteins. Acta Virol 2006, 50: 39-44.
Liebler-Tenorio EM, Greiser-Wilke I, Pohlenz JF: Organ and tissue distribution of the antigen of the cytopathogenic bovine virus diarrhea virus in the early and advanced phase of experimental mucosal disease. Arch Virol 1997, 142: 1613-1634. 10.1007/s007050050184
Liebler-Tenorio EM, Ridpath JF, Neill JD: Lesions and tissue distribution of viral antigen in severe acute versus subclinical acute infection with BVDV2. Biologicals 2003, 31: 119-122. 10.1016/S1045-1056(03)00026-5
Liebler-Tenorio EM, Ridpath JF, Neill JD: Distribution of viral antigen and development of lesions after experimental infection of calves with a BVDV 2 strain of low virulence. J Vet Diagn Invest 2003, 15: 221-232. 10.1177/104063870301500303
Hamers C, Couvreur B, Dehan P, Letellier C, Lewalle P, Pastoret PP, Kerkhofs P: Differences in experimental virulence of bovine viral diarrhoea viral strains isolated from haemorrhagic syndromes. Vet J 2000, 160: 250-258. 10.1053/tvjl.2000.0500
Castrucci G, Osburn BI, Ferrari M, Traldi V: An experimental contribution to the study of the pathogenesis of bovine viral diarrhea virus infection. Comp Immunol Microbiol Infect Dis 1992, 15: 163-169. 10.1016/0147-9571(92)90089-A
I would like to express my gratitude to all the people who supported and assisted us in conducting this study. This work was performed under the auspices of the “Special Fund for Agro-scientific Research in the Public Interest project, No. 201003060” and we are thankful to SINOVET (BEIJING) BIOTHCHNOLOGY CO., LTD for providing the necessary facilities to carry out this work.
The authors declare that they have no competing interests.
WW participated in the design and conducted the majority of the experiments in the study and drafted the manuscript. XS and QT participated in the molecular genetic studies and the sequence alignment. YW and YJ took the animal experiment and collected samples. MX was performed the statistical analysis. WX and HW conceived of the study, and participated in its design. All authors read and approved the final manuscript.