- Open Access
Characterization of one sheep border disease virus in China
Virology Journalvolume 12, Article number: 15 (2015)
Border disease virus (BDV) causes border disease (BD) affecting mainly sheep and goats worldwide. BDV in goat herds suffering diarrhea was recently reported in China, however, infection in sheep was undetermined. Here, BDV infections of sheep herds in Jiangsu, China were screened; a BDV strain was isolated and identified from the sheep flocks in China. The genomic characteristics and pathogenesis of this new isolate were studied.
In 2012, samples from 160 animals in 5 regions of Jiangsu province of China were screened for the presence of BDV genomic RNA and antibody by RT-PCR and ELISA, respectively. 44.4% of the sera were detected positively, and one slowly grown sheep was analyzed to be pestivirus RNA positive and antibody-negative. The sheep kept virus positive and antibody negative in the next 6 months of whole fattening period, and was defined as persistent infection (PI). The virus was isolated in MDBK cells without cytopathic effect (CPE) and named as JSLS12-01. Near-full-length genome sequenced was 12,227 nucleotides (nt). Phylogenetic analysis based on 5'-UTR and Npro fragments showed that the strain belonged to genotype 3, and shared varied homology with the other 3 BDV strains previously isolated from Chinese goats. The genome sequence of JSLS12-01 also had the highest homology with genotype BDV-3 (the strain Gifhorn). Experimental infections of sheep had mild clinical signs as depression and short-period mild fever (5 days). Viremia was detected in 1–7 days post-infection (dpi), and seroconversion began after 14 dpi.
This study reported the genomic and pathogenesis characterizations of one sheep BDV strain, which confirmed the occurrence of BDV infection in Chinese sheep. This sheep derived BDV strain was classified as BDV-3, together with the goat derived strains in China. These results might be helpful for further understanding of BDV infection in China and useful for prevention and control of BDV infections in the future.
Pestivirus is a genus within family Flaviviridae, comprises four recognized species, namely border disease virus (BDV) of small ruminants, bovine viral diarrhea viruses 1 and 2 (BVDV-1, BVDV-2) of cattle and classical swine fever virus (CSFV) of pigs, respectively . Recently, new pestivirus species, such as Giraffe virus, Pronghorn virus, Bungowannah virus and HoBi-like virus have been discovered [2,3].
BDV infects small ruminants (sheep and goats) and causes border disease worldwide. Several BDV strains have been proved to infect pigs and cattle under experimental or nature conditions [4,5]. BD shares some similar characteristics with BVDV infections in cattle and goats, however, with more pronounced emphasis on a wide range of productive diseases, such as abortions, stillbirths or mummified fetuses, barren ewes, malformations and the birth of weak lambs, which leads to a considerable economic effect. The affected lambs show abnormal body conformation and “hairy-shaker” syndrome named due to the hairy fleeces and tremors of the suffered lambs . Small weak lambs can be persistently infected (PI) and the virus is wide spread in all organs in infected animals. However, BDV infections in pregnant goats result exclusively in abortions and malformations in fetuses and neonates, and PI goats are rarely found out via placenta . An experiment infection was carried out in pregnant ewes with BDV-4, besides a high number of stillbirths up to 32%, significantly reduced bodyweight of lambs was also observed . Although severe clinical outbreaks of BD are unusual, several epizootics have been reported in sheep and goats [9-11].
The pestivirus genome consists of a positive single-stranded RNA approximately 12.3 kb in length, encoding a single open reading frame (ORF) flanked by 5'- and 3'- untranslated regions (UTR). The genome codes 4 structural proteins, the capsid (C) and three envelope proteins (Erns, E1 and E2), plus seven or eight non-structural proteins [12,13]. Except the 5'-UTR region, the Npro and E2 genes have also been used for genetic classification of new virus isolates [14,15]. Based on recent reports, BDV isolates have been divided into seven genotypes at least, and widely distributed in different countries, such as many European countries, Australia, New Zealand, Canada, the United States, India, Turkey, and Japan [7,16]. In 2012, BDV infections were first confirmed in several goat herds suffering serious diarrhea, and three BDV strains were isolated in China . This was the first confirmed evidence of BDV genotype 3 circulations in Chinese goats, In order to further investigate the epidemic information of BDV in goats and sheep in the same regions of China, BDV sero-epidemiological survey was carried out in our lab. One BDV strain named JSLS12-01 was isolated from one slow-grown sheep. And the genomic characteristics and pathogenesis of the isolate was determined. The data might be more supplement for BDV epidemiology, pathogenicity as well as its relationship with clinical diseases in China.
Antibody and viral detection
One hundred and sixty sera from sheep and goats in five regions of Jiangsu, China were collected and detected by BDV ELISA kits (SVANOVA). The positive rate was 44.4% (71/160), with 46.3% (63/136) of goats and 33.3% (8/24) of sheep, respectively. The positive rates varied from 16.7% (2/12) to 85.7% (12/14) with the flocks (Table 1). Sero-prevalence was at least 50% in two sheep herds and one goat herd from two regions of Jiangsu (Table 1).
In the sheep herds, an antibody-negative sheep (2 month old) from Nanjing, Jiangsu was detected positive by RT-PCR, which showed the very bright bands of 290 bp and 225 bp amplified by panpesti generic primers and BDV specific primers PBD1/PBD2, respectively. The RT-PCR products were purified and sequenced, which shared high homology with BDV strains by BLAST analysis. The RT-PCR results also showed two samples from goats were tested positive, 5'-UTR sequences were compared to sequences from GenBank, and shared 100% homology with JS12-04 strain reported previously.
The BDV RNA-positive serum was cultured and passaged in MDBK cells. After the third passage, BDV cultures have been positively detected by RT-PCR and confirmed with sequencing. However, no CPE was observed during the process of passage. The new BDV isolate was named as JSLS12-01.
Subsequent clinical observations and serological analysis
Four sampling sheep including the BDV positive were continued to rear for clinical investigation subsequently, with the same feed and management conditions. Compared to the BDV negative sheep, the infected lamb was thinner, weaker, and poorer growth, and the body weight was about 20% less at the end of fattening period (6 month long). However, no other clinical signs were observed. The four sheep were tested once 4 weeks apart, and BDV RNA was positively detected by RT-PCR and the BDV specific antibody was negative for the original BDV positive animal in the whole grown period.
Complete genome sequencing and phylogenetic analysis
RT-PCR products (Figure 1) were purified and cloned to pJET1.2 vectors for sequencing. The near full genome sequence of BDV JSLS12-01 strain was obtained and deposited in GenBank under accession number KC963426 and reported recently , the obtained genome of JSLS12-01 was 12, 227 nucleotides (nt) in length. The isolate shared 80.3% nucleotide homology and 89.9% amino acid homology with goat derived strain Gifhorn (a prototype of BDV-3), respectively (Table 2). It shared 72.2-77.6% homology with other BDV genotypes strains available in GenBank (Table 2). The homology with CSFV and BVDV strains were about 71% and 67%, respectively (Table 2). The whole genome sequence based phylogenetic analysis indicated that JSLS12-01 was classified into the same branch with Gifhorn with 100% bootstrap value, and matched up with the result of the Blast analysis. The isolate clearly differed from other BDV strains and other pestivirus species (CSFVs or BVDVs) (Figure 2).
5'-UTR and Npro sequences of the virus were aligned with the corresponding sequences of BDV reference strains and the evolutionary relationship between these isolates was estimated by phylogenetic analysis to further characterize the JSLS12-01 isolate. The phylogenetic tree based on 5'-UTR region (225 bp) showed that the JSLS12-01 strain was grouped with BDV-3 respective strains (Figure 3). Comparison of 5'-UTR sequences revealed that BDV JSLS12-01 shared nucleotide identities of 78.6%, 83.5% and 94.4% with BDV-3 Chinese strains AH12-02, AH12-01 and JS12-04, respectively. The virus was 87.7%, 87.4% and 89.0% homology with other BDV-3 strains including Gifhorn from Spanish goats, 90-F-6338 and 90-F-6227 from sheep of France, respectively (Table 3); and also 82.5% to 85.7% with BDV-1; 84.1% with BDV-2; 76.4% to 85.1% with BDV-4; 84.2% to 85.0% with BDV-5; 86.7% with BDV- 6; 71.5% to 72.9% with BDV-7; and 70.5% to 72.6% with BDV Tunisian; respectively (Table 3). The highest identity of nucleotide sequence was 94.4% with the strain JS12-04 isolated in 2012 from Chinese goats (Table 3). JSLS12-01 shared the highest homology (73.2% to 80.7%) with other BDV-3 members on the Npro sequences (Table 3). And based on the phylogenetic tree of the Npro sequences also grouped the strain JSLS12-01 into BDV-3 (data not shown).
Pathogenicity of BDV JSLS12-01
All 3 lambs (45 day old) infected with BDV JSLS12-01 cell cultures showed only moderated depression without other clinical signs. The infected animals developed high rectal temperatures (40.0-41.0°C), with peak temperatures appearing on 3–7 day-post-infection (dpi) and returned to the normal level on 8 dpi, while the control sheep kept normal during the same period (Figure 4). RT-PCR and virus isolation revealed that the virus could lead to viremia in the infected animals at 1-7dpi. The infected animals began seroconverted on 14 dpi and kept to increase in the following period (Figure 5).
Goats are major rearing small ruminants, and the numbers of sheep and rearing sheep farms are relative a few in Jiangsu and Anhui provinces of eastern China. BDV had been identified in Chinese diseased goats . In the present study, BDV sero-epidemiological survey was carried out for goat and sheep by BDV ELISA kit and RT-PCR. Serum samples from seven flocks (4 of goats and 3 of sheep) in five regions were tested positive by ELISA (Table 1), which indicated the ubiquity of BDV infection in the small ruminant flocks. Four sheep including one was BDV RNA positive with RT-PCR were kept for the next 6 month rearing. The sheep was BDV Ab negative and viremia for the observation period; and identified as one PI animal according to the reports . The PI sheep mainly showed poor-growing were agreement with BDV-4 virus infections ; however, classical clinical signs of BD such as paralysis and an abnormal fleece known as ‘hairy shaker’ were not appeared. In addition, the previous goat derived BDVs were identified from herds with diarrhea . These uncharacteristic signs emphasized the difficulty of BD diagnosis based on clinical signs and the requirement for routine and accurate laboratory tests.
Phylogenetic segregations of pestiviruses into individual species and subgroups are only identified by the branching order of the phylogenetic tree . The genetic diversity of BDV is greater than that of other pestivirus species. According to the molecular diversity, BDV have been divided to seven groups at least [7,19,20]. The near full sequence of the JSLS12-01 had a closer relationship with BDV-3 strains at gene level comparing with other pestivirus reference strains (BDV-1,BDV-2, BDV-3, BDV-4,BDV-5 and BDV-7). And the strain should be one new BDV-3 strain with a lot of nucleotide or amino acid variation to the respective BDV strains (Figure 2, Table 2). The isolate JSLS12-01 was the highest homology with the strain JS12-04 on the 5'-UTR and Npro gene compositions, and the two viruses were isolated from sheep or goats in different regions of Jiangsu province, however, the virus JSLS12-01 had much lesser identity with the other two viruses AH12-01 and AH12-02 on the gene levels (Table 3 and Figure 3), and the late two were from Anhui province (neighborhood each other). It indicated the epidemic BDV strains in China might have complex endemic situations, and have existed for a relative long period. The BDVs in China have mutated a lot, and may form different sub-branches, such as the two strains JS12-04 and JSLS12-01 isolated in Jiangsu province formed one cluster, and the strains AH12-01 and AH12-02 from Anhui were much closed each other (Figure 3).
In addition, based on sequence alignment and phylogenetic analysis, the 4 Chinese BDV strains showed the high diversity of on the 5'-UTR compositions (Table 3) even though all of them were in BDV-3 genotype (Figure 3) . The strain JSLS12-01 was only 83.5% and 78.6% homology with the strains AH12-01 and AH12-02 on 5'-UTR gene (Table 3), respectively. Furthermore, it was also found that JSLS12-01 had a very high homology of 5'-UTR sequences with other genotype viruses, such as 85.1%, 86.7% and 85% with the strains BU-1CRA22 (BDV-4), 92-F-7119 (BDV-6) and 93-F-7289 (BDV-5) (Table 3), respectively. It was true that the homology of BDV strains in the same genotype might be much less than the viruses being other genotypes. And it appeared that the BDV genotyping based on the 5'-UTR sequences could be still consummated.
Animal experiment was performed upon conventionally rearing sheep to evaluate the pathogenicity of the new BDV strain. JSLS12-01 could induce mild clinical diseases; the infected animals showed viremia at 1–7 dpi and disappeared after 7 dpi, during the same period, mild short-term pyrexia was observed (Figure 4). In addition, the infected animals seroconverted from 14 dpi (Figure 5). No other classical BD symptom was observed. The clinical symptoms of natural infected sheep mainly showed slow-growing and persisted infection. Symptoms of natural and experimental infection were very diverse, this might be related to the immune state and age of the virus infected animals and so on. To date, experimental infections of BDV have not been performed in China, the results would help to understand the pathogenicity of Chinese BDV strains in small ruminants and further systemic studies will be done in the future.
In conclusion, BDV infections of sheep were demonstrated in China, and the BDV JSLS12-01 strain was clustered in BDV-3 subgroups, the predominant genotype in China. The BDV infections of goats and sheep appeared with various clinical signs and epidemic information. More studies should be carried out to make sure of the virulence and biologic characteristics of BDV.
This study was performed in strict accordance with the guidelines of Jiangsu Province Animal Regulations (Government Decree No 45). The protocol was approved by the Committee on the Ethics of Animal Experiments of the Institute of Veterinary Medicine, Jiangsu Academy of Agricultural Sciences (JAAS No 20100604).
Clinical sampling and observation
Serum samples were collected from 3 young sheep farms and 4 goat farms (about 2 months old) from 5 different regions in Jiangsu province of China (Table 1). A total of 24 sera from sheep and 136 sera from goats were collected. The sheep tested for BDV positive were continuously carried out for subsequent clinical observations and serological analysis in next 6 months (the fattening period) after the first sampling.
BDV antibody detection was performed on serum samples of goats and sheep using a commercially available kit (SVANOVA BDV-Ab kit, SVANOVA Biotech, Uppsala, Sweden) according to the manufacturer’s instructions.
Virus isolation was carried out on the positive samples of pestivirus 5'-UTR RT-PCR as reports  and briefly introduced as following: the serum samples were centrifuged at 12,000 rpm for 20 min at 4°C and filtered through 0.22 μm filter, and subsequently inoculated onto confluent monolayers of Madin-Darby bovine kidney (MDBK) cells (obtained from China Institute of Veterinary Drugs Control) to culture for 96 hours with 1% FCS DMEM at 37°C and 5% CO2 conditions. The original cells and FCS were proven to be free of pestivirus antigen and antibodies. The isolation of BDV was checked by RT-PCR as next described using Panpesti generic primers and BDV specific primers PBD1/PBD2 [21,22].
RT-PCR detection and complete genomic sequence analysis
Sera from sheep and goats were tested for the presence of pestivirus by RT-PCR method. To detect the BDV virus, total RNA was extracted from serum samples and cell cultures by using the TRizol reagent (Invitrogen) according to the manufacturer’s instruction. RT-PCR was carried out in a 50 μL reaction mixture containing 1× RT-PCR buffer (TAKARA, Bio, Inc.), 20 pM of each primer (Panpesti generic primers and BDV specific primers PBD1/PBD2 targeting 5'-UTR, with expected product sizes of 290 bp and 225 bp, respectively) [21,22], 2 U of one-step Enzyme Mix (TAKARA, Bio, Inc.) and 4 μL of RNA. The reaction was run in a thermocycler (Mjmini, BIO-RAD) according to the following program: reverse transcription at 50°C for 30 min, denaturation at 95°C for 5 min, 35 cycles at 94°C for 30 s, 54°C for 30 s and 72°C for 45 s, terminated with a final extension of 10 min at 72°C.
For the genomic long RNA RT-PCR, the extension time was 1 kb for 1 min. Amplification products were detected by electrophoresis in 1.2% agarose gels. Positive RT-PCR fragments were purified (Axygen), cloned to pJET1.2 vector (Thermo), and then transformed to E. coli DH5α. Positive clones, as confirmed by PCR and enzyme digestion, were sequenced. Three positive clones of each RT-PCR fragment were sequenced using the appropriate PCR primers for correct check.
Briefly, six pairs of primers were designed to amplify the 6 overlapping fragments covering the virus genome, and summarized as Table 4. The retrieved sequences were edited and assembled with SeqManTM program version 5.03 of the DNASTAR package to obtain the complete sequence of this new BDV strain.
The complete coding sequence of the virus was aligned with some represented BDV, BVDV 1, BVDV 2 and CSFV strain genome sequences. The 5'-UTR and Npro sequences were analyzed with sequences of BDV reference strains using 1.83 and MEGA 4.0.2, the 225 bp 5'-UTR fragments (PBD1/PBD2 product) and 487 bp Npro gene (corresponding to 394-880 bp of Gifhorn genome) sequences were used for analysis, respectively. Phylogenetic analysis was carried out using the neighbor-joining (NJ) method using 1000 replicates for determination the bootstrap values.
Six one-month-old healthy sheep were tested negative for pestivirus (BDV and BVDV) infections by commercial ELISA kit (BDV: SVANOVA and BVDV: BIO-X) and RT-PCR mentioned above. They were further confirmed to be free of micoplasma infections by PCR. The sheep were randomly divided into two groups, with 3 animals in each group. Sheep of the experimental group was infected by intramuscular injection with 105 TCID50 of BDV JSLS12-01 cell cultures, while the sheep in control group were inoculated with PBS buffer. All animals were monitored daily for clinical signs including depression, nasal discharge, diarrhea, coughing and rectal temperature. Serum samples were collected at day −2 to 0 prior to infection and 1, 3, 5, 7, 14, 21, 28, 35 and 42 dpi. Serum samples of days 1, 3, 5, 7, 14 and 21 were tested for viremia by RT-PCR described above. And the procedures to isolate BDV from the sera were described above. Serum samples of days 0, 7, 14, 21, 35 and 42 were tested for BDV specific antibodies using commercial ELISA kit (SVANOVA).
Vilcek S, Leskova V, Meyer D, Postel A, Becher P: Molecular characterization of border disease virus strain Aveyron. Vet Microbiol 2014, 171:87–92.
Liu L, Xia H, Wahlberg N, Belák S, Baule C: Phylogeny, classification and evolutionary insights into pestiviruses. Virology 2009, 385:351–357.
Mao L, Li W, Zhang W, Yang L, Jiang J: Genome sequence of a novel Hobi-like pestivirus in China. J Virol 2012, 86:12444.
Cabezón O, Rosell R, Sibila M, Lavín S, Marco I, Segalés J: Experimental infection of pigs with Border disease virus isolated from Pyrenean chamois (Rupicapra pyrenaica). J Vet Diagn Investig 2010, 22:360–365.
Sharp MW, Rawson BC: The cost of border disease infection in a commercial flock. Vet Rec 1986, 119:128–130.
Nettleton PF, Gilray JA, Russo P, Dlissi E: Border disease of sheep and goats. Vet Res 1998, 29:327–340.
Oguzoglu TC, Tan MT, Toplu N, Demir AB, Bilge-Dagalp S, Karaoglu T, Ozkul A, Alkan F, Burgu I, Haas L, Greiser-Wilke I: Border disease virus (BDV) infections of small ruminants in Turkey: A new BDV subgroup? Vet Microbiol 2009, 135:374–379.
García-Pérez AL, Minguijón E, Estévez L, Barandika JF, Aduriz G, Juste RA, Hurtado A: Clinical and laboratorial findings in pregnant ewes and their progeny infected with Border disease virus (BDV-4 genotype). Res Vet Sci 2009, 86:345–352.
Mainar-Jaime RC, Vazquez-Boland JA: Associations of veterinary services and farmer characteristics with the prevalences of brucellosis and border disease in small ruminants in Spain. Prev Vet Med 1999, 40:193–205.
Valdazo-González B, Alvarez-Martinez M, Greiser-Wilke I: Genetic typing and prevalence of Border disease virus (BDV) in small ruminant flocks in Spain. Vet Microbiol 2006, 117:141–153.
Li W, Mao L, Zhao Y, Sun Y, He K, Jiang J: Detection of border disease virus (BDV) in goat herds suffering diarrhea in eastern China. Virol J 2013, 10:80. doi:10.1186/1743-422X-10-80.
Becher P, Orlich M, Thiel HJ: Complete genomic sequence of border disease virus, a pestivirus from sheep. J Virol 1998, 72:5165–5173.
Collett MS, Larson R, Gold C, Strick D, Anderson DK, Purchio AF: Molecular cloning and nucleotide sequence of the pestivirus bovine viral diarrhea virus. Virology 1988, 165:191–199.
Couvreur B, Letellier C, Collard A, Quenon P, Dehan P, Hamers C, Pastoret PP, Kerkhofs P: Genetic and antigenic variability in bovine viral diarrhea virus (BVDV) isolates from Belgium. Virus Res 2002, 85:17–28.
Hurtado A, García-Pérez AL, Aduriz G, Juste RA: Genetic diversity of ruminant pestiviruses from Spain. Virus Res 2003, 92:67–73.
Strong R, La Rocca SA, Ibata G, Sandvik T: Antigenic and genetic characterisation of border disease viruses isolated from UK cattle. Vet Microbiol 2010, 141:208–215.
Liu X, Mao L, Li W, Yang L, Zhang W, Wei J, Jiang J: Genome Sequence of a Border Disease Virus in China. Genome Annouc 2013, 1(6).
Becher P, Avalos Ramirez R, Orlich M, Cedillo Rosales S, König M, Schweizer M, Stalder H, Schirrmeier H, Thiel HJ: Genetic and antigenic characterization of novel pestivirus genotypes: implication for classification. Virology 2003, 311:96–104.
Dubois E, Russo P, Prigent M, Thiéry R: Genetic characterization of ovine pestiviruses isolated in France, between 1985 and 2006. Vet Microbiol 2008, 130:69–79.
Thabti F, Letellier C, Hammami S, Pépin M, Ribière M, Mesplède A, Kerkhofs P, Russo P: Detection of a novel border disease virus subgroup in Tunisian sheep. Arch Virol 2005, 150:215–229.
Albayrak H, Gumusova SO, Ozan E, Yazici Z: Molecular detection of pestiviruses in aborted foetuses from provinces in northern Turkey. Trop Anim Health Prod 2012, 44:677–680.
Vilcek S, Herring AJ, Herring JA, Nettleton PF, Lowings JP, Paton DJ: Pestiviruses isolated from pigs, cattle and sheep can be allocated into at least three genogroups using polymerase chain reaction and restriction endonuclease analysis. Arch Virol 1994, 136:309–323.
This work was supported by the Special Fund for Independent Innovation of Agricultural Science and Technology in Jiangsu province [CX (14)2090] and Jiangsu Provincial Natural Science Foundation of China [BK20130729].
The authors declare that they have no competing interests.
LM, XL, WL, and JJ participated in the design and conducted the majority of the experiments in the study. LM drafted the manuscript. LY and WZ contributed to the samples detection. JJ revised the manuscript. All authors read and approved the final manuscript.