key: cord-0985294-p0kv1pht
authors: Yip, Chi Wai; Hon, Chung Chau; Shi, Mang; Lam, Tommy Tsan-Yuk; Chow, Ken Yan-Ching; Zeng, Fanya; Leung, Frederick Chi-Ching
title: Phylogenetic perspectives on the epidemiology and origins of SARS and SARS-like coronaviruses
date: 2009-09-30
journal: Infect Genet Evol
DOI: 10.1016/j.meegid.2009.09.015
sha: aef37fccc86113ee96ca502e37f8144b27673b0d
doc_id: 985294
cord_uid: p0kv1pht

Severe Acute Respiratory Syndrome (SARS) is a respiratory disease caused by a zoonotic coronavirus (CoV) named SARS-CoV (SCoV), which rapidly swept the globe after its emergence in rural China during late 2002. The origins of SCoV have been mysterious and controversial, until the recent discovery of SARS-like CoV (SLCoV) in bats and the proposal of bats as the natural reservior of the Coronaviridae family. In this article, we focused on discussing how phylogenetics contributed to our understanding towards the emergence and transmission of SCoV. We first reviewed the epidemiology of SCoV from a phylogenetic perspective and discussed the controversies over its phylogenetic origins. Then, we summarized the phylogenetic findings in relation to its zoonotic origins and the proposed inter-species viral transmission events. Finally, we also discussed how the discoveries of SCoV and SLCoV expanded our knowledge on the evolution of the Coronaviridae family as well as its implications on the possible future re-emergence of SCoV.

A B S T R A C T Severe Acute Respiratory Syndrome (SARS) is a respiratory disease caused by a zoonotic coronavirus (CoV) named SARS-CoV (SCoV), which rapidly swept the globe after its emergence in rural China during late 2002. The origins of SCoV have been mysterious and controversial, until the recent discovery of SARS-like CoV (SLCoV) in bats and the proposal of bats as the natural reservior of the Coronaviridae family. In this article, we focused on discussing how phylogenetics contributed to our understanding towards the emergence and transmission of SCoV. We first reviewed the epidemiology of SCoV from a phylogenetic perspective and discussed the controversies over its phylogenetic origins. Then, we summarized the phylogenetic findings in relation to its zoonotic origins and the proposed inter-species viral transmission events. Finally, we also discussed how the discoveries of SCoV and SLCoV expanded our knowledge on the evolution of the Coronaviridae family as well as its implications on the possible future re-emergence of SCoV.

ß 2009 Published by Elsevier B.V.

1. Introduction

Severe acute respiratory syndrome (SARS) is a highly contagious respiratory disease caused by a previously unknown coronavirus (CoV) named as SARS-CoV (SCoV) Peiris et al., 2003) . The major outbreak started from November 2002 as a rapid wave of epidemic from Guangdong Province of China and spread through 29 regions around the world (WHO, 2004a) . The epidemic was effectively controlled under vigorous quarantine measures and no new case was reported after July 2003. Six months after its disappearance, SCoV re-emerged in December 2003 as four sporadic cases in Guangdong Province, causing no fatality or secondary transmission (WHO, 2004b) . SCoV is probably one of the very few examples of zoonotic viral emergence 'caught-in-the-act' with adequate sequences sampled from different phases of the epidemic (CSMEC, 2004; Wang et al., 2005a) , as well as highly relevant samples from its zoonotic origins Song et al., 2005; Ren et al., 2006) . Taking advantage of the wealth of these sequence data, the evolutionary behaviors of SCoV during the epidemic have been investigated extensively, e.g. chains of transmission, theoretic time of epidemic onset, rate of substitutions and mode of natural selection acting on the viral genome, etc. In the first half of this article, applications of phylogenetics for investigating the emergence and transmission of SCoV were reviewed.

Prior to the discovery of SCoV, members of the Coronaviridae family can be unambiguously classified into three phylogenetic groups. Although initial phylogenetic analyses did not confidently classify SCoV as a member of any of the three existing groups of CoV, futher analyses suggested SCoV might share an ancient and distant ancestry with Group 2 CoVs (Snijder et al., 2003; Gorbalenya et al., 2004) . SCoVs were also isolated from small mammals such as civets in wet markets, suggesting these mammals may have been the direct zoonotic origin of SCoV . Excitingly, wider animal surveys revealed unanticipated high levels of genetic diversity of CoV in bats (Cui et al., 2007) . A group of CoVs named as SARS-like CoVs (SLCoVs) which shares 88-92% nucleotide identity with SCoV, was also identified from bats (Lau et al., 2005; Li et al., 2005b) . These findings lead to the hypothesis that bats are the natural reservoir of SCoV and other members of the Coronaviridae family (Tang et al., 2006; Vijaykrishna et al., 2007) . The biological and evolutionary aspects related to the conjectured inter-species transmission of SLCoV from its natural reservoir (i.e. bats) to the intermediate host (i.e. civets), and finally to human, have been the center of discussion (Shi and Hu, 2008) . In the second half of this article, the phylogenetic origins of SCoV, as well as the diversity of CoVs in its zoonotic sources and the phylogenetic aspects of the speculated inter-species transmission events, were reviewed.

2. Dissecting the epidemiology of SARS outbreaks from a phylogenetic perspective 2.1. Super-spreading 'caught-in-the-act'

Epidemiology of the SARS epidemic has been well-documented and one of its most intriguing characteristics is the concept of ''super-spreading events'' (SSE) (Lloyd-Smith et al., 2005) , which contributed significantly to the rapid spread of the disease locally and globally Shen et al., 2004; Chen et al., 2006; Zhong et al., 2003) . In mid-November 2002, the epidemic started as a series of seemingly independent cases and followed by the first documented SSE in Hospital HSZ-2 in Guangdong at the end of January 2003 (Zhong et al., 2003) . By then, an infected nephrologist traveled from Guangdong to Hong Kong and initiated another SSE in Hotel M at the end of February 2003, leading to the worldwide transmission of SARS and subsequent outbreaks in Hong Kong (Ruan et al., 2003; Guan et al., 2004) . The epidemic has been divided into three phases ( Fig. 1) based on the occurrence of the above two mentioned SSEs (CSMEC, 2004) . The early phase refers to the period prior to the SSE in Hospital HSZ-2 while the late phase refers to the period after to the SSE in Hotel M. The middle phase refers to the period between these two SSEs.

Phylogenetic analyses enabled researchers to better understand the transmission chains and trace the sources of viral epidemics, which provide important information for making public health policy. Broadly speaking, phylogenies of SCoV sequences generally agreed with the documented contact histories, and additionally provided evidence to support some uncorroborated epidemiological speculations (Zhao, 2007) . Phylogenies reconstructed by Neighbor Joining (NJ) (CSMEC, 2004; Guan et al., 2004) , Maximum Likelihood (ML) Tang et al., 2007) and Bayesian methods (Tang et al., 2009) consistently distinguished the late phase strains as a monophyletic cluster, which included the index patient of the SSE in Hotel M and the primary cases in Vietnam, Singapore and Canada (Guan et al., 2004) , supporting the viewpoint that the SARS outbreaks later in Hong Kong, as well as those in other parts of the world, were largely, if not completely, originated from a common source. In fact, analyses of the S gene sequences suggested multiple strains might have been independently introduced into Hong Kong from Guangdong before the SSE in Hotel M, although none of these strains contributed substantially to the subsequent outbreaks Guan et al., 2004) . On the other hand, the strains from early and middle phases are relatively more diverse and appeared as multiple distinct clusters on the phylogenies (CSMEC, 2004; Song et al., 2005; Wang et al., 2005b) , agreeing with the epidemiological investigations that SCoVs have been circulating in Guangdong and caused seemingly independent outbreaks prior to the SSE in Hotel M (Zhong et al., 2003) . This observation fits the model put forward by Antia et al. (2003) and implies the possible occurrence of multiple zoonotic transmissions of phylogenetically distinct, but yet similar, SCoVs to human in the early phase of the epidemic.

The epidemiology of some local outbreaks in the late phase was also investigated using phylogenetics. Firstly, two major subsequent outbreaks in Hong Kong, the Amoy Gardens outbreak (Ng, 2003) and the Hospital P outbreak (Tomlinson and Cockram, 2003) , were phylogenetically demonstrated to be directly linked to the SSE in Hotel M Guan et al., 2004) . Secondly, phylogenies suggested the Taiwan outbreaks were likely to be derived from multiple sources (Lan et al., 2005b; Shih et al., 2005) , including the Amoy Gardens outbreak Lan et al., 2005a) and the SSE in Hotel M (Yeh et al., 2004) . Thirdly, the Singapore outbreak was traced back to two separate initial introductions, but both were directly linked to the SSE in Hotel M (Ruan et al., 2003; Vega et al., 2004) , which later led to several other subsequent outbreaks within Singapore and a case in Germany . Lastly, in Beijing, although the first case was reported approximately a week after the SSE in Hotel M (Bi et al., 2003) , phylogenetic analyses suggested not all the cases in Beijing were related to the SSE in Hotel M and some of the cases were likely originated from the Guangdong prior SSE (Liu et al., 2005c; Zhao, 2007) . 1 . A phylogeny of spike gene nucleotide sequences from SCoV isolated from humans, civets and raccoon dogs. This phylogeny was modified and adopted from Lam et al. (2008a) . Sequences from humans, civets, raccoon dogs and bats were indicated with symbols &, *,~and^, respectively. The tree was constructed using ML method, with confidences of topology summarized from 5000 trees sampled from ML and NJ bootstrap replicates and BMCMC samples. Only confidence values of major clusters were shown (ML/NJ/BMCMC, in the parenthesis). The human epidemic cluster (2002) (2003) was divided into late, early and middle phases according to a previous study (CSMEC, 2004) . Accession numbers of the sequences are shown within round brackets after their strain names (in bold). The distance unit was substitutions/site. Rp3 isolated from bats (^) was used as an out-group to root the tree, and the genetic distance of its branch is not shown.

Robust estimations on the theoretical onset of a viral epidemic provide evidence to testify the hypotheses on viral origins, e.g. AIDS pandemic (Korber et al., 2000; Lemey et al., 2003 Lemey et al., , 2004 Worobey et al., 2004) . Under the assumption of adequate sampling, the theoretical onset of a viral epidemic can be inferred as the time of the most recent common ancestor (tMRCA) of all sampled sequences. tMRCA can be estimated using a range of phylogenetic methods (Drummond et al., 2003) , but all methods are primarily based on the assumption of molecular clocks of different extents (Pybus, 2006) . It is noted that this tMRCA should theoretically be older than the earliest documented case, given the fact that there must be a window period between the emergence and recognition of the epidemic and the length of this window period is often a significant epidemiological interest (Stumpf and Pybus, 2002; Worobey et al., 2004) .

In the case of SCoV, linear regression methods which correlate divergence with sampling dates were commonly used to estimate the tMRCA in the early studies. Based on simple linear regression, Zeng et al. (2003a) Zhao et al. (2004) adopted three strategies to minimize possible errors and inferred the oldest bound of this tMRCA as early as spring 2002. A least square method, evaluated using Monte Carlo simulations, has also been applied to estimate the onset of the epidemic as August to September 2002 (Lu et al., 2004) . While all the above studies relied on the assumption of a constant evolutionary rate over all lineages, i.e. a strict molecular clock (Pybus, 2006) , enforcing a molecular clock to non-clocklike datasets could lead to biased estimations (Yoder and Yang, 2000) . More recently, Salemi et al. (2004) applied a ML method to testify the clocklike behavior of a dataset of ten taxa and placed this tMRCA between September and November of 2002.

Despite the differences in methods and datasets, the above studies generally inferred the oldest bound of the onset of the epidemic at around mid-late 2002. The earliest documented case of SARS was identified on 16th November 2002 (Zhong et al., 2003) , which is remarkably close to the oldest bounds of its theoretical onset when compared with other well-known viral epidemics (Stumpf and Pybus, 2002; Worobey et al., 2004) . This finding suggests the SCoV has been quickly recognized after its emergence in human, partially reflecting the explosive mode of transmission in early phase of the epidemic.

After the discovery of SCoV in civets and sporadic re-emergence of SCoV in December 2003, these SCoV strains have been incorporated into datasets for estimation of their tMCRA (Song et al., 2005; Vijaykrishna et al., 2007; Hon et al., 2008) . This tMRCA represents the upper limit of the time (i.e. oldest bound) of the inter-species transmission of SCoV from civets to human, which is phylogenetically different from the tMRCA of all human SCoVs discussed above and should theoretically be older. The estimated time of this inter-species transmission event will be discussed later in Section 4.

In December 2003, four seemingly independent cases of SCoV infection were reported (Liang et al., 2004; WHO, 2004b) . All the patients have direct or indirect contact history with wild animals (Song et al., 2005; Wang et al., 2005a) . Phylogenetic analyses of human SCoVs from these sporadic cases and civet SCoVs collected from the same period formed a monophyletic group distinct from the SCoVs of the 2002-2003 epidemic ( Fig. 1 ), suggesting these sporadic cases were likely caused by inter-species transmissions that were independent from the previous outbreak (Song et al., 2005; Lam et al., 2008a) . These epidemiological and phylogenetic findings provided convincing evidence for direct transmission of SCoV from civets to human (discussed later in Section 4).

In summary, the above findings demonstrated the vital roles of phylogenetics in understanding the emergence and transmission of SCoV in such a short-lasting but sweeping epidemic. The adequacy of early phase sampling and the relatively short duration of the SARS epidemic make it an excellent textbook example on how phylogenetic analyses can be used as an auxiliary tool in combination with classical epidemiological investigations to study the emergence of a viral epidemic.

Which of the existing groups of CoV is phylogenetically closest to SCoV? This question has been under the spotlight of most discussions because the answer might help to trace the zoonotic origin of SCoV that is fundamentally important to public health. Prior to the availability of its complete genome sequences, initial phylogenetic analyses of a fragment of ORF1 suggested that SCoV may represent a novel group that is independent to the other three existing groups Ksiazek et al., 2003) . Similar conclusions were reached based on phylogenetic analyses on multiple viral proteins after the complete genome sequences were available (Marra et al., 2003; Rota et al., 2003; Zeng et al., 2003b) . Based on the observation that SCoV appears to be genetically equidistant to other known CoV groups, Holmes and Enjuanes (2003) concluded that SCoV is neither a recent host-range mutant of a known CoV nor a recent recombinant between known CoVs, but it probably evolved separately from an ancestor of the known CoVs in an unidentified host for a remarkably long period of time before its emergence in human.

Although the above argument was generally well-received, two follow-up questions were then raised. The first question is: If SCoV represents a lineage that anciently diverged from an ancestor of the existing CoVs, would this ancestor belong to one of the existing CoV groups? It later led to the concept of ''early split-off from Group 2 CoV'' for describing the phylogenetic origin of SCoV (Snijder et al., 2003) . The second question is: Is SCoV an ''ancientrecombinant'' (if not a ''recent-recombinant'') between the ancestors of any existing CoV groups? Due to relatively high divergence between SCoV and other CoVs as well as the lack of a robust out-group (Holmes and Rambaut, 2004) , the conclusions derived from various recombination analyses were yet intricate and inconsistent (Gorbalenya et al., 2004) . A number of reasons were then proposed to explain the observed phylogenetic incongruence might not truly reflect the recombinant origin of SCoV (Holmes and Rambaut, 2004) . The following sessions summarized the findings from numerous studies regarding the controversial phylogenetic origins of SCoV.

Following the initial analyses (Marra et al., 2003; Rota et al., 2003; Zeng et al., 2003b) , the phylogenetic position of SCoV have been vigorously reevaluated with various approaches that were mainly different in their choices of genome regions Zhu and Chen, 2004; Kim et al., 2006) and rooting strategies (Snijder et al., 2003; Gibbs et al., 2004; Lio and Goldman, 2004) . Inspired by the positional conservation of cysteine residues between the S1 region of SCoV and Group 2 CoVs, Snijder et al. (2003) constructed an unrooted NJ phylogeny based on the amino acid sequences of S1 region and concluded SCoV is ''closely related'' to Group 2 CoVs . In two later studies, concatenated amino acid alignments of multiple viral proteins were used to construct unrooted phylogenies using various methods and the results consistently demonstrated the monophyletic relationship between SCoV and Group 2 CoVs Kim et al., 2006) . In addition to these unrooted phylogenies, rooted phylogenies based on ORF1 alignments reached similar conclusions that SCoV and Group 2 CoV shared the last common ancestor (Snijder et al., 2003; Gibbs et al., 2004; Lio and Goldman, 2004) . However, depends on the choices of out-groups, the rooting positions in these studies were inconsistent, which might be related to the possibly compromised accuracy of alignment if the out-groups are too diverged to be aligned (Holmes and Rambaut, 2004) .

Nonetheless, the consensus of above findings is that SCoV could be classified as a ''distant'' member of Group 2 CoVs (Gorbalenya et al., 2004) . This conclusion was reached primarily based on the robustness of the monophyletic relationship between SCoV and Group 2 CoVs. However, considering the actual genetic distance between SCoV and Group 2 CoVs, which is comparable to the intergroup distance among the three CoV groups, the biological significance of classifying SCoV as a distant member of Group 2 CoV is relatively low at best. We believe classification of a novel CoV should not be merely based on the branching orders of phylogenies and quantitative measurement of genetic distance should also be taken into account.

On the other hand, reconstructing the evolutionary relationships between highly divergent taxa has been proven as a difficult phylogenetic task (Philippe and Laurent, 1998) . In particular, the observed branching order of the deep nodes in the highly divergent CoV phylogenies may not reflect their true evolutionary history, due to the possible influences of the long-branch attraction and rate variation among lineages (Gribaldo and Philippe, 2002; Holmes and Rambaut, 2004) . In fact, Holmes and Rambaut reevaluated the tree topologies depicting the three possible phylogenetic positions of SCoV using Kishino-Hasegawa test, which is not biased by rate variations, and that demonstrated the topology with the best likelihood is not necessarily significantly better than the other two topologies, even the deep nodes were well-supported by quartet puzzling support values (Holmes and Rambaut, 2004) . This finding indicated careful interpretations are needed before reaching a firm conclusion from the branching order of deep phylogenies.

Until recently, the addition of several newly discovered CoV lineages to the phylogeny broke up some of those long branches (Holmes and Rambaut, 2004) , which tends to distribute the convergent and parallel mutations more evenly across the tree and hence reduce the problem of long-branch attraction (Hillis, 1996) . In Fig. 2 , the closer phylogenetic relationship between SCoV and Group 2 CoVs suggested from the previous studies is also supported in this relatively well-sampled phylogeny. To this end, the question of whether it is appropriate to classify SCoV as a distant member of Group 2 CoVs seems to be a taxonomic problem more than a virological one, since SCoV and Group 2 CoVs are so divergent that they are likely to possess a substantial number of unique biological characteristics, e.g. host ranges and receptor usages (Haijema et al., 2003; Li et al., 2003) . With the recent discoveries of novel CoV lineages, the classification of coronaviruses has to be revised systematically in both biological and phylogenetic context.

Phylogenetic incongruence, which is often justified by incompatible tree topologies for different genome regions (Holmes et al., 1999) , has been widely used as an indicator for homologous recombination among viral genomes (Posada et al., 2002) . Several studies demonstrated phylogenetic incongruence within the genome of SCoV based on a wide range of methods (Rest and Mindell, 2003; Stanhope et al., 2004; Stavrinides and Guttman, 2004; Zhang et al., 2005b) . Despite the statistical significance of the phylogenetic incongruence observed in these studies, we could not conclude a consensus pattern of recombination since their findings were generally inconsistent. As an example, the potential recombination events within S gene proposed by Zhang et al. (2005b) were different from those proposed by Stavrinides and Guttman (2004) , and were undetectable from the study of Rest and Mindell (2003) . In fact, the considerable divergence between SCoV and the existing CoV groups has already excluded the possibility that SCoV is a recent recombinant from the existing CoV groups. Alternatively speaking, if recombination events had occurred, they have to be the ancient ones (Bosch, 2004) , i.e. both the parents and daughter could have been evolved considerably after the recombination events. A simulation study demonstrated the accuracy for detecting these ancient recombination events can be substantially diminished if they were significantly obscured by the subsequent post-recombination substitutions . Therefore, the intricacy of the findings on the phylogenetic incongruence within the genome of SCoV may be a reflection of the varying sensitivity of the detection methods on the ancient and obscured recombination events, if any.

To further investigate the reported phylogenetic incongruence within the genome of SCoV, Rambaut and Holmes re-evaluated the study of Stavrinides and Guttman (2004) and suggested the patterns cited as evidence for recombination are more probably caused by a variation in substitution rate among lineages (Holmes and Rambaut, 2004) . In addition, the author also stated the effect of the long-branch attraction on the branching order of deep phylogenies may also be a source of artifact, although unproven. Moreover, even if recombination did not occurred, given the stochastic nature of evolution, the authors would expect to observe phylogenetic incongruence among small genome fragments of a set of divergent taxa like SCoV and other CoVs. Therefore, the observed phylogenetic incongruence among the highly divergent genomes of CoV and its possible indication on ancient recombination events should be interpreted with extra cautions. Up to this point, the current phylogenetic evidence supporting the recombinant history of SCoV is weak at best (Holmes and Rambaut, 2004) .

Putting aside the intricate phylogenetic evidence, the presence of stem-loop II motif (S2m) in the genome of SCoV have also been taken as an indication for recombination (Marra et al., 2003) . S2m is a conserved RNA motif present in the genomes of several members of Astroviridae, Coronaviridae, and Picornaviridae family (Jonassen et al., 1998) while SCoV and Group 3 CoVs are the only members of the Coronaviridae that posse the motif. Assuming the motif in SCoV and Group 3 CoV were not acquired independently, the explanation of the co-presence of S2m should be either, (1) SCoV and Group 3 CoVs share the same ancestry, or (2) SCoV have acquired the motif from Group 3 CoVs through recombination (or vice versa). Currently, the phylogenetic data does not provide a better support for any of the above scenarios but a wider survey for the presence of S2m in other unknown CoV lineages will certainly provide insights to the ancient evolutionary history of the Coronaviridae family.

Since the early SARS cases in Guangdong seems to be related to restaurant workers handling wild mammals, Guan et al. (2003) surveyed the wild animals in a local market and isolated CoVs from Himalayan palm civets, which share 99.8% genome sequence identity to the SCoVs in human. Initial phylogenetic analyses suggested SCoVs from human and civets formed two distinct clades , but the SCoVs from civets are phylogenetically closer to the SCoVs of the early-phase epidemic than to those in the late-phase epidemic (Kan et al., 2005) . These findings strongly suggest civets were the immediate sources of the SCoVs leading to the earliest SARS cases in Guangdong. The role of civets as the immediate zoonotic source of the SARS epidemic was further revealed during the sporadic re-emergence of SCoV in December 2003 (Wang et al., 2005a) , based on the observation that the SCoVs isolated from civets and those patients of the same period formed a monophyletic cluster (Song et al., 2005; Lam et al., 2008a) . These findings suggested that the emergence of SCoV in human is likely to be resulted from direct transmissions of SCoVs from civets. In addition, phylogenetic analyses demonstrated that the SCoVs from the 2002-2003 epidemic and the 2003-2004 reemergence are phylogenetically distinct, suggesting the interspecies transmission events from civets to human in the two outbreaks might be independent (Song et al., 2005; Wang et al., 2005a) .

However, a large-scale survey of SCoVs of civets in market in China suggested the lack of widespread infections in civets (Kan et al., 2005) . In addition, the genetic diversity of civet SCoVs was relatively limited and was comparable to that of human SCoVs Fig. 2 . A phylogeny of all known CoVs (n = 40). The phylogeny was constructed based on the amino acid sequences of the RNA-dependent RNA Polymerase region (length = 163 a.a.). The phylogeny was constructed using BEAST (Drummond and Rambaut, 2007 ) under a uncorrelated lognormally relaxed clock model (Drummond et al., 2006) . The number at the nodes indicates the Bayesian posterior probability support (as percentages) summarized from trees sampled at every 1000th step of a BMCMC chain of 10,000,000 steps where values lower than 80% were not shown. The mean substitution rate was fixed at 1.0 and the branch length was expressed in units of substitutions per site. (Song et al., 2005; Lam et al., 2008a) . These findings suggested civets might not be the natural reservoir of SCoV and it only acquired SCoVs shortly before the emergence of SCoV in human. Based on phylogenetic analyses using dynamic homology, Janies et al. (2008) speculated that civets might also acquired SCoV from other species, possibly human, even after the emergence of SCoV in 2002. An experimental evidence suggesting civets may not be the natural reservoir of SCoV is the observed symptoms in civets experimentally infected with SCoV , based on the fact that the natural reservoir hosts usually do not display severe signs of infection (Hudson et al., 2002) . The above observations lead to the speculation of another natural reservoir host, which harbors a diverse group of SCoV-related CoVs and transmitted SCoV to civets prior to the emergence of SCoV in human.

As a result of extensive searches for the natural reservoir of SCoV, two groups of researchers independently identified a diverse group of CoVs from various species of horseshoe bats (Rhinolophus spp.). These CoVs shared 87-92% genome nucleotide identity with SCoVs and formed a distinct monophyletic cluster with SCoVs, therefore they were named SARS-like CoV (Lau et al., 2005; Li et al., 2005b) . The close evolutionary relationship between SCoV and SLCoV is also supported by the presence of the S2m motif in their 3 0 UTR (Tang et al., 2006; Shi and Hu, 2008) . Based on the phylogenetic analyses of the four characterized full genomes of SLCoV in horseshoe bats, the remarkably high genetic diversity of SLCoVs strongly suggested horseshoe bats are the natural reservoir of SCoV and SLCoV (Ren et al., 2006) . This concept is further supported by the relatively high prevalence of SLCoV in R. sinicus (Lau et al., 2005) and the geographically widespread infections of SLCoVs in bats from distinct locations in China . The current hypothesis is that civets might have acquired SLCoVs from horseshoe bats and transferred to human, which is consistent with the observation that the 29-nt deletion in ORF8 are retained in bat SLCoVs, civet SCoVs and early phase human SCoV (Lau et al., 2005) .

Although the hypothesis of horseshoe bats as the natural reservoir for SLCoV and SCoV is relatively sensible, how bats SLCoV were transmitted to civets is still unexplained. If civets acquired SLCoV from bats shortly prior to its emergence in human as proposed, this bat SLCoV strain should be genetically very similar to the SCoVs sampled from civets. However, based on the relatively distant phylogenetic relationship between SCoVs and SLCoVs, none of the currently sampled SLCoVs in bats is the descendant of the direct ancestor of SCoVs in human and civet (Hon et al., 2008) . In particular, Li and coworkers pointed out that substantial genetic changes in the S protein of the currently sampled SLCoV are likely to be necessary for the virus to infect civets or human . Therefore, the direct ancestor of SCoVs in human and civets remains elusive.

More recently, Hon et al. (2008) demonstrated significant phylogenetic discordances among different genome regions of SLCoV strain Rp3 and speculated its potentially recombinant origin. Phylogenetic analysis of the parental regions of Rp3 genome suggested the presence of an uncharacterized bat SLCoV lineage (i.e. HB-SLCoV in Fig. 3 ) that is phylogenetically closer to SCoVs than any of the currently sampled bat SLCoVs. Based on the relatively high genetic diversity among the currently sampled SLCoVs in bats, the existence of a phylogenetically distinct lineage of SLCoV not yet sampled is highly possible. Thus, the authors speculated that the direct ancestor of SCoVs was as a descendant derived from this not yet sampled lineage, which crossed from a horseshoe bat species to civets (Hon et al., 2008) .

Determining the time of inter-species transmission events might help us to comprehend the viral zoonosis of the virus from an evolutionary standpoint. The oldest bound of the inter-species transmission of SCoV from civets to human is theoretically Fig. 3 . A time-scaled phylogeny of SCoV and SLCoV. This phylogeny was modified and adopted from Hon et al. (2008) . The phylogeny was summarized from all MCMC phylogenies of the Orf1 data set analyzed under a Bayesian relaxed clock model. Height of the nodes was represented by the median of its estimates. The window period between the cross-species event and the onset of SARS epidemic was indicated as a dotted line. In the taxa labels, H, C and B represent host of human, civets and bats, respectively. correspondent to the tMRCA of all human and civet SCoVs. Firstly, Song et al. (2005) estimated the synonymous substitution rate of SCoVs using linear regression and placed this tMRCA to be around early December 2002 without providing CI. In a later study, Vijaykrishna et al. (2007) investigated this tMRCA by applying a Bayesian relaxed clock model (Drummond et al., 2006) to a phylogeny of all representatives from the CoV family and estimated it at around 1999 with 95% posterior bounds of 13 years (i.e. 1990-2003) . More recently, Hon et al. (2008) reconsidered this tMRCA using various clock models and estimated it at around September 2002 with 95% CIs between January and December of 2002. This tMRCA is very close to the observed first case of SARS (16th November 2002) as well as the estimated onset of the epidemic (as discussed earlier), suggesting SCoVs might have crossed from civets to human just months before the outbreak, supporting the view that civets are the immediate zoonotic source of SCoVs in human. It also implies the civet SCoVs might have adapted quickly, or alternatively speaking, a 'by-pass' host with only minimal adaptation is needed, to establish a sustainable chain of transmission in human.

On the other hand, the time of the speculated inter-species transmission of SLCoV from bats to civets has also been investigated. In the study described earlier , the oldest bound of the time of this event was first estimated at a mean of 1986 with a 38-year credible interval (i.e. 1964-2002) . Later, Hon et al. (2008) employed the concept of estimating the period between time of divergence (tDIV) and tMCRA proposed by Lam et al. (2008b) , and speculated this inter-species transmission event might have happened with a median of 4.08 years before onset of the epidemic (credible intervals of 1.45-8.84 years) (Fig. 3) . The above two estimates are not contradictory since the later estimate falls within the credible intervals of the former estimate but with an improved precision. Based on this relatively short window period between the inter-species transmission event and onset of the epidemic, Hon et al. (2008) speculated that civets might have acquired the ancestor of SCoV from the host species of SLCoV strain Rp3 directly and the involvement of other intermediate species may be unlikely. Therefore, authors suggested more focused surveillance on the host species of SLCoV strain Rp3 may shed light on the zoonotic origin of the direct ancestor of SCoV in civets.

The host specificity of CoVs is mainly determined by the binding between the spike (S) protein and its cellular receptors (Haijema et al., 2003) . Angiotensin I-converting enzyme 2 (ACE2) is a functional receptor for SCoVs in both human and civets (Li et al., 2005c) , and was demonstrated to interact with the receptor binding domain (RBD) in the S1 subunit of the S protein . Based on the observed mutations within the RBD in strains from different epidemic phases and hosts, Li et al. (2005c) demonstrated that binding of S1 subunit to human and civet ACE2 can be significantly altered by mutating only two residues on the RBD (residue 479 and 487), suggesting these two residues may contribute to the adaptation of SCoV from civet to human. Later, this viewpoint was further supported by the location of these two residues at the binding interface between RBD and ACE2 in crystal structures Li, 2008) .

As observed in other examples of viral host shifts (Parrish and Kawaoka, 2005) , the molecular determinants for host specificity, e.g. RBD residues related to adaptation to human in the case of SCoV, are likely to be subjected to an elevated level of selection pressure during the acquisition of a new host. The ratio of nonsynonymous (dN) to synonymous substitution rate (dS), i.e. v, which is widely used as a measure for selection pressure (Yang and Nielsen, 2002) , has been employed in several studies to detect positive selection on viral genes of SCoV and SLCoV. Firstly, the v values of structural genes were found to be higher than that of the non-structural region, suggesting the structural proteins might have been under a stronger selection pressure Song et al., 2005) . Moreover, the v value of S gene in the strains from early phase is significantly larger than that of those from middle and late phases (CSMEC, 2004) . The above findings provided preliminary phylogenetic evidence for the potential role of S protein in adaptation of SCoV from civets to human.

Additionally, lineage-specific v in the S gene phylogeny (n = 11) was estimated using a codon-based genetic algorithm (Kosakovsky Pond et al., 2006) , and a significantly higher v value was observed along the lineage leading to the human cluster . More recently, Tang et al. (2009) comprehensively analyzed a larger dataset (n = 59) using a ML branch-site codon model On the other hand, positively selected residues on S protein have been identified in a number of similar studies, primarily by applying ML codon models to similar datasets but with different epidemic groupings Shi et al., 2006) . Although the residues identified from these studies are not completely overlapping, probably due to the differences in taxa groupings, these sites are mainly located in the S1 subunit and residue 479 is consistently detected to be under positive selection, supporting its speculated roles in adaptation for new hosts (Song et al., 2005) . It should be noted that positively selected residues were not identified in the S genes from neither the late epidemic phase nor the SLCoV in bats . The mode of selection relevant to the inter-species transmission from bats to civets has not been investigated systematically, primarily due to the lack of bat SLCoV S gene sequences that are phylogenetically close enough for robust analysis.

As a result of the extensive efforts in searching for the zoonotic origins of SCoV, our knowledge on the host range and diversity of CoVs has been expanded rapidly. For example, in the last couple of years, the number of avian species detected to harbor Group 3 CoVs has been doubled (Cavanagh, 2005) . Furthermore, in addition to the two newly discovered human CoVs in the existing Group 1 and 2, i.e. NL63 (Fouchier et al., 2004) and HKU1 (Woo et al., 2005a) , respectively, a number of divergent CoVs have been identified from Asian Leopard Cats and Chinese Ferret Badgers (Dong et al., 2007) , wild birds (Woo et al., 2009 ) and a beluga whale (Mihindukulasuriya et al., 2008) . These findings largely expanded the known diversity of CoVs and the long-established classification of the Coronaviridae family as three distinct groups needs to be revised systematically. On the other hand, a number of diverse CoVs belonging to Group 1 and 2 have been identified from a range of bat species (Woo et al., 2006) . The surprisingly high diversity of bat CoVs rationally leads to the hypothesis that bats are the natural reservoir of all CoVs (Tang et al., 2006) . According to a time-scaled phylogeny of representative CoVs from all groups, Vijaykrishna et al. (2007) concluded that bats are likely to be the host of the ancestor for all presently known CoV lineages. Furthermore, based on the results of their Bayesian coalescent analyses, the authors speculated a diverse group of CoVs may be endemic in various bat species, with repeated introduction to other animals and occasional establishment of new lineages in other species . In fact, other than the speculated inter-species transmission of CoVs from bats to other animal species, hosts shifts of bat CoVs between different Rhinolophus spp. were also proposed based on the incongruence between the phylogenies of CoVs and their host Rhinolophus spp. (Cui et al., 2007) .

Additionally, inter-species transmissions of CoVs between other non-bat species were also proposed. The first animal-human zoonotic pair of CoVs being detailedly analyzed was Bovine CoV and Human CoV-OC43 (Vijgen et al., 2005) . Based on a combination of molecular clock analyses, the authors estimated the tMRCA of these two CoVs at around 1890 and speculated an inter-species transmission event from Bovine CoV to human might have occurred around this period (Vijgen et al., 2005) . Moreover, a number of CoVs are documented to infect multiple closely related host species, e.g. bovine CoVs has been isolated from captive wild ruminants (Alekseev et al., 2008) ; closely related Group 3 CoVs have been isolated from various avian species (Cavanagh et al., 2002; Jonassen et al., 2005; Liu et al., 2005b) . These findings implied the relatively promiscuous nature of the host specificity of CoVs, which was best demonstrated by the evidence of interspecies transmission of highly similar if not identical SCoVs between two distantly related species, i.e. human and civets (Wang et al., 2005a) .

Last but not least, the relatively high rate of homologous recombination has been speculated to facilitate the inter-species transmissions of CoVs (Baric et al., 1995 (Baric et al., , 1997 . Despite the lack of relevant examples to support this hypothesis, naturally occurring recombination between CoV strains of the same species (Jia et al., 1995) , as well as between strains of different CoV species, have been documented (Herrewegh et al., 1998; Hon et al., 2008; Decaro et al., 2009) , suggesting the generation of CoV diversity through recombination has been happening in the field. In fact, the two recently discovered human CoVs were proposed to have a recombinant history based on the observed phylogenetic incongruence between different genome regions (Woo et al., 2005b; Pyrc et al., 2006) . As discussed earlier, due to the relatively high divergence between different groups of CoV, cautions have to be taken to interpret these phylogenetic incongruences as evidence for ancient recombination events (Holmes and Rambaut, 2004; Chan et al., 2006) . In summary, the current data only suggests the occurrence of recombination between closely related CoV strains, but provides no direct evidence to support a role of recombination in the emergence of novel CoVs in novel host species, e.g. emergence of SCoV in human.

The SARS epidemic in 2003 offers a solid lesson on the application of phylogenetics in understanding the epidemiology of a newly established epidemic, as well as the evolutionary basis of a zoonotic viral emergence. Phylogenetics has played an indispensable role in the prompt identification of the transmission chains and its zoonotic origins, which provided important clues for the policy-making in public health, e.g. customs and border control, quarantine measures, culling of civets and the continuous search for the viral natural reservoir. These valuable experiences should help the community to better prepare for the next zoonotic viral epidemic.

More importantly, while most of the attentions have been focused on preparing for the known potential pandemics like avian influenza (Fauci, 2006) , the unanticipated strike of SARS epidemic alarmed public health officials and researchers for the neglected possibility of deadly viral emergence from an unexpected origin. For many years, CoVs have long been regarded as relatively ''mild'' viruses with broad but yet restricted host ranges. However, according to our current understanding to the diversity and evolution of CoVs , at least five divergent species of CoV are known to have zoonotic transmission into the human population and this cross species event will be likely to continue, and the zoonosis of SCoV was just the consequence of one of these inter-species transmission events. Although the prevention of zoonosis from an unexpected origin seems to be impractical, sufficiently flexible and stringent surveillance strategies provide us an opportunity to anticipate the disease in a population and prevent its further spreading by implementing appropriate control measures. The sporadic re-emergence of SCoV in the early 2004 illustrated the importance of stringent surveillance in preventing the further spread of a zoonotic virus in the early phase. One step further, the surveillance strategies must be adopted in a way that we could learn more about the diversity of potentially zoonotic viruses in their reservoirs. With CoV as an example, given the relatively promiscuous nature of its host specificity as discussed above, extensive and regular surveillance of known CoVs in pets and agricultural animals, which are in close contact with human, should be set up. Wellestablished knowledge on the diversity of potentially zoonotic viruses certainly accelerates the identification of its animal origin once it emerges in human. Last but not least, although the immediate zoonotic source of SCoV seems to be eliminated by culling of civets in Mainland China, the uncertainties about the diversity of SLCoV in bats still pose threat on the recurrence SCoV (Hon et al., 2008) . Therefore, we emphasize the importance of continuous surveillance on the genetic diversity of SLCoV in bats.

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

The role of evolution in the emergence of infectious diseases

High recombination and mutation rates in mouse hepatitis virus suggest that coronaviruses may be potentially important emerging viruses

Episodic evolution mediates interspecies transfer of a murine coronavirus

SARS: an amalgam of avian and mammalian viruses?

Coronaviruses in poultry and other birds

Coronaviruses from pheasants (Phasianus colchicus) are genetically closely related to coronaviruses of domestic fowl (infectious bronchitis virus) and turkeys

Detecting recombination in evolving nucleotide sequences

Understanding the super-spreading events of SARS in Singapore

Genomic characterisation of the severe acute respiratory syndrome coronavirus of Amoy Gardens outbreak in Hong Kong

Molecular epidemiology of SARS-from Amoy Gardens to Taiwan

Molecular evolution of the SARS coronavirus during the course of the SARS epidemic in China

Evolutionary relationships between bat coronaviruses and their hosts

Recombinant canine coronaviruses related to transmissible gastroenteritis virus of Swine are circulating in dogs

Detection of a novel and highly divergent coronavirus from asian leopard cats and Chinese ferret badgers in Southern China

Identification of a novel coronavirus in patients with severe acute respiratory syndrome

Inference of viral evolutionary rates from molecular sequences

Relaxed phylogenetics and dating with confidence

BEAST: Bayesian evolutionary analysis by sampling trees

Phylogeny of the SARS coronavirus

Pandemic influenza threat and preparedness

A previously undescribed coronavirus associated with respiratory disease in humans

The phylogeny of SARS coronavirus

Severe acute respiratory syndrome coronavirus phylogeny: toward consensus

Ancient phylogenetic relationships. Theor

Molecular epidemiology of the novel coronavirus that causes severe acute respiratory syndrome

Isolation and characterization of viruses related to the SARS coronavirus from animals in southern China

Switching species tropism: an effective way to manipulate the feline coronavirus genome

Feline coronavirus type II strains 79-1683 and 79-1146 originate from a double recombination between feline coronavirus type I and canine coronavirus

Inferring complex phylogenies

Viral evolution and the emergence of SARS coronavirus

Phylogenetic evidence for recombination in dengue virus

Virology. The SARS coronavirus: a postgenomic era

Evidence of the recombinant origin of a bat severe acute respiratory syndrome (SARS)-like coronavirus and its implications on the direct ancestor of SARS coronavirus

The Ecology of Wildlife Disease

Evolution of genomes, host shifts and the geographic spread of SARS-CoV and related coronaviruses

A novel variant of avian infectious bronchitis virus resulting from recombination among three different strains

A common RNA motif in the 3' end of the genomes of astroviruses, avian infectious bronchitis virus and an equine rhinovirus

Molecular identification and characterization of novel coronaviruses infecting graylag geese (Anser anser), feral pigeons (Columbia livia) and mallards (Anas platyrhynchos)

Molecular evolution analysis and geographic investigation of severe acute respiratory syndrome coronavirus-like virus in palm civets at an animal market and on farms

Close relationship between SARS-coronavirus and group 2 coronavirus

Timing the ancestor of the HIV-1 pandemic strains

GARD: a genetic algorithm for recombination detection

A novel coronavirus associated with severe acute respiratory syndrome

Comments to the predecessor of human SARS coronavirus in 2003-2004 epidemic

Evolutionary analyses of European H1N2 swine influenza A virus by placing timestamps on the multiple reassortment events

Phylogenetic analysis and sequence comparisons of structural and non-structural SARS coronavirus proteins in Taiwan

Molecular epidemiology of severe acute respiratory syndrome-associated coronavirus infections in Taiwan

Severe acute respiratory syndrome coronaviruslike virus in Chinese horseshoe bats

The molecular population genetics of HIV-1 group O

Tracing the origin and history of the HIV-2 epidemic

Structural analysis of major species barriers between humans and palm civets for severe acute respiratory syndrome coronavirus infections

Structure of SARS coronavirus spike receptor-binding domain complexed with receptor

Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus

Bats are natural reservoirs of SARS-like coronaviruses

Animal origins of the severe acute respiratory syndrome coronavirus: insight from ACE2-S-protein interactions

Receptor and viral determinants of SARS-coronavirus adaptation to human ACE2

Predicting super spreading events during the 2003 severe acute respiratory syndrome epidemics in Hong Kong and Singapore

Laboratory diagnosis of four recent sporadic cases of communityacquired SARS

Phylogenomics and bioinformatics of SARS-CoV

SARS transmission pattern in Singapore reassessed by viral sequence variation analysis

Isolation of avian infectious bronchitis coronavirus from domestic peafowl (Pavo cristatus) and teal (Anas)

Molecular epidemiology of SARS-associated coronavirus

Superspreading and the effect of individual variation on disease emergence

Date of origin of the SARS coronavirus strains

Identification of a novel coronavirus from a beluga whale by using a panviral microarray

Possible role of an animal vector in the SARS outbreak at Amoy Gardens

The origins of new pandemic viruses: the acquisition of new host ranges by canine parvovirus and influenza A viruses

Coronavirus as a possible cause of severe acute respiratory syndrome

How good are deep phylogenetic trees?

Recombination in evolutionary genomics

Mosaic structure of human coronavirus NL63, one thousand years of evolution

Full-length genome sequences of two SARS-like coronaviruses in horseshoe bats and genetic variation analysis

SARS associated coronavirus has a recombinant polymerase and coronaviruses have a history of host-shifting

Characterization of a novel coronavirus associated with severe acute respiratory syndrome

Comparative full-length genome sequence analysis of 14 SARS coronavirus isolates and common mutations associated with putative origins of infection

Severe acute respiratory syndrome coronavirus sequence characteristics and evolutionary rate estimate from maximum likelihood analysis

Evolutionary implications of Avian Infectious Bronchitis Virus (AIBV) analysis

A review of studies on animal reservoirs of the SARS coronavirus

SARS-CoV infection was from at least two origins in the Taiwan area

Unique and conserved features of genome and proteome of SARS-coronavirus, an early split-off from the coronavirus group 2 lineage

Cross-host evolution of severe acute respiratory syndrome coronavirus in palm civet and human

Evidence from the evolutionary analysis of nucleotide sequences for a recombinant history of SARS-CoV

Mosaic evolution of the severe acute respiratory syndrome coronavirus

Genetic diversity and models of viral evolution for the hepatitis C virus

Characterizing 56 complete SARS-CoV S-gene sequences from Hong Kong

Differential stepwise evolution of SARS coronavirus functional proteins in different host species

Prevalence and genetic diversity of coronaviruses in bats from China

SARS: experience at Prince of Wales Hospital, Hong Kong

Coronavirus genomic-sequence variations and the epidemiology of the severe acute respiratory syndrome

Mutational dynamics of the SARS coronavirus in cell culture and human populations isolated in 2003

Evolutionary insights into the ecology of coronaviruses

Complete genomic sequence of human coronavirus OC43: molecular clock analysis suggests a relatively recent zoonotic coronavirus transmission event

SARS-CoV infection in a restaurant from palm civet

Molecular evolution and multilocus sequence typing of 145 strains of SARS-CoV

Update 4: Review of Probable and Laboratory-Confirmed SARS Cases in Southern China

A 193-amino acid fragment of the SARS coronavirus S protein efficiently binds angiotensinconverting enzyme 2

Characterization and complete genome sequence of a novel coronavirus, coronavirus HKU1, from patients with pneumonia

Phylogenetic and recombination analysis of coronavirus HKU1, a novel coronavirus from patients with pneumonia

Comparative analysis of complete genome sequences of three avian coronaviruses reveals a novel group 3c coronavirus

Molecular diversity of coronaviruses in bats

Origin of AIDS: contaminated polio vaccine theory refuted

Civets are equally susceptible to experimental infection by two different severe acute respiratory syndrome coronavirus isolates

Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages

Characterization of severe acute respiratory syndrome coronavirus genomes in Taiwan: molecular epidemiology and genome evolution

Estimation of primate speciation dates using local molecular clocks

Estimated timing of the last common ancestor of the SARS coronavirus

The complete genome sequence of severe acute respiratory syndrome coronavirus strain HKU-39849 (HK-39)

Adaptive evolution of the spike gene of SARS coronavirus: changes in positively selected sites in different epidemic groups

Evaluation of an improved branch-site likelihood method for detecting positive selection at the molecular level

Testing the hypothesis of a recombinant origin of the SARS-associated coronavirus

SARS molecular epidemiology: a Chinese fairy tale of controlling an emerging zoonotic disease in the genomics era

Moderate mutation rate in the SARS coronavirus genome and its implications

Epidemiology and cause of severe acute respiratory syndrome (SARS) in Guangdong, People's Republic of China

Monophyletic relationship between severe acute respiratory syndrome coronavirus and group 2 coronaviruses

This work is partially supported by Research Fund for the Control of Infectious Diseases (reference number 06060672) from Hong Kong SAR government and the Strategic Research Theme of Infection and Immunology, The University of Hong Kong.