Genetic Diversity of Bovine Viral Diarrhea Virus Infection in Goats in Southwestern China.
Bovine viral diarrhea virus-1 (BVDV-1) and BVDV-2, together with classical swine fever virus and border disease virus (BDV), are currently classified as the genus Pestivirus and family Flaviviridae . A new unclassified Pestivirus species, HoBi-like virus, was first identified in fetal bovine serum from Brazil and was related to BVDV at the antigenic and genetic levels. Thus, the viruses were referred to as "BVDV-3" . To date, a growing number unassigned atypical Pestivirus species have been detected in Europe, North America, South America, and Asia ; however, such emerging pestiviruses have not been officially recognized. To resolve the taxonomy of pestiviruses, a new taxonomy for the genus Pestivirus (family Flaviviridae) was proposed, where all species were redesignated as Pestivirus A-Pestivirus K .
The natural host of BVDV is cattle though it is also capable of infecting sheep, goats, piglets, and other domestic and wild ruminants by persistently infected (PI) animals and transiently infected individuals . PI cattle are the main source of BVDV; however, persistent infections also occur in nonbovid hosts such as sheep and deer. BVDV infections in goats typically result in reproductive diseases and viable PI goats are rare . There is increasing evidence that BVDV infections occur in a wider range of species, including mountain goats and domestic goats that have been reservoirs for BVDV infection and seriously impact the domestic livestock industry.
BVDV has a single stranded positive sense RNA genome of 12.3-12.5 kb. The full-length viral genome is flanked at both ends by 5' and 3' untranslated regions (UTRs) and contains a single open reading frame that encodes a single polyprotein that is approximately 4000 amino acids. The genetic diversity of novel Pestivirus strains is usually based on the characteristics of partial sequences from the 5 -UTR, [N.sup.pro] or E2 regions of the genome. However, the 5'-UTR of Pestiviruses has the highest degree of sequence conservation and is therefore most frequently used for phylogenetic analysis. According to sequence comparison analyses of this region, there are 21 subtypes of BVDV-1(1a-1u) [7-9] and four subtypes of BVDV-2 (2a-2d) [10-12].
A previous study based on a strain of BVDV-1 that infected goats in eastern China showed that a high proportion (29/236) of BVDV-positive goats came from eastern China . However, studies regarding the genetic diversity of BVDV in diseased goats in southwestern China remain rare. Thus, this study detected and genotyped BVDV from infected goats in this region.
2. Materials and Methods
2.1. Samples. A total of 217 blood samples taken from diseased goats that presented with symptoms, such as diarrhea, respiratory tract infection, and mucositis, were collected between March 2013 and April 2016, from multiple goat farms in southwestern China.
2.2. BVDV Antigen-Capture ELISA (ACE). Sera were screened for BVDV using IDEXX BVD Ag Test (IDEXX laboratories Inc., Shanghai, China), which can be applied for detection of genetically diverse BVDV strains, according to the manufacturer's instructions. The antigen-positive samples were then followed up with RT-PCR detection.
2.3. RNA Isolation and RT-PCR Detection of BVDV-1, BVDV-2, BVDV-3, and BDV. Viral RNA was extracted from 200 [micro]L aliquots of goat serum using the QIAamp Viral RNA Mini Kit (QIAGEN China, Shanghai, China), following the manufacturer's recommendations. Four different PCR assays for BVDV-1 , BVDV-2 , BVDV-3 , and BDV  were performed to amplify the respective targets. The PCR products were purified and sequenced by the JIELI Biology Company (Shanghai, China).
2.4. Phylogenetic Analysis. The 5'-UTR genomic sequences obtained were compared with the relevant sequences published in GenBank (http://www.ncbi.nlm.nih.gov/genbank) using the Basic Local Alignment Search Tool (BLAST). The sequences in this study were trimmed and analyzed with MEGA7.0 software  to obtain 224 nucleotide sequences, corresponding region of BVDV. The 224 nucleotide sequences were then subjected to alignment using the ClustalW program using MegAlign in the Lasergene 12.0 software (DNASTAR Inc. Madison, WI, USA), and genetic diversity analysis was performed using MEGA7.0 software. The evolutionary distances in the phylogenetic tree were computed using the Tamura 3-parameter model with 1000 bootstrap replicates, each constructed by the Neighbor-Joining method .
The results of ACE showed that out of the 217 diseased goat sera samples tested using ACE, 30 were antigen-positive for BVDV. Antigen-positive samples were confirmed as BVDV-1 using BVDV-1 specific nested RT-PCR. Fifteen additional samples that were close to the cutoff value of ACE were also tested with nested RT-PCR, eight of which were positive for BVDV-1. Taken together, 38 sequences from 217 total samples corresponded to BVDV-1 by nested RT-PCR, which resulted in a 197-bp amplicon from the BVDV-1 5'-UTR. However, BVDV-2, BVDV-3, and BDV were not detected in any of these samples using RT-PCR. These results indicated that the average infection rate of BVDV-1 in goats was 17.51% (38/217) in southwestern China.
BVDV-1 was detected in 12 of 38 BVDV-positive samples by first round RT-PCR, which produced a 305 bp product that was not detected in the other 26 BVDV-1-positive samples. The sequences detected with the first round RT-PCR in the 12 BVDV-1 positive samples were deposited in GenBank under accession numbers MG323516-MG323527 (Table 1).
Twelve of those sequences were selected for the construction of a phylogenetic tree and were divided into five different subtypes. Five isolates (XC1522, HL1407, DC1301, SC1601, and MN1319) were clustered within BVDV-1b, one (DS1324) was typed as BVDV-1a, one (FS1573) as BVDV-1c, one (FS1426) as BVDV-1p, and four (FS1605, CF1374, DF1312, and DB1503) were characterized as BVDV-1m. No BVDV-2, BVDV-3, or BDV was detected in field isolates by phylogenetic analysis.
In this study, we identified 12 selected BVDV-1-positive samples that represented different flocks or geographical origins of blood samples collected between 2013 and 2016.
Experimental or natural infection in goats has been previously confirmed . Goats infected with BVDV represent a great risk for disease transmission, thus making the identification of BVDV infection in goats an important endeavor. This investigation described the genetic diversity of BVDV isolated from goats in southwestern China and demonstrated the occurrence of at least five subtypes of BVDV, including BVDV-1a, 1b, 1c, 1m, and 1p.
Among the BVDV-1 types, the most frequent subtype was BVDV-1b (n=5) (Table 1 and Figure 1). Five BVDV-1b-positive samples from different counties of southwestern China shared 98.1%-99.0% identity with each other and had 91.6%-98.0% identity with the reference strain, BVDV-1b 24-15 (AF298060), and Osloss (M96687). Those five strains were also clustered with the two BVDV-1b reference strains (Figure 1). In cattle, BVDV-1b was first detected in the tissue of a bovine fetus in 1980 in China and continually isolated and identified in Tianjin, Hebei, Gansu, Xinjiang, and Heilongjiang province from then onward [8, 20, 21].
The 1b subtype is widely considered to be a major contributor to infection in Chinese cattle [20, 22]. Similarly, BVDV-1b is also a main subtype in Chinese goat herds . BVDV transmission among different animal species has been proven and, therefore, increased surveillance of BVDV infection in goats should be an important factor to consider in the control and understanding of BVDV infection in cattle.
BVDV-1m is another major subtype in this region (Figure 1). Four BVDV-1-positive samples identified in this study were classified as BVDV-1m. The isolates, FS1605, CF1374, DF1312, and DB1503, had 94.5%-98.7% homology with each other and shared 95.4%-97.0% homology with the BVDV-1m reference strain ZM-95 (AF526381) (Table 1). Therefore, those four positive samples were clustered into BVDV-1m. BVDV-1m was first isolated from a diseased pig in 1995 in China . According to a survey conducted between 2005 and 2013, this subtype has been widely distributed in China and is the most prevalent among BVDV-1 subtypes .
Our positive samples, DS1324, FS1573, and FS1426, respectively, shared 87.9%, 97.5%, and 98.4% sequence identity with the reference strains BVDV-1a NADL (M31182), BVDV-1c Manasi (EU159702), Bega (AF049221), BVDV-1p TJ06 (GU120246), and TJ07 (GU120247) (Table 1). These positive samples were classified into BVDV-1a, 1c, and 1p, which are not commonly found.
In summary, this study demonstrated that BVDV-1 was a common causative agent for infection of goats in southwestern China and that 17.51% of the detected samples from ill goats were BVDV-1 positive. However, we did not assess samples from the healthy goats and cannot determine the prevalence of BVDV in healthy goats. To further our understanding of the genetic diversity of BVDV strains detected, all antigen-positive samples were confirmed and characterized with RT-PCR as being infected by BVDV-1, but not BVDV-2, BVDV-3, or BDV. Twelve BVDV-1-positive samples were identified from these RT-PCR positive samples and were clustered into the six subtypes BVDV-1a, 1b, 1c, 1m, and 1p. Due to the possibility of cross infections of BVDV that might occur among cattle, piglets, and goats, much more emphasis should be placed on epidemiology surveys in these animals. This is the first description of the genetic diversity of BVDV-1-positive samples collected from goats in southwestern China.
The data used to support the findings of this study are included within the article.
No studies involving human participants or animals performed by any of authors are described in this article.
Conflicts of Interest
The authors declare that there are no conflicts of interest in this study.
Yu Deng and Guiying Hao designed and performed the experiments; Silu Wang contributed to the clinical diagnosis and collected samples; Yu Deng and Runxia Liu analyzed data and wrote the manuscript. All authors read and approved the final manuscript.
This study was supported by the Applied Foundation Project of the Sichuan Provincial Science and Technology Department (2018JY0243), the Key Foundation Project of the Sichuan Provincial Education Department (18ZD0442), the "Double high talent" Project of Xichang College (LGZ201712), and the postdoctoral scholarship of China Scholarship council (CSC201408515154).
 L. Liu, H. Xia, N. Wahlberg, S. Belak, and C. Baule, "Phylogeny, classification and evolutionary insights into pestiviruses," Virology, vol. 385, no. 2, pp. 351-357, 2009.
 F. V. Bauermann, J. F. Ridpath, R. Weiblen, and E. F. Flores, "HoBi-like viruses: An emerging group of pestiviruses," Journal of Veterinary Diagnostic Investigation, vol. 25, no. 1, pp. 6-15, 2013.
 J. F. Ridpath, "Emerging pestiviruses infecting domestic and wildlife hosts," Animal Health Research Reviews, vol. 16, no. 1, pp. 55-59, 2015.
 D. B. Smith, G. Meyers, J. Bukh et al., "Proposed revision to the taxonomy of the genus Pestivirus, family Flaviviridae," Journal of General Virology, vol. 98, no. 8, pp. 2106-2112, 2017.
 J. F. Ridpath, "Bovine viral diarrhea virus: global status," Veterinary Clinics of North America: Food Animal Practice, vol. 26, no. 1, pp. 105-121, 2010.
 T. Passler, K. P. Riddell, M. A. Edmondson et al., "Experimental infection of pregnant goats with bovine viral diarrhea virus (BVDV) 1 or 2," Veterinary Research, vol. 45, no. 1, p. 38, 2014.
 Y. Abe, T. Tamura, S. Torii et al., "Genetic and antigenic characterization of bovine viral diarrhea viruses isolated from cattle in Hokkaido, Japan," Journal of Veterinary Medical Science, vol. 78, no. 1, pp. 61-70, 2016.
 M. Deng, S. Ji, W. Fei et al., "Prevalence study and genetic typing of Bovine Viral Diarrhea Virus (BVDV) in four bovine species in China," PLoS ONE, vol. 10, no. 4, 2015.
 S. Silveira, M. N. Weber, A. C. S. Masena et al., "Genetic diversity of brazilian bovine pestiviruses detected between 1995 and 2014," Transboundary and Emerging Diseases, vol. 64, no. 2, pp. 613-623, 2017.
 L. Elvira Partida, M. Fern Ndez, J. Guti Rrez et al., "Detection of Bovine Viral Diarrhoea Virus 2 as the Cause of Abortion Outbreaks on Commercial Sheep Flocks," Transboundary and Emerging Diseases, vol. 64, pp. 19-26, 2017.
 M. Jenckel, D. Hoper, H. Schirrmeier et al., "Mixed triple: Allied viruses in unique recent isolates of highly virulent type 2 bovine viral diarrhea virus detected by deep sequencing," Journal of Virology, vol. 88, no. 12, pp. 6983-6992, 2014.
 K. Yesilbag, G. Alpay, and P. Becher, "Variability and global distribution of subgenotypes of bovine viral diarrhea virus," Viruses, vol. 9, p. 128, 2017.
 L. Mao, W. Li, L. Yang et al., "Primary surveys on molecular epidemiology of bovine viral diarrhea virus 1 infecting goats in Jiangsu province, China," BMC Veterinary Research, vol. 12, p. 181, 2016.
 Y. Deng, C.-Q. Sun, S.-J. Cao et al., "High prevalence of bovine viral diarrhea virus 1 in Chinese swine herds," Veterinary Microbiology, vol. 159, no. 3-4, pp. 490-493, 2012.
 M. N. Weber, S. Silveira, G. Machado et al., "High frequency of bovine viral diarrhea virus type 2 in Southern Brazil," Virus Research, vol. 191, no. 1, pp. 117-124, 2014.
 V. Mari, M. Losurdo, M. S. Lucente et al., "Multiplex real-time RT-PCR assay for bovine viral diarrhea virus type 1, type 2 and HoBi-like pestivirus," Journal of Virological Methods, vol. 229, pp. 1-7, 2016.
 H. Albayrak, S. O. Gumusova, E. Ozan, and Z. Yazici, "Molecular detection of pestiviruses in aborted foetuses from provinces in northern Turkey," Tropical Animal Health and Production, vol. 44, no. 4, pp. 677-680, 2012.
 S. Kumar, G. Stecher, and K. Tamura, "MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets," Molecular Biology and Evolution, vol. 33, no. 7, pp. 1870-1874, 2016.
 D. D. Nelson, M. J. Dark, D. S. Bradway et al., "Evidence for persistent Bovine viral diarrhea virus infection in a captive mountain goat (Oreamnos americanus)," Journal of Veterinary Diagnostic Investigation, vol. 20, no. 6, pp. 752-759, 2008.
 F. Xue, Y.-M. Zhu, J. Li et al., "Genotyping of bovine viral diarrhea viruses from cattle in China between 2005 and 2008," Veterinary Microbiology, vol. 143, pp. 379-383, 2010.
 F. Zhong, N. Li, X. Huang et al., "Genetic typing and epidemiologic observation of bovine viral diarrhea virus in Western China," Virus Genes, vol. 42, no. 2, pp. 204-207, 2011.
 S. Gao, J. Luo, J. Du et al., "Serological and molecular evidence for natural infection of Bactrian camels with multiple subgenotypes of bovine viral diarrhea virus in Western China," Veterinary Microbiology, vol. 163, no. 1-2, pp. 172-176, 2013.
 X. Xu, Q. Zhang, X. Yu et al., "Sequencing and comparative analysis of a pig bovine viral diarrhea virus genome," Virus Research, vol. 122, no. 1-2, pp. 164-170, 2006.
Yu Deng [ID], (1) Silu Wang, (1) Runxia Liu, (2) and Guiying Hao (1)
(1) School of Animal Science, Xichang College, Xichang 615013, China
(2) South Dakota State University, Brookings, SD 57007, USA
Correspondence should be addressed to Yu Deng; firstname.lastname@example.org
Received 9 August 2018; Accepted 11 October 2018; Published 18 November 2018
Academic Editor: William Alberto Canon-Franco
Caption: FIGURE 1: Phylogenetic tree of a 224 bp portion of the 5'-UTR nucleotide sequences. The BVDV strains isolated from goats in southwestern China, with Pestivirus reference strains of BVDV-1, BVDV-2, BVDV-3, and BDV registered in GenBank database. The tree was constructed using the Neighbor-Joining method, bootstrapping (1,000 replications), and the Kimura three parameter statistical model of MEGA 7.0 software. The 12 isolated sequences from goat serum in this study are designated with "blue circle".
TABLE 1: Description of the 12 isolates used in this study. Name of the Year Symptoms Region of the positive samples samples collected DF1312 2013 Respiratory DF symptom DC1301 2013 Diarrhea DC DS1324 2013 Weak DS MN1319 2013 Diarrhea MN CF1374 2013 Diarrhea CF HL1407 2014 Abortion HL FS1426 2014 Weak FS XC1522 2015 Mucosal XC disease DB1503 2015 Abortion DB FS1573 2015 Respiratory FS symptom FS1605 2016 Respiratory FS symptom SC1601 2016 Respiratory SC symptom Name of the GenBank subtype positive Accession samples number DF1312 MG323517 1m DC1301 MG323518 1b DS1324 MG323520 1a MN1319 MG323522 1b CF1374 MG323527 1m HL1407 MG323523 1b FS1426 MG323525 1p XC1522 MG323516 1b DB1503 MG323519 1m FS1573 MG323524 1c FS1605 MG323521 1m SC1601 MG323526 1b Name of the Identity (%) positive samples 1a 1b 1c NADL Osloss Mannasi DF1312 DC1301 98.0 DS1324 87.9 MN1319 98.0 CF1374 HL1407 91.6 FS1426 XC1522 97.0 DB1503 FS1573 97.5 FS1605 SC1601 97.0 Name of the Identity (%) positive samples 1m 1p ZM95 TJ06 DF1312 95.4 DC1301 DS1324 MN1319 CF1374 95.7 HL1407 FS1426 98.4 XC1522 DB1503 97.0 FS1573 FS1605 95.7 SC1601
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Deng, Yu; Wang, Silu; Liu, Runxia; Hao, Guiying|
|Publication:||Journal of Veterinary Medicine|
|Date:||Jan 1, 2018|
|Previous Article:||Comparative Seroepidemiological Study of Brucellosis in Sheep under Smallholder Farming and Governmental Breeding Ranches of Central and North East...|
|Next Article:||Phylogenetic Groups and Antimicrobial Susceptibility Patterns of Escherichia coli from Healthy Chicken in Eastern and Central Uganda.|