Discovery of seven novel Mammalian and avian coronaviruses in the genus deltacoronavirus supports bat coronaviruses as the gene source of alphacoronavirus and betacoronavirus and avian coronaviruses as the gene source of gammacoronavirus and deltacoronavirus.
Journal: 2012/June - Journal of Virology
ISSN: 1098-5514
Abstract:
Recently, we reported the discovery of three novel coronaviruses, bulbul coronavirus HKU11, thrush coronavirus HKU12, and munia coronavirus HKU13, which were identified as representatives of a novel genus, Deltacoronavirus, in the subfamily Coronavirinae. In this territory-wide molecular epidemiology study involving 3,137 mammals and 3,298 birds, we discovered seven additional novel deltacoronaviruses in pigs and birds, which we named porcine coronavirus HKU15, white-eye coronavirus HKU16, sparrow coronavirus HKU17, magpie robin coronavirus HKU18, night heron coronavirus HKU19, wigeon coronavirus HKU20, and common moorhen coronavirus HKU21. Complete genome sequencing and comparative genome analysis showed that the avian and mammalian deltacoronaviruses have similar genome characteristics and structures. They all have relatively small genomes (25.421 to 26.674 kb), the smallest among all coronaviruses. They all have a single papain-like protease domain in the nsp3 gene; an accessory gene, NS6 open reading frame (ORF), located between the M and N genes; and a variable number of accessory genes (up to four) downstream of the N gene. Moreover, they all have the same putative transcription regulatory sequence of ACACCA. Molecular clock analysis showed that the most recent common ancestor of all coronaviruses was estimated at approximately 8100 BC, and those of Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Deltacoronavirus were at approximately 2400 BC, 3300 BC, 2800 BC, and 3000 BC, respectively. From our studies, it appears that bats and birds, the warm blooded flying vertebrates, are ideal hosts for the coronavirus gene source, bats for Alphacoronavirus and Betacoronavirus and birds for Gammacoronavirus and Deltacoronavirus, to fuel coronavirus evolution and dissemination.
Relations:
Content
Citations
(275)
References
(52)
Diseases
(2)
Chemicals
(1)
Genes
(11)
Organisms
(12)
Processes
(4)
Affiliates
(1)
Similar articles
Articles by the same authors
Discussion board
J Virol 86(7): 3995-4008

Discovery of Seven Novel Mammalian and Avian Coronaviruses in the Genus <em class="genus-species">Deltacoronavirus</em> Supports Bat Coronaviruses as the Gene Source of <em class="genus-species">Alphacoronavirus</em> and <em class="genus-species">Betacoronavirus</em> and Avian Coronaviruses as the Gene Source of <em class="genus-species">Gammacoronavirus</em> and <em class="genus-species">Deltacoronavirus</em>

+4 authors

INTRODUCTION

Coronaviruses (CoVs) are found in a wide variety of animals, in which they can cause respiratory, enteric, hepatic, and neurological diseases of varying severity. Based on genotypic and serological characterization, CoVs were traditionally divided into three distinct groups (3, 22, 54). Recently, the Coronavirus Study Group of the International Committee for Taxonomy of Viruses has proposed three genera, Alphacoronavirus, Betacoronavirus, and Gammacoronavirus, to replace the traditional CoV groups 1, 2, and 3. As a result of the unique mechanism of viral replication, CoVs have a high frequency of recombination (22). Their tendency for recombination and the inherently high mutation rates in RNA virus may allow them to adapt to new hosts and ecological niches (18, 47).

The recent severe acute respiratory syndrome (SARS) epidemic, the discovery of SARS coronavirus (SARS-CoV), and the identification of SARS-CoV-like viruses from Himalayan palm civets and a raccoon dog from wild live markets in China have boosted interest in the discovery of novel CoVs in both humans and animals (5, 16, 33, 36, 39, 40, 46). A novel human CoV (HCoV) of the genus Alphacoronavirus, human coronavirus NL63 (HCoV-NL63), was reported independently by two groups in 2004 (12, 44). In 2005, we also described the discovery, complete genome sequence, clinical features, and molecular epidemiology of another novel HCoV, human coronavirus HKU1 (HCoV-HKU1), in the genus Betacoronavirus (24, 48, 50). As for animal CoVs, we and others have described the discovery of SARS-CoV-like viruses in horseshoe bats in Hong Kong Special Administrative Region (HKSAR) and other provinces of China (25, 30). Based on these findings, we conducted molecular surveillance studies to examine the diversity of CoVs in bats of our locality as well as of the Guangdong province of southern China, where the SARS epidemic originated and wet markets and game food restaurants serving bat dishes are commonly found. In these studies, at least nine other novel CoVs were discovered, including two novel subgroups in Betacoronavirus, subgroups C and D (26, 37, 45, 51). Other groups have also conducted molecular surveillance studies in bats and other animals, and additional novel CoVs were discovered and complete genomes sequenced (4, 6, 7, 9, 10, 1315, 17, 21, 31, 32, 34, 43, 53).

Birds are the reservoir of major emerging viruses, most notably, avian influenza viruses (29). Due to their flocking behavior and abilities to fly over long distances, birds have the potential to disseminate these emerging viruses efficiently among themselves and to other animals and humans. As for CoVs, the number of known CoVs in birds is relatively small compared to that in bats. Recently, we described the discovery of three novel CoVs in three families of birds, named bulbul coronavirus HKU11 (BuCoV HKU11), thrush coronavirus HKU12 (ThCoV HKU12), and munia coronavirus HKU13 (MunCoV HKU13) (49). These three CoVs formed a unique group of CoV, which probably represented a novel genus of CoV, Deltacoronavirus (8). We hypothesize that there are other previously unrecognized CoVs in this novel genus from mammals and other families of birds. To test this hypothesis, we carried out a territory-wide molecular epidemiology study in 3,137 mammals and 3,519 birds in HKSAR. Based on the results of comparative genome and phylogenetic analysis in the present study, we propose seven novel CoVs in Deltacoronavirus. Our model of bats and birds as the gene source of the four genera of coronaviruses is also discussed.

MATERIALS AND METHODS

Animal surveillance and sample collection.

All specimens of bats, cats, dogs, wild rodents, monkeys, and birds were collected with the assistance of the Department of Agriculture, Fisheries and Conservation, HKSAR, and those of pigs, cattle, chickens, and street rodents were collected with the assistance of the Department of Food, Environmental and Hygiene, HKSAR, from various locations in HKSAR over a 53-month period (February 2007 to June 2011). All specimens of Asian leopard cats were collected in the Guangdong province of southern China over an 8-month period (August 2010 to March 2011). Tracheal, rectal, and cloacal swabs were collected using procedures described previously (47, 49). Nasopharyngeal aspirates from humans were collected from patients in Queen Mary Hospital over a 13-month period (February 2010 to February 2011) (24, 47, 50). A total of 7,140 samples from 11 species of bats, 169 pigs, 230 cats, 231 dogs, 47 cattle, 221 chickens, 389 rodents, 235 monkeys, 1,397 humans, 15 Asian leopard cats, and 3,298 dead wild birds of 134 different species in 38 families had been tested.

RNA extraction.

Viral RNA was extracted from the tracheal, rectal, and cloacal swabs and nasopharyngeal aspirates using RNeasy Mini Spin column (Qiagen, Hilden, Germany) (27, 45, 47, 50). The RNA was eluted in 50 μl of RNase-free water and was used as the template for reverse transcription-PCR (RT-PCR).

RT-PCR of RdRp gene of CoVs using Deltacoronavirus conserved primers and DNA sequencing.

Initial CoV screening was performed by amplifying a 440-bp fragment of the RNA-dependent RNA polymerase (RdRp) gene of CoVs using Deltacoronavirus conserved primers (5′-GTGGVTGTMTTAATGCACAGTC-3′ and 5′-TACTGYCTGTTRGTCATRGTG-3′) designed by multiple alignments of the nucleotide sequences of available RdRp genes of BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13 (49). Reverse transcription was performed using the SuperScript III kit (Invitrogen, San Diego, CA). The PCR mixture (25 μl) contained cDNA, PCR buffer (10 mM Tris-HCl, pH 8.3, 50 mM KCl, 3 mM MgCl2, and 0.01% gelatin), 200 μM each deoxynucleoside triphosphate (dNTP), and 1.0 U Taq polymerase (Applied Biosystems, Foster City, CA). The mixtures were amplified with 60 cycles of 94°C for 1 min, 48°C for 1 min, and 72°C for 1 min and a final extension at 72°C for 10 min in an automated thermal cycler (Applied Biosystems, Foster City, CA). Standard precautions were taken to avoid PCR contamination, and no false positive was observed in negative controls.

The PCR products were gel purified using the QIAquick gel extraction kit (Qiagen, Hilden, Germany). Both strands of the PCR products were sequenced twice with an ABI Prism 3700 DNA analyzer (Applied Biosystems, Foster City, CA), using the two PCR primers. The sequences of the PCR products were compared with known sequences of the RdRp genes of CoVs in the GenBank database.

Complete genome sequencing.

Two complete genomes of porcine coronavirus HKU15 (PorCoV HKU15) and one complete genome each of white-eye coronavirus HKU16 (WECoV HKU16), sparrow coronavirus HKU17 (SpCoV HKU17), magpie robin coronavirus HKU18 (MRCoV HKU18), night heron coronavirus HKU19 (NHCoV HKU19), wigeon coronavirus HKU20 (WiCoV HKU20), and common moorhen coronavirus HKU21 (CMCoV HKU21) were amplified and sequenced using the RNA extracted from the original swab specimens as templates. The RNA was converted to cDNA by a combined random-priming and oligo(dT)-priming strategy. The cDNA was amplified by degenerate primers designed by multiple alignments of the genomes of other CoVs with complete genomes available, using strategies described in our previous publications (28, 45, 48, 49) and the CoV database CoVDB (20) for sequence retrieval. Additional primers were designed from the results of the first and subsequent rounds of sequencing. These primer sequences are available on request. The 5′ ends of the viral genomes were confirmed by rapid amplification of cDNA ends (RACE) using the 5′/3′ RACE kit (Roche, Germany). Sequences were assembled and manually edited to produce final sequences of the viral genomes.

Genome analysis.

The nucleotide sequences of the genomes and the deduced amino acid sequences of the open reading frames (ORFs) were compared to those of other CoVs using EMBOSS needle (http://www.ebi.ac.uk). Phylogenetic tree construction was performed using the neighbor joining method with ClustalX 1.83. Protein family analysis was performed using PFAM and InterProScan (1, 2). Prediction of transmembrane domains was performed using TMpred and TMHMM (19, 41).

Estimation of divergence dates.

Divergence times for the four genera of CoVs, Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Deltacoronavirus, were calculated using a Bayesian Markov chain Monte Carlo (MCMC) approach as implemented in BEAST (Version 1.6.1) as described previously (11, 23, 27, 47). One parametric model (Constant Size) and one nonparametric model (Bayesian Skyline) tree priors were used for the inference. Analyses were performed under the GTR+I+G substitution model for RdRp gene sequence data and using both a strict and an unrelaxed log-normal-distributed (Ucld) relaxed molecular clock. The MCMC run was 5 × 10 steps long, with sampling every 1,000 steps. Convergence was assessed on the basis of the effective sampling size after a 10% burn-in using Tracer software version 1.5 (11). The mean time of the most recent common ancestor (tMRCA) and the highest posterior density regions at 95% (HPD) (i.e., a credible set that contains 95% of the sampled values) were calculated, and the best-fitting model was selected by a Bayes factor, using marginal likelihoods implemented in Tracer (see Table S1 in the supplemental material) (42). Bayesian Skyline under a relaxed-clock model with Ucld was adopted for making inferences, as Bayes factor analysis indicated that this model fitted the data better than other models tested (see Table S1). The trees were summarized in a target tree by the Tree Annotator program included in the BEAST package by choosing the tree with the maximum sum of posterior probabilities (maximum clade credibility) after a 10% burn-in.

Nucleotide sequence accession numbers.

The nucleotide sequences of the eight genomes of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 have been lodged within the GenBank sequence database under accession no. {"type":"entrez-nucleotide","attrs":{"text":"JQ065042","term_id":"1027948168","term_text":"JQ065042"}}JQ065042 to {"type":"entrez-nucleotide","attrs":{"text":"JQ065049","term_id":"380005522","term_text":"JQ065049"}}JQ065049.

Animal surveillance and sample collection.

All specimens of bats, cats, dogs, wild rodents, monkeys, and birds were collected with the assistance of the Department of Agriculture, Fisheries and Conservation, HKSAR, and those of pigs, cattle, chickens, and street rodents were collected with the assistance of the Department of Food, Environmental and Hygiene, HKSAR, from various locations in HKSAR over a 53-month period (February 2007 to June 2011). All specimens of Asian leopard cats were collected in the Guangdong province of southern China over an 8-month period (August 2010 to March 2011). Tracheal, rectal, and cloacal swabs were collected using procedures described previously (47, 49). Nasopharyngeal aspirates from humans were collected from patients in Queen Mary Hospital over a 13-month period (February 2010 to February 2011) (24, 47, 50). A total of 7,140 samples from 11 species of bats, 169 pigs, 230 cats, 231 dogs, 47 cattle, 221 chickens, 389 rodents, 235 monkeys, 1,397 humans, 15 Asian leopard cats, and 3,298 dead wild birds of 134 different species in 38 families had been tested.

RNA extraction.

Viral RNA was extracted from the tracheal, rectal, and cloacal swabs and nasopharyngeal aspirates using RNeasy Mini Spin column (Qiagen, Hilden, Germany) (27, 45, 47, 50). The RNA was eluted in 50 μl of RNase-free water and was used as the template for reverse transcription-PCR (RT-PCR).

RT-PCR of RdRp gene of CoVs using Deltacoronavirus conserved primers and DNA sequencing.

Initial CoV screening was performed by amplifying a 440-bp fragment of the RNA-dependent RNA polymerase (RdRp) gene of CoVs using Deltacoronavirus conserved primers (5′-GTGGVTGTMTTAATGCACAGTC-3′ and 5′-TACTGYCTGTTRGTCATRGTG-3′) designed by multiple alignments of the nucleotide sequences of available RdRp genes of BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13 (49). Reverse transcription was performed using the SuperScript III kit (Invitrogen, San Diego, CA). The PCR mixture (25 μl) contained cDNA, PCR buffer (10 mM Tris-HCl, pH 8.3, 50 mM KCl, 3 mM MgCl2, and 0.01% gelatin), 200 μM each deoxynucleoside triphosphate (dNTP), and 1.0 U Taq polymerase (Applied Biosystems, Foster City, CA). The mixtures were amplified with 60 cycles of 94°C for 1 min, 48°C for 1 min, and 72°C for 1 min and a final extension at 72°C for 10 min in an automated thermal cycler (Applied Biosystems, Foster City, CA). Standard precautions were taken to avoid PCR contamination, and no false positive was observed in negative controls.

The PCR products were gel purified using the QIAquick gel extraction kit (Qiagen, Hilden, Germany). Both strands of the PCR products were sequenced twice with an ABI Prism 3700 DNA analyzer (Applied Biosystems, Foster City, CA), using the two PCR primers. The sequences of the PCR products were compared with known sequences of the RdRp genes of CoVs in the GenBank database.

Complete genome sequencing.

Two complete genomes of porcine coronavirus HKU15 (PorCoV HKU15) and one complete genome each of white-eye coronavirus HKU16 (WECoV HKU16), sparrow coronavirus HKU17 (SpCoV HKU17), magpie robin coronavirus HKU18 (MRCoV HKU18), night heron coronavirus HKU19 (NHCoV HKU19), wigeon coronavirus HKU20 (WiCoV HKU20), and common moorhen coronavirus HKU21 (CMCoV HKU21) were amplified and sequenced using the RNA extracted from the original swab specimens as templates. The RNA was converted to cDNA by a combined random-priming and oligo(dT)-priming strategy. The cDNA was amplified by degenerate primers designed by multiple alignments of the genomes of other CoVs with complete genomes available, using strategies described in our previous publications (28, 45, 48, 49) and the CoV database CoVDB (20) for sequence retrieval. Additional primers were designed from the results of the first and subsequent rounds of sequencing. These primer sequences are available on request. The 5′ ends of the viral genomes were confirmed by rapid amplification of cDNA ends (RACE) using the 5′/3′ RACE kit (Roche, Germany). Sequences were assembled and manually edited to produce final sequences of the viral genomes.

Genome analysis.

The nucleotide sequences of the genomes and the deduced amino acid sequences of the open reading frames (ORFs) were compared to those of other CoVs using EMBOSS needle (http://www.ebi.ac.uk). Phylogenetic tree construction was performed using the neighbor joining method with ClustalX 1.83. Protein family analysis was performed using PFAM and InterProScan (1, 2). Prediction of transmembrane domains was performed using TMpred and TMHMM (19, 41).

Estimation of divergence dates.

Divergence times for the four genera of CoVs, Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Deltacoronavirus, were calculated using a Bayesian Markov chain Monte Carlo (MCMC) approach as implemented in BEAST (Version 1.6.1) as described previously (11, 23, 27, 47). One parametric model (Constant Size) and one nonparametric model (Bayesian Skyline) tree priors were used for the inference. Analyses were performed under the GTR+I+G substitution model for RdRp gene sequence data and using both a strict and an unrelaxed log-normal-distributed (Ucld) relaxed molecular clock. The MCMC run was 5 × 10 steps long, with sampling every 1,000 steps. Convergence was assessed on the basis of the effective sampling size after a 10% burn-in using Tracer software version 1.5 (11). The mean time of the most recent common ancestor (tMRCA) and the highest posterior density regions at 95% (HPD) (i.e., a credible set that contains 95% of the sampled values) were calculated, and the best-fitting model was selected by a Bayes factor, using marginal likelihoods implemented in Tracer (see Table S1 in the supplemental material) (42). Bayesian Skyline under a relaxed-clock model with Ucld was adopted for making inferences, as Bayes factor analysis indicated that this model fitted the data better than other models tested (see Table S1). The trees were summarized in a target tree by the Tree Annotator program included in the BEAST package by choosing the tree with the maximum sum of posterior probabilities (maximum clade credibility) after a 10% burn-in.

Nucleotide sequence accession numbers.

The nucleotide sequences of the eight genomes of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 have been lodged within the GenBank sequence database under accession no. {"type":"entrez-nucleotide","attrs":{"text":"JQ065042","term_id":"1027948168","term_text":"JQ065042"}}JQ065042 to {"type":"entrez-nucleotide","attrs":{"text":"JQ065049","term_id":"380005522","term_text":"JQ065049"}}JQ065049.

RESULTS

Animal surveillance and identification of seven novel mammalian and avian CoVs.

A total of 7,140 respiratory and alimentary specimens from 3,298 dead wild birds, 221 chickens, and 3,137 mammals were obtained (Table 1). RT-PCR for a 440-bp fragment in the RdRp genes of CoVs was positive in specimens from 17 pigs and 35 dead wild birds. Sequencing results suggested the presence of seven novel CoVs (Fig. 1 and Table 1). These seven novel CoVs were most closely related to our recently described BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13, sharing <66% nucleotide identity with all other known CoVs (Fig. 1). No positive results were obtained from any of the 15 Asian leopard cats, 434 bats, 230 cats, 47 cattle, 221 chickens, 231 dogs, 1,387 humans, 235 monkeys, and 389 rodents tested (Table 1).

Table 1

Animals screened and associated CoVs in the present surveillance study

AnimalSample typeNo. of specimens testedNo. (%) of specimens positive for CoVCoV
Asian leopard catRectal swab and tracheal swab300
BatRectal swab4340
BirdaRectal swab3,30635 (1.1%)WECoV HKU16 (n = 3), SpCoV HKU17 (n = 7), MRCoV HKU18 (n = 1), NHCoV HKU19 (n = 5), WiCoV HKU20 (n = 1), CMCoV HKU21 (n = 1), BuCoV HKU11 (n = 10), ThCoV HKU12 (n = 1), MunCoV HKU13 (n = 6)
CatRectal swab and tracheal swab4600
CattleRectal swab470
ChickenCloacal swab2210
DogRectal swab and tracheal swab4620
HumanNPAb1,3870
MonkeyRectal swab2350
PigRectal swab16917 (10.1%)PorCoV HKU15
RodentRectal swab3890
No. of birds tested for individual species and their associated CoVs are listed in Table S2 in the supplemental material.
NPA, nasopharyngeal aspirate.
An external file that holds a picture, illustration, etc.
Object name is zjv9990958200001.jpg

Phylogenetic analysis of amino acid sequences of the 228-bp fragment (excluding primer sequences) of RNA-dependent RNA polymerase (RdRp) of CoVs identified from dead wild birds and pigs in the present study. The tree was constructed by the neighbor joining method using Kimura correction and bootstrap values calculated from 1,000 trees. The scale bar indicates the estimated number of substitutions per 20 amino acids. The eight genomes completely sequenced are shown in bold. PEDV, porcine epidemic diarrhea virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_003436","term_id":"19387576","term_text":"NC_003436"}}NC_003436); Sc-BatCoV-512, Scotophilus bat coronavirus 512 ({"type":"entrez-nucleotide","attrs":{"text":"NC_009657","term_id":"152994036","term_text":"NC_009657"}}NC_009657); TGEV, transmissible gastroenteritis virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_002306","term_id":"315192962","term_text":"NC_002306"}}NC_002306); FIPV, feline infectious peritonitis virus ({"type":"entrez-nucleotide","attrs":{"text":"AY994055","term_id":"62836705","term_text":"AY994055"}}AY994055); CCoV, canine coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"GQ477367","term_id":"283771347","term_text":"GQ477367"}}GQ477367); PRCV, porcine respiratory coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"DQ811787","term_id":"110746812","term_text":"DQ811787"}}DQ811787); Rh-BatCoV-HKU2, Rhinolophus bat coronavirus HKU2 ({"type":"entrez-nucleotide","attrs":{"text":"EF203064","term_id":"148283139","term_text":"EF203064"}}EF203064); Mi-BatCoV 1A, Miniopterus bat coronavirus 1A ({"type":"entrez-nucleotide","attrs":{"text":"NC_010437","term_id":"169822550","term_text":"NC_010437"}}NC_010437); Mi-BatCoV 1B, Miniopterus bat coronavirus 1B ({"type":"entrez-nucleotide","attrs":{"text":"NC_010436","term_id":"169822542","term_text":"NC_010436"}}NC_010436); Mi-BatCoV-HKU8, Miniopterus bat coronavirus HKU8 ({"type":"entrez-nucleotide","attrs":{"text":"NC_010438","term_id":"169822558","term_text":"NC_010438"}}NC_010438); HCoV-229E, human coronavirus 229E ({"type":"entrez-nucleotide","attrs":{"text":"NC_002645","term_id":"12175745","term_text":"NC_002645"}}NC_002645); HCoV-NL63, human coronavirus NL63 ({"type":"entrez-nucleotide","attrs":{"text":"NC_005831","term_id":"49169782","term_text":"NC_005831"}}NC_005831); HCoV OC43, human coronavirus OC43 ({"type":"entrez-nucleotide","attrs":{"text":"NC_005147","term_id":"38018022","term_text":"NC_005147"}}NC_005147); BCoV, bovine coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_003045","term_id":"15081544","term_text":"NC_003045"}}NC_003045); AntelopeCoV, sable antelope CoV ({"type":"entrez-nucleotide","attrs":{"text":"EF424621","term_id":"145208956","term_text":"EF424621"}}EF424621); GiCoV, giraffe coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"EF424622","term_id":"145208968","term_text":"EF424622"}}EF424622); ECoV, equine coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_010327","term_id":"167600353","term_text":"NC_010327"}}NC_010327); PHEV, porcine hemagglutinating encephalomyelitis virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_007732","term_id":"85718614","term_text":"NC_007732"}}NC_007732); MHV, murine hepatitis virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_001846","term_id":"9629812","term_text":"NC_001846"}}NC_001846); RCoV, rat coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_012936","term_id":"253750530","term_text":"NC_012936"}}NC_012936); HCoV-HKU1, human coronaivurs HKU1 ({"type":"entrez-nucleotide","attrs":{"text":"NC_006577","term_id":"85667876","term_text":"NC_006577"}}NC_006577); Ty-BatCoV-HKU4, Tylonycteris bat coronavirus HKU4 ({"type":"entrez-nucleotide","attrs":{"text":"NC_009019","term_id":"126030112","term_text":"NC_009019"}}NC_009019); Pi-BatCoV-HKU5, Pipistrellus bat coronavirus HKU5 ({"type":"entrez-nucleotide","attrs":{"text":"NC_009020","term_id":"126030122","term_text":"NC_009020"}}NC_009020); SARS CoV, SARS-related human coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_004718","term_id":"30271926","term_text":"NC_004718"}}NC_004718); SARSr-Rh-BatCoV HKU3, SARS-related Rhinolophus bat coronavirus HKU3 ({"type":"entrez-nucleotide","attrs":{"text":"DQ022305","term_id":"76160337","term_text":"DQ022305"}}DQ022305); SARSr CoV CFB, SARS-related Chinese ferret badger coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"AY545919","term_id":"42563758","term_text":"AY545919"}}AY545919); SARSr-CiCoV, SARS-related palm civet coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"AY304488","term_id":"34482139","term_text":"AY304488"}}AY304488); Ro-BatCoV-HKU9, Rousettus bat coronavirus HKU9 ({"type":"entrez-nucleotide","attrs":{"text":"NC_009021","term_id":"126030132","term_text":"NC_009021"}}NC_009021); IBV, infectious bronchitis virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_001451","term_id":"9626535","term_text":"NC_001451"}}NC_001451); IBV-partridge, partridge coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"AY646283","term_id":"50235229","term_text":"AY646283"}}AY646283); TCoV, turkey coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_010800","term_id":"189313868","term_text":"NC_010800"}}NC_010800); IBV-peafowl, peafowl coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"AY641576","term_id":"50082761","term_text":"AY641576"}}AY641576); BWCoV-SW1, beluga whale coronavirus SW1 ({"type":"entrez-nucleotide","attrs":{"text":"NC_010646","term_id":"187251953","term_text":"NC_010646"}}NC_010646); ALCCoV, Asian leopard cat coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"EF584908","term_id":"146411511","term_text":"EF584908"}}EF584908); BuCoV HKU11, bulbul coronavirus HKU11({"type":"entrez-nucleotide","attrs":{"text":"FJ376619","term_id":"212377306","term_text":"FJ376619"}}FJ376619); ThCoV HKU12, thrush coronavirus HKU12 ({"type":"entrez-nucleotide","attrs":{"text":"FJ376621","term_id":"211907050","term_text":"FJ376621"}}FJ376621); MunCoV HKU13, munia coronavirus HKU13 ({"type":"entrez-nucleotide","attrs":{"text":"FJ376622","term_id":"211907060","term_text":"FJ376622"}}FJ376622); PorCoV HKU15, porcine coronavirus HKU15; WECoV HKU16, white-eye coronavirus HKU16; SpCoV HKU17 (TrSp, tree sparrow), sparrow coronavirus HKU17; MRCoV HKU18 (OMR, oriental magpie robin), magpie robin coronavirus HKU18; NHCoV HKU19 (BlCrNH, black-crowned night heron), night heron coronavirus HKU19; WiCoV HKU20 (EuWi, Eurasian wigeon), wigeon coronavirus HKU20; CMCoV HKU21, common moorhen (CM) coronavirus HKU21. Mu, munia; ChMu, chestnut munia; GbTh, gray-backed thrush; ShBu, sooty-headed bulbul; RwBu, red-whiskered bulbul; ChBu, chestnut bulbul.

Genome organization and coding potential of the seven novel mammalian and avian CoVs.

Complete genome sequence data of two strains of PorCoV HKU15 and one complete genome each of WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 were obtained by assembly of the sequences of the RT-PCR products from the RNA extracted from the corresponding individual specimens.

The size of the genomes of the seven novel CoVs ranged from 25,416 bases (PorCoV HKU15) to 26,674 (MRCoV HKU18) and their G+C contents ranged from 35% (CMCoV HKU21) to 47% (MRCoV HKU18) (Table 2). Their genome organizations are similar to those of other CoVs, with the characteristic gene order 5′-replicase ORF1ab, spike (S), envelope (E), membrane (M), nucleocapsid (N)-3′ (Fig. 2 and Table 3). Both 5′ and 3′ ends contain short untranslated regions. The replicase ORF1ab occupies 18.620 to 18.887 kb of the genomes (Table 3). This ORF encodes a number of putative proteins, including nsp3 [which contains the putative papain-like protease (PL)], nsp5 [putative chymotrypsin-like protease (3CL)], nsp12 (putative RdRp), nsp13 (putative helicase), and other proteins of unknown functions. Notably, the amino acids upstream to the putative cleavage sites at nsp2/nsp3, nsp3/nsp4, and nsp4/nsp5 are all AG, AG, and LQ for PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, and CMCoV HKU21; however, those at nsp2/nsp3 are VG and DG, those at nsp3/nsp4 are TG and GG, and those at nsp4/nsp5 are VQ for NHCoV HKU19 and WiCoV HKU20 (see Table S3 in the supplemental material).

Table 2

Comparison of genomic features and amino acid identities among CoVs with complete genome sequences availablea

CoVGenome features
Pairwise amino acid identity (%)
Pairwise amino acid identity (%)
Size (bases)G+C contentPorCoV HKU15
WECoV HKU16
SpCoV HKU17
MRCoV HKU18
NHCoV HKU19
WiCoV HKU20
CMCoV HKU21
3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN
Alphacoronavirus
    PEDV28,0330.4235.848.749.338.023.437.449.347.838.722.436.548.948.939.224.136.848.949.140.822.337.248.847.936.320.639.750.048.537.921.138.149.646.739.421.7
    TGEV28,5860.3834.949.651.635.523.234.549.449.636.124.535.349.851.439.423.534.649.650.739.923.732.950.850.837.122.337.350.250.236.723.733.749.349.436.422.9
    FIPV29,3550.3835.749.751.235.124.535.549.649.336.325.135.749.951.138.525.035.550.150.438.922.232.451.250.736.920.837.950.150.036.522.534.549.749.136.022.1
    CCoV29,3630.3835.649.751.634.923.335.249.449.635.624.035.949.851.438.823.335.349.650.739.323.032.950.850.736.423.637.349.950.736.323.733.349.149.435.223.4
    PRCV27,5500.3734.949.551.640.323.234.549.349.640.523.235.349.751.444.823.534.649.550.744.422.832.950.750.841.422.037.350.350.241.923.733.749.149.440.121.8
    HCoV-229E27,3170.3834.449.350.642.521.635.449.048.342.422.534.249.550.245.523.034.849.349.844.022.534.649.249.039.920.436.349.649.644.121.534.848.647.543.222.8
    HCoV-NL6327,5530.3435.948.849.938.222.138.149.248.140.123.035.649.249.639.322.636.949.249.539.324.634.848.949.036.221.539.150.648.638.823.536.949.547.238.820.6
    Rh-BatCoV-HKU227,1650.3934.450.151.425.020.834.350.049.125.222.334.450.251.126.220.934.449.850.325.921.734.350.549.725.120.935.550.649.126.325.234.150.149.127.322.5
    Mi-BatCoV 1A28,3260.3833.549.051.435.824.435.049.450.135.723.334.249.451.139.425.234.249.051.438.224.532.650.548.735.723.834.749.948.936.422.533.447.949.138.622.4
    Mi-BatCoV 1B28,4760.3934.248.551.135.624.635.448.849.436.122.134.848.850.739.124.933.548.851.138.223.831.949.848.235.722.335.749.648.935.724.132.847.748.338.222.8
    Mi-BatCoV-HKU828,7730.4233.149.349.835.919.436.049.947.536.018.833.449.649.338.920.434.449.349.340.419.834.650.148.137.222.636.750.548.837.622.336.048.547.437.021.6
    Sc-BatCoV-51228,1790.4033.848.649.139.024.836.049.247.538.723.734.148.748.841.325.235.448.348.841.123.434.949.047.836.822.337.348.949.138.123.134.748.747.839.624.3
Betacoronavirus
    Subgroup A
        HCoV-OC4330,7380.3738.151.648.326.022.238.951.548.625.923.237.851.848.326.922.437.551.348.326.421.034.154.548.425.922.538.751.848.425.720.437.851.549.025.624.2
        BCoV31,0280.3738.551.848.425.722.938.851.748.625.821.738.551.848.426.722.837.851.548.426.924.034.454.548.525.623.438.351.748.526.021.538.251.649.025.723.4
        PHEV30,4800.3738.551.748.326.922.138.151.648.626.123.138.551.648.327.222.337.851.448.326.922.234.454.548.526.122.738.751.748.527.121.638.251.649.025.424.1
        AntelopeCoV30,9950.3738.551.848.425.822.938.851.748.525.621.738.551.848.427.022.137.851.448.427.024.034.454.448.526.123.938.351.748.526.221.538.251.649.225.823.4
        GiCoV30,9790.3738.851.848.425.922.938.851.748.525.721.738.851.848.427.222.137.851.448.427.024.034.454.448.525.723.938.351.748.526.521.538.551.649.225.923.4
        ECoV30,9920.3738.551.749.826.023.938.851.649.026.422.638.551.749.926.524.037.851.449.827.523.534.454.648.525.324.638.351.548.526.922.038.251.649.725.624.9
        MHV31,3570.4238.351.948.126.324.339.051.348.526.124.038.351.848.326.525.337.651.948.326.324.235.053.647.525.324.639.650.847.927.124.039.251.248.626.024.6
        HCoV-HKU129,9260.3238.151.249.326.125.238.051.448.226.424.837.951.349.425.726.036.451.448.826.426.036.354.447.425.424.738.150.948.525.822.738.351.247.825.025.4
        RCoV31,2500.4138.751.847.927.224.539.551.448.327.024.338.551.748.125.525.138.251.848.225.825.035.053.647.424.325.239.950.747.727.424.138.351.048.426.423.5
    Subgroup B
        SARS CoV29,7510.4134.550.751.426.126.536.150.350.627.924.734.251.151.625.325.634.250.851.425.426.232.150.550.226.322.734.849.850.326.924.332.950.851.027.324.8
        SARSr-CiCoV29,7280.4134.550.751.426.226.536.150.350.628.024.734.251.151.625.225.634.250.851.425.526.232.150.550.226.222.734.849.850.327.024.332.950.851.027.124.8
        SARSr-Rh-BatCoV HKU329,7040.4134.250.551.426.425.235.850.350.826.224.333.951.151.625.624.933.950.951.426.025.732.150.450.625.623.034.849.750.526.023.532.650.651.327.224.1
        SARSr CoV CFB29,7340.4134.550.651.426.126.536.150.250.628.024.734.251.051.625.525.634.250.751.425.226.232.150.450.225.922.734.849.950.326.824.332.950.751.027.224.8
    Subgroup C
        Ty-BatCoV-HKU430,2860.3836.951.249.826.625.136.651.049.426.124.736.951.549.727.026.235.751.049.727.325.932.751.349.927.324.435.850.949.926.424.435.651.948.926.825.1
        Pi-BatCoV-HKU530,4880.4335.751.150.026.025.637.850.349.025.525.735.451.449.827.225.335.051.149.726.326.233.750.949.626.225.634.651.249.825.324.736.050.949.025.726.1
    Subgroup D
        Ro-Bat-CoV HKU929,1140.4136.451.651.228.425.139.252.650.126.523.136.451.750.926.624.935.851.950.927.725.033.852.250.927.022.835.051.449.227.222.936.952.351.427.823.3
Gammacoronavirus
    IBV27,6080.3843.954.856.630.330.042.654.654.529.930.844.254.956.627.628.943.354.356.228.430.343.653.654.828.729.747.152.455.429.429.241.753.955.129.428.1
    TCoV27,6570.3843.654.957.130.129.243.354.555.330.329.643.955.057.129.529.542.354.457.130.329.244.353.855.330.330.046.252.855.430.329.941.353.655.830.628.2
    BWCoV-SW131,6860.3938.852.952.827.132.139.552.951.628.331.138.852.952.828.531.937.252.352.227.331.641.152.854.227.430.242.052.251.127.230.038.852.152.728.231.5
Deltacoronavirus
    BuCoV HKU1126,4760.3981.188.289.469.874.882.490.996.062.573.280.888.289.643.575.179.588.390.444.571.957.072.376.841.150.658.370.875.443.351.477.584.891.051.860.7
    ThCoV HKU1226,3960.3882.188.289.747.979.783.189.594.747.881.081.888.289.946.779.481.486.889.945.876.757.371.976.443.649.457.771.374.543.649.678.284.490.546.263.3
    MunCoV HKU1326,5520.4382.790.195.871.276.876.587.989.161.374.683.490.196.043.878.894.594.698.046.187.553.172.978.041.453.455.471.775.744.053.272.084.785.452.264.4
    PorCoV HKU1525,4210.4376.988.188.461.975.897.097.899.244.896.884.390.696.144.477.954.072.578.241.852.257.771.074.943.852.473.684.584.650.162.0
    WECoV HKU1626,0270.4076.988.188.461.975.877.288.388.646.476.477.287.389.145.575.455.071.776.642.149.659.671.574.343.450.676.984.890.551.564.3
    SpCoV HKU1726,0670.4597.097.899.244.896.877.288.388.646.476.484.991.096.368.179.052.872.378.447.251.958.071.075.145.853.273.084.784.946.063.5
    MRCoV HKU1826,6740.4784.390.696.144.477.977.287.389.145.575.484.991.096.368.179.054.072.577.746.453.856.471.275.146.353.173.385.184.845.763.9
    NHCoV HKU1926,0640.3854.072.578.241.852.255.071.776.642.149.652.872.378.447.251.954.072.577.746.453.858.369.375.441.054.555.571.977.643.654.5
    WiCoV HKU2026,2110.3957.771.074.943.852.459.671.574.343.450.658.071.075.145.853.256.471.275.146.353.158.369.375.441.054.558.370.876.444.157.0
    CMCoV HKU2126,2160.3573.684.584.650.162.076.984.890.551.564.373.084.784.946.063.573.385.184.845.763.955.571.977.643.654.558.370.876.444.157.0
Comparison of genomic features of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 and other CoVs with complete genome sequences available and of amino acid identities between the predicted 3CL, RNA-dependent RNA (RdRp), helicase (Hel), S, and N proteins of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 and the corresponding proteins of other CoVs. PEDV, porcine epidemic diarrhea virus; TGEV, porcine transmissible gastroenteritis virus; FIPV, feline infectious peritonitis virus; CCoV, canine coronavirus; PRCV, porcine respiratory coronavirus; HCoV-229E, human coronavirus 229E; HCoV-NL63, human coronavirus NL63; Rh-BatCoV-HKU2, Rhinolophus bat coronavirus HKU2; Mi-BatCoV 1A, Miniopterus bat coronavirus 1A; Mi-BatCoV 1B, Miniopterus bat coronavirus 1B; Mi-BatCoV-HKU8, Miniopterus bat coronavirus HKU8; Sc-BatCoV-512, Scotophilus bat coronavirus 512; HCoV OC43, human coronavirus OC43; BCoV, bovine coronavirus; PHEV, porcine hemagglutinating encephalomyelitis virus; AntelopeCoV, sable antelope coronavirus; GiCoV, giraffe coronavirus; ECoV, equine coronavirus; MHV, murine hepatitis virus; HCoV-HKU1, human coronavirus HKU1; RCoV, rat coronavirus; SARS CoV, SARS-related human coronavirus; SARSr-CiCoV, SARS-related palm civet coronavirus; SARSr-Rh-BatCoV HKU3, SARS-related Rhinolophus bat coronavirus HKU3; SARSr CoV CFB, SARS-related Chinese ferret badger coronavirus; Ty-BatCoV-HKU4, Tylonycteris bat coronavirus HKU4; Pi-BatCoV-HKU5, Pipistrellus bat coronavirus HKU5; Ro-BatCoV-HKU9, Rousettus bat coronavirus HKU9; IBV, infectious bronchitis virus; TCoV, turkey coronavirus; BWCoV-SW1, Beluga whale coronavirus SW1; BuCoV HKU11, bulbul coronavirus HKU11; ThCoV HKU12, thrush coronavirus HKU12; MunCoV HKU13, munia coronavirus HKU13.
An external file that holds a picture, illustration, etc.
Object name is zjv9990958200002.jpg

Genome organization of members in Deltacoronavirus. ORFs downstream of S gene are magnified to show the differences among the genomes of the 10 CoVs. Papain-like protease (PL), chymotrypsin-like protease (3CL), and RNA-dependent RNA polymerase (RdRp) are represented by orange boxes. Spike (S), envelope (E), membrane (M), and nucleocapsid (N) are represented by green boxes. Putative accessory proteins are represented by blue boxes. The seven CoVs discovered in this study are shown in bold.

Table 3

Coding potential and putative transcription regulatory sequences of CoV genomesa

CoVORFLocation (nt)Length (nt)Length (aa)FramePutative TRS
TRS location (nt)TRS sequence(s) (distance in bases to AUG)b
PorCoV HKU151ab540–1934218,8036,268+3, +275ACACCA(459)AUG
S19324–228063,4831,161+119178ACACCA(145)AUG
E22800–2305125284+322777ACACCG(17)AUG
M23044–23697654218+123018ACACCA(20)AUG
NS623697–2398128595+323645ACACCA(46)AUG
N24002–250301,029343+223989ACACCA(7)AUG
NS724096–24698603201+324008GCACCA(82)AUG
WECoV HKU161ab511–1939718,8876,296+1, +366ACACCA(439)AUG
S19379–229183,5401,180+219233ACACCA(140)AUG
E22912–2316024983+122886ACACCA(20)AUG
M23153–23809657219+223130ACACCA(17)AUG
NS623809–2409028294+123768ACAUCA(35)AUG
N24115–251581,044348+124101ACACCA(8)AUG
NS7a24143–24811669223+224101ACACCA(36)AUG
NS7b25139–2527013244+225039AAACCA(94)AUG
SpCoV HKU171ab520–1935218,8336,278+1, +357ACACCA(452)AUG
S19334–229543,6211,207+219188ACACCA(140)AUG
E22948–2319624983+122925ACACCG(17)AUG
M23189–23842654218+223166ACACCA(17)AUG
NS623842–2412928896+123790ACACCA(46)AUG
N24150–251781,029343+324137ACACCA(7)AUG
NS7a25189–25623435145+125179ACACCA(4)AUG
NS7b25539–2575121371+325523ACUCCA(10)AUG
MRCoV HKU181ab596–1935618,7616,254+2, +164ACACCA(526)AUG
S19338–229913,6541,218+319192ACACCA(140)AUG
E22985–2323324983+222945ACACCG(34)AUG
M23226–23882657219+323203ACACCA(17)AUG
NS623882–2417229197+223857ACGCCA(19)AUG
N24355–253951,041347+124340ACACCA(9)AUG
NS7a25407–2558017458+325396ACACCA(5)AUG
NS7b25561–25932372124+1
NS7c25941–2619525585+325910ACACCA(25)AUG
NHCoV HKU191ab482–1932318,8426,281+2, +167ACACCG(409)AUG
S19305–230693,7651,255+319156ACACCG(143)AUG
E23069–2331724983+223013ACACCA(50)AUG
M23310–23960651217+323211ACACCG(93)AUG
NS623960–2423827993+223951ACACCU(3)AUG
N24248–252761,029343+224231ACACCU(8)AUG
NS7a25277–2557329799+225248ACACCG(23)AUG
NS7b25583–2587629498+225560ACACCA(17)AUG
WiCoV HKU201ab219–1883818,6206,207+3, +260ACACCA(153)AUG
S18817–224553,6391,213+118731ACACCU(80)AUG
E22455–2271526187+322380ACACCA(69)AUG
M22708–23358651217+122597ACACCG(105)AUG
NS623358–2363027391+3
N23646–246981,053351+323631ACACCA(9)AUG
NS7a24695–2492823478+224609AAACCA(80)AUG
NS7b25218–2546624983+325177ACACCG(35)AUG
NS7c25450–2571626789+125444ACACCGAUG
NS7d25752–2595220167+325735AAACCU(11)AUG
CMCoV HKU211ab478–1910318,6266,209+1, +363ACACCA(409)AUG
S19085–227293,6451,215+218939ACACCA(140)AUG
E22723–2297124983+122697ACACCA(20)AUG
M22973–23779807269+222938ACACCA(29)AUG
NS623779–2402424682+123727ACACCA(46)AUG
N24052–251071,056352+124039ACACCG(7)AUG
NS7a25107–2537927391+325036ACACCU(65)AUG
NS7b25391–2557618662+225379ACACCU(6)AUG
NS7c25500–25916417139+2
PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21. aa, amino acid; nt, nucleotide.
Boldface indicates putative TRS sequences. The nucleotide variations are in italic.

The seven novel CoVs display similar genome organizations and differ only in the number of ORFs downstream of N (Fig. 2). Their transcription regulatory sequences (TRSs) conform to the consensus motif 5′-ACACCA-3′ (Table 3), which appears to be unique to members of the genus Deltacoronavirus. Interestingly, similar to BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13, the perfect TRSs of S in the genomes of the seven novel CoVs were separated from the corresponding AUG by 80 to 145 bases (Table 3). This is in contrast to the relatively small number of bases between the TRSs for S and the corresponding AUG (range: from 0 bases in HCoV-NL63, Rhinolophus bat coronavirus HKU2 [Rh-BatCoV-HKU2], HCoV-HKU1, bovine coronavirus [BCoV], HCoV-OC43, mouse hepatitis virus [MHV], porcine hemagglutinating encephalomyelitis virus, SARS-CoV, and SARS-related Rhinolophus bat coronavirus HKU3 [SARSr-Rh-batCoV HKU3] to 52 bases in infectious bronchitis virus [IBV]) in members of Alphacoronavirus, Betacoronavirus, and Gammacoronavirus. Similar to BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13, the genomes of the seven novel CoVs have putative PL, which are homologous to PL2 of Alphacoronavirus and Betacoronavirus subgroup A and PL of Betacoronavirus subgroups B, C, and D and Gammacoronavirus (Fig. 2). Similar to BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13, one ORF (NS6) is found between M and N of the genomes of the seven novel CoVs. On the other hand, one ORF (NS7) is present overlapping with N in PorCoV HKU15, two ORFs (NS7a and 7b) are present overlapping or downstream of N in WECoV HKU16, SpCoV HKU17, and NHCoV HKU19, three ORFs (NS7a, 7b, and 7c) are present downstream of N in MRCoV HKU18 and CMCoV HKU21, and four ORFs (NS7a, 7b, 7c, and 7d) are present overlapping or downstream of N in WiCoV HKU20. For NS7 of PorCoV, the presence of an imperfect TRS (GCACCA) and its relatively high Ka/Ks ratio (number of nonsynonymous substitutions per nonsynonymous site/number of synonymous substitutions per synonymous site) of 1.046 (data not shown) implied that this ORF may not be expressed. BLAST search revealed no amino acid similarities between these putative nonstructural proteins and other known proteins, and no functional domain was identified by PFAM and InterProScan, except that NS7a of NHCoV HKU19 was found to be homologous to the NS7a of BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13. NS7b of WiCoV HKU20 and CMCoV HKU21, and NS7d of WiCoV HKU20, were also found to be homologous to the NS3b of IBV and hypothetical protein of goose coronavirus, respectively. Transmembrane helices, predicted by TMHMM and TMpred, in putative accessory proteins downstream to the N genes in the genomes of SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 are listed in Table S4 in the supplemental material. Each of the genomes of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, and CMCoV HKU21 contains a stem-loop II motif (s2m) (residues 25,220 to 25,251, 25,825 to 25,856, 25,865 to 25,896, 26,472 to 26,503, and 26,013 to 26,044, respectively), a conserved RNA element downstream of N and upstream of the poly(A) tail, similar to those in IBV, TCoV, SARSr-Rh-BatCoV, and SARS-CoV, as well as other CoVs discovered in Asian leopard cat, graylag geese, feral pigeons, and mallards, for which complete genomes are not available (Fig. 3) (14, 21, 38).

An external file that holds a picture, illustration, etc.
Object name is zjv9990958200003.jpg

Multiple alignments of conserved s2m of infectious bronchitis virus (IBV), SARS-related human coronavirus (SARS CoV), SARS-related Rhinolophus bat coronavirus HKU3 (SARSr-Rh-BatCoV HKU3), BuCoV HKU11, ThCoV HKU12, MunCoV HKU13, Asian leopard cat coronavirus (ALCCoV), PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, and CMCoV HKU21. Identical nucleotides are marked by asterisks. Acc. No., accession no.

Comparison of the amino acid identities of the seven conserved replicase domains for species demarcation (ADRP, nsp5 [3CL], nsp12 [RdRp], nsp13 [Hel], nsp14 [ExoN], nsp15 [NendoU], and nsp16 [O-MT]) (8) among the 10 deltacoronaviruses is shown in Table S5 in the supplemental material. In all the seven domains, the amino acid sequences of PorCoV HKU15 and SpCoV HKU17 showed more than 90% identity, indicating that these two coronaviruses should be subspecies of the same species.

Phylogenetic analyses.

The phylogenetic trees constructed using the nucleotide sequences of the 3CL, RdRp, Hel, S, and N of the seven novel CoVs and other CoVs are shown in Fig. 4 and the corresponding pairwise amino acid identities are shown in Table 2. For all five genes, the seven novel CoVs possessed higher amino acid identities to each other and BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13 than to any other known CoVs with complete genomes available (Table 2). In all five trees, the seven novel CoVs were clustered with BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13 (Fig. 4). For Hel, S, and N, PorCoVs were also clustered with a CoV found in Asian leopard cat (10), for which the sequences of these genes were available (Fig. 4). There were <2% base differences between the Hel, S, and N genes of PorCoV and those of the Asian leopard cat coronavirus. Based on both phylogenetic tree analyses and amino acid differences, the seven novel CoVs as well as BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13 should belong to the same genus, Deltacoronavirus.

An external file that holds a picture, illustration, etc.
Object name is zjv999095820004a.jpg
An external file that holds a picture, illustration, etc.
Object name is zjv999095820004b.jpg
An external file that holds a picture, illustration, etc.
Object name is zjv999095820004c.jpg

Phylogenetic analyses of 3CL, RdRp, helicase (Hel), S, and N proteins of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21. The trees were constructed by using the neighbor joining method using Kimura correction and bootstrap values calculated from 1,000 trees. Two hundred ninety-five, 892, 590, 802, and 249 amino acid positions in 3CL, RdRp, Hel, S, and N, respectively, were included in the analyses. The trees were midpoint rooted. For 3CL and S, the scale bar indicates the estimated number of substitutions per 10 amino acids. For RdRp and Hel, the scale bar indicates the estimated number of substitutions per 20 amino acids. For N, the scale bar indicates the estimated number of substitutions per 5 amino acids. Viruses characterized in this study are in bold. Virus name abbreviations are the same as those in the Fig. 1 legend.

Estimation of divergence dates.

Using the Bayesian Skyline under a relaxed-clock model with an uncorrelated log-normal distribution, the mean evolutionary rate of CoVs was estimated at 1.3 × 10 nucleotide substitutions per site per year for the RdRp gene. Molecular clock analysis using the RdRp gene showed that the tMRCA of all CoVs was estimated at ∼8100 BC (HPDs, 20607 to 974 BC), that of Alphacoronavirus at ∼2400 BC (HPDs, 7659 to 722 BC), that of Betacoronavirus at ∼3300 BC (HPDs, 9713 to 447 BC), that of Gammacoronavirus at ∼2800 BC (HPDs, 8840 to 700 BC), and that of Deltacoronavirus at ∼3000 BC (HPDs, 9073 to 555 BC) (Fig. 5).

An external file that holds a picture, illustration, etc.
Object name is zjv9990958200005.jpg

Estimation of the time to the most recent common ancestor for Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Deltacoronavirus. The time-scaled phylogeny was summarized from all MCMC phylogenies of the RdRp gene data set analyzed under the relaxed-clock model with an uncorrelated log-normal distribution in BEAST version 1.6.1. Viruses characterized in this study are in bold. The numbers indicate number of years ago. This is shown in the scale bar. Virus name abbreviations are the same as those in the legends of Fig. 1.

Animal surveillance and identification of seven novel mammalian and avian CoVs.

A total of 7,140 respiratory and alimentary specimens from 3,298 dead wild birds, 221 chickens, and 3,137 mammals were obtained (Table 1). RT-PCR for a 440-bp fragment in the RdRp genes of CoVs was positive in specimens from 17 pigs and 35 dead wild birds. Sequencing results suggested the presence of seven novel CoVs (Fig. 1 and Table 1). These seven novel CoVs were most closely related to our recently described BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13, sharing <66% nucleotide identity with all other known CoVs (Fig. 1). No positive results were obtained from any of the 15 Asian leopard cats, 434 bats, 230 cats, 47 cattle, 221 chickens, 231 dogs, 1,387 humans, 235 monkeys, and 389 rodents tested (Table 1).

Table 1

Animals screened and associated CoVs in the present surveillance study

AnimalSample typeNo. of specimens testedNo. (%) of specimens positive for CoVCoV
Asian leopard catRectal swab and tracheal swab300
BatRectal swab4340
BirdaRectal swab3,30635 (1.1%)WECoV HKU16 (n = 3), SpCoV HKU17 (n = 7), MRCoV HKU18 (n = 1), NHCoV HKU19 (n = 5), WiCoV HKU20 (n = 1), CMCoV HKU21 (n = 1), BuCoV HKU11 (n = 10), ThCoV HKU12 (n = 1), MunCoV HKU13 (n = 6)
CatRectal swab and tracheal swab4600
CattleRectal swab470
ChickenCloacal swab2210
DogRectal swab and tracheal swab4620
HumanNPAb1,3870
MonkeyRectal swab2350
PigRectal swab16917 (10.1%)PorCoV HKU15
RodentRectal swab3890
No. of birds tested for individual species and their associated CoVs are listed in Table S2 in the supplemental material.
NPA, nasopharyngeal aspirate.
An external file that holds a picture, illustration, etc.
Object name is zjv9990958200001.jpg

Phylogenetic analysis of amino acid sequences of the 228-bp fragment (excluding primer sequences) of RNA-dependent RNA polymerase (RdRp) of CoVs identified from dead wild birds and pigs in the present study. The tree was constructed by the neighbor joining method using Kimura correction and bootstrap values calculated from 1,000 trees. The scale bar indicates the estimated number of substitutions per 20 amino acids. The eight genomes completely sequenced are shown in bold. PEDV, porcine epidemic diarrhea virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_003436","term_id":"19387576","term_text":"NC_003436"}}NC_003436); Sc-BatCoV-512, Scotophilus bat coronavirus 512 ({"type":"entrez-nucleotide","attrs":{"text":"NC_009657","term_id":"152994036","term_text":"NC_009657"}}NC_009657); TGEV, transmissible gastroenteritis virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_002306","term_id":"315192962","term_text":"NC_002306"}}NC_002306); FIPV, feline infectious peritonitis virus ({"type":"entrez-nucleotide","attrs":{"text":"AY994055","term_id":"62836705","term_text":"AY994055"}}AY994055); CCoV, canine coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"GQ477367","term_id":"283771347","term_text":"GQ477367"}}GQ477367); PRCV, porcine respiratory coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"DQ811787","term_id":"110746812","term_text":"DQ811787"}}DQ811787); Rh-BatCoV-HKU2, Rhinolophus bat coronavirus HKU2 ({"type":"entrez-nucleotide","attrs":{"text":"EF203064","term_id":"148283139","term_text":"EF203064"}}EF203064); Mi-BatCoV 1A, Miniopterus bat coronavirus 1A ({"type":"entrez-nucleotide","attrs":{"text":"NC_010437","term_id":"169822550","term_text":"NC_010437"}}NC_010437); Mi-BatCoV 1B, Miniopterus bat coronavirus 1B ({"type":"entrez-nucleotide","attrs":{"text":"NC_010436","term_id":"169822542","term_text":"NC_010436"}}NC_010436); Mi-BatCoV-HKU8, Miniopterus bat coronavirus HKU8 ({"type":"entrez-nucleotide","attrs":{"text":"NC_010438","term_id":"169822558","term_text":"NC_010438"}}NC_010438); HCoV-229E, human coronavirus 229E ({"type":"entrez-nucleotide","attrs":{"text":"NC_002645","term_id":"12175745","term_text":"NC_002645"}}NC_002645); HCoV-NL63, human coronavirus NL63 ({"type":"entrez-nucleotide","attrs":{"text":"NC_005831","term_id":"49169782","term_text":"NC_005831"}}NC_005831); HCoV OC43, human coronavirus OC43 ({"type":"entrez-nucleotide","attrs":{"text":"NC_005147","term_id":"38018022","term_text":"NC_005147"}}NC_005147); BCoV, bovine coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_003045","term_id":"15081544","term_text":"NC_003045"}}NC_003045); AntelopeCoV, sable antelope CoV ({"type":"entrez-nucleotide","attrs":{"text":"EF424621","term_id":"145208956","term_text":"EF424621"}}EF424621); GiCoV, giraffe coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"EF424622","term_id":"145208968","term_text":"EF424622"}}EF424622); ECoV, equine coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_010327","term_id":"167600353","term_text":"NC_010327"}}NC_010327); PHEV, porcine hemagglutinating encephalomyelitis virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_007732","term_id":"85718614","term_text":"NC_007732"}}NC_007732); MHV, murine hepatitis virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_001846","term_id":"9629812","term_text":"NC_001846"}}NC_001846); RCoV, rat coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_012936","term_id":"253750530","term_text":"NC_012936"}}NC_012936); HCoV-HKU1, human coronaivurs HKU1 ({"type":"entrez-nucleotide","attrs":{"text":"NC_006577","term_id":"85667876","term_text":"NC_006577"}}NC_006577); Ty-BatCoV-HKU4, Tylonycteris bat coronavirus HKU4 ({"type":"entrez-nucleotide","attrs":{"text":"NC_009019","term_id":"126030112","term_text":"NC_009019"}}NC_009019); Pi-BatCoV-HKU5, Pipistrellus bat coronavirus HKU5 ({"type":"entrez-nucleotide","attrs":{"text":"NC_009020","term_id":"126030122","term_text":"NC_009020"}}NC_009020); SARS CoV, SARS-related human coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_004718","term_id":"30271926","term_text":"NC_004718"}}NC_004718); SARSr-Rh-BatCoV HKU3, SARS-related Rhinolophus bat coronavirus HKU3 ({"type":"entrez-nucleotide","attrs":{"text":"DQ022305","term_id":"76160337","term_text":"DQ022305"}}DQ022305); SARSr CoV CFB, SARS-related Chinese ferret badger coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"AY545919","term_id":"42563758","term_text":"AY545919"}}AY545919); SARSr-CiCoV, SARS-related palm civet coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"AY304488","term_id":"34482139","term_text":"AY304488"}}AY304488); Ro-BatCoV-HKU9, Rousettus bat coronavirus HKU9 ({"type":"entrez-nucleotide","attrs":{"text":"NC_009021","term_id":"126030132","term_text":"NC_009021"}}NC_009021); IBV, infectious bronchitis virus ({"type":"entrez-nucleotide","attrs":{"text":"NC_001451","term_id":"9626535","term_text":"NC_001451"}}NC_001451); IBV-partridge, partridge coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"AY646283","term_id":"50235229","term_text":"AY646283"}}AY646283); TCoV, turkey coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"NC_010800","term_id":"189313868","term_text":"NC_010800"}}NC_010800); IBV-peafowl, peafowl coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"AY641576","term_id":"50082761","term_text":"AY641576"}}AY641576); BWCoV-SW1, beluga whale coronavirus SW1 ({"type":"entrez-nucleotide","attrs":{"text":"NC_010646","term_id":"187251953","term_text":"NC_010646"}}NC_010646); ALCCoV, Asian leopard cat coronavirus ({"type":"entrez-nucleotide","attrs":{"text":"EF584908","term_id":"146411511","term_text":"EF584908"}}EF584908); BuCoV HKU11, bulbul coronavirus HKU11({"type":"entrez-nucleotide","attrs":{"text":"FJ376619","term_id":"212377306","term_text":"FJ376619"}}FJ376619); ThCoV HKU12, thrush coronavirus HKU12 ({"type":"entrez-nucleotide","attrs":{"text":"FJ376621","term_id":"211907050","term_text":"FJ376621"}}FJ376621); MunCoV HKU13, munia coronavirus HKU13 ({"type":"entrez-nucleotide","attrs":{"text":"FJ376622","term_id":"211907060","term_text":"FJ376622"}}FJ376622); PorCoV HKU15, porcine coronavirus HKU15; WECoV HKU16, white-eye coronavirus HKU16; SpCoV HKU17 (TrSp, tree sparrow), sparrow coronavirus HKU17; MRCoV HKU18 (OMR, oriental magpie robin), magpie robin coronavirus HKU18; NHCoV HKU19 (BlCrNH, black-crowned night heron), night heron coronavirus HKU19; WiCoV HKU20 (EuWi, Eurasian wigeon), wigeon coronavirus HKU20; CMCoV HKU21, common moorhen (CM) coronavirus HKU21. Mu, munia; ChMu, chestnut munia; GbTh, gray-backed thrush; ShBu, sooty-headed bulbul; RwBu, red-whiskered bulbul; ChBu, chestnut bulbul.

Genome organization and coding potential of the seven novel mammalian and avian CoVs.

Complete genome sequence data of two strains of PorCoV HKU15 and one complete genome each of WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 were obtained by assembly of the sequences of the RT-PCR products from the RNA extracted from the corresponding individual specimens.

The size of the genomes of the seven novel CoVs ranged from 25,416 bases (PorCoV HKU15) to 26,674 (MRCoV HKU18) and their G+C contents ranged from 35% (CMCoV HKU21) to 47% (MRCoV HKU18) (Table 2). Their genome organizations are similar to those of other CoVs, with the characteristic gene order 5′-replicase ORF1ab, spike (S), envelope (E), membrane (M), nucleocapsid (N)-3′ (Fig. 2 and Table 3). Both 5′ and 3′ ends contain short untranslated regions. The replicase ORF1ab occupies 18.620 to 18.887 kb of the genomes (Table 3). This ORF encodes a number of putative proteins, including nsp3 [which contains the putative papain-like protease (PL)], nsp5 [putative chymotrypsin-like protease (3CL)], nsp12 (putative RdRp), nsp13 (putative helicase), and other proteins of unknown functions. Notably, the amino acids upstream to the putative cleavage sites at nsp2/nsp3, nsp3/nsp4, and nsp4/nsp5 are all AG, AG, and LQ for PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, and CMCoV HKU21; however, those at nsp2/nsp3 are VG and DG, those at nsp3/nsp4 are TG and GG, and those at nsp4/nsp5 are VQ for NHCoV HKU19 and WiCoV HKU20 (see Table S3 in the supplemental material).

Table 2

Comparison of genomic features and amino acid identities among CoVs with complete genome sequences availablea

CoVGenome features
Pairwise amino acid identity (%)
Pairwise amino acid identity (%)
Size (bases)G+C contentPorCoV HKU15
WECoV HKU16
SpCoV HKU17
MRCoV HKU18
NHCoV HKU19
WiCoV HKU20
CMCoV HKU21
3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN3CLproRdRpHelSN
Alphacoronavirus
    PEDV28,0330.4235.848.749.338.023.437.449.347.838.722.436.548.948.939.224.136.848.949.140.822.337.248.847.936.320.639.750.048.537.921.138.149.646.739.421.7
    TGEV28,5860.3834.949.651.635.523.234.549.449.636.124.535.349.851.439.423.534.649.650.739.923.732.950.850.837.122.337.350.250.236.723.733.749.349.436.422.9
    FIPV29,3550.3835.749.751.235.124.535.549.649.336.325.135.749.951.138.525.035.550.150.438.922.232.451.250.736.920.837.950.150.036.522.534.549.749.136.022.1
    CCoV29,3630.3835.649.751.634.923.335.249.449.635.624.035.949.851.438.823.335.349.650.739.323.032.950.850.736.423.637.349.950.736.323.733.349.149.435.223.4
    PRCV27,5500.3734.949.551.640.323.234.549.349.640.523.235.349.751.444.823.534.649.550.744.422.832.950.750.841.422.037.350.350.241.923.733.749.149.440.121.8
    HCoV-229E27,3170.3834.449.350.642.521.635.449.048.342.422.534.249.550.245.523.034.849.349.844.022.534.649.249.039.920.436.349.649.644.121.534.848.647.543.222.8
    HCoV-NL6327,5530.3435.948.849.938.222.138.149.248.140.123.035.649.249.639.322.636.949.249.539.324.634.848.949.036.221.539.150.648.638.823.536.949.547.238.820.6
    Rh-BatCoV-HKU227,1650.3934.450.151.425.020.834.350.049.125.222.334.450.251.126.220.934.449.850.325.921.734.350.549.725.120.935.550.649.126.325.234.150.149.127.322.5
    Mi-BatCoV 1A28,3260.3833.549.051.435.824.435.049.450.135.723.334.249.451.139.425.234.249.051.438.224.532.650.548.735.723.834.749.948.936.422.533.447.949.138.622.4
    Mi-BatCoV 1B28,4760.3934.248.551.135.624.635.448.849.436.122.134.848.850.739.124.933.548.851.138.223.831.949.848.235.722.335.749.648.935.724.132.847.748.338.222.8
    Mi-BatCoV-HKU828,7730.4233.149.349.835.919.436.049.947.536.018.833.449.649.338.920.434.449.349.340.419.834.650.148.137.222.636.750.548.837.622.336.048.547.437.021.6
    Sc-BatCoV-51228,1790.4033.848.649.139.024.836.049.247.538.723.734.148.748.841.325.235.448.348.841.123.434.949.047.836.822.337.348.949.138.123.134.748.747.839.624.3
Betacoronavirus
    Subgroup A
        HCoV-OC4330,7380.3738.151.648.326.022.238.951.548.625.923.237.851.848.326.922.437.551.348.326.421.034.154.548.425.922.538.751.848.425.720.437.851.549.025.624.2
        BCoV31,0280.3738.551.848.425.722.938.851.748.625.821.738.551.848.426.722.837.851.548.426.924.034.454.548.525.623.438.351.748.526.021.538.251.649.025.723.4
        PHEV30,4800.3738.551.748.326.922.138.151.648.626.123.138.551.648.327.222.337.851.448.326.922.234.454.548.526.122.738.751.748.527.121.638.251.649.025.424.1
        AntelopeCoV30,9950.3738.551.848.425.822.938.851.748.525.621.738.551.848.427.022.137.851.448.427.024.034.454.448.526.123.938.351.748.526.221.538.251.649.225.823.4
        GiCoV30,9790.3738.851.848.425.922.938.851.748.525.721.738.851.848.427.222.137.851.448.427.024.034.454.448.525.723.938.351.748.526.521.538.551.649.225.923.4
        ECoV30,9920.3738.551.749.826.023.938.851.649.026.422.638.551.749.926.524.037.851.449.827.523.534.454.648.525.324.638.351.548.526.922.038.251.649.725.624.9
        MHV31,3570.4238.351.948.126.324.339.051.348.526.124.038.351.848.326.525.337.651.948.326.324.235.053.647.525.324.639.650.847.927.124.039.251.248.626.024.6
        HCoV-HKU129,9260.3238.151.249.326.125.238.051.448.226.424.837.951.349.425.726.036.451.448.826.426.036.354.447.425.424.738.150.948.525.822.738.351.247.825.025.4
        RCoV31,2500.4138.751.847.927.224.539.551.448.327.024.338.551.748.125.525.138.251.848.225.825.035.053.647.424.325.239.950.747.727.424.138.351.048.426.423.5
    Subgroup B
        SARS CoV29,7510.4134.550.751.426.126.536.150.350.627.924.734.251.151.625.325.634.250.851.425.426.232.150.550.226.322.734.849.850.326.924.332.950.851.027.324.8
        SARSr-CiCoV29,7280.4134.550.751.426.226.536.150.350.628.024.734.251.151.625.225.634.250.851.425.526.232.150.550.226.222.734.849.850.327.024.332.950.851.027.124.8
        SARSr-Rh-BatCoV HKU329,7040.4134.250.551.426.425.235.850.350.826.224.333.951.151.625.624.933.950.951.426.025.732.150.450.625.623.034.849.750.526.023.532.650.651.327.224.1
        SARSr CoV CFB29,7340.4134.550.651.426.126.536.150.250.628.024.734.251.051.625.525.634.250.751.425.226.232.150.450.225.922.734.849.950.326.824.332.950.751.027.224.8
    Subgroup C
        Ty-BatCoV-HKU430,2860.3836.951.249.826.625.136.651.049.426.124.736.951.549.727.026.235.751.049.727.325.932.751.349.927.324.435.850.949.926.424.435.651.948.926.825.1
        Pi-BatCoV-HKU530,4880.4335.751.150.026.025.637.850.349.025.525.735.451.449.827.225.335.051.149.726.326.233.750.949.626.225.634.651.249.825.324.736.050.949.025.726.1
    Subgroup D
        Ro-Bat-CoV HKU929,1140.4136.451.651.228.425.139.252.650.126.523.136.451.750.926.624.935.851.950.927.725.033.852.250.927.022.835.051.449.227.222.936.952.351.427.823.3
Gammacoronavirus
    IBV27,6080.3843.954.856.630.330.042.654.654.529.930.844.254.956.627.628.943.354.356.228.430.343.653.654.828.729.747.152.455.429.429.241.753.955.129.428.1
    TCoV27,6570.3843.654.957.130.129.243.354.555.330.329.643.955.057.129.529.542.354.457.130.329.244.353.855.330.330.046.252.855.430.329.941.353.655.830.628.2
    BWCoV-SW131,6860.3938.852.952.827.132.139.552.951.628.331.138.852.952.828.531.937.252.352.227.331.641.152.854.227.430.242.052.251.127.230.038.852.152.728.231.5
Deltacoronavirus
    BuCoV HKU1126,4760.3981.188.289.469.874.882.490.996.062.573.280.888.289.643.575.179.588.390.444.571.957.072.376.841.150.658.370.875.443.351.477.584.891.051.860.7
    ThCoV HKU1226,3960.3882.188.289.747.979.783.189.594.747.881.081.888.289.946.779.481.486.889.945.876.757.371.976.443.649.457.771.374.543.649.678.284.490.546.263.3
    MunCoV HKU1326,5520.4382.790.195.871.276.876.587.989.161.374.683.490.196.043.878.894.594.698.046.187.553.172.978.041.453.455.471.775.744.053.272.084.785.452.264.4
    PorCoV HKU1525,4210.4376.988.188.461.975.897.097.899.244.896.884.390.696.144.477.954.072.578.241.852.257.771.074.943.852.473.684.584.650.162.0
    WECoV HKU1626,0270.4076.988.188.461.975.877.288.388.646.476.477.287.389.145.575.455.071.776.642.149.659.671.574.343.450.676.984.890.551.564.3
    SpCoV HKU1726,0670.4597.097.899.244.896.877.288.388.646.476.484.991.096.368.179.052.872.378.447.251.958.071.075.145.853.273.084.784.946.063.5
    MRCoV HKU1826,6740.4784.390.696.144.477.977.287.389.145.575.484.991.096.368.179.054.072.577.746.453.856.471.275.146.353.173.385.184.845.763.9
    NHCoV HKU1926,0640.3854.072.578.241.852.255.071.776.642.149.652.872.378.447.251.954.072.577.746.453.858.369.375.441.054.555.571.977.643.654.5
    WiCoV HKU2026,2110.3957.771.074.943.852.459.671.574.343.450.658.071.075.145.853.256.471.275.146.353.158.369.375.441.054.558.370.876.444.157.0
    CMCoV HKU2126,2160.3573.684.584.650.162.076.984.890.551.564.373.084.784.946.063.573.385.184.845.763.955.571.977.643.654.558.370.876.444.157.0
Comparison of genomic features of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 and other CoVs with complete genome sequences available and of amino acid identities between the predicted 3CL, RNA-dependent RNA (RdRp), helicase (Hel), S, and N proteins of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 and the corresponding proteins of other CoVs. PEDV, porcine epidemic diarrhea virus; TGEV, porcine transmissible gastroenteritis virus; FIPV, feline infectious peritonitis virus; CCoV, canine coronavirus; PRCV, porcine respiratory coronavirus; HCoV-229E, human coronavirus 229E; HCoV-NL63, human coronavirus NL63; Rh-BatCoV-HKU2, Rhinolophus bat coronavirus HKU2; Mi-BatCoV 1A, Miniopterus bat coronavirus 1A; Mi-BatCoV 1B, Miniopterus bat coronavirus 1B; Mi-BatCoV-HKU8, Miniopterus bat coronavirus HKU8; Sc-BatCoV-512, Scotophilus bat coronavirus 512; HCoV OC43, human coronavirus OC43; BCoV, bovine coronavirus; PHEV, porcine hemagglutinating encephalomyelitis virus; AntelopeCoV, sable antelope coronavirus; GiCoV, giraffe coronavirus; ECoV, equine coronavirus; MHV, murine hepatitis virus; HCoV-HKU1, human coronavirus HKU1; RCoV, rat coronavirus; SARS CoV, SARS-related human coronavirus; SARSr-CiCoV, SARS-related palm civet coronavirus; SARSr-Rh-BatCoV HKU3, SARS-related Rhinolophus bat coronavirus HKU3; SARSr CoV CFB, SARS-related Chinese ferret badger coronavirus; Ty-BatCoV-HKU4, Tylonycteris bat coronavirus HKU4; Pi-BatCoV-HKU5, Pipistrellus bat coronavirus HKU5; Ro-BatCoV-HKU9, Rousettus bat coronavirus HKU9; IBV, infectious bronchitis virus; TCoV, turkey coronavirus; BWCoV-SW1, Beluga whale coronavirus SW1; BuCoV HKU11, bulbul coronavirus HKU11; ThCoV HKU12, thrush coronavirus HKU12; MunCoV HKU13, munia coronavirus HKU13.
An external file that holds a picture, illustration, etc.
Object name is zjv9990958200002.jpg

Genome organization of members in Deltacoronavirus. ORFs downstream of S gene are magnified to show the differences among the genomes of the 10 CoVs. Papain-like protease (PL), chymotrypsin-like protease (3CL), and RNA-dependent RNA polymerase (RdRp) are represented by orange boxes. Spike (S), envelope (E), membrane (M), and nucleocapsid (N) are represented by green boxes. Putative accessory proteins are represented by blue boxes. The seven CoVs discovered in this study are shown in bold.

Table 3

Coding potential and putative transcription regulatory sequences of CoV genomesa

CoVORFLocation (nt)Length (nt)Length (aa)FramePutative TRS
TRS location (nt)TRS sequence(s) (distance in bases to AUG)b
PorCoV HKU151ab540–1934218,8036,268+3, +275ACACCA(459)AUG
S19324–228063,4831,161+119178ACACCA(145)AUG
E22800–2305125284+322777ACACCG(17)AUG
M23044–23697654218+123018ACACCA(20)AUG
NS623697–2398128595+323645ACACCA(46)AUG
N24002–250301,029343+223989ACACCA(7)AUG
NS724096–24698603201+324008GCACCA(82)AUG
WECoV HKU161ab511–1939718,8876,296+1, +366ACACCA(439)AUG
S19379–229183,5401,180+219233ACACCA(140)AUG
E22912–2316024983+122886ACACCA(20)AUG
M23153–23809657219+223130ACACCA(17)AUG
NS623809–2409028294+123768ACAUCA(35)AUG
N24115–251581,044348+124101ACACCA(8)AUG
NS7a24143–24811669223+224101ACACCA(36)AUG
NS7b25139–2527013244+225039AAACCA(94)AUG
SpCoV HKU171ab520–1935218,8336,278+1, +357ACACCA(452)AUG
S19334–229543,6211,207+219188ACACCA(140)AUG
E22948–2319624983+122925ACACCG(17)AUG
M23189–23842654218+223166ACACCA(17)AUG
NS623842–2412928896+123790ACACCA(46)AUG
N24150–251781,029343+324137ACACCA(7)AUG
NS7a25189–25623435145+125179ACACCA(4)AUG
NS7b25539–2575121371+325523ACUCCA(10)AUG
MRCoV HKU181ab596–1935618,7616,254+2, +164ACACCA(526)AUG
S19338–229913,6541,218+319192ACACCA(140)AUG
E22985–2323324983+222945ACACCG(34)AUG
M23226–23882657219+323203ACACCA(17)AUG
NS623882–2417229197+223857ACGCCA(19)AUG
N24355–253951,041347+124340ACACCA(9)AUG
NS7a25407–2558017458+325396ACACCA(5)AUG
NS7b25561–25932372124+1
NS7c25941–2619525585+325910ACACCA(25)AUG
NHCoV HKU191ab482–1932318,8426,281+2, +167ACACCG(409)AUG
S19305–230693,7651,255+319156ACACCG(143)AUG
E23069–2331724983+223013ACACCA(50)AUG
M23310–23960651217+323211ACACCG(93)AUG
NS623960–2423827993+223951ACACCU(3)AUG
N24248–252761,029343+224231ACACCU(8)AUG
NS7a25277–2557329799+225248ACACCG(23)AUG
NS7b25583–2587629498+225560ACACCA(17)AUG
WiCoV HKU201ab219–1883818,6206,207+3, +260ACACCA(153)AUG
S18817–224553,6391,213+118731ACACCU(80)AUG
E22455–2271526187+322380ACACCA(69)AUG
M22708–23358651217+122597ACACCG(105)AUG
NS623358–2363027391+3
N23646–246981,053351+323631ACACCA(9)AUG
NS7a24695–2492823478+224609AAACCA(80)AUG
NS7b25218–2546624983+325177ACACCG(35)AUG
NS7c25450–2571626789+125444ACACCGAUG
NS7d25752–2595220167+325735AAACCU(11)AUG
CMCoV HKU211ab478–1910318,6266,209+1, +363ACACCA(409)AUG
S19085–227293,6451,215+218939ACACCA(140)AUG
E22723–2297124983+122697ACACCA(20)AUG
M22973–23779807269+222938ACACCA(29)AUG
NS623779–2402424682+123727ACACCA(46)AUG
N24052–251071,056352+124039ACACCG(7)AUG
NS7a25107–2537927391+325036ACACCU(65)AUG
NS7b25391–2557618662+225379ACACCU(6)AUG
NS7c25500–25916417139+2
PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21. aa, amino acid; nt, nucleotide.
Boldface indicates putative TRS sequences. The nucleotide variations are in italic.

The seven novel CoVs display similar genome organizations and differ only in the number of ORFs downstream of N (Fig. 2). Their transcription regulatory sequences (TRSs) conform to the consensus motif 5′-ACACCA-3′ (Table 3), which appears to be unique to members of the genus Deltacoronavirus. Interestingly, similar to BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13, the perfect TRSs of S in the genomes of the seven novel CoVs were separated from the corresponding AUG by 80 to 145 bases (Table 3). This is in contrast to the relatively small number of bases between the TRSs for S and the corresponding AUG (range: from 0 bases in HCoV-NL63, Rhinolophus bat coronavirus HKU2 [Rh-BatCoV-HKU2], HCoV-HKU1, bovine coronavirus [BCoV], HCoV-OC43, mouse hepatitis virus [MHV], porcine hemagglutinating encephalomyelitis virus, SARS-CoV, and SARS-related Rhinolophus bat coronavirus HKU3 [SARSr-Rh-batCoV HKU3] to 52 bases in infectious bronchitis virus [IBV]) in members of Alphacoronavirus, Betacoronavirus, and Gammacoronavirus. Similar to BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13, the genomes of the seven novel CoVs have putative PL, which are homologous to PL2 of Alphacoronavirus and Betacoronavirus subgroup A and PL of Betacoronavirus subgroups B, C, and D and Gammacoronavirus (Fig. 2). Similar to BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13, one ORF (NS6) is found between M and N of the genomes of the seven novel CoVs. On the other hand, one ORF (NS7) is present overlapping with N in PorCoV HKU15, two ORFs (NS7a and 7b) are present overlapping or downstream of N in WECoV HKU16, SpCoV HKU17, and NHCoV HKU19, three ORFs (NS7a, 7b, and 7c) are present downstream of N in MRCoV HKU18 and CMCoV HKU21, and four ORFs (NS7a, 7b, 7c, and 7d) are present overlapping or downstream of N in WiCoV HKU20. For NS7 of PorCoV, the presence of an imperfect TRS (GCACCA) and its relatively high Ka/Ks ratio (number of nonsynonymous substitutions per nonsynonymous site/number of synonymous substitutions per synonymous site) of 1.046 (data not shown) implied that this ORF may not be expressed. BLAST search revealed no amino acid similarities between these putative nonstructural proteins and other known proteins, and no functional domain was identified by PFAM and InterProScan, except that NS7a of NHCoV HKU19 was found to be homologous to the NS7a of BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13. NS7b of WiCoV HKU20 and CMCoV HKU21, and NS7d of WiCoV HKU20, were also found to be homologous to the NS3b of IBV and hypothetical protein of goose coronavirus, respectively. Transmembrane helices, predicted by TMHMM and TMpred, in putative accessory proteins downstream to the N genes in the genomes of SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 are listed in Table S4 in the supplemental material. Each of the genomes of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, and CMCoV HKU21 contains a stem-loop II motif (s2m) (residues 25,220 to 25,251, 25,825 to 25,856, 25,865 to 25,896, 26,472 to 26,503, and 26,013 to 26,044, respectively), a conserved RNA element downstream of N and upstream of the poly(A) tail, similar to those in IBV, TCoV, SARSr-Rh-BatCoV, and SARS-CoV, as well as other CoVs discovered in Asian leopard cat, graylag geese, feral pigeons, and mallards, for which complete genomes are not available (Fig. 3) (14, 21, 38).

An external file that holds a picture, illustration, etc.
Object name is zjv9990958200003.jpg

Multiple alignments of conserved s2m of infectious bronchitis virus (IBV), SARS-related human coronavirus (SARS CoV), SARS-related Rhinolophus bat coronavirus HKU3 (SARSr-Rh-BatCoV HKU3), BuCoV HKU11, ThCoV HKU12, MunCoV HKU13, Asian leopard cat coronavirus (ALCCoV), PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, and CMCoV HKU21. Identical nucleotides are marked by asterisks. Acc. No., accession no.

Comparison of the amino acid identities of the seven conserved replicase domains for species demarcation (ADRP, nsp5 [3CL], nsp12 [RdRp], nsp13 [Hel], nsp14 [ExoN], nsp15 [NendoU], and nsp16 [O-MT]) (8) among the 10 deltacoronaviruses is shown in Table S5 in the supplemental material. In all the seven domains, the amino acid sequences of PorCoV HKU15 and SpCoV HKU17 showed more than 90% identity, indicating that these two coronaviruses should be subspecies of the same species.

Phylogenetic analyses.

The phylogenetic trees constructed using the nucleotide sequences of the 3CL, RdRp, Hel, S, and N of the seven novel CoVs and other CoVs are shown in Fig. 4 and the corresponding pairwise amino acid identities are shown in Table 2. For all five genes, the seven novel CoVs possessed higher amino acid identities to each other and BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13 than to any other known CoVs with complete genomes available (Table 2). In all five trees, the seven novel CoVs were clustered with BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13 (Fig. 4). For Hel, S, and N, PorCoVs were also clustered with a CoV found in Asian leopard cat (10), for which the sequences of these genes were available (Fig. 4). There were <2% base differences between the Hel, S, and N genes of PorCoV and those of the Asian leopard cat coronavirus. Based on both phylogenetic tree analyses and amino acid differences, the seven novel CoVs as well as BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13 should belong to the same genus, Deltacoronavirus.

An external file that holds a picture, illustration, etc.
Object name is zjv999095820004a.jpg
An external file that holds a picture, illustration, etc.
Object name is zjv999095820004b.jpg
An external file that holds a picture, illustration, etc.
Object name is zjv999095820004c.jpg

Phylogenetic analyses of 3CL, RdRp, helicase (Hel), S, and N proteins of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21. The trees were constructed by using the neighbor joining method using Kimura correction and bootstrap values calculated from 1,000 trees. Two hundred ninety-five, 892, 590, 802, and 249 amino acid positions in 3CL, RdRp, Hel, S, and N, respectively, were included in the analyses. The trees were midpoint rooted. For 3CL and S, the scale bar indicates the estimated number of substitutions per 10 amino acids. For RdRp and Hel, the scale bar indicates the estimated number of substitutions per 20 amino acids. For N, the scale bar indicates the estimated number of substitutions per 5 amino acids. Viruses characterized in this study are in bold. Virus name abbreviations are the same as those in the Fig. 1 legend.

Estimation of divergence dates.

Using the Bayesian Skyline under a relaxed-clock model with an uncorrelated log-normal distribution, the mean evolutionary rate of CoVs was estimated at 1.3 × 10 nucleotide substitutions per site per year for the RdRp gene. Molecular clock analysis using the RdRp gene showed that the tMRCA of all CoVs was estimated at ∼8100 BC (HPDs, 20607 to 974 BC), that of Alphacoronavirus at ∼2400 BC (HPDs, 7659 to 722 BC), that of Betacoronavirus at ∼3300 BC (HPDs, 9713 to 447 BC), that of Gammacoronavirus at ∼2800 BC (HPDs, 8840 to 700 BC), and that of Deltacoronavirus at ∼3000 BC (HPDs, 9073 to 555 BC) (Fig. 5).

An external file that holds a picture, illustration, etc.
Object name is zjv9990958200005.jpg

Estimation of the time to the most recent common ancestor for Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Deltacoronavirus. The time-scaled phylogeny was summarized from all MCMC phylogenies of the RdRp gene data set analyzed under the relaxed-clock model with an uncorrelated log-normal distribution in BEAST version 1.6.1. Viruses characterized in this study are in bold. The numbers indicate number of years ago. This is shown in the scale bar. Virus name abbreviations are the same as those in the legends of Fig. 1.

DISCUSSION

The diversity of CoVs in birds is comparable to that observed in bats. In the last 7 years, we and others have demonstrated a previously unrecognized diversity of CoVs in bats (4, 6, 23, 25, 26, 28, 43). More than 10 CoVs were discovered in bats, with at least nine present in our locality, and complete genome sequences are available for eight, which includes SARSr-Rh-BatCoV HKU3, Rh-BatCoV-HKU2, Miniopterus bat coronavirus 1, Miniopterus bat coronavirus HKU8, Scotophilus bat coronavirus 512, Tylonycteris bat coronavirus HKU4, Pipistrellus bat coronavirus HKU5, and Rousettus bat coronavirus HKU9 (4, 6, 25, 26, 43). Due to the similarities between bats and birds, such as their abilities to fly and high species diversity, we hypothesized that there should be previously unrecognized CoVs in birds. In our previous study and the present one, we demonstrated that there are at least nine CoVs, in addition to IBV and its close relatives, in birds (49). Potentially novel CoVs in Gammacoronavirus were also observed in another study, although complete genome sequences are not available and therefore detailed genomic and phylogenetic analysis are not possible (35). The nine CoVs discovered in the present and previous studies were found in birds of nine different families, showing host specificity. This phenomenon of host specificity is similar to that observed in bats, in which different genera are hosts of different CoVs (26, 45, 51, 52). We speculate that this diversity and host specificity of bat and bird CoVs is due to the large variety of species in bats and birds, giving rise to a large variety of cell types and receptors for the different CoVs to attach and replicate.

The presence of a huge diversity of bat CoVs in Alphacoronavirus and Betacoronavirus but not Gammacoronavirus and Deltacoronavirus and a huge diversity of bird CoVs in Gammacoronavirus and Deltacoronavirus but not Alphacoronavirus and Betacoronavirus supports our model of CoV evolution, in which bats are the gene source of Alphacoronavirus and Betacoronavirus and birds the gene source of Gammacoronavirus and Deltacoronavirus (Fig. 6) (52). It is not known whether the first CoVs occurred in bats and jumped to birds or vice versa. In the bat CoV lineage, the bat CoV jumped to another species of bat, giving rise to Alphacoronavirus and Betacoronavirus. These bat CoVs in turn jumped to other bat species and other mammals, including humans, with each interspecies jumping evolving dichotomously. As for the bird CoV lineage, the bird CoV jumped to another species of bird, giving rise to Gammacoronavirus and Deltacoronavirus. These bird CoVs in turn jumped to other bird species and occasionally to some mammalian species, such as whale and pig, with each interspecies jumping evolving dichotomously. Although PorCoV HKU15 was closely related to a CoV previously found in Asian leopard cats and Chinese ferret badgers, further experiments are warranted to confirm whether these viruses really replicate in the corresponding animals. Of note is that the estimation of divergence time was based on a relaxed-clock assumption with no recombination among the genomes. Since CoVs have a tendency to recombine, the estimated divergence time gives only a rough approximation of the actual divergence time. When more complete genomes of CoVs in the four different genera at different time points are available, such divergence time estimation can be performed using multiple gene loci to achieve more accurate estimation.

An external file that holds a picture, illustration, etc.
Object name is zjv9990958200006.jpg

A model of CoV evolution. CoVs in bats are the gene source of Alphacoronavirus and Betacoronavirus, and CoVs in birds are the gene source of Gammacoronavirus and Deltacoronavirus.

Both avian and mammalian CoVs are members of Deltacoronavirus, with similar genome characteristics and structures. In all the 10 members of Deltacoronavirus with complete genome sequences available, all have a very small genome size, from 25.421 (PorCoV HKU15) to 26.674 (MRCoV HKU18) kb, the smallest among all CoVs. Only one papain-like protease domain is observed in the nsp3 gene of their genomes. As for their gene contents, ORF NS6 was present between the M and N genes, and one to four ORFs were also observed downstream to the N gene. As for the TRSs, they all have the same putative TRS of ACACCA and separation of the TRS from the AUG of the S gene by a long stretch of nucleotides. Despite these similar genome characteristics among members of Deltacoronavirus, NHCoV HKU19 and WiCoV HKU20 possessed genomic features distinct from the other members of Deltacoronavirus, including the amino acids upstream of the putative cleavage sites at the junction of nsp2/nsp3, nsp3/nsp4, and nsp4/nsp5. It is also notable that NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 occupied the first three branches in the phylogenetic trees constructed using 3CL, Hel, RdRp, and N, indicating that they could be more ancestral than the other members. Furthermore, these three CoVs were found in large birds, including black-crowned night heron, Eurasian wigeon, and common moorhen, in contrast to BuCoV HKU11, ThCoV HKU12, MunCoV HKU13, WECoV HKU16, SpCoV HKU17, and MRCoV HKU18, which were found in small birds, including bulbuls, blackbird, gray-backed thrush, munias, Japanese white-eye, Eurasian tree sparrow, and oriental magpie robin. We speculate that the change in genome characteristics (e.g., acquisition of s2m) could have occurred during interspecies jumping of the CoV within the large birds before the jump to the small birds. Interestingly, the fact that PorCoV HKU15 and SpCoV HKU17 are the same species implies that interspecies jumping from birds to pigs may have occurred relatively recently. It is possible that a deletion of 3′ Ns7a and Ns7b had occurred during interspecies jumping from birds to pigs, which is similar to the observation of interspecies jumping of SARS-CoV from civets to humans, with the deletion of 29 bp in ORF 8 (25). As for the Asian leopard cat coronavirus, with only the Hel, S, E, M, and N gene sequences available, the sequences of these gene fragments differ from the corresponding ones in PorCoV by less than 2.1% nucleotides or 1.7% amino acids, including that for the S gene, which is responsible for receptor binding. BEAST analysis showed that the CoV jumped from birds to mammals around 523 years ago (Fig. 5). The mixing of birds, pigs, and other mammals in domestic environments and wildlife markets as well as their close contacts with humans may provide the correct environment for interspecies jumping and could subsequently pose risks of further genetic changes for adapting to human host as in the case of SARS (5). More extensive epidemiological studies in different varieties of mammalian species in other parts of the world for members of Deltacoronavirus would further improve our understanding on the diversity of this genus as well as its evolutionary history.

Supplementary Material

Supplemental material:
Department of Microbiology
State Key Laboratory of Emerging Infectious Diseases
Research Centre of Infection and Immunology
the Carol Yu Centre for Infection, The University of Hong Kong, Hong Kong
Guangzhou Center for Disease Control and Prevention, Guangzhou, China
Corresponding author.
Address correspondence to Kwok-Yung Yuen, kh.ukh.ccukh@neuyyk.
P.C.Y.W. and S.K.P.L. contributed equally to this article.
Address correspondence to Kwok-Yung Yuen, kh.ukh.ccukh@neuyyk.
P.C.Y.W. and S.K.P.L. contributed equally to this article.
Received 2011 Oct 12; Accepted 2012 Jan 17.

Abstract

Recently, we reported the discovery of three novel coronaviruses, bulbul coronavirus HKU11, thrush coronavirus HKU12, and munia coronavirus HKU13, which were identified as representatives of a novel genus, Deltacoronavirus, in the subfamily Coronavirinae. In this territory-wide molecular epidemiology study involving 3,137 mammals and 3,298 birds, we discovered seven additional novel deltacoronaviruses in pigs and birds, which we named porcine coronavirus HKU15, white-eye coronavirus HKU16, sparrow coronavirus HKU17, magpie robin coronavirus HKU18, night heron coronavirus HKU19, wigeon coronavirus HKU20, and common moorhen coronavirus HKU21. Complete genome sequencing and comparative genome analysis showed that the avian and mammalian deltacoronaviruses have similar genome characteristics and structures. They all have relatively small genomes (25.421 to 26.674 kb), the smallest among all coronaviruses. They all have a single papain-like protease domain in the nsp3 gene; an accessory gene, NS6 open reading frame (ORF), located between the M and N genes; and a variable number of accessory genes (up to four) downstream of the N gene. Moreover, they all have the same putative transcription regulatory sequence of ACACCA. Molecular clock analysis showed that the most recent common ancestor of all coronaviruses was estimated at approximately 8100 BC, and those of Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Deltacoronavirus were at approximately 2400 BC, 3300 BC, 2800 BC, and 3000 BC, respectively. From our studies, it appears that bats and birds, the warm blooded flying vertebrates, are ideal hosts for the coronavirus gene source, bats for Alphacoronavirus and Betacoronavirus and birds for Gammacoronavirus and Deltacoronavirus, to fuel coronavirus evolution and dissemination.

Abstract
Click here to view.

ACKNOWLEDGMENTS

We thank York Y. N. Chow, Health, Welfare and Food, HKSAR, The Peoples' Republic of China; Alan Chi-Kong Wong, Siu Fai Leung, Chik Chuen Lay, Thomas Sit, Elaine Lee, and Geraldine Luk of the Agriculture, Fisheries, and Conservation Department; and Clement Leung, Constance Chan, and Wing Ka Au of the Food, Environmental and Hygiene Department of the HKSAR.

We are grateful to the generous support of Hui Hoy and Hui Ming in the genomic sequencing platform and Eunice Lam for her generous donation to emerging infectious disease research. This work is partly supported by Research Grant Council grant HKU 780709 M; University Development Fund and Outstanding Young Researcher Award, The University of Hong Kong; The Tung Wah Group of Hospitals Fund for Research in Infectious Diseases; the HKSAR Research Fund for the Control of Infectious Diseases of the Health, Welfare and Food Bureau; and the Shaw Foundation.

ACKNOWLEDGMENTS

Footnotes

Published ahead of print 25 January 2012

Supplemental material for this article may be found at http://jvi.asm.org/.

Footnotes

REFERENCES

REFERENCES

References

  • 1. Apweiler R, et al. 2001. The InterPro database, an integrated documentation resource for protein families, domains and functional sites. Nucleic Acids Res.29:37–40
  • 2. Bateman A, et al. 2002. The Pfam protein families database. Nucleic Acids Res.30:276–280
  • 3. Brian DA, Baric RS. 2005 Coronavirus genome structure and replication. Curr. Top. Microbiol. Immunol.287:1–30 [[PubMed][Google Scholar]
  • 4. Cao J, Wu CC, Lin TL. 2008 Complete nucleotide sequence of polyprotein gene 1 and genome organization of turkey coronavirus. Virus Res.136:43–49 [[PubMed][Google Scholar]
  • 5. Cheng VC, Lau SK, Woo PC, Yuen KY. 2007 Severe acute respiratory syndrome coronavirus as an agent of emerging and reemerging infection. Clin. Microbiol. Rev.20:660–694 [Google Scholar]
  • 6. Chu DK, Peiris JS, Chen H, Guan Y, Poon LL. 2008 Genomic characterizations of bat coronaviruses (1A, 1B and HKU8) and evidence for co-infections in Miniopterus bats. J. Gen. Virol.89:1282–1287 [[PubMed][Google Scholar]
  • 7. Circella E, et al. 2007. Coronavirus associated with an enteric syndrome on a quail farm. Avian Pathol.36:251–258 [[PubMed]
  • 8. de Groot RJ, et al. 2011. Coronaviridae, p 806–828 In King AMQ, Adams MJ, Carstens EB, Lefkowitz EJ, editors. (ed), Virus taxonomy: ninth report of the International Committee on Taxonomy of Viruses, International Union of Microbiological Societies, Virology Division. Elsevier Academic Press, London, United Kingdom [PubMed]
  • 9. Dominguez SR, O'Shea TJ, Oko LM, Holmes KV. 2007 Detection of group 1 coronaviruses in bats in North America. Emerg. Infect. Dis.13:1295–1300 [Google Scholar]
  • 10. Dong BQ, et al. 2007. Detection of a novel and highly divergent coronavirus from Asian leopard cats and Chinese ferret badgers in Southern China. J. Virol.81:6920–6926
  • 11. Drummond AJ, Rambaut A. 2007 BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol.7:214. [Google Scholar]
  • 12. Fouchier RA, et al. 2004. A previously undescribed coronavirus associated with respiratory disease in humans. Proc. Natl. Acad. Sci. U. S. A.101:6212–6216
  • 13. Gloza-Rausch F, et al. 2008. Detection and prevalence patterns of group I coronaviruses in bats, northern Germany. Emerg. Infect. Dis.14:626–631
  • 14. Gomaa MH, Barta JR, Ojkic D, Yoo D. 2008 Complete genomic sequence of turkey coronavirus. Virus Res.135:237–246 [[PubMed][Google Scholar]
  • 15. Gough RE, Drury SE, Culver F, Britton P, Cavanagh D. 2006 Isolation of a coronavirus from a green-cheeked Amazon parrot (Amazon viridigenalis Cassin). Avian Pathol.35:122–126 [[PubMed][Google Scholar]
  • 16. Guan Y, et al. 2003. Isolation and characterization of viruses related to the SARS coronavirus from animals in southern China. Science302:276–278 [[PubMed]
  • 17. Hasoksuz M, et al. 2007. Biologic, antigenic, and full-length genomic characterization of a bovine-like coronavirus isolated from a giraffe. J. Virol.81:4981–4990
  • 18. Herrewegh AA, Smeenk I, Horzinek MC, Rottier PJ, de Groot RJ. 1998 Feline coronavirus type II strains 79-1683 and 79-1146 originate from a double recombination between feline coronavirus type I and canine coronavirus. J. Virol.72:4508–4514 [Google Scholar]
  • 19. Hofmann K, Stoffel W. 1993 TMBASE - a database of membrane spanning protein segments. Biol. Chem. Hoppe-Seyler374:166 [PubMed][Google Scholar]
  • 20. Huang Y, Lau SK, Woo PC, Yuen KY. 2008 CoVDB: a comprehensive database for comparative analysis of coronavirus genes and genomes. Nucleic Acids Res.36:D504–D511 [Google Scholar]
  • 21. Jonassen CM, et al. 2005. Molecular identification and characterization of novel coronaviruses infecting graylag geese (Anser anser), feral pigeons (Columbia livia) and mallards (Anas platyrhynchos). J. Gen. Virol.86:1597–1607 [[PubMed]
  • 22. Lai MM, Cavanagh D. 1997 The molecular biology of coronaviruses. Adv. Virus Res.48:1–100 [[PubMed][Google Scholar]
  • 23. Lau SK, et al. 2010. Ecoepidemiology and complete genome comparison of different strains of severe acute respiratory syndrome-related Rhinolophus bat coronavirus in China reveal bats as a reservoir for acute, self-limiting infection that allows recombination events. J. Virol.84:2808–2819
  • 24. Lau SK, et al. 2006. Coronavirus HKU1 and other coronavirus infections in Hong Kong. J. Clin. Microbiol.44:2063–2071
  • 25. Lau SK, et al. 2005. Severe acute respiratory syndrome coronavirus-like virus in Chinese horseshoe bats. Proc. Natl. Acad. Sci. U. S. A.102:14040–14045
  • 26. Lau SK, et al. 2007. Complete genome sequence of bat coronavirus HKU2 from Chinese horseshoe bats revealed a much smaller spike gene with a different evolutionary lineage from the rest of the genome. Virology367:428–439 [[PubMed]
  • 27. Lau SK, et al. 2011. Molecular epidemiology of human coronavirus OC43 reveals evolution of different genotypes over time and recent emergence of a novel genotype due to natural recombination. J. Virol.85:11325–11337
  • 28. Lau SK, et al. 2010. Coexistence of different genotypes in the same bat and serological characterization of Rousettus bat coronavirus HKU9 belonging to a novel Betacoronavirus subgroup. J. Virol.84:11385–11394
  • 29. Li KS, et al. 2004. Genesis of a highly pathogenic and potentially pandemic H5N1 influenza virus in eastern Asia. Nature430:209–213 [[PubMed]
  • 30. Li W, et al. 2005. Bats are natural reservoirs of SARS-like coronaviruses. Science310:676–679 [[PubMed]
  • 31. Liu S, et al. 2005. Isolation of avian infectious bronchitis coronavirus from domestic peafowl (Pavo cristatus) and teal (Anas). J. Gen. Virol.86:719–725 [[PubMed]
  • 32. Mardani K, Noormohammadi AH, Hooper P, Ignjatovic J, Browning GF. 2008 Infectious bronchitis viruses with a novel genomic organization. J. Virol.82:2013–2024 [Google Scholar]
  • 33. Marra MA, et al. 2003. The genome sequence of the SARS-associated coronavirus. Science300:1399–1404 [[PubMed]
  • 34. Mihindukulasuriya KA, Wu G, St. Leger J, Nordhausen RW, Wang D. 2008. Identification of a novel coronavirus from a beluga whale by using a panviral microarray. J. Virol.82:5084–5088
  • 35. Muradrasoli S, et al. 2010. Prevalence and phylogeny of coronaviruses in wild birds from the Bering Strait area (Beringia). PLoS One5:e13640.
  • 36. Peiris JS, et al. 2003. Coronavirus as a possible cause of severe acute respiratory syndrome. Lancet361:1319–1325 [[PubMed]
  • 37. Poon LL, et al. 2005. Identification of a novel coronavirus in bats. J. Virol.79:2001–2009
  • 38. Robertson MP, Igel H, Baertsch R, Haussler D, Ares M, Jr, Scott WG. 2005 The structure of a rigorously conserved RNA element within the SARS virus genome. PLoS Biol.3:e5. [Google Scholar]
  • 39. Rota PA, et al. 2003. Characterization of a novel coronavirus associated with severe acute respiratory syndrome. Science300:1394–1399 [[PubMed]
  • 40. Snijder EJ, et al. 2003. Unique and conserved features of genome and proteome of SARS-coronavirus, an early split-off from the coronavirus group 2 lineage. J. Mol. Biol.331:991–1004 [[PubMed]
  • 41. Sonnhammer EL, von Heijne G, Krogh A. 1998 A hidden Markov model for predicting transmembrane helices in protein sequences. Proc. Int. Conf. Intell. Syst. Mol. Biol.6:175–182 [[PubMed][Google Scholar]
  • 42. Suchard MA, Weiss RE, Sinsheimer JS. 2001 Bayesian selection of continuous-time Markov chain evolutionary models. Mol. Biol. Evol.18:1001–1013 [[PubMed][Google Scholar]
  • 43. Tang XC, et al. 2006. Prevalence and genetic diversity of coronaviruses in bats from China. J. Virol.80:7481–7490
  • 44. van der Hoek L, et al. 2004. Identification of a new human coronavirus. Nat. Med.10:368–373 [[PubMed]
  • 45. Woo PC, et al. 2007. Comparative analysis of twelve genomes of three novel group 2c and group 2d coronaviruses reveals unique group and subgroup features. J. Virol.81:1574–1585
  • 46. Woo PC, et al. 2004. Relative rates of non-pneumonic SARS coronavirus infection and SARS coronavirus pneumonia. Lancet363:841–845 [[PubMed]
  • 47. Woo PC, et al. 2006. Comparative analysis of 22 coronavirus HKU1 genomes reveals a novel genotype and evidence of natural recombination in coronavirus HKU1. J. Virol.80:7136–7145
  • 48. Woo PC, et al. 2005. Characterization and complete genome sequence of a novel coronavirus, coronavirus HKU1, from patients with pneumonia. J. Virol.79:884–895
  • 49. Woo PC, et al. 2009. Comparative analysis of complete genome sequences of three avian coronaviruses reveals a novel group 3c coronavirus. J. Virol.83:908–917
  • 50. Woo PC, et al. 2005. Clinical and molecular epidemiological features of coronavirus HKU1-associated community-acquired pneumonia. J. Infect. Dis.192:1898–1907 [[PubMed]
  • 51. Woo PC, et al. 2006. Molecular diversity of coronaviruses in bats. Virology351:180–187 [[PubMed]
  • 52. Woo PC, Lau SK, Huang Y, Yuen KY. 2009 Coronavirus diversity, phylogeny and interspecies jumping. Exp. Biol. Med. (Maywood)234:1117–1127 [[PubMed][Google Scholar]
  • 53. Zhang J, et al. 2007. Genomic characterization of equine coronavirus. Virology369:92–104 [[PubMed]
  • 54. Ziebuhr J. 2004 Molecular biology of severe acute respiratory syndrome coronavirus. Curr. Opin. Microbiol.7:412–419 [[PubMed][Google Scholar]
Collaboration tool especially designed for Life Science professionals.Drag-and-drop any entity to your messages.