Skip to main content

Serological evidence of bluetongue virus infection and serotype distribution in dairy cattle in South Korea

Abstract

Background

Bluetongue is a vector-borne viral disease, and bluetongue virus (BTV) outbreaks can cause substantial economic losses. Even subclinical infection may carry significant associated costs, including a loss of condition, reduced milk yield, and infertility and abortion, and indirect costs, largely due to the export restrictions and surveillance requirements imposed to limit the spread of the virus. However, the BTV epidemiology in the Far East remains incompletely understood, especially in the cattle population in South Korea. In this study, the seroprevalence of BTV antibodies and distribution of BTV serotypes in dairy cattle in South Korea were evaluated to improve the understanding of the BTV epidemiological situation in the Asia-Pacific region.

Results

Between 2012 and 2013, a total of 37 out of 171 dairy cattle herds (21.6%) and 85 out of 466 dairy cattle heads (18.2%) showed antibodies against BTV. Neutralizing antibodies to BTV-1, − 2, − 3, − 4, − 7, − 15, and − 16 serotypes were identified, and the RNAs of the BTV-1, − 2, − 3, − 15, and − 16 serotypes were detected, indicating that BTV was circulating in the dairy cattle population in South Korea.

Conclusions

These findings indicate that BTV is widespread and has circulated in dairy cattle in South Korea. This is the first report presenting evidence of circulating antibodies against BTV and the serotype distribution in bovine populations in South Korea.

Background

Bluetongue is a vector-borne viral disease that affects wild and domestic ruminant species and reduces herd productivity [1]. Even subclinical infection may carry significant associated costs, including a loss of condition, reduced milk yield, and infertility and abortion [2], and indirect costs, which largely arise from the export restrictions and surveillance requirements imposed to limit the spread of the bluetongue virus (BTV) [3]. Because it can spread rapidly and cause serious economic consequences in affected countries, bluetongue has been classified as a notifiable disease by the World Organization for Animal Health (OIE) [4].

Bluetongue in sheep is characterized by fever, facial edema, hyperemia, congestion, and erosion of mucous membranes. Bluetongue disease results from vascular injury, likely through a process analogous to that of human hemorrhagic viral fevers, in which the production of vasoactive mediators from virus-infected macrophages and dendritic cells results in enhanced endothelial paracellular permeability with subsequent vascular leakage and hypovolemic shock [5]. BTV infections in ruminants can be subclinical, though fatal diseases predominantly occur in sheep, deer, and wild ruminants [6, 7].

Although bluetongue was recognized and described more than 100 years ago in southern Africa, BTV infection of ruminants, as an enzootic disease, and its vector, Culicoides (Diptera: Ceratopogonidae) genus insects are traditionally found in tropical and temperate regions of the world between the latitudes of approximately 40° North and 35° South. Exceptions include regions of Asia and western North America, where BTV infection of ruminants occurs as far as 50° North [5, 8,9,10]. However, there have been drastic recent regional alterations in the global distribution of BTV infection, particularly in Europe, since 1998 [5]. BTV is considered an emerging disease in Europe, where the disease has spread with unprecedented speed and virulence. There is a substantial body of evidence linking this emergence to climate change [3], and it has been proposed, but certainly not proven, that global climate change is responsible for these events [5].

The epidemiology of BTV infection is poorly defined in much of the world, including extensive portions of Asia. In Asia, bluetongue was first recognized in Israel in 1949 and was later recognized in the People’s Republic of China in 1979 [11], Indonesia in 1981 [12], Japan in 1985 [13], and Malaysia in 1987 [14]. Although the importance of BTV as a transboundary and emerging disease in the world has been suggested, the epidemiology of BTV in South Korea and its neighboring countries in the Asia-Pacific region remains unclear. In the Far East, evidence of BTV has been found in countries neighboring South Korea. It has been accepted that BTV may be endemically established in Japan and China [13, 15,16,17,18,19,20,21,22,23,24,25,26,27,28], and epidemiological studies of the wildlife in Japan have noted the current existence of interconnected domestic and wild cycles that could account for the maintenance of BTV [29]. There was no BTV outbreak or case report in South Korea until 2015, when BTV was isolated from whole-blood samples taken from cattle at abattoirs [30]. Because there was no information on BTV in cattle in South Korea, a study was designed to evaluate the seroprevalence of BTV antibodies and distribution of BTV serotypes in dairy cattle to improve the understanding of the epidemiological situation of BTV in the Asia-Pacific region.

Results

Between 2012 and 2013, at the national level, 37 of 171 dairy cattle herds (21.6, 95% confidence interval: 16.1–28.4%) analyzed and 85 of 466 dairy cattle heads (18.2, 95% confidence interval: 15.0–22.0%) analyzed showed antibodies against BTV using cELISA as shown in Table 1. The agreement of the serologic status between dairy cattle sampled within the same herd in our prevalence study, as measured by the intra-class correlation coefficient, was 0.21. Our results identified the population sizes of herds to be a protective factor. An increased number of animals inside the farm led to a decreased risk of being positive. In the univariable analysis (Table 2), no significant risk differences in land use, adult/calf ratios, and experience of reproductive problems were found. Our results identified herd size as a protective factor. An increased number of animals inside the farm led to risk of being positive (OR = 2.762, p = 0.003 in a herd of ≤5 animals; OR = 4.174, p = < 0.001 in a herd of 6–30), while older age was shown to be a significant risk factor (OR = 0.294, p = < 0.001 in juveniles; OR = 0.443, p = 0.004 in a herd of subadults). The cattle density for the class of 11–20 animals was another significant risk factor (OR = 1.737, p = 0.070), while cattle densities for the classes of ≤10 or ≥ 21 were not significantly associated. We observed a significant difference in the individual likelihood of being positive in southern provinces with respect to northern provinces (OR = 1.757, p = 0.019). A significant difference was also observed in the individual likelihood of being positive in western provinces with respect to eastern provinces (OR = 0.531, p = 0.026). The risk factors identified in the multivariable logistic model (p < 0.001) were as follows: older animals (adults), southern, and northern provinces. The herd size was confirmed to be a protective factor (Table 3). Of the positive dairy cattle herds, 82.4% (28/34) were clustered in either < 30% or > 81%, suggesting a bimodal frequency distribution (Fig. 2). Additionally, there were substantial regional differences in the seroprevalence within South Korea (Fig. 1). Of the 85 ELISA-positive samples, only 59 were positive by SNT (Table 1) and neutralized one or more BTV serotypes: 1 (26 serum samples), 2 (11 serum samples), 3 (15 serum samples), 4 (13 serum samples), 7 (6 serum samples), 15 (9 serum samples), and 16 (13 serum samples). These findings support the notion that BTV-1 is the most prevalent serotype in South Korea. By contrast, 26 dairy cattle serum samples failed to neutralize any known BTV serotype. Additionally, the RNA of the BTV-1, − 2, − 3, − 15, and − 16 serotypes was detected in 13 serologically positive blood samples by RT-PCR, indicating that several BTV serotypes were actively circulating in the dairy cattle populations in the studied area. Further phylogenetic analysis and virus isolation for these blood samples could not be performed using these samples because of an insufficient blood volume and quality of the positive samples.

Table 1 Seroprevalence of BTV infection in dairy cattle in the Republic of Korea (2012–2013)
Table 2 Univariable analysis of the exposure variables relative to the animal seropositivity outcome
Table 3 Multivariable analysis of the exposure variables relative to the animal seropositivity outcome
Fig. 1
figure 1

a Geographical location of the provinces in South Korea. b Density of dairy cattle in South Korea per km2. The herds sampled in this study are mapped by a dot. c Geographical distribution patterns of the seroprevalence to BTV. Percentages of tested animals seropositive for bluetongue virus in 11 provinces in South Korea. The serotypes identified by serum neutralization tests or the BTV serotype identified by RT-PCR of blood sample are shown. Differences in the intensity of the darkness on the map represent different percentages of seroprevalence. The map depicted in figure 1 is our own

Discussion

Our study provides evidence for the serological prevalence of circulating antibodies against BTV and RNA of the BTV-1, − 2, − 3, − 15, and − 16 serotypes in the dairy cattle population in South Korea. The present study demonstrated that BTV infection was prevalent in the dairy cattle populations analyzed, in which approximately one in five dairy cattle herds and one in six dairy cattle were infected (Table 1). Further studies might include virological and serological investigations of BTV-4 and -7 circulating in South Korea because neutralizing antibodies against BTV-4 and -7 have been detected in dairy cattle serum samples. However, none of the blood samples in the present study showed virological evidence of BTV-4 and -7. To understand the results of the present study, it is important to know that South Korea has no vaccination program for bluetongue. Thus, the high seroprevalence of BTV infection in South Korea can be assumed to reflect a natural infection of the dairy cattle evaluated.

Briefly, with respect to Asian countries and the Middle East, the overall seroprevalence obtained was lower than that previously reported in this species in Taiwan, Nepal, India, and Japan [31,32,33,34], but higher than that previously reported in central and south-east Iran [35, 36]. The seroprevalence of BTV infection described in our study (18.2%) is comparable to that described among sentinel ruminants, such as cattle, buffalo, and goats, in China (17.1% seroprevalence) [37] and Indonesian ruminants (2–23% seroprevalence) [38]. The seroprevalence determined in this study is lower than that described among ruminants in Taiwan (32.7% in cattle and 8.2% in goats) [31], sentinel cattle in Japan (5–71%) [13], and sheep and cattle in Israel (16.7% in sheep and 63.2% cattle) [39] but higher than that reported in domestic yaks in China (2–5%) [23].

It was noted that 81.1% (28/34) of the seropositive dairy cattle herds clustered at either < 30% or > 81%, suggesting a “bimodal” distribution (Fig. 2), similar to that described by Taiwanese and Dutch research groups [31, 40]. The total dairy cattle population of South Korea is mostly composed of the Holstein breed. In 2015, the dairy cattle population in South Korea was composed of 402,405 bovines, which were kept on 5,407 holdings (approximately 74 dairy cattle per holding) [41]. The highest density of dairy cattle was in Gyeonggi Province, with 15.7 dairy cattle per km2 (Fig. 2). The lack of clinical signs of BTV might also be because the development of bluetongue in sheep has attracted little attention in South Korea because of the very small number of domestic sheep. Consequently, most BTV episodes throughout the world may be completely silent. Therefore, the results of the seroepidemiological surveillance performed in this study suggest that subclinical or mild BTV infection of ruminants in South Korea is prevalent almost every year and that there are repeated and recurrent infections among domestic cattle and other ruminants in the affected regions.

Fig. 2
figure 2

Frequency distribution of the seroprevalence proportions of bluetongue virus infection in dairy cattle in South Korea (2012–2013). Note the bimodal distribution; one cluster is < 30% and the other cluster > 81% among individual cattle herds

Despite the serological evidence of BTV infection and the serotype distribution in dairy cattle, there is no clinical report of bluetongue in any species in South Korea. BTV is apparently still infecting dairy cattle but causing subclinical disease. The factors related to the lack of a routine monitoring system, difficulties in making a clinical diagnosis, and misdiagnoses of similar viral diseases should also not to be ignored. For example, Aino virus, Akabane virus, Chuzan virus, and Ibaraki virus present the greatest difficulty for differential diagnosis because these viral infections cause similar symptoms and have high seropositive rates in sentinel cattle in South Korea, 33.2% for Aino virus, 40.2% for Akabene virus, 29.1% for Chuzan virus, and 7.5% for Ibaraki virus [42]. In the present study, an intra-class correlation coefficient of 0.21 was found, indicating that the correlation between two animals within a herd with respect to the BTV result was relatively high. The intra-class correlation coefficient of livestock infectious disease is usually < 0.2 and ranges from 0.04 to 0.42 [43]. This finding indicates that in any particular dairy cattle herd, it is likely that either most dairy cattle in that herd will be serologically positive or most will be negative. Our results revealed a significant association between age and BTV seroprevalence in South Korea, and a similar association was previously reported by several groups studying other cattle populations [44,45,46]. This association between age class and BTV seropositivity can be explained by a longer time of exposure for adults [46, 47]. The proportion of seropositive animals increased with age, probably resulting more from prolonged exposure of the adults to the vector than to any resistant status of juveniles. BTV infection significantly decreased with the increase in the number of dairy cattle inside a farm. The virus pressure in the vectors and ruminants would decrease due to the dilution effect caused by the increase in the density of susceptible dairy cattle on the farm. Nevertheless, this interesting finding requires further investigation. The results obtained from this study indicate that northern and eastern regions are areas of relatively low BTV seroprevalence. Relative to the factors that affect the likelihood of infection at an individual animal level, living in southern and western regions in South Korea increased the risk of being BTV positive compared with living in northern and eastern areas. This result could partially be explained by the spatial distribution of the Culicoides group, but the association between the BTV epidemiology and Culicoides vector distribution in South Korea is not fully understood.

Conclusions

Collectively, in this study, the significance of the serological evidence and the findings of BTV serotypes distributed in dairy cattle South Korea are discussed and further studies to isolate local BTV strains from sentinel animals, including dairy cattle, are proposed. Entomological surveillance of biting Culicoides midges has been consistently conducted in South Korea since 2010 [48,49,50,51,52], and an epidemiological analysis is being assessed based on the vector species abundance, diversity, and competence of Culicoides midges distributed in South Korea. In addition, future surveillance programs for BTV should be extended to include other susceptible animals, such as sheep, cattle, and other ruminants. To our knowledge, this is the first report of evidence of circulating antibodies against BTV and the serotype distribution in the bovine population in South Korea. These findings indicate that BTV is widespread and actively circulating in dairy cattle in South Korea. The results also support the idea that serosurveillance of this species would be useful for detecting BTV circulation [34, 53,54,55,56]. Thorough research on the hosts and vectors involved in BTV circulation in natural ecosystems in the Far East is also needed, as has previously been suggested elsewhere [30].

Methods

Field samples for survey of BTV serology and the serotype distribution

A serological survey was conducted between 2012 and 2013 on dairy cattle for specific antibodies to BTV. Dairy cattle were selected as the target population because seroconversions usually occur earlier in bovine species (they are the preferred host [57]) and BTV antibody prevalence is considered to be higher in cattle (49–60%) than in sheep or goats (13.7–20%) [44, 58,59,60,61,62]. For this study, serum from as many animals as possible was collected from each herd based on the estimated seroprevalence of 10% and a confidence limit of 95% [63]. Herds and animals within the herds were selected by a simple random sampling method in each province based on the government’s national statistics [64]. The sampling frame was established using the dairy cattle farm ID and herd size obtained from the KAHIS (Korea Animal Health Integrated System, Animal and Plant quarantine agency, Anyang, South Korea). Animals younger than 6 months were excluded from the sampling frame to avoid the detection of antibodies due to maternal immunity. The apparent prevalence rates were considered to be the animal-level prevalence, defined as the proportion of cELISA-positive animals out of the total number of animals tested in the study area, and the herd prevalence, which was the proportion of cELISA-positive herds out of the total number of tested herds in the area. A herd was classified as positive if at least one animal was found to be cELISA positive. The intra-class correlation coefficient μ was calculated to measure the serologic status agreement between dairy cattle sampled within the same herd. Samples from dairy cattle were obtained from the blood and serum bank of the National Surveillance Program maintained by the Foreign Animal Diseases Division of the National Veterinary Research and Quarantine Service (Anyang, South Korea) in close collaboration with local veterinary practitioners and/or government veterinary officers. The number of samples from each province is shown in Table 1. The serum separated from the blood samples was stored at − 20 °C until further analysis. The seroprevalence rates were estimated at the herd and animal levels in dairy cattle, and the study was conducted in South Korea (33°06′ N - 39° 25′ N, 124°36′ E - 131°52′ W) from 2012 through 2013.

Determination of BTV antibodies

It was previously reported that specific serodiagnostic techniques, such as competitive enzyme-linked immunosorbent assay (cELISA) or BTV neutralization tests, should be used for bluetongue surveillance in Ibaraki virus (IBV)-endemic areas because IBV-positive serum samples may result in false-positive bluetongue agar gel immunodiffusion (AGID) test reactions [65]. Because South Korea is an IBV-endemic region, the BTV VP7 Antibody Test Kit (IDEXX cELISA; IDEXX Laboratories, Inc., Institute Pourquier, Montpellier, France) was used to identify the presence of antibodies against BTV in the serum samples. BTV VP7 is the major structural protein of the inner core and contains group-specific antigenic determinants [66]. Therefore, most diagnostic methods for BTV-specific antibody detection have been based on BTV VP7 [67,68,69,70,71,72,73,74,75]. The assays were performed according to the manufacturer’s instructions. The specificity and sensitivity of the IDEXX cELISA kit are 100 and 82.8%, as stated by the manufacturer, respectively [76]. A herd was defined as “positive” when at least one seropositive sample was present. The seroprevalence rates and 95% confidence intervals were calculated using the program ‘Survey Toolbox for Livestock Diseases’ (Ausvet, Australia). The serotype-specific serum neutralization test (SNT) was used to confirm the positive ELISA results. Positive and negative controls for the SNT were obtained from the Institute for Animal Health, Pirbright, United Kingdom. The microtiter neutralization method was used in this study according to the Manual of Diagnostic Tests and Vaccines for Terrestrial Animals [77]. The BTV strains were used in this study as follows: RSArrr/01 (serotype 1), RSArrr/02 (serotype 2), RSArrr/03 (serotype 3), RSArrr/04 (serotype 4), RSArrr/05 (serotype 5), RSArrr/06 (serotype 6), RSArrr/07 (serotype 7), RSArrr/08 (serotype 8), RSArrr/09 (serotype 9), RSArrr/10 (serotype 10), RSArrr/11 (serotype 11), RSArrr/12 (serotype 12), RSArrr/13 (serotype 13), RSArrr/14 (serotype 14), RSArrr/15 (serotype 15), RSArrr/16 (serotype 16), RSArrr/17 (serotype 17), RSArrr/18 (serotype 18), RSArrr/19 (serotype 19), RSArrr/20 (serotype 20), RSArrr/21 (serotype 21), RSArrr/22 (serotype 22), RSArrr/23 (serotype 23), and RSArrr/24 (serotype 24). Briefly, approximately 100 TCID50 (50% tissue culture infective dose) of the standard or untyped virus was added to a volume of 50 μl to test the wells of a flat-bottomed microtiter plate and was mixed with an equal volume of standard antiserum that had been serially diluted in tissue culture medium. Approximately 104 Vero cells (CCL-81, American Type Culture Collection, Manassas, VA, USA) were added per well in a volume of 100 μl and were assessed after incubation for 4–6 days using an inverted microscope (Olympus, Tokyo, Japan). The wells were scored for the degree of cytopathic effects observed.

RT-PCR for BTV RNA detection

A total of 466 blood samples was tested for the presence of BTV RNA using non-serotype-specific reverse transcription PCR targeting a conserved region within the VP3 gene of the BTV genome according to the method described by Yeh et al. [78]. It is known that RNA segment 3 of BTV encodes a serogroup-reactive protein, VP3, and that VP3 sequences are often used for the genetic characterization of BTV serotypes worldwide [79,80,81]. Briefly, total nucleic acids were extracted from 400 μl of whole blood. Automated extraction was performed using a BioRobot M48 workstation apparatus (Qiagen, GmbH, Hilden, Germany) with a MagAttract Virus Mini M48 kit (Qiagen). Nucleic acids were recovered in 50 μl of elution buffer. RT-PCR was performed using a one-step RT-PCR kit (Qiagen). The reactions were prepared in a volume of 25 μl containing 2 μl of RNA, 1× buffer [Tris-Cl, KCl, (NH4)2SO4], 2.5 mM MgCl2, 0.2 mM deoxynucleoside triphosphates, 0.4 μM aliquots of each of the specific primers, 5 U of RNase inhibitor (Intron Biotechnology, Korea), and 1 μl of enzyme mix (Omniscript and Sensiscript RTs, HotStartTaq DNA polymerase; Qiagen). Reverse transcription amplification was accomplished in one step using the following optimized incubation program: 30 min at 50 °C; 15 min at 95 °C; 40 cycles of 94 °C for 30 s, 61 °C for 90 s, and 72 °C for 50 s; and 1 min at 72 °C. RT-PCR amplifications were performed using an Eppendorf MasterCycler gradient thermal cycler (Eppendorf, Germany). The RT-PCR amplification products (5 μl) were analyzed by gel electrophoresis on a 3% agarose gel containing 0.5 μg of ethidium bromide/ml. For the identification and differentiation of bluetongue virus serotypes, previously published serotype-specific RT-PCR [82] was performed using BTV RNA-positive samples identified by RT-PCR to detect pan-BTV. The positive control for the BTV serotype RNAs was prepared according to a method previously described by Yeh et al. [78].

Statistical analysis and risk factor analysis

The true prevalence (TP) was calculated from the apparent prevalence (AP) using the Rogan and Gladen equation. The formula for the calculation is as follows: [83]

$$ \mathrm{TP}=\left(\mathrm{AP}+\mathrm{Sp}-1\right)/\left(\mathrm{Se}+\mathrm{Sp}-1\right). $$

The TP and Blaker’s 95% confidence intervals (CI) [84] were calculated using the Epitools epidemiological calculator [85]. Clustering of disease was explored using the analysis of variance (ANOVA) estimator of the intra-class correlation. A measure for agreement in serologic status between animals within a herd is given by the intra-class correlation coefficient μ. The intra-class correlation coefficient (minimum 0, maximum 1) was estimated using analysis of variance, with herd as the independent variable and the serologic status of individual animals (seropositive or seronegative) as the dependent variable [86]. In this study, univariable and multivariable analyses were performed. The following individual exposure variables were considered for the univariable and multivariable analyses: land use (agricultural, woodland and semi-natural and urban areas, according to KAHIS), the population sizes of the herds, and the cattle population size inside a 1 km-buffer around the sampling farm (hereafter called cattle density). The radius size was chosen considering the most likely Culicoides flying range < 1 km [87]. The herd size and number of animals were divided into classes following the tertile classification method: herd size (≤5, 6–30, ≥31), cattle density (≤10, 11–20, ≥21), and adult/calf ratios ≤1.0, 1.1–2.0, ≥2.1). The animals were also classified into three age groups based on tooth replacement and livestock owner questionnaires: juveniles (between 6 months and 1 year old), sub-adults (between 1 and 2 years old), and adults (> 2 years old). Regional risk factors, such as geographic location and local seropositivity risk factors based on livestock owner questionnaires—e.g., reproductive problems, including abortion—were also investigated. A logistic regression model was used to check the association of the animal seropositivity outcome with potential risk factors. The effect of the exposure variables on individual seropositivity was analyzed using univariable logistic regression models, and the variables in the univariable analysis were screened for pair-wise collinearity or associations using Pearson’s correlation coefficient or the chi-squared test for continuous or categorized variables, respectively. The strength of association was calculated using odds ratios at 95% CI. A p value < 0.05 was considered to be statistically significant. All statistical analyses were performed using the statistical software SPSS Statistics version 25 (IBM Corp., Armonk, NY, USA).

Availability of data and materials

All data generated or analyzed during the study are included in this published article. The datasets used and/or analyzed in the current study are available from the corresponding author on reasonable request. Individual farm information may be protected or not provided by privacy and security concerns.

Abbreviations

BTV:

Bluetongue virus

cELISA:

Competitive enzyme-linked immunosorbent assay

IBV:

Ibaraki virus

OIE:

World Organization for Animal Health

References

  1. Ruiz-Fons F, Sanchez-Matamoros A, Gortazar C, Sanchez-Vizcaino JM. The role of wildlife in bluetongue virus maintenance in Europe: lessons learned after the natural infection in Spain. Virus Res. 2014;182:50–8.

    Article  CAS  PubMed  Google Scholar 

  2. Osburn BI. The impact of bluetongue virus on reproduction. Comp Immunol Microbiol Infect Dis. 1994;17(3–4):189–96.

    Article  CAS  PubMed  Google Scholar 

  3. Wilson A, Mellor P. Bluetongue in Europe: vectors, epidemiology and climate change. Parasitol Res. 2008;103(Suppl 1):S69–77.

    Article  PubMed  Google Scholar 

  4. Pascual-Linaza AV, Martinez-Lopez B, Pfeiffer DU, Moreno JC, Sanz C, Sanchez-Vizcaino JM. Evaluation of the spatial and temporal distribution of and risk factors for bluetongue serotype 1 epidemics in sheep Extremadura (Spain), 2007-2011. Prev Vet Med. 2014;116(3):279–95.

    Article  PubMed  Google Scholar 

  5. Maclachlan NJ. Bluetongue: history, global epidemiology, and pathogenesis. Prev Vet Med. 2011;102(2):107–11.

    Article  PubMed  Google Scholar 

  6. Maclachlan NJ, Drew CP, Darpel KE, Worwa G. The pathology and pathogenesis of bluetongue. J Comp Pathol. 2009;141(1):1–16.

    Article  CAS  PubMed  Google Scholar 

  7. Maclachlan NJ, Mayo CE, Daniels PW, Savini G, Zientara S, Gibbs EP. Bluetongue. Rev Sci Tech. 2015;34(2):329–40.

    Article  CAS  PubMed  Google Scholar 

  8. Clavijo A, Heckert RA, Dulac GC, Afshar A. Isolation and identification of bluetongue virus. J Virol Methods. 2000;87(1–2):13–23.

    Article  CAS  PubMed  Google Scholar 

  9. Lundervold M, Milner-Gulland EJ, O'Callaghan CJ, Hamblin C. First evidence of bluetongue virus in Kazakhstan. Vet Microbiol. 2003;92(3):281–7.

    Article  CAS  PubMed  Google Scholar 

  10. Tabachnick WJ. Culicoides and the global epidemiology of bluetongue virus infection. Vet Ital. 2004;40(3):144–50.

    CAS  PubMed  Google Scholar 

  11. Zhang N, MacLachlan NJ, Bonneau KR, Zhu J, Li Z, Zhang K, et al. Identification of seven serotypes of bluetongue virus from the People's Republic of China. Vet Rec. 1999;145(15):427–9.

    Article  CAS  PubMed  Google Scholar 

  12. Daniels PW, Sendow I, Soleha E, Sukarsih, Hunt NT, Bahri S. Australian-Indonesian collaboration in veterinary arbovirology--a review. Vet Microbiol. 1995;46(1–3):151–74.

    Article  CAS  PubMed  Google Scholar 

  13. Goto Y, Yamaguchi O, Kubo M. Epidemiological observations on bluetongue in sheep and cattle in Japan. Vet Ital. 2004;40(3):78–82.

    CAS  PubMed  Google Scholar 

  14. Sharifah SH, Ali MA, Gard GP, Polkinghorne IG. Isolation of multiple serotypes of bluetongue virus from sentinel livestock in Malaysia. Trop Anim Health Prod. 1995;27(1):37–42.

    Article  CAS  PubMed  Google Scholar 

  15. Kato T, Shirafuji H, Tanaka S, Sato M, Yamakawa M, Tsuda T, et al. Bovine arboviruses in Culicoides biting midges and sentinel cattle in southern Japan from 2003 to 2013. Transbound Emerg Dis. 2016;63(6):e160–72.

    Article  CAS  PubMed  Google Scholar 

  16. Shirafuji H, Yanase T, Kato T, Yamakawa M. Genetic and phylogenetic characterization of genome segments 2 and 6 of bluetongue virus isolates in Japan from 1985 to 2008. J Gen Virol. 2012;93(Pt 7):1465–73.

    Article  CAS  PubMed  Google Scholar 

  17. Miura Y, Inaba Y, Tsuda T, Tokuhisa S, Sato K, Akashi H. Seroepizootiological survey on bluetongue virus infection in cattle in Japan. Natl Inst Anim Health Q (Tokyo). 1982;22(4):154–8.

    CAS  Google Scholar 

  18. Ma JG, Zhang XX, Zheng WB, Xu YT, Zhu XQ, Hu GX, et al. Seroprevalence and risk factors of bluetongue virus infection in Tibetan sheep and yaks in Tibetan plateau, China. Biomed Res Int. 2017;2017:5139703.

    PubMed  PubMed Central  Google Scholar 

  19. Gao X, Qin H, Xiao J, Wang H. Meteorological conditions and land cover as predictors for the prevalence of bluetongue virus in the Inner Mongolia Autonomous Region of mainland China. Prev Vet Med. 2017;138:88–93.

    Article  PubMed  Google Scholar 

  20. Yang H, Xiao L, Wang J, Meng J, Lv M, Liao D, et al. Phylogenetic Characterization Genome Segment 2 of Bluetongue Virus Strains Belonging to Serotypes 5, 7 and 24 Isolated for the First Time in China During 2012 to 2014. Transbound Emerg Dis. 2017;64(4):1317–21.

    Article  PubMed  Google Scholar 

  21. Yang H, Lv M, Sun M, Lin L, Kou M, Gao L, et al. Complete genome sequence of the first bluetongue virus serotype 7 isolate from China: evidence for entry of African-lineage strains and reassortment between the introduced and native strains. Arch Virol. 2016;161(1):223–7.

    Article  CAS  PubMed  Google Scholar 

  22. Sun EC, Huang LP, Xu QY, Wang HX, Xue XM, Lu P, et al. Emergence of a novel bluetongue virus serotype, China 2014. Transbound Emerg Dis. 2016;63(6):585–9.

    Article  CAS  PubMed  Google Scholar 

  23. Li J, Li K, Shahzad M, Han Z, Nabi F, Gao J, et al. Seroprevalence of bluetongue virus in domestic yaks (Bos grunniens) in Tibetan regions of China based on circulating antibodies. Trop Anim Health Prod. 2015;47(6):1221–3.

    Article  PubMed  Google Scholar 

  24. Zhu J, Yang H, Li H, Xiao L, Wang J, Li N, Zhang N. 2013. Full-genome sequence of bluetongue virus serotype 1 (BTV-1) strain Y863, the first BTV-1 isolate of Eastern origin found in China. Genome Announc. 1(4):e00403-13. https://0-doi-org.brum.beds.ac.uk/10.1128/genomeA.00403-13.

  25. Yang T, Liu N, Xu Q, Sun E, Qin Y, Zhao J, et al. Complete genomic sequence of bluetongue virus serotype 1 from China. J Virol. 2012;86(2):1288–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Yang H, Zhu J, Li H, Xiao L, Wang J, Li N, et al. Full genome sequence of bluetongue virus serotype 4 from China. J Virol. 2012;86(23):13122–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Yang T, Liu N, Xu Q, Sun E, Qin Y, Zhao J, et al. Complete genomic sequence of bluetongue virus serotype 16 from China. J Virol. 2011;85(24):13472.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Zhang N, Li Z, Zhang F, Zhu J. Studies on bluetongue disease in the People’s Republic of China. Vet Ital. 2004;40(3):51–6.

    PubMed  Google Scholar 

  29. Sugiyama I, Shimizu E, Nogami S, Suzuki K, Miura Y, Sentsui H. Serological survey of arthropod-borne viruses among wild boars in Japan. J Vet Med Sci. 2009;71(8):1059–61.

    Article  PubMed  Google Scholar 

  30. Seo HJ, Park JY, Cho YS, Cho IS, Yeh JY. First report of bluetongue virus isolation in the Republic of Korea and analysis of the complete coding sequence of the segment 2 gene. Virus Genes. 2015;50(1):156–9.

    Article  CAS  PubMed  Google Scholar 

  31. Lee F, Ting LJ, Jong MH, Chang WM, Wang FI. Subclinical bluetongue virus infection in domestic ruminants in Taiwan. Vet Microbiol. 2010;142(3–4):225–31.

    Article  PubMed  Google Scholar 

  32. Khanal DR, Prajapati M, Shrestha P, Acharya MP, Paudyal N, Bowen R, et al. Detection of antibodies against bluetongue virus among domestic ruminants in the highlands of Nepal. Vet Ital. 2016;52(3–4):313–7.

    PubMed  Google Scholar 

  33. Ayanur A, Singh KP, Cherian S, Singh V, Nandi S, Saminathan M, et al. Sero-epidemiology and molecular detection of bluetongue virus in Indian ruminants. Vet Ital. 2016;52(3–4):305–11.

    PubMed  Google Scholar 

  34. Gaire TN, Karki S, Dhakal IP, Khanal DR, Bowen RA. Serosurveillance and factors associated with the presence of antibodies against bluetongue virus in dairy cattle in two eco-zones of Nepal. Rev Sci Tech. 2016;35(3):779–85.

    Article  CAS  PubMed  Google Scholar 

  35. Noaman V, Shirvani E, Hosseini SM, Shahmoradied AH, Heidari MR, Raiszadeh H, et al. Serological surveillance of bluetongue virus in cattle in Central Iran. Vet Ital. 2013;49(2):141–4.

    PubMed  Google Scholar 

  36. Mozaffari AA, Khalili M, Yahyazadeh F. A serological investigation of bluetongue virus in cattle of south-East Iran. Vet Ital. 2012;48(1):41–4.

    PubMed  Google Scholar 

  37. Kirkland PD, Zhang N, Hawkes RA, Li Z, Zhang F, Davis RJ, et al. Studies on the epidemiology of bluetongue virus in China. Epidemiol Infect. 2002;128(2):257–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Sendow I, Daniels PW, Cybinski DH, Young PL, Ronohardjo P. Antibodies against certain bluetongue and epizootic haemorrhagic disease viral serotypes in Indonesian ruminants. Vet Microbiol. 1991;28(1):111–8.

    Article  CAS  PubMed  Google Scholar 

  39. Kedmi M, Levi S, Galon N, Bomborov V, Yadin H, Batten C, et al. No evidence for involvement of sheep in the epidemiology of cattle virulent epizootic hemorrhagic disease virus. Vet Microbiol. 2011;148(2–4):408–12.

    Article  CAS  PubMed  Google Scholar 

  40. Elbers AR, Popma J, Oosterwolde S, van Rijn PA, Vellema P, van Rooij EM. A cross-sectional study to determine the seroprevalence of bluetongue virus serotype 8 in sheep and goats in 2006 and 2007 in the Netherlands. BMC Vet Res. 2008;4:33.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Ministry of Agriculture FaRA. Agriculture, forestry and livestock products : Key statistics 2015; 2016.

    Google Scholar 

  42. Kim YH, Oem JK, Lee EY, Lee KK, Kim SH, Lee MH, et al. Seroprevalence of five arboviruses in sentinel cattle as part of nationwide surveillance in South Korea, 2009-2012. J Vet Med Sci. 2015;77(2):247–50.

    Article  PubMed  Google Scholar 

  43. Otte MJ, Gumm ID. Intra-cluster correlation coefficients of 20 infections calculated from the results of cluster-sample surveys. Prev Vet Med. 1997;31(1–2):147–50.

    Article  CAS  PubMed  Google Scholar 

  44. Ward MP, Carpenter TE, Osburn BI. Host factors affecting seroprevalence of bluetongue virus infections of cattle. Am J Vet Res. 1994;55(7):916–20.

    CAS  PubMed  Google Scholar 

  45. Linden A, Gregoire F, Nahayo A, Hanrez D, Mousset B, Massart AL, et al. Bluetongue virus in wild deer, Belgium, 2005-2008. Emerg Infect Dis. 2010;16(5):833–6.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Garcia-Bocanegra I, Arenas-Montes A, Lorca-Oro C, Pujols J, Gonzalez MA, Napp S, et al. Role of wild ruminants in the epidemiology of bluetongue virus serotypes 1, 4 and 8 in Spain. Vet Res. 2011;42:88.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Falconi C, Lopez-Olvera JR, Gortazar C. BTV infection in wild ruminants, with emphasis on red deer: a review. Vet Microbiol. 2011;151(3–4):209–19.

    Article  PubMed  Google Scholar 

  48. Kim HC, Bellis GA, Kim MS, Chong ST, Lee DK, Park JY, et al. Seasonal abundance of biting midges, Culicoides spp. (Diptera: Ceratopogonidae), collected at cowsheds in the southern part of the Republic of Korea. Korean J Parasitol. 2012;50(2):127–31.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Bellis G, Kim HC, Kim MS, Klein TA, Lee DK, Gopurenko D. Three species of Culicoides Latreille (Diptera: Ceratopogonidae) newly recorded from the Republic of Korea. Zootaxa. 2013;3718:171–82.

    Article  PubMed  Google Scholar 

  50. Oem JK, Chung JY, Kwon MS, Kim TK, Lee TU, Bae YC. Abundance of biting midge species (Diptera: Ceratopogonidae, Culicoides spp.) on cattle farms in Korea. J Vet Sci. 2013;14(1):91–4.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Kim HC, Bellis GA, Kim MS, Klein TA, Chong ST, Park JY. Seasonal abundance of Culicoides (Diptera: Ceratopogonidae) collected by mosquito magnet(R) in northern Gyeonggi-do (province), Korea. Korean J Parasitol. 2014;52(1):57–62.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Kim HC, Bellis GA, Kim MS, Klein TA, Gopurenko D, Cai DC, et al. Species diversity and seasonal distribution of Culicoides spp. (Diptera: Ceratopogonidae) in Jeju-do, Republic of Korea. Korean J Parasitol. 2015;53(4):501–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Nusinovici S, Madouasse A, Fourichon C. Quantification of the increase in the frequency of early calving associated with late exposure to bluetongue virus serotype 8 in dairy cows: implications for syndromic surveillance. Vet Res. 2016;47:18.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Mayo CE, Mullens BA, Reisen WK, Osborne CJ, Gibbs EP, Gardner IA, et al. Seasonal and interseasonal dynamics of bluetongue virus infection of dairy cattle and Culicoides sonorensis midges in northern California--implications for virus overwintering in temperate zones. PLoS One. 2014;9(9):e106975.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Mayo CE, Mullens BA, Gerry AC, Barker CM, Mertens PP, Maan S, et al. The combination of abundance and infection rates of Culicoides sonorensis estimates risk of subsequent bluetongue virus infection of sentinel cattle on California dairy farms. Vet Parasitol. 2012;187(1–2):295–301.

    Article  PubMed  Google Scholar 

  56. Mayo CE, Gardner IA, Mullens BA, Barker CM, Gerry AC, Guthrie AJ, et al. Anthropogenic and meteorological factors influence vector abundance and prevalence of bluetongue virus infection of dairy cattle in California. Vet Microbiol. 2012;155(2–4):158–64.

    Article  PubMed  Google Scholar 

  57. Afshar A, Thomas FC, Wright PF, Shapiro JL, Anderson J. Comparison of competitive ELISA, indirect ELISA and standard AGID tests for detecting blue-tongue virus antibodies in cattle and sheep. Vet Rec. 1989;124(6):136–41.

    Article  CAS  PubMed  Google Scholar 

  58. Behymer DE, Riemann HP, Utterback W, D-Elmi C, Franti CE. Mass screening of cattle sera against 14 infectious disease agents, using an ELISA system for monitoring health in livestock. Am J Vet Res. 1991;52(10):1699–705.

    CAS  PubMed  Google Scholar 

  59. Bonneau KR, DeMaula CD, Mullens BA, MacLachlan NJ. Duration of viraemia infectious to Culicoides sonorensis in bluetongue virus-infected cattle and sheep. Vet Microbiol. 2002;88(2):115–25.

    Article  CAS  PubMed  Google Scholar 

  60. Lelli R, Di Ventura M, Mercante MT, Tittarelli M, Mangana-Vougiouka O, Nomikou K, et al. Bluetongue laboratory diagnosis: a ring test to evaluate serological results using a competitive ELISA kit. Vet Ital. 2004;40(4):577–80.

    CAS  PubMed  Google Scholar 

  61. Nevill EM. The use of cattle to protect sheep from bluetongue infection. J S Afr Vet Assoc. 1978;49(2):129–30.

    CAS  PubMed  Google Scholar 

  62. Singer RS, Boyce WM, Gardner IA, Johnson WO, Fisher AS. Evaluation of bluetongue virus diagnostic tests in free-ranging bighorn sheep. Prev Vet Med. 1998;35(4):265–82.

    Article  CAS  PubMed  Google Scholar 

  63. Thrusfield M. Veterinary Epidemiology. Oxford: Blackwell Science Publications; 1995.

    Google Scholar 

  64. Livestock_Management_Division. Domestic animals statistics. In: Animal Farm Policy Bureau, Ministry of Agriculture, Food and Rural Affairs; 2012.

    Google Scholar 

  65. Shimizu S, Toyota I, Arishima T, Goto Y. Frequency of serological cross-reactions between Ibaraki and bluetongue viruses using the agar gel immunodiffusion test. Vet Ital. 2004;40(4):583–6.

    CAS  PubMed  Google Scholar 

  66. Mertens PP, Brown F, Sangar DV. Assignment of the genome segments of bluetongue virus type 1 to the proteins which they encode. Virology. 1984;135(1):207–17.

    Article  CAS  PubMed  Google Scholar 

  67. Chand K, Biswas SK, Pandey AB, Saxena A, Tewari N, Mondal B. A competitive ELISA for detection of group specific antibody to bluetongue virus using anti-core antibody. Biologicals. 2017;46:168–71.

    Article  CAS  PubMed  Google Scholar 

  68. Venkatesan G, Biswas SK, Bhanuprakash V, Singh RK, Mondal B. Evaluation of thermo-stability of bluetongue virus recombinant VP7 antigen in indirect ELISA. Virusdisease. 2015;26(1–2):19–26.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Xu QY, Sun EC, Yang T, Feng YF, Geng HW, Wei T, et al. Isolation of a bluetongue virus group-specific monoclonal antibody and application to a diagnostic competitive ELISA. Appl Microbiol Biotechnol. 2015;99(2):729–39.

    Article  CAS  PubMed  Google Scholar 

  70. Yin HQ, Jia MX, Shi LJ, Yang S, Zhang LY, Zhang QM, et al. Nanoparticle-based bio-barcode assay for the detection of bluetongue virus. J Virol Methods. 2011;178(1–2):225–8.

    Article  CAS  PubMed  Google Scholar 

  71. Chand K, Biswas SK, Sing B, De A, Mondal B. A sandwich ELISA for the detection of bluetongue virus in cell culture using antiserum against the recombinant VP7 protein. Vet Ital. 2009;45(3):443–8.

    PubMed  Google Scholar 

  72. Yang J, Hua Q, Chen H, Lv J, Qin Z, Jin M, et al. Development and evaluation of an immunochromatographic strip for the detection of serum antibodies against bluetongue virus. J Virol Methods. 2010;163(1):68–73.

    Article  CAS  PubMed  Google Scholar 

  73. Mecham JO, Wilson WC. Antigen capture competitive enzyme-linked immunosorbent assays using baculovirus-expressed antigens for diagnosis of bluetongue virus and epizootic hemorrhagic disease virus. J Clin Microbiol. 2004;42(2):518–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Zhou EM, Ridd D, Riva J, Fernando L, Clavijo A. Development and evaluation of an IgM-capture ELISA for detection of recent infection with bluetongue viruses in cattle. J Virol Methods. 2001;91(2):175–82.

    Article  CAS  PubMed  Google Scholar 

  75. Cloete M, du Plessis DH, van Dijk AA, Huismans H, Viljoen GJ. Vaccinia virus expression of the VP7 protein of south African bluetongue virus serotype 4 and its use as an antigen in a capture ELISA. Arch Virol. 1994;135(3–4):405–18.

    Article  CAS  PubMed  Google Scholar 

  76. Niedbalski W. Evaluation of commercial ELISA kits for the detection of antibodies against bluetongue virus. Pol J Vet Sci. 2011;14(4):615–9.

    Article  CAS  PubMed  Google Scholar 

  77. Baetza HJ. Eradication of bluetongue disease in Germany by vaccination. Vet Immunol Immunopathol. 2014;158(1–2):116–9.

    Article  PubMed  Google Scholar 

  78. Yeh JY, Lee JH, Seo HJ, Park JY, Moon JS, Cho IS, et al. Simultaneous detection of Rift Valley fever, bluetongue, rinderpest, and Peste des petits ruminants viruses by a single-tube multiplex reverse transcriptase-PCR assay using a dual-priming oligonucleotide system. J Clin Microbiol. 2011;49(4):1389–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. McColl KA, Gould AR. Detection and characterisation of bluetongue virus using the polymerase chain reaction. Virus Res. 1991;21(1):19–34.

    Article  CAS  PubMed  Google Scholar 

  80. Huismans H, Erasmus BJ. Identification of the serotype-specific and group-specific antigens of bluetongue virus. Onderstepoort J Vet Res. 1981;48(2):51–8.

    CAS  PubMed  Google Scholar 

  81. Gould AR, Pritchard LI. Phylogenetic analyses of the complete nucleotide sequence of the capsid protein (VP3) of Australian epizootic haemorrhagic disease of deer virus (serotype 2) and cognate genes from other orbiviruses. Virus Res. 1991;21(1):1–18.

    Article  CAS  PubMed  Google Scholar 

  82. Maan NS, Maan S, Belaganahalli MN, Ostlund EN, Johnson DJ, Nomikou K, et al. Identification and differentiation of the twenty six bluetongue virus serotypes by RT-PCR amplification of the serotype-specific genome segment 2. PLoS One. 2012;7(2):e32601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Rogan WJ, Gladen B. Estimating prevalence from the results of a screening test. Am J Epidemiol. 1978;107(1):71–6.

    Article  CAS  PubMed  Google Scholar 

  84. Blaker H. Confidence curves and improved exact confidence intervals for discrete distributions. Can J Stat. 2000;28:783–98.

    Article  Google Scholar 

  85. Sergeant ES. Epitools epidemiological calculators. Canberra: Ausvet Pty Ltd; 2016.

  86. Donald A, Donner A. Adjustments to the mantel-Haenszel chi-square statistic and odds ratio variance estimator when the data are clustered. Stat Med. 1987;6(4):491–9.

    Article  CAS  PubMed  Google Scholar 

  87. Sanders CJ, Harrup LE, Tugwell LA, Brugman VA, England M, Carpenter S. Quantification of within- and between-farm dispersal of Culicoides biting midges using an immunomarking technique. J Appl Ecol. 2017;54:1429–39.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

The authors thank all of the farmers who volunteered to participate in the study. The authors would also like to acknowledge the individuals who helped to recruit farms to the study and the farmers who allowed us to collect serum samples at their farms. The authors are very grateful to the veterinary officers and veterinary practitioners of the local government for their help in collecting blood samples from dairy cattle throughout the country. Bluetongue-specific antibody testing, SNT, and RT-PCR were mainly performed at the National Veterinary Research and Quarantine Service (NVRQS, Anyang, South Korea). Additional laboratory experiments, the literature review, data analysis, statistical analysis, and the writing of the paper were conducted at Incheon National University (Incheon, South Korea).

Funding

This work was supported by an Incheon National University Research Grant and Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT and Future Planning (NRF-2017R1C1B5017328).

Author information

Authors and Affiliations

Authors

Contributions

Conceived and designed the study: JYY. Laboratory testing: JYY and JMH. Drafted the preliminary and revised manuscript: JYY and JGK. The authors read and approved the manuscript.

Corresponding author

Correspondence to Jung-Yong Yeh.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Non-applicable.

Competing interests

None of the authors of this paper had a financial or personal relationship with other people or organizations that could inappropriately influence or bias the content of the paper. The opinions expressed by authors contributing to this journal do not necessarily reflect the opinions of the institutions with which the authors are affiliated.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hwang, JM., Kim, J.G. & Yeh, JY. Serological evidence of bluetongue virus infection and serotype distribution in dairy cattle in South Korea. BMC Vet Res 15, 255 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s12917-019-2000-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12917-019-2000-z

Keywords