key: cord-265768-hwki5lk2
authors: Abi, Keha-mo; Zhang, Qi; Zhang, Bin; Zhou, Long; Yue, Hua; Tang, Cheng
title: An emerging novel bovine coronavirus with a 4-amino-acid insertion in the receptor-binding domain of the hemagglutinin-esterase gene
date: 2020-10-06
journal: Arch Virol
DOI: 10.1007/s00705-020-04840-y
sha: 
doc_id: 265768
cord_uid: hwki5lk2

The hemagglutinin-esterase (HE) protein of betacoronavirus lineage A is a secondary receptor in the infection process and is involved in the emergence of new betacoronavirus genotypes with altered host specificity and tissue tropism. We previously reported a novel recombinant bovine coronavirus (BCoV) strain that was circulating in dairy cattle in China, but this virus was not successfully isolated, and the genetic characteristics of BCoV are still largely unknown. In this study, 20 diarrheic faecal samples were collected from a farm in Liaoning province that had an outbreak of calf diarrhea (≤ 3 months of age) in November 2018, and all of the samples tested positive for BCoV by RT-PCR. In addition, a BCoV strain with a recombinant HE (designated as SWUN/A1/2018) and another BCoV strain with a recombinant HE containing an insertion (designated as SWUN/A10/2018) were successfully isolated in cell culture (TCID50: 10(4.25)/mL and 10(4.73)/mL, respectively). Unexpectedly, we identified the emergence of a novel BCoV variant characterized by a 12-nt bovine gene insertion in the receptor-binding domain in a natural recombinant HE gene, suggesting a novel evolutionary pattern in BCoV. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00705-020-04840-y) contains supplementary material, which is available to authorized users.

Betacoronavirus lineage A includes human coronavirus HKU1, human coronavirus OC43 (HCoV-OC43), equine coronavirus, mouse hepatitis virus (MHV), canine respiratory coronavirus (CRCoV), porcine hemagglutinating encephalomyelitis virus, and bovine coronavirus (BCoV). CRCoV and HCoV-OC43 are closely genetically related to BCoV [1, 2] . Members of this lineage have a unique gene encoding a surface hemagglutinin-esterase (HE) protein that is not found in other coronaviruses The HE protein serves as a secondary receptor in the infection process [3] and is involved in the emergence of novel genotypes and in shifts in host specificity and tissue tropism [4] . In addition, HE protein can induce the production of neutralizing antibodies against the BCoV [5] .

BCoV causes gastrointestinal and respiratory diseases in cattle, leading to serious economic losses all over the word [6] [7] [8] . No consistent genetic or antigenic markers have been identified for discriminating BCoV in different clinical syndromes [9, 10] , and the virus is thought to only has one genotype [10, 11] . BCoV can infect wild ruminant animals, including alpacas, sambar deer, and giraffes [8, 12, 13] and probably can be transmitted to dogs [1, 14, 15] . Interestingly, a coronavirus isolated from human diarrheic samples has been shown to be related to BCoV [16] , indicating its public health significance.

Recently, we reported the wide circulation of a novel BCoV strain with a recombinant HE in dairy cattle in China [17] . In this study, in order to monitor the epidemic of the recombinant strain, a total of 20 diarrheic faecal samples were collected from a farm that had an outbreak of diarrhea in calves three months of age or younger in Liaoning province in November, 2018, and all of the samples tested positive for BCoV by RT-PCR as described previously [17] . Full-length HE genes (1275 nt) were successfully cloned from 18 of the 20 BCoV-positive samples as described previously [18] [19] [20] , and PCR products were purified and cloned into the pMD19-T simple vector prior to sequencing. A Handling Editor: T. K. Frey.

The online version of this article (https ://doi.org/10.1007/s0070 5-020-04840 -y) contains supplementary material, which is available to authorized users. [17] clustered in a large independent branch (Fig. 1) , and three of the 18 HE gene sequences (GenBank accession numbers: MN982186, MN982190, and MN982199) clustered in an independent small branch. Further analysis showed that all 18 HE genes identified in this study had undergone the same recombination event identified in our previous report [17] . This suggests an increase in the prevalence of this recombinant, since it was found in only two out of five strains detected in the earlier study (samples collected in January 2018 [17] ) but was present in 18 of the 18 isolates from the same farm characterized in this study (samples collected in November 2018). This indicates that the HE recombinant BCoV strain is still present in Liaoning, China, and may be increasing in prevalence. The BCoV HE receptor-binding domain is composed of six surface loops, five of which are grafted onto the beta-sandwich core of the lectin domain, the R1-loop, R2-loop, R3-loop, R4-loop, and the RBS-hairpin, and the other is present on the esterase domain, the E-loop ( Fig. 2A) . The 18 isolates with a recombinant HE were sequenced and found to contain the same amino acid (aa) variation (F181V) in the R2-loop that was described in our previous report [17] . The R2-loop in the lectin domain acts in ligand binding, and aa substitutions in this domain may affect receptor binding. For example, HCoV-OC43 may have evolved during adaptation of BCoV to a human host by gradually losing its lectin activity due to progressive accumulation of aa mutations in the R2-loop of the lectin-binding domain in the HE protein [2] . Unexpectedly, three isolates that clustered in an independent small branch all contained an identical insertion of 12 nucleotides (AAG GCT ACT GTT ), resulting in four additional amino acids in the R3-loop (Fig. 2B ). This sequence was not present in any of the available BCoV HE sequences in the GenBank database. Interestingly, the inserted nucleotides were identified as originating from the bovine host, matching the sequence of a transcript variant of ORF159 on Bos taurus chromosome 16 C1 (GenBank accession no. XM027565276). It is worth noting that the R3-loop is composed of 13 aa (aa 207 -219) in BCoV HE, and residues F 211, L 212, S 213, and N 214 are essential for receptor-ligand interaction [3] . A 3D model constructed based on the crystal structure of bovine coronavirus hemagglutinin-esterase (SMTL ID: 3cl4.1) using the online software SWISSMODEL (https ://www.swiss model .expas y.org/inter activ e) showed that the position of the 4-aa insertion in the R3-loop between F 211 and F 212 in HE in the recombinant BCoV strains would alter the spatial conformation of the receptor-binding site relative to that in the HE recombinant strains and prototype strains lacking this insertion ( Fig. 2B and C) . MHV is also a member of betacoronavirus lineage A, and a variant has been reported with an aa insertion in the R3-loop that causes a change in its receptor specificity to 4-O-Ac-Sia from 9-O-Ac-Sia [4] . It is reasonable to speculate that the 4-aa insertion in the R3-loop in HE may affect HE receptor binding in the BCoV strain. The betacoronavirus HE protein is involved in receptor binding, host spectrum and tissue tropism [2, 4, 21] , and it would be valuable to further investigate the biological significance of this HE variant.

Betacoronavirus HE may have arisen from an influenza-C-like HE fusion protein (HEF), with transformation of the betacoronavirus HE from a trimer into a dimer for increased plasticity of the receptor-binding site (RBS) compared to that of HEF [3] . The integration of a bovine gene into other viruses has been reported previously. For example, insertion of the bovine SMT3B gene into the genome of bovine viral diarrhea virus correlates with the cytopathogenicity of the virus [22] . Interestingly, the inserted sequence in HE in this study may have originated from bovine chromosomal DNA, which may provide an evolutionary advantage for BCoV, and its identification should enhance our current understanding of the genetic evolution of BCoV.

The BCoV strain with a recombinant HE (designated as SWUN/A1/2018) and BCoV strain with a recombinant HE containing an insertion (designated as SWUN/A10/2018) successfully isolated in Vero cell culture. A cytopathic effect (CPE) was observed after three blind passages, and typical CPE characterized by cell rounding, cell lysis, and detachment from the culture flask was observed after 48 h at passage 6. After plaque purification, the TCID 50 values of strains SWUN/A1/2018 and SWUN/A10/2018 were calculated as 10 4.25 /mL and 10 4.73 /mL, respectively, by the Reed-Muench method [23] . To better understand the genetic evolution of SWUN/A10/2018, genomic RNA was successfully obtained from strains SWUN/A1/2018 and SWUN/ A10/2018 (GenBank accession number: MN982198 and MN982199). The two complete genomes shared 99.4% nt sequence identity with each other, and 97.0%-98.7% nt sequence identity with all 47 other BCoV genomes with sequences in the GenBank database. Phylogenetic trees constructed by the neighbor-joining method based on the full-length genome sequence and the ORF1a, ORF1b and S genes showed that the two strains have a close genetic relationship (see appendix). Thus, the BCoV SWUN/A10/2018 strain might have evolved from the SWUN/A1/2018 strain.

In summary, we report a novel bovine coronavirus variant with an insertion of 12 nucleotides in the HE gene that results in the addition of four amino acids in the receptorbinding portion of the HE protein. This disruption could alter the receptor-binding site by changing its spatial conformation, and it may reflect an emerging pattern in the evolution of BCoV. The betacoronavirus HE gene is involved in the emergence of new betacoronavirus genotypes, host specificity, and tissue tropism [2, 4, 21] . Thus, further investigation of the biological and epidemiological characteristics of this variant BCoV strain is warranted.

Isolation and sequence analysis of canine respiratory coronavirus

Betacoronavirus adaptation to humans involved progressive loss of hemagglutinin-esterase lectin activity

Structure of coronavirus hemagglutinin-esterase offers insight into corona and influenza virus evolution

The murine coronavirus hemagglutinin-esterase receptorbinding site: a major shift in ligand specificity through modest changes in architecture

Monoclonal antibodies to bovine coronavirus: characteristics and topographical mapping of neutralizing epitopes on the E2 and E3 glycoproteins

Factors affecting calf mortality in Iranian Holstein dairy herds

Molecular and antigenic characterization of bovine Coronavirus circulating in Argentinean cattle during 1994-2010

Market impacts of reducing the prevalence of bovine respiratory disease in United States beef cattle feedlots

Bovine respiratory coronavirus

Coronaviruses in cattle

Serological characterization of genotypically distinct enteric and respiratory bovine coronaviruses

Bovine-like coronaviruses isolated from four species of captive wild ruminants are homologous to bovine coronaviruses, based on complete genomic sequences

Detection and characterization of bovine-like coronaviruses from four species of zoo ruminants

The infectivity and pathogenicity of a group 2 bovine coronavirus in pups

Detection of a group 2 coronavirus in dogs with canine infectious respiratory disease

Biological and genetic characterization of a hemagglutinating coronavirus isolated from a diarrhoeic child

Prevalence of a novel bovine coronavirus strain with a recombinant hemagglutinin/esterase gene in dairy calves in China

Bovine coronaviruses associated with enteric and respiratory diseases in Canadian dairy cattle display different reactivities to anti-HE monoclonal antibodies and distinct amino acid changes in their HE, S and ns49 protein

Molecular epidemiology of human coronavirus OC43 reveals evolution of different genotypes over time and recent emergence of a novel genotype due to natural recombination

Molecular and phylogenetic analysis of bovine coronavirus based on the spike glycoprotein gene

Comparative analysis of 22 coronavirus HKU1 genomes reveals a novel genotype and evidence of natural recombination in coronavirus HKU1

Insertion of a bovine SMT3B gene in NS4B and duplication of NS3 in a bovine viral diarrhea virus genome correlate with the cytopathogenicity of the virus

A simple method of estimating fifty per cent endpoints

Acknowledgements This work was funded by the 13th Five-Year Plan National Science and Technology Support Program (Grant number