Genetic relatedness of the novel human group C betacoronavirus to Tylonycteris bat coronavirus HKU4 and Pipistrellus bat coronavirus HKU5

The recent outbreak of severe respiratory infections associated with a novel group C betacoronavirus (HCoV-EMC) from Saudi Arabia has drawn global attention to another highly probable “SARS-like” animal-to-human interspecies jumping event in coronavirus (CoV). The genome of HCoV-EMC is most closely related to Tylonycteris bat coronavirus HKU4 (Ty-BatCoV HKU4) and Pipistrellus bat coronavirus HKU5 (Pi-BatCoV HKU5) we discovered in 2006. Phylogenetically, HCoV-EMC is clustered with Ty-BatCoV HKU4/Pi-BatCoV HKU5 with high bootstrap supports, indicating that HCoV-EMC is a group C betaCoV. The major difference between HCoV-EMC and Ty-BatCoV HKU4/Pi-BatCoV HKU5 is in the region between S and E, where HCoV-EMC possesses five ORFs (NS3a-NS3e) instead of four, with low (31%–62%) amino acid identities to Ty-BatCoV HKU4/Pi-BatCoV HKU5. Comparison of the seven conserved replicase domains for species demarcation shows that HCoV-EMC is a novel CoV species. More intensive surveillance studies in bats and other animals may reveal the natural host of HCoV-EMC.

he recent outbreak of severe respiratory infections associated with a novel group C betacoronavirus (HCoV-EMC) from Saudi Arabia has drawn global attention to another highly probable ''SARSlike'' animal-to-human interspecies jumping event in coronavirus (CoV). The genome of HCoV-EMC is most closely related to Tylonycteris bat coronavirus HKU4 (Ty-BatCoV HKU4) and Pipistrellus bat coronavirus HKU5 (Pi-BatCoV HKU5) we discovered in 2006. Phylogenetically, HCoV-EMC is clustered with Ty-BatCoV HKU4/Pi-BatCoV HKU5 with high bootstrap supports, indicating that HCoV-EMC is a group C betaCoV. The major difference between HCoV-EMC and Ty-BatCoV HKU4/Pi-BatCoV HKU5 is in the region between S and E, where HCoV-EMC possesses five ORFs (NS3a-NS3e) instead of four, with low (31%-62%) amino acid identities to Ty-BatCoV HKU4/Pi-BatCoV HKU5. Comparison of the seven conserved replicase domains for species demarcation shows that HCoV-EMC is a novel CoV species. More intensive surveillance studies in bats and other animals may reveal the natural host of HCoV-EMC.
The recent outbreak of severe respiratory tract infections associated with a novel human group C betacoronavirus originating from Saudi Arabia has drawn global attention to another highly probable ''SARS-like'' interspecies jumping event of coronavirus (CoV) from animal to human. In June 2012, a novel CoV was isolated using Vero cells from the lung tissue of a 60-year-old resident of Saudi Arabia with fatal acute pneumonia and renal failure. In September 2012, another 49-year-old male resident of Qatar with severe acute pneumonia and renal failure and recent travel history to Saudi Arabia was admitted to an intensive care unit in Qatar. RT-PCR and sequencing of a short fragment of RNA-dependent RNA polymerase (RdRp) confirmed the presence of the same CoV as detected in the first Saudi Arabian case. 1 Complete genome sequencing of the virus isolated from the first patient was performed by Fouchier et al. at the Erasmus University Medical Centre, the Netherlands, and the sequence was released on September 28, 2012 (GenBank accession NO JX869059 and named as human betacoronavirus 2c EMC/2012). So far, there is no evidence of human-to-human transmission. The source of the virus remains obscure. In this article, this novel human group C betaCoV is abbreviated as HCoV-EMC.
After the SARS epidemic, we started to focus on CoV biodiversity, genomics and phylogeny and built up an evolutionary map of CoV evolution. Before 2003, there were less than 10 CoVs with complete genomes available, which include two human CoVs, human coronavirus 229E (HCoV-229E) and human coronavirus OC43 (HCoV-OC43). By September 2012, the number of CoVs with complete genomes sequenced had tripled. It includes two additional human CoVs, human coronavirus NL63 (HCoV-NL63) and human coronavirus HKU1 (HCoV-HKU1). 2,3 Traditionally, CoVs were classified into groups 1, 2 and 3. In 2011, the Coronavirus Study Group of the International Committee for Taxonomy of Viruses has re-classified these three groups of CoVs as three genera, Alphacoronavirus, Betacoronavirus and Gammacoronavirus; and we have discovered a fourth genus of CoV, Deltacoronavirus, which includes at least nine avian CoVs and a porcine coronavirus HKU15. 4,5 Within the betaCoVs, they are further subclassified into group A, including HCoV-HKU1, HCoV-OC43, bovine coronavirus (BCoV), sable antelope coronavirus, giraffe coronavirus, equine coronavirus, porcine hemagglutinating encephalomyelitis virus, murine hepatitis virus, rat coronavirus and rabbit coronavirus HKU14 (RbCoV HKU14); 6 group B, including the human and civet SARS-related CoVs (SARSr-CoV) and SARS-related Rhinolophus bat coronavirus (SARSr-Rh-BatCoV); 7,8 group C, including Tylonycteris bat coronavirus HKU4 (Ty-BatCoV HKU4) and Pipistrellus bat coronavirus HKU5 (Pi-BatCoV HKU5) we discovered in 2006; 9,10 and group D, including Rousettus bat coronavirus HKU9 (Ro-BatCoV HKU9). 10,11 In addition to Ty-BatCoV HKU4 and Pi-BatCoV HKU5, other group C bat betaCoVs should also be present, but their complete genome sequences are not available. 12,13 Based on the CoVs discovered, we have constructed a model of CoV evolution, with evidence supporting that bat CoVs are the gene source of alphaCoVs and betaCoVs and avian CoVs are the gene source of gammaCoVs and deltaCoVs. 5 All these works have laid down an evolutionary map for rapid phylogenetic and bioinformatics analyses of HCoV-EMC. The diversity of CoVs is a result T of the infidelity of RdRp which make CoV genomes especially plastic, a high frequency of homologous RNA recombination due to their unique random template switching during RNA replication, and their large genomes. In addition to biodiversity, a number of natural recombination and possible interspecies jumping events has also been documented in betaCoVs. 6,11,[14][15][16][17][18] For group A betaCoVs, molecular clock analysis has shown that HCoV-OC43 is a relatively recent zoonotic virus of bovine origin that emerged in around 1890 likely from bovine-to-human transmission. 17 We have also recently discovered RbCoV HKU14, closely related to other members of the species Betacoronavirus 1 including HCoV-OC43 and BCoV, with recombination events that may have played a role in interspecies transmission of these HCoV-OC43related viruses between human, cattle, rabbits, swine and horses. 6 Despite having circulated in humans for more than a century, HCoV-OC43 is also found to be continuously evolving, with the recent emergence of a novel genotype due to natural recombination. 15 For group B betaCoVs, SARSr-CoV is believed to be transmitted from civet to humans, although it is the horseshoe bat that was likely the primary host. 7,8 Civet SARSr-CoV was also likely a recombinant virus arising from different strains of SARSr-Rh-BatCoV from different geographical locations in China. 14,16 Although no interspecies transmission events have been documented in group D betaCoVs, we have also identified recombination events between different Ro-BatCoV HKU9 strains from different bat individuals, which may have allowed for the generation of different genotypes. 11 While these findings supported that betaCoVs have the propensity to recombine and cause interspecies transmission, such events were unknown in group C betaCoVs. As HCoV-EMC is most closely related to Ty-BatCoV HKU4 and Pi-BatCoV HKU5, it would be important to study their genetic relatedness, which may provide clues on whether bats are the possible origin as in SARSr-CoV.
The genome characteristics and organization of HCoV-EMC are similar to those of Ty-BatCoV HKU4 and Pi-BatCoV HKU5. Ty-BatCoV HKU4 was discovered from lesser bamboo bats (Tylonycteris pachypus) and Pi-BatCoV HKU5 was discovered from Japanese pipistrelles (Pipistrellus abramus) in Hong Kong. 9 Both lesser bamboo bats and Japanese pipistrelles are insectivorous microbats found in China and some other parts of Asia. The size of the genome of HCoV-EMC is 30 106 bases, slightly smaller than those of Ty-BatCoV HKU4 (30 286 to 30 316 bases) and Pi-BatCoV HKU5 (30 482 to 30 488 bases); and the G1C content is 41%, in between those of Ty-BatCoV HKU4 (38%) and Pi-BatCoV HKU5 (43%). The replicase ORF1ab occupies 21.5 kb of the genome. This ORF encodes 16 putative non-structural proteins, including nsp3 (which contains the putative papain-like protease (PL pro )), nsp5 (putative chymotrypsin-like protease (3CL pro )), nsp12 (putative RdRp), nsp13 (putative helicase (Hel)) and other proteins of unknown functions. These proteins are produced by proteolytic cleavage of the large replicase polyprotein by PL pro and 3CL pro at specific sites which are conserved with those in Ty-BatCoV HKU4 and/or Pi-BatCoV HKU5 (Table 1).
Novel human group C betacoronavious PCY Woo et al 2 103 of NS3d in Ty-BatCoV HKU4 and Pi-BatCoV HKU5 (35%-49% amino acid identities), with a stop codon UAG present at nucleotide position 27 160, leading to premature termination. NS3e of HCoV-EMC is homologous to amino acids 116/122 to 223/227 of NS3d in Ty-BatCoV HKU4 and Pi-BatCoV HKU5 (60%-62% amino acid identities). NS3c and NS3e do not possess any TRS or internal ribosomal entry site. BLAST search revealed no amino acid similarities between these putative non-structural proteins and other known proteins and no functional domains were identified by PFAM and Inter-ProScan. TMHMM and TMpred analyses show one and two putative transmembrane domains in NS3a (residues 9 to 29) and NS3d (residues 36 to 56 and 71 to 91), respectively. Similar to Ty-BatCoV HKU4 and Pi-BatCoV HKU5, the 3' untranslated region of the genome of HCoV-EMC contains predicted bulged stem-loop structures 16 to 76 nucleotides downstream of the N genes. Downstream to the bulged stem-loop structure, 97 to 121 nucleotides downstream of the N genes, a pseudoknot structure is present. Bootscan analysis did not show any recombination between HCoV-EMC, Ty-BatCoV HKU4 and Pi-BatCoV HKU5.
Using the sequences available at the moment and Yule process speciation under a relaxed clock model with an uncorrelated lognormal distribution, the mean evolutionary rate of betaCoVs was estimated at 2.37310 24 nucleotide substitutions per site per year for the RdRp gene. Molecular clock analysis using the RdRp gene showed that HCoV-EMC diverged from the most recent common ancestor of group C betaCoVs at ,year 941 (HPDs, 529 BC to 1878). Compared to the human and civet SARSr-CoV and SARSr-Rh-BatCoV cluster, the human/civet SARSr-CoV diverged from the Novel human group C betacoronavious PCY Woo et al 3

Emerging Microbes and Infections
Novel human group C betacoronavious PCY Woo et al 4 most recent common ancestor of the human/ civet SARSr-CoV and SARSr-Rh-BatCoV at ,year 1653 (HPDs, 1150 to 1968). By definition, the human and civet SARSr-CoV and SARSr-Rh-BatCoV are the same CoV species. These observations suggest that there should be one or more intermediate hosts between Ty-BatCoV HKU4, Pi-BatCoV HKU5 and HCoV-EMC. Sequencing more strains of Ty-BatCoV HKU4, Pi-BatCoV HKU5 and HCoV-EMC, as well as other group C betaCoVs collected at different time points, should be performed to achieve a more accurate estimation of the divergence time.
In the last decade, we have already witnessed the discovery of two novel human CoVs and an animal-to-human CoV interspecies jumping event on SARSr-CoVs. In contrast to HCoV-229E, HCoV-OC43, HCoV-NL63 and HCoV-HKU1, which are notoriously difficult to culture, HCoV-EMC and human SARS-CoV are both readily cultivable using primate cell lines. This may suggest a possible correlation between cultivability and virulence/recent interspecies jumping. Sequencing more genomes and performing evolutionary analysis will help us understand whether HCoV-EMC represent another recent interspecies jumping event from animal to human or another human CoV that has stably infected human. Our most recent findings showed that CoVs can be transmitted between two bat species of different suborders, suggesting that different degrees of interspecies jumping can occur in nature. 19 More intensive surveillance studies for group C betaCoVs in bats and other animals may reveal the natural host of this novel human group C betaCoV. As coronaviruses are prone to recombination and mutation and it has been documented that different levels of interspecies jumping can indeed occur in nature, we should not underestimate the potential of coronaviruses being the cause of another major ''SARS-like'' pandemic.