Carol S. F. Lam

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.

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.

Comparative Analysis of Twelve Genomes of Three Novel Group 2c and Group 2d Coronaviruses Reveals Unique Group and Subgroup Features

ABSTRACT Twelve complete genomes of three novel coronaviruses—bat coronavirus HKU4 (bat-CoV HKU4), bat-CoV HKU5 (putative group 2c), and bat-CoV HKU9 (putative group 2d)—were sequenced. Comparative genome analysis showed that the various open reading frames (ORFs) of the genomes of the three coronaviruses had significantly higher amino acid identities to those of other group 2 coronaviruses than group 1 and 3 coronaviruses. Phylogenetic trees constructed using chymotrypsin-like protease, RNA-dependent RNA polymerase, helicase, spike, and nucleocapsid all showed that the group 2a and 2b and putative group 2c and 2d coronaviruses are more closely related to each other than to group 1 and 3 coronaviruses. Unique genomic features distinguishing between these four subgroups, including the number of papain-like proteases, the presence or absence of hemagglutinin esterase, small ORFs between the membrane and nucleocapsid genes and ORFs (NS7a and NS7b), bulged stem-loop and pseudoknot structures downstream of the nucleocapsid gene, transcription regulatory sequence, and ribosomal recognition signal for the envelope gene, were also observed. This is the first time that NS7a and NS7b downstream of the nucleocapsid gene has been found in a group 2 coronavirus. The high Ka/Ks ratio of NS7a and NS7b in bat-CoV HKU9 implies that these two group 2d-specific genes are under high selective pressure and hence are rapidly evolving. The four subgroups of group 2 coronaviruses probably originated from a common ancestor. Further molecular epidemiological studies on coronaviruses in the bats of other countries, as well as in other animals, and complete genome sequencing will shed more light on coronavirus diversity and their evolutionary histories.

Comparative Analysis of Complete Genome Sequences of Three Avian Coronaviruses Reveals a Novel Group 3c Coronavirus

ABSTRACT In this territory-wide molecular epidemiology study of coronaviruses (CoVs) in Hong Kong involving 1,541 dead wild birds, three novel CoVs were identified in three different bird families (bulbul CoV HKU11 [BuCoV HKU11], thrush CoV HKU12 [ThCoV HKU12], and munia CoV HKU13 [MuCoV HKU13]). Four complete genomes of the three novel CoVs were sequenced. Their genomes (26,396 to 26,552 bases) represent the smallest known CoV genomes. In phylogenetic trees constructed using chymotrypsin-like protease (3CLpro), RNA-dependent RNA polymerase (Pol), helicase, spike, and nucleocapsid proteins, BuCoV HKU11, ThCoV HKU12, and MuCoV HKU13 formed a cluster distantly related to infectious bronchitis virus and turkey CoV (group 3a CoVs). For helicase, spike, and nucleocapsid, they were also clustered with a CoV recently discovered in Asian leopard cats, for which the complete genome sequence was not available. The 3CLpro, Pol, helicase, and nucleocapsid of the three CoVs possessed higher amino acid identities to those of group 3a CoVs than to those of group 1 and group 2 CoVs. Unique genomic features distinguishing them from other group 3 CoVs include a distinct transcription regulatory sequence and coding potential for small open reading frames. Based on these results, we propose a novel CoV subgroup, group 3c, to describe this distinct subgroup of CoVs under the group 3 CoVs. Avian CoVs are genetically more diverse than previously thought and may be closely related to some newly identified mammalian CoVs. Further studies would be important to delineate whether the Asian leopard cat CoV was a result of interspecies jumping from birds, a situation analogous to that of bat and civet severe acute respiratory syndrome CoVs.

Genetic Characterization of Betacoronavirus Lineage C Viruses in Bats Reveals Marked Sequence Divergence in the Spike Protein of Pipistrellus Bat Coronavirus HKU5 in Japanese Pipistrelle: Implications for the Origin of the Novel Middle East Respiratory Syndrome Coronavirus

ABSTRACT While the novel Middle East respiratory syndrome coronavirus (MERS-CoV) is closely related to Tylonycteris bat CoV HKU4 (Ty-BatCoV HKU4) and Pipistrellus bat CoV HKU5 (Pi-BatCoV HKU5) in bats from Hong Kong, and other potential lineage C betacoronaviruses in bats from Africa, Europe, and America, its animal origin remains obscure. To better understand the role of bats in its origin, we examined the molecular epidemiology and evolution of lineage C betacoronaviruses among bats. Ty-BatCoV HKU4 and Pi-BatCoV HKU5 were detected in 29% and 25% of alimentary samples from lesser bamboo bat (Tylonycteris pachypus) and Japanese pipistrelle (Pipistrellus abramus), respectively. Sequencing of their RNA polymerase (RdRp), spike (S), and nucleocapsid (N) genes revealed that MERS-CoV is more closely related to Pi-BatCoV HKU5 in RdRp (92.1% to 92.3% amino acid [aa] identity) but is more closely related to Ty-BatCoV HKU4 in S (66.8% to 67.4% aa identity) and N (71.9% to 72.3% aa identity). Although both viruses were under purifying selection, the S of Pi-BatCoV HKU5 displayed marked sequence polymorphisms and more positively selected sites than that of Ty-BatCoV HKU4, suggesting that Pi-BatCoV HKU5 may generate variants to occupy new ecological niches along with its host in diverse habitats. Molecular clock analysis showed that they diverged from a common ancestor with MERS-CoV at least several centuries ago. Although MERS-CoV may have diverged from potential lineage C betacoronaviruses in European bats more recently, these bat viruses were unlikely to be the direct ancestor of MERS-CoV. Intensive surveillance for lineage C betaCoVs in Pipistrellus and related bats with diverse habitats and other animals in the Middle East may fill the evolutionary gap.

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

ABSTRACT Despite the identification of severe acute respiratory syndrome-related coronavirus (SARSr-CoV) in Rhinolophus Chinese horseshoe bats (SARSr-Rh-BatCoV) in China, the evolutionary and possible recombination origin of SARSr-CoV remains undetermined. We carried out the first study to investigate the migration pattern and SARSr-Rh-BatCoV genome epidemiology in Chinese horseshoe bats during a 4-year period. Of 1,401 Chinese horseshoe bats from Hong Kong and Guangdong, China, that were sampled, SARSr-Rh-BatCoV was detected in alimentary specimens from 130 (9.3%) bats, with peak activity during spring. A tagging exercise of 511 bats showed migration distances from 1.86 to 17 km. Bats carrying SARSr-Rh-BatCoV appeared healthy, with viral clearance occurring between 2 weeks and 4 months. However, lower body weights were observed in bats positive for SARSr-Rh-BatCoV, but not Rh-BatCoV HKU2. Complete genome sequencing of 10 SARSr-Rh-BatCoV strains showed frequent recombination between different strains. Moreover, recombination was detected between SARSr-Rh-BatCoV Rp3 from Guangxi, China, and Rf1 from Hubei, China, in the possible generation of civet SARSr-CoV SZ3, with a breakpoint at the nsp16/spike region. Molecular clock analysis showed that SARSr-CoVs were newly emerged viruses with the time of the most recent common ancestor (tMRCA) at 1972, which diverged between civet and bat strains in 1995. The present data suggest that SARSr-Rh-BatCoV causes acute, self-limiting infection in horseshoe bats, which serve as a reservoir for recombination between strains from different geographical locations within reachable foraging range. Civet SARSr-CoV is likely a recombinant virus arising from SARSr-CoV strains closely related to SARSr-Rh-BatCoV Rp3 and Rf1. Such frequent recombination, coupled with rapid evolution especially in ORF7b/ORF8 region, in these animals may have accounted for the cross-species transmission and emergence of SARS.

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.

Abstract Apart from bat-SARS-CoV, we have identified a novel group 1 coronavirus, bat-CoV HKU2, in Rhinolophus sinicus (Chinese horseshoe bats). Since it has been suggested that the receptor-binding motif (RBM) of SARS-CoV may have been acquired from a group 1 coronavirus, we conducted a surveillance study and identified bat-SARS-CoV and bat-CoV HKU2 in 8.7% and 7.5% respectively of R. sinicus in Hong Kong and Guangdong. Complete genome sequencing of four strains of bat-CoV HKU2 revealed the smallest coronavirus genome (27164 nucleotides) and a unique spike protein evolutionarily distinct from the rest of the genome. This spike protein, sharing similar deletions with other group 2 coronaviruses in its C-terminus, also contained a 15-amino acid peptide homologous to a corresponding peptide within the RBM of spike protein of SARS-CoV, which was absent in other coronaviruses except bat-SARS-CoV. These suggest a common evolutionary origin in the spike protein of bat-CoV HKU2, bat-SARS-CoV, and SARS-CoV.

Complete genome analysis of three novel picornaviruses from diverse bat species

ABSTRACT Although bats are important reservoirs of diverse viruses that can cause human epidemics, little is known about the presence of picornaviruses in these flying mammals. Among 1,108 bats of 18 species studied, three novel picornaviruses (groups 1, 2, and 3) were identified from alimentary specimens of 12 bats from five species and four genera. Two complete genomes, each from the three picornaviruses, were sequenced. Phylogenetic analysis showed that they fell into three distinct clusters in the Picornaviridae family, with low homologies to known picornaviruses, especially in leader and 2A proteins. Moreover, group 1 and 2 viruses are more closely related to each other than to group 3 viruses, which exhibit genome features distinct from those of the former two virus groups. In particular, the group 3 virus genome contains the shortest leader protein within Picornaviridae, a putative type I internal ribosome entry site (IRES) in the 5′-untranslated region instead of the type IV IRES found in group 1 and 2 viruses, one instead of two GXCG motifs in 2A, an L→V substitution in the DDLXQ motif in 2C helicase, and a conserved GXH motif in 3C protease. Group 1 and 2 viruses are unique among picornaviruses in having AMH instead of the GXH motif in 3Cpro. These findings suggest that the three picornaviruses belong to two novel genera in the Picornaviridae family. This report describes the discovery and complete genome analysis of three picornaviruses in bats, and their presence in diverse bat genera/species suggests the ability to cross the species barrier.

Comparative analysis of six genome sequences of three novel picornaviruses, turdiviruses 1, 2 and 3, in dead wild birds, and proposal of two novel genera, Orthoturdivirus and Paraturdivirus, in the family Picornaviridae

In this territory-wide molecular epidemiology study of picornaviruses, involving 6765 dead wild birds from 201 species in 50 families over a 12 month period, three novel picornaviruses, turdiviruses 1, 2 and 3 (TV1, TV2 and TV3), were identified from birds of different genera in the family Turdidae. In contrast to many other viruses in birds of the family Turdidae or viruses of the family Picornaviridae, TV1, TV2 and TV3 were found exclusively in the autumn and winter months. Two genomes each of TV1, TV2 and TV3 were sequenced. Regions P1, P2 and P3 of the three turdiviruses possessed, respectively, <40, <40 and <50 % amino acid identities with those of other picornaviruses. Moreover, P1, P2 and P3 of TV1 also possessed, respectively, <40, <40 and <50 % amino acid identities with those of TV2 and TV3. Phylogenetic analysis revealed that TV1, TV2 and TV3 were distantly related to members of the genus Kobuvirus. Among the three turdiviruses, TV2 and TV3 were always clustered together, with high bootstrap supports of 1000. The genomic features of TV2 and TV3 were also distinct from TV1, including lower G+C contents, shorter leader protein and a preference for codon sequence NNT rather than NNC for amino acids that can use either NNT or NNC as codons (P<0.001 by χ(2)-test). Based on our results we propose two novel genera, Orthoturdivirus for TV1, and Paraturdivirus for TV2 and TV3, in the family Picornaviridae. The type of internal ribosomal entry site for TV1, TV2 and TV3 remains to be determined.

Isolation and Characterization of a Novel Betacoronavirus Subgroup A Coronavirus, Rabbit Coronavirus HKU14, from Domestic Rabbits

ABSTRACT We describe the isolation and characterization of a novel Betacoronavirus subgroup A coronavirus, rabbit coronavirus HKU14 (RbCoV HKU14), from domestic rabbits. The virus was detected in 11 (8.1%) of 136 rabbit fecal samples by reverse transcriptase PCR (RT-PCR), with a viral load of up to 108 copies/ml. RbCoV HKU14 was able to replicate in HRT-18G and RK13 cells with cytopathic effects. Northern blotting confirmed the production of subgenomic mRNAs coding for the HE, S, NS5a, E, M, and N proteins. Subgenomic mRNA analysis revealed a transcription regulatory sequence, 5′-UCUAAAC-3′. Phylogenetic analysis showed that RbCoV HKU14 formed a distinct branch among Betacoronavirus subgroup A coronaviruses, being most closely related to but separate from the species Betacoronavirus 1. A comparison of the conserved replicase domains showed that RbCoV HKU14 possessed <90% amino acid identities to most members of Betacoronavirus 1 in ADP-ribose 1″-phosphatase (ADRP) and nidoviral uridylate-specific endoribonuclease (NendoU), indicating that RbCoV HKU14 should represent a separate species. RbCoV HKU14 also possessed genomic features distinct from those of other Betacoronavirus subgroup A coronaviruses, including a unique NS2a region with a variable number of small open reading frames (ORFs). Recombination analysis revealed possible recombination events during the evolution of RbCoV HKU14 and members of Betacoronavirus 1, which may have occurred during cross-species transmission. Molecular clock analysis using RNA-dependent RNA polymerase (RdRp) genes dated the most recent common ancestor of RbCoV HKU14 to around 2002, suggesting that this virus has emerged relatively recently. Antibody against RbCoV was detected in 20 (67%) of 30 rabbit sera tested by an N-protein-based Western blot assay, whereas neutralizing antibody was detected in 1 of these 20 rabbits.

Recent Transmission of a Novel Alphacoronavirus, Bat Coronavirus HKU10, from Leschenault's Rousettes to Pomona Leaf-Nosed Bats: First Evidence of Interspecies Transmission of Coronavirus between Bats of Different Suborders

ABSTRACT Although coronaviruses are known to infect various animals by adapting to new hosts, interspecies transmission events are still poorly understood. During a surveillance study from 2005 to 2010, a novel alphacoronavirus, BatCoV HKU10, was detected in two very different bat species, Ro-BatCoV HKU10 in Leschenaults rousettes (Rousettus leschenaulti) (fruit bats in the suborder Megachiroptera) in Guangdong and Hi-BatCoV HKU10 in Pomona leaf-nosed bats (Hipposideros pomona) (insectivorous bats in the suborder Microchiroptera) in Hong Kong. Although infected bats appeared to be healthy, Pomona leaf-nosed bats carrying Hi-BatCoV HKU10 had lower body weights than uninfected bats. To investigate possible interspecies transmission between the two bat species, the complete genomes of two Ro-BatCoV HKU10 and six Hi-BatCoV HKU10 strains were sequenced. Genome and phylogenetic analyses showed that Ro-BatCoV HKU10 and Hi-BatCoV HKU10 represented a novel alphacoronavirus species, sharing highly similar genomes except in the genes encoding spike proteins, which had only 60.5% amino acid identities. Evolution of the spike protein was also rapid in Hi-BatCoV HKU10 strains from 2005 to 2006 but stabilized thereafter. Molecular-clock analysis dated the most recent common ancestor of all BatCoV HKU10 strains to 1959 (highest posterior density regions at 95% [HPDs], 1886 to 2002) and that of Hi-BatCoV HKU10 to 1986 (HPDs, 1956 to 2004). The data suggested recent interspecies transmission from Leschenaults rousettes to Pomona leaf-nosed bats in southern China. Notably, the rapid adaptive genetic change in BatCoV HKU10 spike protein by ∼40% amino acid divergence after recent interspecies transmission was even greater than the ∼20% amino acid divergence between spike proteins of severe acute respiratory syndrome-related Rhinolophus bat coronavirus (SARSr-CoV) in bats and civets. This study provided the first evidence for interspecies transmission of coronavirus between bats of different suborders.