Primary surveys on molecular epidemiology of bovine viral diarrhea virus 1 infecting goats in Jiangsu province, China
© The Author(s). 2016
Received: 8 January 2016
Accepted: 1 September 2016
Published: 5 September 2016
Bovine viral diarrhea virus (BVDV) is a pathogen of domestic and wildlife animals worldwide and is associated with several diseases. In China, there are many reports about genotyping of BVDV strains originated from cattle and pigs, and some of them focused on the geographical distributions of BVDV. Currently, the goat industry in Jiangsu province of China is under going a rapid expansion. Most of these goat farms are backyard enterprises and in close proximity to pig and cattle farms. However, there was very limited information about BVDV infections in goats. The objective of this study was to assess the frequency of BVDV infections of goats, the relationship of these infections to clinical signs and determine what BVDV genotypes are circulating in Jiangsu province.
From 236 goat sera collected from six regions in Jiangsu province between 2011 and 2013, BVDV-1 was identified in 29 samples from the five regions by RT-PCR. The BVDV-1 infections occurred with/without clinical signs. Eight different BVDV-1 strains were identified from these positive samples based on the 5′-untranslated region (5′-UTR) sequences, and further clustered into four BVDV-1 subtypes on the phylogenetic analysis. Three were BVDV-1b, two BVDV-1m, two BVDV-1o, and one BVDV-1p, respectively.
To our knowledge, this is the first report to investigate the occurrence of BVDV and the genotypes of BVDV infecting goats in China. The results indicated that BVDV-1 infections were indeed present and the viruses were with genetic variations in Chinese goat herds. The information would be very useful for prevention and control of BVDV-1 infections in China.
KeywordsBVDV-1 Goats Genotyping
Bovine viral diarrhea viruses (BVDVs) are members of the genus Pestivirus (family Flaviviridae), and cause considerable economic losses in cattle industry worldwide [1–3]. Traditionally, two different BVDV species have been recognized as BVDV-1 and BVDV-2, and co-circulated in cattle, sheep, goats and pigs [4–7]. In recent years, several atypical bovine pestivirus strains named HoBi-like pestivirus were detected from contaminated fetal bovine sera and infected animals . Phylogenetic analysis of these HoBi-like viruses showed a closer relationship to the recognized BVDV-1 and BVDV-2, and these viruses have been proposed as a new pestivirus species, BVDV-3 or HoBi-like pestivirus [8–11]. BVDV-1 and -2 infections are associated with several cattle diseases, including subclinical infections, immunosuppression, acute diarrhea, respiratory diseases, reproductive failures, and mucosal diseases in persistently infected calves . Reproductive disorders caused by BVDVs appear with various clinical signs including infertility, abortion, stillbirth and persistently infected calves depending on the stage of pregnancy a cow is infected at [7, 12].
Genetic typing of the pestiviruses is based on different genomic regions, such as 5′-UTR, Npro, and E2 genes [13–15]. The 5′-UTR is the most frequently analyzed portion of the genome , and provides meaningful phylogenetic inferences as this region has the highest degree of sequence conservation and efficient amplification by RT-PCR [13, 17]. Based on the 5′-UTR sequences, BVDV-1 strains have further been divided into 21 subtypes (1a-1u) [13, 17–20], BVDV-2 into four subtypes (2a-2d) [21, 22], and BVDV-3 species has also been divided into two genotypes of Brazilian and Thai origin .
Investigating on the molecular epidemiology of BVDV can provide invaluable information about the diversity of viruses in a population and, in turn, inform control programs, drive vaccine development and determine likely infection sources. So far, epidemiologic surveys on BVDV infections have been carried out in cattle, pig and yak populations in several provinces of China [19, 23–26]. These findings highlighted that corresponding measures should be taken for control and prevention against BVDV infections.
Recently, the goat industry has developed increasingly in Jiangsu province, China. There were over five million goats reared in the province during 2013 according to the official statistics. However, there were different sizes and kinds of the goat farms because the farmers had various economic conditions, most of them were backyard reared with various sizes and numbers, and some goat flocks were adjoined to pigs and cattle. The risk of cross infection might occur between cattle, pigs and goats, so we investigated the BVDV infections in the goat herds and further determined the subtypes of BVDVs on the basis of 5′-UTR genomic region amplification and sequencing.
A total of 236 goat serum samples in the laboratory were from 31 farms in six regions of Jiangsu province, China during 2011 to 2013. Of the samples examined in this study, 71 were from goats with clinical signs of diarrhea or abortion, while the remaining 165 samples from the animals with no overt clinical sign. No BVDV vaccines have been used for animal immunization in Jiangsu province yet.
RNA extraction and RT-PCR detection
Total RNA was extracted from the serum samples using TRIzol Reagent (Invitrogen), and suspended in 20 μL of ultrapure water, following the manufacturer’s recommendations. RT-PCR was carried out in a 50 μL reaction mixture containing 1× RT-PCR buffer (TAKARA, Bio, Inc.), 20 pM of Panpesti generic primers targeting 5′-UTR, 2 U of one-step Enzyme Mix (TAKARA, Bio, Inc.) and 4 μL of RNA for the expected product sizes of 290 bp , 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, and termination with a final extension of 10 min at 72 °C. The amplification products were electrophoresed in 2 % agarose gels. Part of the interesting Npro genes were amplified as the same method and the primer sequences were as followed: N1F: ATGGAGTTGATTTCAAATGAACT; N504R: GCAGCTTGAAACCCATAGAG, the expected product was 504 bp.
The expected bands were excised and recovered from the agarose gels using the Axygen gel extraction kit (Axygen, Hangzhou, China), and cloned into the pMD-18-T vectors (TAKARA, Bio, Inc.). The recombinant plasmids of the three positive bacteria were sequenced. The resultant gene sequences were compared with the sequences in the GenBank database using the Basic Local Alignment Search Tool (BLAST) and submitted to GenBank for confirmation of their identities.
Sequence data were assembled and analyzed using the DNAStar version 7.0 package (DNAStar Inc., USA). Multiple sequences were aligned with the corresponding regions of BVDV-1, BVDV-2 and atypical BVDV reference sequences retrieved from GenBank by the ClustalW program. Phylogenetic reconstructions for genotyping were compiled using the 224 bp fragment of the 5′-UTR region and 504 bp of Npro gene, phylogenetic trees were constructed by MEGA version 4.1. The reliability of the neighbor-joining tree was estimated by boostrap analysis with 1000 replicates.
The prevalence of pestivirus in Jiangsu province
Year of collection
BVDV prevalence (%)
BDV prevalence (%)
Details of the BVDV positive samples identified
Year of collection
Clinical signs (Quantity)
Diarrhea (6)/None (4)
Based on the 5′-UTR sequences, 29 BVDV samples were classified as eight BVDV-1 strains (Table 2): the strains JS12/02 and JR1-2 were identified from two abortion goats, BH789 from animals with diarrhea, HA2-12 from six diarrhea and four healthy goats, XY-3 from two weakness kids, HA6-5 from four healthy goats, the strains XZ5-9 and XZ5-8 from healthy goats, respectively. The 5′-UTR sequences of these strains were deposited in GenBank under following accession number: KP749794- KP749800, KP749802 (Table 2), the interested Npro sequences of HA2-12, JS12-02 and JR1-2 were as KX218370-218372.
Identity of 5′-UTR between the new strains and reference strains
(VEDEVAC, 289 bp)
(ZM-95, 290 bp)
(IS25CP/01, 256 bp a)
(BJ0701, 245 bp a)
The 5′-UTR phylogenetic tree was created using their nucleotide sequences and 44 BVDV reference strains covering all subtypes retrieved from GenBank, and the eight detected BVDV-1 strains were further subdivided into four subtypes (Fig. 2): BVDV-1b (XZ5-9, XZ5-8, and JR1-2), BVDV-1p (XY-3), BVDV-1m (BH789 and HA6-5) and BVDV-1o (JS12/02 and HA2-12), respectively. The same result showed in Npro phylogenetic tree including the strains HA2-12, JS12-02 and JR1-2 (Fig. 3).
Cattle are generally considered as the main host of BVDV infections, however, BVDV has also been found to infect an extensive range of animals including sheep, goats, swine, yaks, deer and members of the Camelidae family [4–7, 25, 29]. In China, BVDV infections were first detected in cattle in 1980 while the pig origin BVDV strain ZM-95 was isolated in 1995 . BVDV infections in cattle, pigs, camels and Sika deer (Cervus nippon) have continually been reported in different regions of China [6, 18, 19, 23–26, 31–34]. The BVDV prevalence of pig population has significantly increased recent years [6, 26, 35]. Until now, the BVDV infections of goats in China were not examined, and the viral epidemic situation of BVDV in the goat flocks was not clear either. In this study, the serum samples from 31 goat herds of six regions in Jiangsu were detected and the results showed that BVDV infections were indeed present as 12.3 % of the prevalence in the samples. Eight different BVDV-1 strains were identified from five regions, and the infections might occur in the goat herds no matter the animals were with clinical signs or normal situations.
Currently, ten BVDV subtypes have been detected to be circulating among multiple hosts in China: BVDV-1a, BVDV-1b, BVDV-1c, BVDV-1d, BVDV-1m, BVDV-1o, BVDV-1p, BVDV-1q, BVDV-1u and BVDV-2a [18, 19, 23, 24, 26, 31, 33]. The four BVDV-1 subtypes determined in goats of Jiangsu were consistent with the previous subtypes reported from pigs and cattle on BVDV phylogenetic groups. At the same time, the different subtypes infecting goat herds might indicate that the epidemic situation of BVDV infections in the animal population of China was very complicated. The same subtype viruses found in cattle, pigs and goats meant the possibility of interspecies transmission of the viruses.
BVDV-1b was first isolated in Jilin province, China in 1980, and continually detected in Hebei, Xinjiang, Heilongjiang, Tianjin and Qinghai [23, 24], it is considered to be the major predominant subtype in Chinese cattle [23, 24]. The BVDV-1b strains JR1-2, XZ5-8 and XZ5-9 from two different regions were 96.2–99.7 % of identity with each other, and had very high identity (>96 %) with BVDV vaccine strain VEDEVAC , especially the JR1-2 had the 100 % identity in 5′-UTR and 98.2 % in Npro gene sequences with the strain VEDEVAC, respectively, and these viruses also showed closed relationship based on the evolutionary trees. VEDEVAC was the main strain in batches of BVDV vaccine Oregon C24V in Hungary, however, no BVDV vaccines have been licensed for cattle and goats in China yet. Additionally, the current 1b strains shared 95.5–100.0 % identity with BVDV-1b isolates YL07 and TJ0802 originated from cattle in 5′-UTR gene, which showed that cross infections of the BVDV-1b viruses might occur between cattle and goats.
Chinese BVDV-1 strain ZM-95 of pig origin was first isolated in 1995 and classified into subtype BVDV-1m [30, 35]. There were no further reports until 2010, new BVDV-1m strains in Chinese cattle were detected . According to the surveys, BVDV-1m infections appeared to be more popular in China recently . As Table 3, two BVDV-1m strains BH789 and HA6-5 showed as high as 95.5 and 94.8 % identity with ZM-95, respectively. BVDV-1m infections in goats might be closely related with its prevalence in pigs since there were many pig farms in the regions. Some goats infected with BVDV-1m had diarrhea, however, the pathogenicity of BVDV-1m infections in goats was not clear yet.
It should be noticed that four isolates detected from cattle in Beijing and Tianjing of China in 2010, formed a separate branch as BVDV-1p in the 5′-UTR phylogenetic tree, three of them were detected in the aborted bovine fetuses . Now, the strain XY-3 detected from weak kids was in the same cluster with the BVDV-1p strains. It might indicate that the subtype viruses could infect animals to cause productive disorders.
This is the first report that BVDV-1 infections occurred in goats in Jiangsu province, China, and the genetic diversity of the viruses was determined as four BVDV-1 subtypes. Although numbers of the infected goats with clinical signs were limited, these infections might cause economic loss to the developing goat industry. More epidemiological investigations should be performed to better understand BVDV infections in goats as well as other farm animals, and design and develop effective control strategies.
The study showed that BVDV were circulating in Jiangsu, and present as 12.3 % of the prevalence in the goats, eight different BVDV-1 strains were identified from the samples, which were further divided into four subtypes: 1b, 1p, 1m, 1o. BVDV cross infections might be occur between goats, cattle and pigs, for control and prevention against BVDV infections, epidemiologic surveys should be carried out in these animals.
The Authors acknowledged the precious collaboration of all the farmers involved during the study. We thanked the technical staff at institute of agricultural ecomomics and information, Jiangsu academy of agricultural sciences for their support on the mapping. We thanked Dr. Kevin Coombs (University of Manitoba in Canada) read and revised the manuscript kindly and carefully.
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).
Availability of data and materials
All the data supporting the findings of the present study are included within the manuscript. Sequence data are deposited in GenBank under following accession number: KP749794- KP749800, KP749802, KX218370-KX218372.
LM and JJ designed and initiated the study. LM, WZ and LY carried out RT-PCR work; LM, WL, and FH collected samples and analyzed the sequence information; JW, SC, YW, QW, YD and YS joined in collecting the goat serum samples at the different regions; ML drafted the manuscript, JJ revised and all authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Sample collections were performed in strict accordance with the guidelines of Jiangsu Province Animal Regulations (Government Decree No 45). The owners consented to the use and disclosure of questionnaire data for the current study.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Peterhans E, Schweizer M. BVDV: a pestivirus inducing tolerance of the innate immune response. Biologicals. 2013;41(1):39–51.View ArticlePubMedGoogle Scholar
- Baker JC. Bovine viral diarrhea virus: a review. J Am Vet Med Assoc. 1987;190(11):1449–58.PubMedGoogle Scholar
- Houe H. Epidemiological features and economical importance of bovine virus diarrhoea virus (BVDV) infections. Vet Microbiol. 1999;64(2–3):89–107.View ArticlePubMedGoogle Scholar
- Julia S, Craig MI, Jimenez LS, Pinto GB, Weber EL. First report of BVDV circulation in sheep in Argentina. Prev Vet Med. 2009;90(3–4):274–7.View ArticlePubMedGoogle Scholar
- Kim IJ, Hyun BH, Shin JH, Lee KK, Lee KW, Cho KO, Kang MI. Identification of bovine viral diarrhea virus type 2 in Korean native goat (Capra hircus). Virus Res. 2006;121(1):103–6.View ArticlePubMedGoogle Scholar
- Tao J, Liao J, Wang Y, Zhang X, Wang J, Zhu G. Bovine viral diarrhea virus (BVDV) infections in pigs. Vet Microbiol. 2013;165(3–4):185–9.View ArticlePubMedGoogle Scholar
- Walz PH, Grooms DL, Passler T, Ridpath JF, Tremblay R, Step DL, Callan RJ, Givens MD. Control of bovine viral diarrhea virus in ruminants. J Vet Intern Med. 2010;24(3):476–86.View ArticlePubMedGoogle Scholar
- Schirrmeier H, Strebelow G, Depner K, Hoffmann B, Beer M. Genetic and antigenic characterization of an atypical pestivirus isolate, a putative member of a novel pestivirus species. J Gen Virol. 2004;85(Pt 12):3647–52.View ArticlePubMedGoogle Scholar
- Decaro N, Lucente MS, Mari V, Sciarretta R, Pinto P, Buonavoglia D, Martella V, Buonavoglia C. Hobi-like pestivirus in aborted bovine fetuses. J Clin Microbiol. 2012;50(2):509–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu L, Xia H, Wahlberg N, Belak S, Baule C. Phylogeny, classification and evolutionary insights into pestiviruses. Virology. 2009;385(2):351–7.View ArticlePubMedGoogle Scholar
- Stahl K, Beer M, Schirrmeier H, Hoffmann B, Belak S, Alenius S. Atypical ‘HoBi’-like pestiviruses–recent findings and implications thereof. Vet Microbiol. 2010;142(1–2):90–3.View ArticlePubMedGoogle Scholar
- Brigstocke T. Cattle Health Certification Standards and BVDV control. Vet Rec. 2010;167(21):834–5.View ArticlePubMedGoogle Scholar
- Vilcek S, Paton DJ, Durkovic B, Strojny L, Ibata G, Moussa A, Loitsch A, Rossmanith W, Vega S, Scicluna MT, et al. Bovine viral diarrhoea virus genotype 1 can be separated into at least eleven genetic groups. Arch Virol. 2001;146(1):99–115.View ArticlePubMedGoogle Scholar
- Vijayaraghavan B, Xia H, Harimoorthy R, Liu L, Belak S. Evaluation of envelope glycoprotein E(rns) of an atypical bovine pestivirus as antigen in a microsphere immunoassay for the detection of antibodies against bovine viral diarrhea virus 1 and atypical bovine pestivirus. J Virol Methods. 2012;185(2):193–8.View ArticlePubMedGoogle Scholar
- Ridpath JF, Bolin SR, Dubovi EJ. Segregation of bovine viral diarrhea virus into genotypes. Virology. 1994;205(1):66–74.View ArticlePubMedGoogle Scholar
- 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(3–4):356–60.View ArticlePubMedGoogle Scholar
- 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(3–4):386–91.View ArticlePubMedGoogle Scholar
- Deng M, Ji S, Fei W, Raza S, He C, Chen Y, Chen H, Guo A. Prevalence study and genetic typing of bovine viral diarrhea virus (BVDV) in four bovine species in China. PLoS One. 2015;10(4):e0121718.View ArticlePubMedPubMed CentralGoogle Scholar
- Gong X, Cao X, Zheng F, Chen Q, Zhou J, Yin H, Liu L, Cai X. Identification and characterization of a novel subgenotype of bovine viral diarrhea virus isolated from dairy cattle in Northwestern China. Virus Genes. 2013;46(2):375–6.View ArticlePubMedGoogle Scholar
- Tajima M, Frey HR, Yamato O, Maede Y, Moennig V, Scholz H, Greiser-Wilke I. Prevalence of genotypes 1 and 2 of bovine viral diarrhea virus in Lower Saxony, Germany. Virus Res. 2001;76(1):31–42.View ArticlePubMedGoogle Scholar
- 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(1):51–60.View ArticlePubMedGoogle Scholar
- Barros SC, Ramos F, Pauperio S, Thompson G, Fevereiro M. Phylogenetic analysis of Portuguese bovine viral diarrhoea virus. Virus Res. 2006;118(1–2):192–5.View ArticlePubMedGoogle Scholar
- 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(2):204–7.View ArticlePubMedGoogle Scholar
- 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(2–4):379–83.View ArticlePubMedGoogle Scholar
- Gong X, Liu L, Zheng F, Chen Q, Li Z, Cao X, Yin H, Zhou J, Cai X. Molecular investigation of bovine viral diarrhea virus infection in yaks (Bos gruniens) from Qinghai, China. Virol J. 2014;11:29.View ArticlePubMedPubMed CentralGoogle Scholar
- 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(3–4):490–3.View ArticlePubMedGoogle Scholar
- 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(3–4):309–23.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Rodriguez-Prieto V, Kukielka D, Rivera-Arroyo B, Martinez-Lopez B, de las Heras AI, Sanchez-Vizcaino JM, Vicente J. Evidence of shared bovine viral diarrhea infections between red deer and extensively raised cattle in south-central Spain. BMC Vet Res. 2016;12:11.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang X, Tu C, Li H, Jin K, Xuan H, Chang G, Sun H, Zhu W, Fei E, Yin Z. Detection and isolation of bovine viral diarrhea virus from classical swine fever suspected pigs. Chin J Vet Sci. 1996;16(4):341–5.Google Scholar
- 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(3):197–202.View ArticlePubMedGoogle Scholar
- Lang Y, Gao S, Du J, Shao J, Cong G, Lin T, Zhao F, Liu L, Chang H. Polymorphic genetic characterization of E2 gene of bovine viral diarrhea virus in China. Vet Microbiol. 2014;174(3–4):554–9.View ArticlePubMedGoogle Scholar
- Gao S, Luo J, Du J, Lang Y, Cong G, Shao J, Lin T, Zhao F, Belak S, Liu L, et al. Serological and molecular evidence for natural infection of Bactrian camels with multiple subgenotypes of bovine viral diarrhea virus in Western China. Vet Microbiol. 2013;163(1–2):172–6.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu X, Zhang Q, Yu X, Liang L, Xiao C, Xiang H, Tu C. Sequencing and comparative analysis of a pig bovine viral diarrhea virus genome. Virus Res. 2006;122(1–2):164–70.View ArticlePubMedGoogle Scholar