Skip to main content

Phylogenetic investigation of enteric bovine coronavirus in Ireland reveals partitioning between European and global strains



Bovine coronavirus is a primary cause of neonatal calf diarrhea worldwide, and is also associated with acute diarrhea in adult cattle during the winter season. There are no reports on molecular characterization of bovine coronavirus in Ireland, and little data exists apart from serological studies.


In this study, 11 neonatal (mean age 9 days) calf BCoV strains from the south of Ireland were collected over a one year period and characterized using molecular methods. The spike gene which encodes a protein involved in viral entry, infectivity and immune response shows the most variability amongst the isolates and was subsequently selected for in depth analysis. Phylogenetic analysis of the spike gene revealed that the Irish strains clustered with novel BCoV strains from Europe in a unique clade, possibly indicating lineage partitioning. Direct analysis of alignments identified amino acid changes in the spike protein unique to the Irish clade.


Thus, monitoring of bovine coronavirus in Ireland is important as the current isolates in circulation in the south of Ireland may be diverging from the available vaccine strain, which may have implications regarding future BCoV vaccine efficacy.


Coronaviruses (CoV, family Coronaviridae) are large enveloped viral particles containing a positive sense single stranded RNA genome (26–30 kb), coding for several structural proteins, including polymerase (Pol), nucleocapsid (N), membrane (M), hemagglutinin-esterase (HE), spike (S) proteins and several non-structural proteins (NSPs). Coronaviruses have been associated with respiratory and enteric infections in humans and ruminants [1].

Enteric Bovine coronavirus (BCoV) replicates in the epithelial cells of the gut, destroying villi, resulting in severe, often bloody diarrhea in calves, which can be life threatening, due to loss of electrolytes and malnutrition [2]. Disease in calves usually occurs within the first month of life, with respiratory and enteric infections being the most common conditions diagnosed. In adult cows, as a result of close confinement during transport or housing, BCoV is associated with winter dysentery and shipping fever [2]. The spike proteins of BCoV play an important role in immune response, eliciting both cellular immune responses and neutralizing antibodies [3].

BCoV has been detected in Ireland using molecular or immunological techniques [46], but it has not been characterized or compared to other global BCoV strains. Enteric pathogens frequently isolated from neonatal calves with enteritis in Ireland are rotavirus, cryptosporidium and much less frequently, coronavirus [7]. Currently in Ireland, a trivalent vaccine is licensed for the immunization of pregnant cows against rotavirus, coronavirus and Escherichia coli, confering passive immunity to calves, the coronavirus aspect of the vaccine is based on an inactivated Mebus strain. [8] In this study we aimed to: (i) characterize bovine coronavirus in the South of Ireland via analysis of the Spike gene, and (ii) compare Irish BCoV to global and vaccine isolates to identify variations in the hyper-variable region of the spike gene.


Faecal samples were collected from the Cork Regional Veterinary Laboratory (CRVL) after they had tested positive for coronavirus using an immunochromatographic commercial kit (ICK), Corona Vet (Serosep, Ireland). Faecal samples were also screened for rotavirus, cryptosporidium and Salmonella. A total of 11 coronavirus positive samples were collected from neonatal calves, mean age 9 days, presenting with diarrhea between 2010 and 2011. Samples were stored at −80 °C prior to analysis.

Prior to extraction, faecal samples were homogenized in an equal volume of 0.89 % NaCl, centrifuged and filtered using a 0.20 μm pore size. The RNA was then extracted from the cell free fluid using Qiagen Viral RNA mini kit (Qiagen), following the manufacturers’ instructions. Extracted RNA was stored at −20 °C prior to analysis. Extracted RNA was tested for the presence of Coronavirus using degenerate oligonucleotide primers described previously [9], targeting a 250 bp region of the polymerase gene. A nested PCR was used to amplify the spike (S) gene [10], specifically the hypervariable region (HVR) [11].

Following analysis of this region, the most variable isolate was selected for complete characterization of the S gene using primers previously described [12, 13]. Reactions were carried out using Enhanced Avian Reverse Transcriptase kit (Sigma-Aldrich), following the manufacturers’ instructions, all reactions were carried out using a Biometra T3000 thermocycler. Amplified products were run on 1.5 % agarose gels, stained with ethidium bromide and visualized using a UV light trans-illuminator. Bands containing positive samples were cleaned using Roche High Pure PCR clean kit (Roche) and sequenced using a commercial service (MWG Eurofins, Germany).

Resulting sequence data was then analysed and edited using Bioedit v7.0.9.0 [14] and online BLAST tool (, to identify homologous strains. Sequence alignment was carried out using Clustal W in Bioedit [14], and a sequence alignment profile generated (Fig. 1). For analysis of the whole S gene, contigs were assembled using DNAstar program Seqman. The phylogenetic trees for the S gene were constructed using Maximum likelihood (ML) in MEGA5.1 [15] with a GTR model, plus gamma distribution and invariant sites with 1000 bootstrap replicates (Figs. 2 and 3). In ML trees shown, all strains are displayed with accession numbers.

Fig. 1

Amino acid sequence alignment profile of the hypervariable region of the spike protein. The profile was produced using Clustal W alignment and Bioedit. Amino acids in columns of the alignment with no alterations are shown as dots. A unique amino acid change occurs at position 21 in isolate RVLC4, novel alterations were detected at position 60 in RVLC7 and position 149 for all isolates

Fig. 2

Maximum Likelihood tree based on partial length of the spike gene specifically the HVR (464 bp), nt positions 1324 to 1787. All Irish isolates from this study are represented by a filled circle ()

Fig. 3

Maximum likelihood tree based on the complete spike gene (4092 bp), including the S1 and S2 subunits. Irish isolates from this study are indicated by a filled circle () in the tree

Partial nucleotide sequences are available for Polymerase gene (JN204179, JN204180), Spike (KF272908, KF272910, KF272912, KF272914). The complete gene sequence of the spike protein of Irish BCoV strain RVLC7 was registered in Genbank under the accession number-KF272919.


Eleven coronavirus ICK positive specimens were collected from the CRVL and then subjected to further analysis at the Virology Unit, Department of Biological Sciences, Cork Institute of Technology (CIT). Of the 11 positive faecal samples tested using an ICK, 8/11 (72.7 %) tested positive for coronavirus using molecular techniques (Table 1). Initial analysis was carried out using the online BLAST tool, Bioedit [14], and using MEGA5.1 [15] for phylogenetic analysis. Initially, a partial sequence representing the HVR of the spike was characterized, 7/8 strains were successfully amplified, and 4 representative strains were selected and sequenced.

Table 1 Summary of the 11 isolates used in this study, number as per main text

Sequence alignment of the spike HVR showed that Irish strains had novel substitutions compared to the Mebus strain sharing between 92.2–93.5 % amino acid identity. The hypervariable region of the spike gene (S) had 7 single nucleotide polymorphisms (SNPs) in a 464 bp sequence, resulting in 3 amino acid changes (Fig. 1); the amino acid changes occurred at positions 21, 60 and 149. Position 21 had a substitution from proline to leucine, found only in isolate RVLC4. Position 60 shows two different residue changes, the most common being proline to serine (RVLC4, 9, 10 However, only 3/4 of the Irish strains contained this proline to serine substitution, while the fourth contained a change from proline to phenylalanine (RVLC 7), which has a different polarity and lacks the ability to be phosphorylated. The position 149 shift was found in all strains and resulted in an aspartic acid residue being changed to glutamic acid here.

From alignment analysis of the HVR sequences, one isolate was selected for complete sequencing of the whole S gene, using primers previously described [12, 13]. Isolate RVLC7 was selected due to the novel polymorphism identified in position 60 of the HVR (proline to phenylalanine, corresponding to position 501 in the complete S protein). As depicted in Table 2, the complete spike protein contained a total of 14 residue changes when compared to the Mebus strain; 9 of these changes were shared with other strains found in Europe (light grey regions in Table 2), while 5 were unique to Irish isolate RVLC7 (dark grey regions in Table 2).

Table 2 Sub-sampling of the complete alignment of the spike gene


In 2012, BCoV was the most commonly isolated virus from live clinical cases of pneumonia and is also an important cause of enteric disease in young Irish calves [7]. In this study, we aimed to characterize bovine coronavirus in Ireland, and compare Irish BCoV strains with field and vaccine global strains. The partial gene tree for the hypervariable region of the S1 subunit of the spike gene showed that the Irish isolates clustered together with other European isolates (Italian, Swedish and Dutch) (Fig. 2). The phylogenetic tree of the complete spike gene shows that the European isolates form a distinct clade in comparison to American, Canadian and Korean isolates these different clades may represent differences in antigenicity [16] and possibly different lineages (Fig. 3). The Irish isolates group together, clustering with European isolates, in a distinct clade. It is known that viral lineages can form natural groups based on geographical location as can be seen with human Group A Rotavirus[17]. Sequence alignment profiles of the HVR revealed a polymorphism in all Irish isolates (position 149 in Fig. 1, position 590 in the whole protein), a change from aspartic acid to glutamic acid. This change is also present in Human Coronavirus 4408, a BCoV-like isolate from a child in USA [18]. A previously identified polymorphism was identified in Irish isolate RVLC7 [11], (position 60 in Fig. 1, residue 501 in whole protein), a substitution from proline to phenylalanine was detected, this substitution has been associated immunological escape mutants through changes in protein secondary structure [12, 19]. The vaccine licensed for use in Ireland is based on the Mebus strain, which is in a different clade. This phylogenetic difference between wild type and vaccine strains has been reported previously in other jurisdictions [20]. The data presented here suggests partitioning of Irish BCoV wild type strains away from the vaccine strain, and we identified the presence of novel mutation in spike HVR in wild type strain (RVLC 7, phenylalanine 501).


Monitoring of bovine coronavirus in Ireland is important as the data suggests that the current isolates in circulation maybe diverging. In addition, phylogenetic analysis demonstrated clear differences between European and global BCoV strains, therefore, continuous monitoring of BCoV is essential in order to detect new BCoV strains which may emerge and provide data which can be used to inform future BCoV vaccine design.


  1. 1.

    Lai MMC, Perlman S, Anderson LJ. Coronaviridae. In: Straus SE, editor. Fields virology. Philadelphia: Lippincott, Williams and Wilkins; 2007. p. 1917–74.

    Google Scholar 

  2. 2.

    Clark MA. Bovine coronavirus. Br Vet J. 1993;149:51–70.

    PubMed  CAS  Article  Google Scholar 

  3. 3.

    Deregt D, Babiuk LA. Monoclonal antibodies to bovine coronavirus: characteristics and topographical mapping of neutralizing epitopes on the E2 and E3 glycoproteins. Virology. 1987;161:410–420.

  4. 4.

    Gunn H. Elimination of transmissible gastroenteritis virus from a pig farm by culling and serological surveillance. Vet Rec. 1996;9:196–8.

    Article  Google Scholar 

  5. 5.

    Priestnall SL, Brownlie J, Dubovi EJ, Erles K. Serological prevalence of canine respiratory coronavirus. Vet Microbiol. 2006;115:43–53.

    PubMed  Article  Google Scholar 

  6. 6.

    Thompson HP, Dooley JSG, Kenny J, McCoy M, Lowery CJ, Moore JE, et al. Genotypes and subtypes of Cryptosporidium spp. in neonatal calves in Northern Ireland. Parasitol Res. 2007;100:619–24.

    PubMed  Article  Google Scholar 

  7. 7.

    AFBI / DAFF Veterinary Laboratories Publication. All-Island Animal Disease Surveillance Report. 2013.

  8. 8.

    Mebus CA, White RG, Stair EL, Rhodes MB, Twiehaus MJ. Neonatal calf diarrhea: results of a field trial using a reo-like virus vaccine. Vet Med Small Anim Clin. 1972;67:173–4.

    PubMed  CAS  Google Scholar 

  9. 9.

    Stephensen CB, Casebolt DB, Gangopadhyay NN. Phylogenetic analysis of a highly conserved region of the polymerase gene from 11 coronaviruses and development of a consensus polymerase chain reaction assay. Virus Res. 1999;60:181–9.

    PubMed  CAS  Article  Google Scholar 

  10. 10.

    Brandão PE, Gregori F, Monteleone GS, Soares RM, Rosales CAR, Jerez JA. Nested pcr assay for detection of bovine coronavirus S1 gene. Arq Inst Biol. 2003;70:1–3.

    Google Scholar 

  11. 11.

    Rekik MR, Dea S. Comparative sequence analysis of a polymorphic region of the spike glycoprotein S1 subunit of enteric bovine coronavirus isolates. Arch Virol. 1994;135:319–31.

    PubMed  CAS  Article  Google Scholar 

  12. 12.

    Hasoksuz M, Sreevatsan S, Cho KO, Hoet AE, Saif LJ. Molecular analysis of the S1 subunit of the spike glycoprotein of respiratory and enteric bovine coronavirus isolates. Virus Res. 2002;84:101–9.

    PubMed  CAS  Article  Google Scholar 

  13. 13.

    Martínez N, Brandão PE, de Souza SP, Barrera M, Santana N, de Arce HD, et al. Molecular and phylogenetic analysis of bovine coronavirus based on the spike glycoprotein gene. Infect Genet Evol. 2012;12:1870–8.

    PubMed  Article  Google Scholar 

  14. 14.

    Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999;41:95–8.

    CAS  Google Scholar 

  15. 15.

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  16. 16.

    Kanno T, Ishihara R, Hatama S, Uchida I. Antigenic variation among recent Japanese isolates of bovine coronaviruses belonging to phylogenetically distinct genetic groups. Arch Virol. 2013;158:1047–53.

    PubMed  CAS  Article  Google Scholar 

  17. 17.

    Dennis AF, McDonald SM, Payne DC, Mijatovic-Rustempasic S, Esona MD, Edwards KM, et al. Molecular epidemiology of contemporary G2P[4] human rotaviruses cocirculating in a single U.S. community: footprints of a globally transitioning genotype. J Virol. 2014;88:3789–801.

    PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Zhang XM, Herbst W, Kousoulas KG, Storz J. Biological and genetic characterization of a hemagglutinating coronavirus isolated from a diarrhoeic child. J Med Virol. 1994;44:152–61.

    PubMed  CAS  Article  Google Scholar 

  19. 19.

    Yoo D, Deregt D. A single amino acid change within antigenic domain II of the spike protein of bovine coronavirus confers resistance to virus neutralization. Clin Diagn Lab Immunol. 2001;8:297–302.

    PubMed  CAS  PubMed Central  Google Scholar 

  20. 20.

    Fulton RW, Ridpath JF, Burge LJ. Bovine coronaviruses from the respiratory tract: antigenic and genetic diversity. Vaccine. 2013;31:886–92.

    PubMed  CAS  Article  Google Scholar 

Download references


The authors wish to thank our colleagues at the CRVL, particularly Maria Colgan for assistance with collection of specimens.


Funding for this research was provided under the National Development Plan, through the Food Institutional Research Measure (grant no. 05/R&D/CIT/365, awarded to HOS) administered by the Department of Agriculture, Food and Marine, Ireland).

Author information



Corresponding author

Correspondence to H. O’Shea.

Additional information

Competing interests

None of the authors of this paper has a financial or personal relationship with other people or organisations that could inappropriately influence or bias the content of the paper. The authors declare that they have no competing interests.

Authors’ contributions

LG carried out the analysis and drafted the manuscript. PC, MJOC and HOS helped draft the manuscript and gave final approval of the version to be published. All authors read and approved the final manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gunn, L., Collins, P.J., O’Connell, M.J. et al. Phylogenetic investigation of enteric bovine coronavirus in Ireland reveals partitioning between European and global strains. Ir Vet J 68, 31 (2015).

Download citation


  • Bovine coronavirus
  • Spike protein
  • Sequence analysis