Next Article in Journal
Patients with Hematological Malignancies Treated with T-Cell or B-Cell Immunotherapy Remain at High Risk of Severe Forms of COVID-19 in the Omicron Era
Next Article in Special Issue
Comparison of the Virulence of Korean African Swine Fever Isolates from Pig Farms during 2019–2021
Previous Article in Journal
Genomic Analysis of SARS-CoV-2 Alpha, Beta and Delta Variants of Concern Uncovers Signatures of Neutral and Non-Neutral Evolution
Previous Article in Special Issue
Development of a Next-Generation Vaccine Platform for Porcine Epidemic Diarrhea Virus Using a Reverse Genetics System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Isolation and Genetic Characterization of a Bovine Coronavirus KBR-1 Strain from Calf Feces in South Korea

1
Virus Disease Division, Animal and Plant Quarantine Agency, Gimcheon 39660, Korea
2
Department of Veterinary Medicine Virology Laboratory, College of Veterinary Medicine and Research Institute for Veterinary Science, Seoul National University, GwanAk-Ro 1, GwanAk-Gu, Seoul 08826, Korea
*
Author to whom correspondence should be addressed.
Viruses 2022, 14(11), 2376; https://doi.org/10.3390/v14112376
Submission received: 1 September 2022 / Revised: 25 October 2022 / Accepted: 26 October 2022 / Published: 27 October 2022
(This article belongs to the Special Issue State-of-the-Art Animal Virus Research in South Korea)

Abstract

:
Bovine coronavirus (BCoV) causes severe diarrhea in neonatal calves, winter dysentery in adult cattle, and respiratory disease in feedlot cattle, resulting in economic losses. A total of 16/140 calf diarrheic feces samples collected in South Korea between 2017 and 2018 were positive for BCoV. Phylogenetic analysis of the complete spike and hemagglutinin/esterase genes revealed that the 16 Korean BCoV strains belonged to group GIIa along with Korean strains isolated after 2000, whereas Korean BCoV strains isolated before 2000 belonged to group GI. Mice and goats inoculated with an inactivated KBR-1 strain (isolated from this study) generated higher antibody titers (96 ± 13.49 and 73 ± 13.49, respectively) when mixed with the Montanide01 adjuvant than when mixed with the Carbopol or IMS1313 adjuvants. Viral antigens were detected in the large intestine, jejunum, and ileum of calves inoculated with inactivated KBR-1 vaccine (104.0 TCID50/mL) at 14 days of post-challenge (DPC). However, no viral antigens were detected in calves vaccinated with a higher dose of inactivated KBR-1 strain (106.0 TCID50/mL) at 14 DPC, and they had high antibody titers and stable diarrhea scores. Currently, the group GIIa is prevalent in cows in South Korea, and although further research is needed in the future, the recently isolated KBR-1 strain has potential value as a new vaccine candidate.

1. Introduction

Bovine coronavirus (BCoV) is a member of the species Betacoronavirus within the family Coronaviridae [1]. BCoV possesses a single-stranded, positive-sense, non-segmented RNA genome of approximately 31 kb, which encodes five major structural proteins: nucleocapsid (N), transmembrane (M), spike (S), small membrane (E), and hemagglutinin/esterase (HE) [2]. Of these, the S protein is a particularly important mediator of host cell attachment, viral entry, and pathogenesis [3]. The S protein comprises the S1 subunit, which contains the dominant neutralizing epitopes, and the S2 subunit, which mediates viral membrane fusion. The HE protein exhibits an esterase receptor-destroying activity that may also be important for viral entry and serves as a second viral attachment protein to initiate infection [4].
BCoV is shed in the feces and nasal secretions, and is associated with three distinct clinical syndromes in cattle [5]: (neonatal) calf diarrhea [(N)CD)] [6,7]; winter dysentery (WD), characterized by hemorrhagic diarrhea in adults [7,8,9,10]; and respiratory infections in cattle of different ages (commonly observed as a part of the bovine respiratory disease complex, or shipping fever in feedlot cattle) [6,11,12]. In South Korea, WD and CD outbreaks were confirmed by surveys conducted during 2002–2004 [13] and 2004–2005 [14], respectively. BCoV isolated from Korean WD cases suggests that the virus has dual tropism and induces pathological changes in both the digestive and respiratory tracts of calves [14]. Another study showed that 1% of Korean goats have antibodies against BCoV [15], whereas another detected a BCoV-like virus in nasal swabs from non-captive Korean wild water deer [16]. A recent study performed genetic characterization of BCoV isolated from diarrheic pre-weaned native Korean calves from 2019–2021 [17]. Because BCoV is common in Korean cows that have recently been inoculated with the BCoV vaccine, cow farm owners and clinical veterinarians are proposing development of a new vaccine.
Here, we compared the sequences of the S and HE genes of BCoV strains detected in this study and compared them with BCoV strains isolate around the world. A vaccine candidate strain KBR-1 was used to inoculate goats and mice, and immunogenicity and optimal adjuvant efficacy were investigated. In addition, we also investigated the potential of the KBR-1 strain as a vaccine candidate for calves to protect them against infection by BCoV.

2. Materials and Methods

2.1. Diarrhea Samples and RNA Extraction

A total of 140 diarrhea samples were obtained from calves under 2 months-of-age. Samples were collected from farms by three clinical veterinarians between February 2017 and May 2018. Samples were collected from eight provinces: nine samples (three farms) from Gyenggi (GG), 11 samples (five farms) from Gangwon (GW), 38 samples (13 farms) from Chungnam (CN), six samples (three farms) from Chungbuk (CB), 24 samples (nine farms) from Gyeongnam (GN), eight samples (five farms) from Gyeongbuk (GB), nine samples (four farms) from Jeonbuk (JB), and 35 samples (14 farms) from Jeonnam (JN). A total of 57 and 83 calves were either vaccinated or unvaccinated against BCoV, respectively. Total RNA was extracted from 40–60 mg of fecal sample using a Qiagen RNeasy kit, and cDNA was prepared using a Helix CriptTM Easy cDNA synthesis kit.

2.2. RT-PCR and Sequences Analysis of BCoV

Total RNA extracted from fecal samples was subjected to a nested RT-PCR assay targeting the BCoV nucleocapsid gene; the primers and methods have been described previously [18]. Samples identified as BCoV-positive by the nested RT-PCR were subjected to complete gene amplification using six primer pairs targeting the S gene, and two primer pairs targeting the HE gene, as described previously [18]. The amplified PCR products were extracted and subcloned into the pGEM-T vector system II (Promega, Madison, WI, USA). Cloned plasmids corresponding to the BCoV S and HE genes were sequenced with T7 and SP6 sequencing primers and an ABI Prism 3730xl DNA sequencer (Cosmo Genetech Co., Seoul, Korea). After the genes from each strain were sequenced, multiple sequence alignment of the complete S and HE genes, as well as phylogenetic analysis, were carried out using Bio Edit Sequence Alignment Editor version 7.2 and MEGA 7.0 software [19], respectively. Bootstrap values based on 1000 replicates were used to assess the confidence limits of the branching. Phylogenetic trees were constructed based on the complete S and HE gene sequences of BCoV reference strains isolated in Asia, America, and Europe, including the 1994–2021 strains detected in South Korea and the OC43 strain (human coronavirus) was used as an outgroup. Complete S and HE gene sequences of BCoV reference strains were obtained from the National Center for Biotechnology Information (NCBI) GenBank database. The accession numbers of the reference sequences reported worldwide including the Korean strains are displayed in Figure 1 and Figure 2.

2.3. Isolation of BCoV

Monolayers of human rectal tumor (HRT-18) cells were used to isolate BCoV. The BCoV isolate was exposed to antibiotics from the initial to the 3rd passage, and the absence of major bovine enteric viruses (bovine rotavirus and bovine viral diarrhea virus) was confirmed at every passage using RT-PCR (LiliF-BD-Multi RT-PCR kit, iNtRON Biotechnology, Seongnam, Republic of Korea). In brief, supernatants from fecal samples were collected by centrifugation and passed through 0.45 µm filters. The fecal supernatants were treated for 30 min with 10 μg/mL trypsin (Gibco, Grand Island, NY, USA), with intermittent rocking. The cultured cells were then washed twice with plain medium and inoculated with fecal supernatant for 1 h. After the inoculum was discarded, cells were cultured for 3–4 days at 37 °C/5% CO2 in RPMI 1640 (Gibco, Grand Island, NY, USA) containing 5 μg/mL trypsin. The cells were observed daily for evidence of cytopathic effects (CPE); if no CPE were observed, isolation of BCoV was confirmed in an indirect immunofluorescence assay (IFA). IFA was performed with a mouse monoclonal antibody specific for BCoV (Invitrogen, Carlsbad, CA, USA), followed by staining with a fluorescence-conjugated goat anti-mouse antibody (anti-FITC, green; SeraCare, Gaithersburg, MD, USA). Nuclei were stained with DAPI (4′, 6-diamidino-2-phenylindole, blue; Sigma-Aldrich, St. Louis, MO, USA). Merged images (FITC and DAPI) were used for analysis. Virus titers (TCID50) were determined using the Reed and Muench method [20]. In brief, BCoV was serially diluted 10-fold in RPMI 1640 medium containing 5 μg/mL trypsin, and then inoculated into HRT-18 cells pre-seeded in 96-well plates (eight replicate wells per dilution). After incubation for 60 min, the medium was replaced with RPMI 1640. The plates were incubated at 37 °C/5% CO2 for 3 days. The medium was discarded from the wells and the plates were immunoreacted with an anti-BCoV monoclonal antibody for IFA. Every well containing FITC-stained cells was counted, and titers were expressed as TCID50/mL. The KBR-1 strain isolated at passage 40 was used for adjuvant selection in mice and goats, and to evaluate vaccine efficacy in calves.

2.4. Adjuvant Selection Using Mice and Goats

The KBR-1 (105.0 TCID50/mL) strain was obtained from HRT-18 cells and inactivated using binary ethyleneimine (BEI). The residual BEI in the inactivated viral fluid was neutralized by sodium thiosulfate. Three adjuvants (Carbopol (0.02% w/v, Lubrizol, USA) [21]; Montanide01 (10% v/v, SEPPIC, France); and IMS1313 (30% v/v, SEPPIC, France)) were mixed with the inactivated KBR-1, in accordance with the manufacturer’s instructions (SEPPIC, France). Three groups of 12 goats (10 months-old) and three groups of 15 mice (15–20 g) were inoculated intramuscularly (IM) with three adjuvants to identify the best one. For each adjuvant, three goats/group were inoculated IM (1 mL/dose) twice, with an interval of 3 weeks; in addition, five mice received a dose of 0.2 mL. Three goats (10 months-old) and five mice (15–20 g) were used as non-vaccine mock controls. Blood samples were collected three times (before vaccination, at 3 weeks after the first vaccination, and at 4 weeks after the second vaccination) to detect BCoV antibodies.

2.5. Immunization of Calves and Examination of Vaccine Efficacy

All six of the calves used for vaccine efficacy evaluation were negative for BCoV antigen and antibodies. The RT-PCR kit (LiliF-BD-Multi RT-PCR kit, iNtRON Biotechnology, Seongnam, Republic of Korea) for cow enteric viruses (BCoV, bovine rotavirus A, and bovine viral diarrhea virus) was used to screen the antigen-negative calves, and an hemagglutination inhibition (HI) test was performed to screen BCoV antibody-negative calves. Four of six calves (2 months-old) per group were inoculated with a high (106.0 TCID50/mL; n = 2) or low (104.0 TCID50/mL; n = 2) dose of KBR-1 vaccine. The remaining two calves comprised the mock group (without vaccine). Calves were vaccinated twice (1 mL per dose), with a 3-week interval in-between. At 2 weeks after the second vaccination, all calves were inoculated orally with 5 mL (106.0 TCID50/mL) of virulent BCoV. Virulent BCoV was obtained from KVCC (Korea Veterinary Culture Collection, Republic of Korea, KVCC number: VR1800062). Blood, nasal, and fecal samples were collected at 0, 21, and 35 DPV (days post-vaccination), and at 1, 2, 3, 6, 9, 12, and 14 DPC (days post-challenge) to measure antibody and antigen levels. Fecal consistency was scored as follows: 0, solid; 1, pasty; 2, semi-liquid; and 3, severe liquid. Autopsy was performed at 14 DPC, and organs were tested for BCoV.

2.6. Serum Neutralizing Antibody and HI Tests

To measure serum neutralizing antibody titers, mouse and goat sera were diluted serially 2-fold in serum-free RPMI 1640 medium (Gibco) containing trypsin, and then mixed with 102.0 TCID50/0.1 mL BCoV. The virus-serum mixture was incubated for 1 h to neutralize the virus and then used to infect HRT-18 cells. After incubation for 3 days, plates containing cells not expressing serum neutralizing BCoV antigen were assessed by IFA. Neutralizing antibody titers were expressed as the reciprocal of the highest dilution that yielded 50% neutralization. Calf antiserum was mixed with a receptor-destroying enzyme to remove non-specific inhibitors. Next, 4–8 HA units of BCoV were added to serially diluted sera. Finally, 1% mouse erythrocytes (RBC) were added to generate RBC-BCoV antigen conjugates. An HI titer ≥ 1:20 was considered positive.

2.7. Real-Time RT-PCR

Real-time RT-PCR was used to quantify BCoV antigens from nasal and fecal samples. Total RNA was extracted using a micro column-based RNeasy Mini kit (Qiagen, CA, USA), and the real-time RT-PCR targeting the N gene [22,23] was performed using the Path-IDTM Multiplex One-Step RT-PCR kit (Life Technologies, Carlsbad, CA, USA). BCoV PCR primer sets (F primer: CTA GTA ACC AGG CTG ATG TCA ATA CC; R primer: GGC GGA AAC CTA GTC GGA ATA) and probe (FAM-CGC CTG ACA TTC TCG ATC–MGB) were used for the real-time RT-PCR [22,23]. The cycling conditions of real-time RT-PCR were as follows: reverse transcription (RT) for 30 min at 45 °C; activation of DNA polymerase for 10 min at 95 °C; 40 cycles of denaturation at 94 °C for 15 s, and annealing/extension at 60 °C for 60 s.

2.8. Statistical Analysis

All statistical analyses were performed using GraphPad Prism software, version 6.0, for Windows. Data are expressed as the mean ± SD.

3. Results

3.1. Prevalence of BCoV

The prevalence of BCoV on calf farms reporting diarrhea was 11.4% (16/140). In terms of age, six calves were between 1 and 20 days old, six were between 21 and 30 days old, and four were between 31 and 60 days old. With respect to regions, BCoV was detected in GW (1/11, 9.1%), CN (6/38, 15.8%), GN (2/24, 8.3%), JB (2/9; 22.2%), and JN (5/35, 14.3%). No cases were detected in GG (0/9), CB (0/6), or GB (0/8).

3.2. Phylogenetic Analysis of the S and HE Genes

Phylogenetic analysis revealed that the complete S gene was split into two groups (GI and GII). The GI group includes old strains (European, American, and Canadian), including classical strains (Figure 1). GII also split into two subgroups: GIIa (America and Asia) and GIIb (European). The GIIa subgroup contains Korean BCoV strains isolated between 2002 and 2021, including all of the strains detected in this study (Figure 1); it also included strains from Japan, China, Vietnam, Cuba, and the USA (Figure 1). The GIIb subgroup included independent strains from European countries, including France, Sweden, and Denmark (Figure 1). Phylogenetic analysis of the HE gene also revealed a split into two groups. Korean strains (SUN5, A3, and BC94) detected in 1994 were included in the GI subgroup, and all Korean strains detected between 2002 and 2021 were included in the GIIa subgroup (Figure 2). The 16 Korean strains detected in 2017 and 2018 were closely related to strains detected in Korea between 2019 and 2021 (Figure 2).

3.3. Comparison of the S and HE Gene Sequences

Comparative analysis of nucleotide (nt) sequence homology revealed that the S and HE genes of the three Korean strains (SUN5, A3, and BC94) detected in 1994 showed high similarity (99.5–99.6% and 99.8%, respectively) with the Mebus strain, which is virulent (Table 1).
However, the nt sequences of the S genes of Korean strains detected in 2017 and 2018 showed somewhat lower similarity (97.3–97.8% and 97.5–97.6%, respectively) with the Mebus strain (Table 1). In addition, the genes (S and HE genes) of strains isolated in Korea after 2000 showed low homology with the Mebus strain (Table 1). Comparison of the deduced amino acid sequences for the S protein of Korean BCoVs showed that the S protein of the Mebus strain has high similarity with that of SNU5 and BC94 strains isolated in 1994, although it harbors mutations in 34 amino acids when compared with Korean strains from 2017–2018 (Table 2). The HE gene of the Mebus strain harbored mutations at six amino acid positions when compared with the Korean BCoV strains from 2002–2021 (Table 3).

3.4. Non-CPE of the BC94 and KBR-1 Strains

Of the 16 Korean BCoVs examined herein, only the KBR-1 strain, detected in the feces of a calf in the CN region in 2017, could be passaged in human rectal tissue 18 (HRT-18) cells. The virus yielded 103.0 TCID50/mL at passage 5, 104.2 TCID50/mL at passage 20, and 106.2 TCID50/mL at passage 40 (Figure 3A). The highest titers were detected between 48 (106.2 TCID50/mL) and 72 h (106.0 TCID50/mL) of culture (Figure 3B). The non-CPE of the KBR-1 strain was confirmed by IFA, staining with DAPI, on examination of merged images (FITC and DAPI), or on bright field microscopy (Figure 3C).

3.5. Antibody Titers and Protection of Mice and Goats Inoculated with the Inactivated KBR-1 Vaccine

Anti-BCoV antibody titers in mice at 2 weeks after the second inoculation with the inactivated KBR-1 vaccine were as follows: 25.3 log2 when mixed with Carbopol (10% v/v, Lubrizol, USA); 26.4 log2 when mixed with Montanide01 (50% v/v, SEPPIC, France); and 24.8 log2 when mixed with IMS1313 (10% v/v, SEPPIC, France) (Figure 4). The BCoV antibody titers in goats at 2 weeks after the second injection of the inactivated KBR-1 vaccine were as follows: 26.4 log2 when mixed with Carbopol, and 27.1 log2 when mixed with Montanide01, or 27.0 log2 when mixed with IMS1313 (Figure 4). After two inoculations with the inactivated KBR-1 vaccine, neither mice nor goats showed specific clinical signs.

3.6. Diarrhea Score, and Antigen and Antibody Levels, after Challenge of Calves with the Inactivated KBR-1 Vaccine

Two calves receiving the inactivated KBR-1 (106.0 TCID50/dose) vaccine mixed with the Montanide01 adjuvant showed a titer of 40–80 HI units at 2 weeks after the second vaccine, and 160 HI units at DPC 14; no antigen was detected in nasal, fecal, and blood samples during the observation period after vaccination and challenge (Table 4). During the observation period, calves vaccinated with the inactivated KBR-1 (106.0 TCID50/dose) had normal feces (Table 4). However, calves vaccinated with a lower dose of the inactivated KBR-1 (104.0 TCID50/dose) vaccine had low antibody titers (20 HI units at DPV 35 and 40–80 HI units at DPC 14) (Table 4).
BCoV was detected in feces at DPC 1–3 and at DPC 2–6, respectively, with a diarrhea score of 1 and 2 at DPC 3 and 6, respectively (Table 4). For all calves in the mock group exposed to virulent BCoV, virus was detected in rectal swabs for 8–10 days (at DPC 3–12 and at DPC 2–9), which was longer than in the vaccination groups. The mock group also developed diarrhea, with a diarrhea score of 1 at DPC 2–14 and 2 at DPC 2–12 (Table 4). Real-time RT-PCR showed that all nasal samples were negative. However, in calves vaccinated with a lower dose of the inactivated KBR-1 (104.0 TCID50/dose) vaccine, virus was detected in feces for 1–3 DPC (ct value: 28.7–33.1) and 2–6 DPC (ct: 29.1–32.9) post-challenge (Table 4). In addition, virus was detected in feces from the mock group for 3–12 DPC (ct: 20.7–31.1) and 2–9 DPC (ct: 22.5–29.2) (Table 4). None of the calves died during the experimental period. All calves were euthanized and autopsied at DPC 14. No BCoV was detected in 13 organs (including brain) harvested from calves vaccinated with the high dose of inactivated KBR-1 (106.0 TCID50/dose) vaccine (Table 5).
However, calves vaccinated with the lower dose of inactivated KBR-1 (104.0 TCID50/dose) vaccine harbored virus in the organs, particularly the large intestine (ct: 32.6 and 33.2), jejunum (ct: 28.9 and 30.5), and ileum (ct: 26.1 and 29.8) (Table 5). In addition, the mock group harbored virus in the organs, particularly the duodenum (ct: 31.8), jejunum (ct: 23.6 and 24.7), and ileum (ct: 26.2 and 25.5) (Table 5).

4. Discussion

Between 2010 and 2019, many countries reported positive BCoV infections in calf feces [24], including China (12.2–19.0%), Thailand (12%), India (16%), Vietnam (6.9%), Australia (21.6%), Iran (7.2%), and Algeria (20.7%). In Korea, the prevalence of BCoV in 2002–2003 was 33.0% [25] and remained at 5.4–15.6% until 2021 [13,17]. Although the overall prevalence of BCoV in Korea has tended to decrease, continuous and damaging BCoV infections on cow farms might be possible. The present study shows that the prevalence of BCoV in diarrheic calves on cow farms located nationwide is 11.4%. Prevalence rates are particularly high in Jeonbuk (22.2%), Chungnam (15.8%), and Jeonam (14.3%) provinces. This result shows that continuous epidemiological monitoring is required to reduce future BCoV outbreaks.
In this study, phylogenetic tree analysis of the S and HE genes identified three branches: GI (classical strains), GIIa (America/Asia strains), and GIIb (European strains). The phylogenetic tree of S and HE genes shows that Korean strains detected in 1994 belonged to GI, but Korean BCoV strains detected from 2002 to 2021 belonged to GIIa. All HE genes of Korean BCoV strains detected after 2002 were closely related to the R-AH187 strain (isolated in the USA in 2000) and all S genes were closely related to the LSU-94LSS-051-2 strain (isolated in the USA in 1994). Alignment of the putative amino acid sequences derived from the S genes of the Korean BCoV strains identified from 1994, 2002–2018, and 2019–2021 in this study highlighted differences in the dominant amino acids at positions 25, 34, and 35 of the S protein (compared with the Mebus strain), respectively. A previous study suggests that Korean BCoVs show more genetic variation than those isolated in other countries due to continuous evolution [17]. The most recent common ancestor (TMRCA) results published in a recent evolutionary study also suggest that Korean BCoVs have been influenced by viruses originating from the USA, which have adapted rapidly to their environment [17]. The BCoV S and HE gene sequences identified in South Korea have formed independent clusters since 2002; in particular, the 16 BCoVs from 2017–2018 are closely associated with the cluster of 15 BCoVs from 2019–2021. The fact that BCoV strains emerging recently in South Korea are included in a single cluster is thought to be due to continuous evolution within the same cluster since 2002. In addition, it is estimated that evolution within the same cluster will continue if BCoV does not flow from other countries.
Although further investigations are needed, ongoing outbreaks have led cow disease experts to suggest that the efficacy of current vaccines requires improvement. Of the samples used in the present study, 40.7% (57/140) were collected from calves born to 13 vaccinated cow farms and 59.3% (83/140) from calves born to 25 unvaccinated cow farms. However, most of the 16 BCoV-positive diarrhea cases identified in this study were calves born to vaccinated cows (62.5%; 10/16) in four farms, whereas 37.5% (6/16) were in calves born to unvaccinated cows in three farms. Given the high prevalence in BCoV-vaccinated farms, future studies should investigate whether the current BCoV vaccine protects against the most recent virulent BCoV strain.
The BCoV antibody titers in goats at 2 weeks after the second injection of the inactivated KBR-1 vaccine (105.0 TCID50/dose) in this study were high (27.1 log2) when mixed with Montanide01 and when mixed with IMS1313 (27.0 log2); however, it was low (26.4 log2) when mixed with Carbopol. IMS1313 and Montanide01 are commercial adjuvants comprising nanoparticles and an immunostimulatory compound, and a synthetic polymer, respectively (SEPPIC). According to the manufacturer, the two adjuvants enhance the strength and duration of the immune response. Montanide01 is an aqueous-based gel adjuvant, which facilitates easier vaccine delivery than IMS1313, which is emulsified in oil. In addition, we selected Montanide01 for the cattle experiments because it is recommended by the manufacturer as a vaccine adjuvant for cattle and small ruminants.
Here, we found that the candidate KBR-1 (106.0 TCID50/dose) vaccine mixed with the Montanide01 adjuvant protected calves against diarrhea and virus shedding, with little antigen detected in organs following challenge with a virulent virus strain; however, the inactivated KBR-1 vaccine at a lower dose (104.0 TCID50/dose) failed to protect against diarrhea and virus shedding, with viral antigen detectable in organs. The ct values (28.7–33.1 and 29.1–32.9) obtained for real-time RT-PCR of feces samples collected after challenge from calves vaccinated with inactivated KBR-1 (104.0 TCID50/dose) were relatively higher than those (20.7–31.1 and 22.2–29.2) from feces collected from mock controls after challenge. In addition, after euthanasia at DPC 14, the ct values were higher for organs in mock group calves than in those from inactivated KBR-1 group calves. These results may imply a weak protective effect, even in calves vaccinated with the inactivated KBR-1 (104.0 TCID50/dose) strain. A previous study suggested that an inactivated vaccine may be useful for preventing WD caused by BCoV [26]. Thus, using the candidate KBR-1 vaccine at a high concentration may be necessary for effective protection of newborn calves from CD or WD.

5. Conclusions

In conclusion, phylogenetic analysis revealed that all Korean BCoV strains isolated after 2002 belong to the GIIa group. Animal studies showed that the recently isolated KBR-1 vaccine strain (group GIIa) candidate was effective and safe and protects cows from infection by virulent BCoV. Although additional experiments with a large number of animals are required, the KBR-1 strain shows potential as a vaccine candidate.

Author Contributions

Conceptualization, J.S. and S.C.; Methodology, G.-N.P., S.S., K.-S.K. and B.-H.A.; Writing—review and editing, B.-H.H. and D.-J.A. All authors read and agreed to the published version of the manuscript.

Funding

This project was supported by grants (Project Code no. B-1543083-2022-23-03) from the Animal and Plant Quarantine Agency, Republic of Korea.

Institutional Review Board Statement

The study was approved by the Animal Care and Use Committee of the Animal and Plant Quarantine Agency (APQA) (Approval Number: 2021-598). All animal experiments were performed in accordance with the guidelines of the APQA Animal Care and Use Committee.

Informed Consent Statement

Not applicable.

Data Availability Statement

The nucleotide sequences of the 2017–2018 viruses obtained in this study were submitted to the GenBank database under accession numbers OP186313-OP186328 for the S gene and OP186329-OP186344 for the HE gene.

Acknowledgments

We would like to thank Sungkyun Kwon, Applevet Animal Hospital, Anseong City, Gyeonggi-do, South Korea, for providing animal diarrhea samples. We also thank Ra Mi Cha for providing virulent BCoV from KVCC (Korea Veterinary Culture Collection).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Groot, R.J.; Baker, S.C.; Baric, R.; Enjuanes, L.; Gorbalenya, A.E.; Holmes, K.V.; Perlman, S.; Poon, L.; Rottier, P.J.M.; Talbot, P.J.; et al. Family Coronaviridae. In Virus Taxonomy, Classification and Nomenclature of Viruses. Ninth Report of the International Committee on Taxonomy of Viruses, 1st ed.; King, A.M.Q., Adams, M.J., Carstens, E.B., Lefkowitz, E.J., Eds.; Elsevier Academic Press: San Diego, CA, USA, 2011; pp. 806–828. [Google Scholar]
  2. Lai, M.M.C.; Cavanagh, D. The molecular biology of coronaviruses. Adv. Virus Res. 1997, 48, 1–100. [Google Scholar]
  3. He, Q.; Guo, Z.; Zhang, B.; Yue, H.; Tang, C. First detection of bovine coronavirus in Yak (Bos grunniens) and a bovine coronavirus genome with a recombinant HE gene. J. Gen. Virol. 2019, 100, 793–803. [Google Scholar] [CrossRef] [PubMed]
  4. Desforges, M.; Desjardins, J.; Zhang, C.; Talbot, P.J. The acetyl-esterase activity of the hemagglutinin-esterase protein of human coronavirus OC43 strongly enhances the production of infectious virus. J. Virol. 2013, 87, 3097–3107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Saif, L.J.; Jung, K. Comparative pathogenesis of bovine and porcine respiratory coronaviruses in the animal host species and SARS-CoV-2 in humans. J. Clin. Microbiol. 2020, 58, 3291–3298. [Google Scholar] [CrossRef] [PubMed]
  6. Clark, M.A. Bovine coronavirus. Br. Vet. J. 1993, 149, 51–70. [Google Scholar] [CrossRef]
  7. Tsunemitsu, H.; Saif, L.J. Antigenic and biological comparisons of bovine coronaviruses derived from neonatal calf diarrhea and winter dysentery of adult cattle. Arch. Virol. 1995, 140, 1303–1311. [Google Scholar] [CrossRef] [PubMed]
  8. Saif, L.J. A review of evidence implicating bovine coronavirus in the etiology of winter dysentery in cows: An enigma resolved? Cornell Vet. 1990, 80, 303–311. [Google Scholar]
  9. Cho, K.O.; Halbur, P.G.; Bruna, J.D.; Sorden, S.D.; Yoon, K.J.; Janke, B.H.; Chang, K.O.; Saif, L.J. Detection and isolation of coronavirus from feces of three herds of feedlot cattle during outbreaks of winter dysentery-like disease. J. Am. Vet. Med. Assoc. 2000, 217, 1191–1194. [Google Scholar] [CrossRef]
  10. Tråvén, M.; Näslund, K.; Linde, N.; Linde, B.; Silván, A.; Fossum, C.; Hedlund, K.O.; Larsson, B. Experimental reproduction of winter dysentery in lactating cows using BCV—Comparison with BCV infection in milk-fed calves. Vet. Microbiol. 2001, 81, 127–151. [Google Scholar] [CrossRef]
  11. Lathrop, S.L.; Wittum, T.E.; Brock, K.V.; Loerch, S.C.; Perino, L.J.; Bingham, H.R.; McCollum, F.T.; Saif, L.J. Association between infection of the respiratory tract attributable to bovine coronavirus and health and growth performance of cattle in feedlots. Am. J. Vet. Res. 2000, 61, 1062–1066. [Google Scholar] [CrossRef]
  12. Ellis, J. What is the evidence that bovine coronavirus is a biologically significant respiratory pathogen in cattle? Can. Vet. J. 2019, 60, 147–152. [Google Scholar] [PubMed]
  13. Jeong, J.H.; Kim, G.Y.; Yoon, S.S.; Park, S.J.; Kim, Y.J.; Sung, C.M.; Jang, O.J.; Shin, S.S.; Koh, H.B.; Lee, B.J.; et al. Detection and isolation of winter dysentery bovine coronavirus circulated in Korea during 2002–2004. J. Vet. Med. Sci. 2005, 67, 187–189. [Google Scholar] [CrossRef] [PubMed]
  14. Park, S.J.; Kim, G.Y.; Choy, H.E.; Hong, Y.J.; Saif, L.J.; Jeong, J.H.; Park, S.I.; Kim, H.H.; Kim, S.K.; Shin, S.S.; et al. Dual enteric and respiratory tropisms of winter dysentery bovine coronavirus in calves. Arch. Virol. 2007, 152, 1885–1900. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Yang, D.K.; Hwang, I.J.; Kim, B.H.; Kweon, C.H.; Lee, K.W.; Kang, M.I.; Lee, C.S.; Cho, K.O. Serosurveillance of viral diseases in Korean native goats (Capra hircus). J. Vet. Med. Sci. 2008, 70, 977–979. [Google Scholar] [CrossRef] [Green Version]
  16. Kim, J.H.; Jang, J.H.; Yoon, S.W.; Noh, J.Y.; Ahn, M.J.; Kim, Y.; Jeong, D.G.; Kim, H.K. Detection of bovine coronavirus in nasal swab of non-captive wild water deer, Korea. Transbound. Emerg. Dis. 2018, 65, 627–631. [Google Scholar] [CrossRef] [Green Version]
  17. Kim, E.M.; Cho, H.C.; Shin, S.U.; Park, J.; Choi, K.S. Prevalence and genetic characterization of bovine coronavirus identified from diarrheic pre-weaned native Korean calves from 2019 to 2021. Infect. Genet. Evol. 2022, 100, 105263. [Google Scholar] [CrossRef]
  18. Park, S.J.; Jeong, C.; Yoon, S.S.; Choy, H.E.; Saif, L.J.; Park, S.H.; Kim, Y.J.; Jeong, J.H.; Park, S.I.; Kim, H.H.; et al. Detection and Characterization of Bovine Coronaviruses in Fecal Specimens of Adult Cattle with Diarrhea during the Warmer Seasons. J. Clin. Microbiol. 2006, 44, 3178–3188. [Google Scholar] [CrossRef] [Green Version]
  19. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  20. Reed, L.J.; Muench, H. A simple method of estimating fifty percent endpoints. Am. J. Hyg. 1938, 27, 493–497. [Google Scholar]
  21. Dey, A.K.; Burke, B.; Sun, Y.; Hartog, K.; Heeney, J.L.; Montefiori, D.; Srivastava, I.K.; Barnett, S.W. Use of a polyanionic carbomer, Carbopol971P, in combination with MF59, improves antibody responses to HIV-1 envelope glycoprotein. Vaccine 2012, 30, 2749–2759. [Google Scholar] [CrossRef] [Green Version]
  22. Cho, Y.I.; Kim, W.I.; Liu, S.; Kinyon, J.M.; Yoon, K.J. Development of a panel of multiplex real-time polymerase chain reaction assays for simultaneous detection of major agents causing calf diarrhea in feces. J. Vet. Diagn. Investig. 2010, 22, 509–517. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Chae, J.B.; Kim, H.C.; Kang, J.G.; Choi, K.S.; Chae, J.S.; Yu, D.H.; Park, B.K.; Oh, Y.S.; Choi, H.J.; Park, J. The prevalence of causative agents of calf diarrhea in Korean native calves. J. Anim. Sci. Technol. 2021, 63, 864–871. [Google Scholar] [CrossRef] [PubMed]
  24. Zhu, Q.; Li, B.; Sun, D. Advances in Bovine Coronavirus Epidemiology. Viruses 2022, 14, 1109. [Google Scholar] [CrossRef] [PubMed]
  25. Jeong, J.H.; Kim, G.Y.; Yoon, S.S.; Park, S.J.; Kim, Y.J.; Sung, C.M.; Shin, S.S.; Lee, B.J.; Kang, M.I.; Park, N.Y.; et al. Molecular analysis of S gene of spike glycoprotein of winter dysentery bovine coronavirus circulated in Korea during 2002–2003. Virus Res. 2005, 108, 207–212. [Google Scholar] [CrossRef]
  26. Takamura, K.; Matsumoto, Y.; Shimizu, Y. Field study of bovine coronavirus vaccine enriched with hemagglutinating antigen for winter dysentery in dairy cows. Can. J. Vet. Res. 2002, 66, 278–281. [Google Scholar]
Figure 1. Phylogenetic tree based on the complete S gene sequences of BCoV reference strains isolated worldwide, including the 2017–2018 strains detected in this study. A maximum-likelihood phylogenetic (LogL = −15690.97) tree was constructed using the Tamura-Nei model and the Nearest-Neighbor-Interchange method. Bootstrap values (1000 replicates) ≥ 70% are shown at the nodes, and the sequences identified in this study are marked by a red box.
Figure 1. Phylogenetic tree based on the complete S gene sequences of BCoV reference strains isolated worldwide, including the 2017–2018 strains detected in this study. A maximum-likelihood phylogenetic (LogL = −15690.97) tree was constructed using the Tamura-Nei model and the Nearest-Neighbor-Interchange method. Bootstrap values (1000 replicates) ≥ 70% are shown at the nodes, and the sequences identified in this study are marked by a red box.
Viruses 14 02376 g001
Figure 2. Phylogenetic tree based on the complete HE gene sequences of BCoV reference strains isolated worldwide, including the 2017–2018 strains detected in this study. The phylogenetic (LogL = −3940.11) tree was constructed using the maximum-likelihood method (Tamura-Nei model), with bootstrap values (1000 replicates) ≥ 70%, in MEGA 7.0 software. The sequences identified in this study are marked by a red box.
Figure 2. Phylogenetic tree based on the complete HE gene sequences of BCoV reference strains isolated worldwide, including the 2017–2018 strains detected in this study. The phylogenetic (LogL = −3940.11) tree was constructed using the maximum-likelihood method (Tamura-Nei model), with bootstrap values (1000 replicates) ≥ 70%, in MEGA 7.0 software. The sequences identified in this study are marked by a red box.
Viruses 14 02376 g002
Figure 3. Proliferation and isolation of the KBR-1 strain using HRT-18 cells. Virus titer (according to passage number of the KBR-1 strain) detected in diarrhea samples from calves (A); virus titer at different times after inoculation with the KBR-1 strain (B); and virus detection and CPE (IFA, DAPI, merged images, and bright field microscopy) of the KBR-1 strain (C).
Figure 3. Proliferation and isolation of the KBR-1 strain using HRT-18 cells. Virus titer (according to passage number of the KBR-1 strain) detected in diarrhea samples from calves (A); virus titer at different times after inoculation with the KBR-1 strain (B); and virus detection and CPE (IFA, DAPI, merged images, and bright field microscopy) of the KBR-1 strain (C).
Viruses 14 02376 g003
Figure 4. Antibody titers in mice and goats inoculated with the inactivated KBR-1 vaccine. Anti-BCoV antibody titers at 2 weeks post-injection of the inactivated KBR-1 vaccine mixed with Carbopol, Montanide01, or IMS1313.
Figure 4. Antibody titers in mice and goats inoculated with the inactivated KBR-1 vaccine. Anti-BCoV antibody titers at 2 weeks post-injection of the inactivated KBR-1 vaccine mixed with Carbopol, Montanide01, or IMS1313.
Viruses 14 02376 g004
Table 1. Comparison of the nucleotide and deduced amino acid sequences of the BCoV spike and hemagglutinin esterase genes according to year.
Table 1. Comparison of the nucleotide and deduced amino acid sequences of the BCoV spike and hemagglutinin esterase genes according to year.
StrainGeneKorean BCoV Strains a
199420022004201520172018201920202021
MebusSpike99.5–99.6 b97.5–98.297.9–98.1-97.3–97.897.5–97.697.6–97.797.5–97.697.5–97.6
(98.8–98.9) c(96.1–97.6)(97.1–97.5)-(96.4–97.6)(97.1–97.4)(97.0)(97.2–97.4)(97.1–97.4)
HE99.897.5–98.097.7–98.097.9–98.097.6–98.097.6–97.896.9–97.697.3–98.097.5–97.8
(99.3–99.5)(97.4–98.6)(98.6–99.1)(98.8)(98.1–98.8)(98.4–98.6)(97.1–97.9)(97.4–98.8)(98.1–98.6)
a The annual detection lists of Korean BCoV strains are shown in Figure 1 (spike genes) and in Figure 2 (HE genes). b Nucleotide sequence homology (%). c Amino acid sequence homology (%).
Table 2. Comparison of the deduced amino acid sequences for the Spike genes of Korean BCoVs.
Table 2. Comparison of the deduced amino acid sequences for the Spike genes of Korean BCoVs.
Ref.
Strain
Korean BCoV Strains a
Year1971199420022004200520172018201920202021
Spike aaMebusBCoV
(n = 2)
BCoV
(n = 19)
BCoV
(n = 10)
BCoV
(n = 0)
BCoV
(n = 13)
BCoV
(n = 3)
BCoV
(n = 4)
BCoV
(n = 5)
BCoV
(n = 6)
11M**/T*/T-TTTTT
16I***-***/V*/V*/V
33A*VV-VVVVV
40I*/TTT-TTTTT
88R*T*/T-TTTTT
100I*TT-TTTTT
115K*DD-DDDDD
147L*FF-FFFFF
149N*SS-SSSSS
169H*NN-NNNNN
173H*N*/N-NNNNN
179KR/QR/QQ-QQQQQ
248L*MM-MMMMM
253S*/YNN-NNNNN
256M**LL-LLLLL
458FSSS-SSSSS
465V*A*/A-AAAAA
470H*DD-*/DDDDD
484STTT-TTTTT
499N**/S/T/KS/T-TTTTT
501P**/SS-S/YSSF/S*/S
510S*TT-TTTTT
531ND*/D*/D-DDDDD
543S*AA-AAAAA
578T*SS-SSSSS
769A*SS-SSSSS
779ANNN-NNNNN
965V*EE-EEEEE
984LWWW-WWWWW
988VAAA-AAAAA
1026D*GG-GGGGG
1030E**/DD-DDDDD
1100VAAA-AAAAA
1241H*P*/P-PPPPP
1341I**/KK-KKKKK
a The annual detection lists of Korean BCoV strains are shown in Figure 1 (spike genes). Asterisks (*) indicate no substitution.
Table 3. Comparison of the deduced amino acid sequences for the HE genes of Korean BCoVs.
Table 3. Comparison of the deduced amino acid sequences for the HE genes of Korean BCoVs.
Ref.
Strain
Korean BCoV Strains a
Year1971199420022004200520172018201920202021
HE aaMebusBCoV
(n = 3)
BCoV
(n = 19)
BCoV
(n = 10)
BCoV
(n = 2)
BCoV
(n = 13)
BCoV
(n = 3)
BCoV
(n = 4)
BCoV
(n = 5)
BCoV
(n = 6)
5L*P*/PPPP*/PP/HP/H
29N*TTTTTTTT
46D*GGGGGG*/GG
339S*NNSNN*/N*/N*/N
367S*PPDPPPPP
392L*IIIIIIII
a The annual detection lists of Korean BCoV strains are shown in Figure 2 (HE genes). Asterisks (*) indicate no substitution.
Table 4. Detection of BCoV antigen and antibodies, and diarrhea scores after vaccination and challenge with BCoV KBR-1.
Table 4. Detection of BCoV antigen and antibodies, and diarrhea scores after vaccination and challenge with BCoV KBR-1.
GroupNo. CowVaccinationChallenge (after Vaccination) a
DPV
0
DPV
21
DPV
35
DPC 1DPC 2DPC 3DPC 6DPC 9DPC 12DPC 14
KBR-1 (106.0TCID50/dose)A1b/<10 c/0 d/40/0/80/0-/NT/0-/NT/0-/80/0-/160/0-/160/0-/160/0-/160/0
A2-/<10/0/20/0/40/0-/NT/0-/NT/0-/160/0-/80/0-/160/0-/160/0-/160/0
KBR-1 (104.0TCID50/dose)B1-/<20/0/20/0/20/030.3/NT/028.7/NT/033.1/40/2-/40/1-/80/0-/40/0-/80/0
B2-/<10/0/10/0/20/0-/NT/029.1/NT/032.9/40/230.8/20/2-/40/0-/40/0-/40/0
MockC1-/<20/0/10/0/10/0-/NT/0-/NT/123.4/10/120.7/10/226.8/20/231.1/10/2-/20/1
C2-/<10/0/10/0/10/0-/NT/028.1/NT/122.5/20/223.1/10/229.2/20/1-/20/2-/40/0
a Challenge (virulent BCoV: 106.0 TCID50/mL); b Real-time RT-PCR (ct value); c HI (hemagglutinin inhibition) titer; d diarrhea score (0: normal, 1: mild, 2: moderate, 3: watery).
Table 5. Detection of BCoV in organs of vaccinated calves challenged with BCV KBR-1.
Table 5. Detection of BCoV in organs of vaccinated calves challenged with BCV KBR-1.
GroupNo. CowOrgans
Br *LvSpLnHtKdSITsLiDdJiIlMl
KBR-1 (106.0TCID50/mL)A1-------------
A2-------------
KBR-1 (104.0TCID50/mL)B1--------32.6-28.926.1-
B2--------33.2-30.529.8-
MockC1----------23.626.2-
C2---------31.824.725.5-
* Br., brain; Lv., liver; Sp., spleen; Ln., lung; Ht., heart; Kd., kidney; Sl., submandibular lymph node; Ts., tonsil; Li., large intestine; Dd., duodenum; Jj., jejunum; Il., ileum; Ml., mesenteric lymph nodes.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shin, J.; Choe, S.; Park, G.-N.; Song, S.; Kim, K.-S.; An, B.-H.; Hyun, B.-H.; An, D.-J. Isolation and Genetic Characterization of a Bovine Coronavirus KBR-1 Strain from Calf Feces in South Korea. Viruses 2022, 14, 2376. https://doi.org/10.3390/v14112376

AMA Style

Shin J, Choe S, Park G-N, Song S, Kim K-S, An B-H, Hyun B-H, An D-J. Isolation and Genetic Characterization of a Bovine Coronavirus KBR-1 Strain from Calf Feces in South Korea. Viruses. 2022; 14(11):2376. https://doi.org/10.3390/v14112376

Chicago/Turabian Style

Shin, Jihye, SeEun Choe, Gyu-Nam Park, Sok Song, Ki-Sun Kim, Byung-Hyun An, Bang-Hun Hyun, and Dong-Jun An. 2022. "Isolation and Genetic Characterization of a Bovine Coronavirus KBR-1 Strain from Calf Feces in South Korea" Viruses 14, no. 11: 2376. https://doi.org/10.3390/v14112376

APA Style

Shin, J., Choe, S., Park, G. -N., Song, S., Kim, K. -S., An, B. -H., Hyun, B. -H., & An, D. -J. (2022). Isolation and Genetic Characterization of a Bovine Coronavirus KBR-1 Strain from Calf Feces in South Korea. Viruses, 14(11), 2376. https://doi.org/10.3390/v14112376

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop