Research article

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

Vinsensius B Vega¹ , Yijun Ruan¹ , Jianjun Liu¹ , Wah Heng Lee¹ , Chia Lin Wei¹ , Su Yun Se-Thoe² , Kin Fai Tang³ , Tao Zhang¹ , Prasanna R Kolatkar¹ , Eng Eong Ooi³ , Ai Ee Ling² , Lawrence W Stanton¹ , Philip M Long⁴ and Edison T Liu¹

¹Genome Institute of Singapore, 60 Biopolis Street, Singapore 138672

²Virology Section, Department of Pathology, Singapore General Hospital, Singapore

³Environmental Health Institute, 41 Science Park Road, Singapore Science Park II, Singapore 117610

⁴Center for Computational Learning Systems, Columbia University, New York, NY 10027 USA

author email corresponding author email

BMC Infectious Diseases 2004, 4:32doi:10.1186/1471-2334-4-32

Published:	6 September 2004

Abstract

Background

The SARS coronavirus is the etiologic agent for the epidemic of the Severe Acute Respiratory Syndrome. The recent emergence of this new pathogen, the careful tracing of its transmission patterns, and the ability to propagate in culture allows the exploration of the mutational dynamics of the SARS-CoV in human populations.

Methods

We sequenced complete SARS-CoV genomes taken from primary human tissues (SIN3408, SIN3725V, SIN3765V), cultured isolates (SIN848, SIN846, SIN842, SIN845, SIN847, SIN849, SIN850, SIN852, SIN3408L), and five consecutive Vero cell passages (SIN2774_P1, SIN2774_P2, SIN2774_P3, SIN2774_P4, SIN2774_P5) arising from SIN2774 isolate. These represented individual patient samples, serial in vitro passages in cell culture, and paired human and cell culture isolates. Employing a refined mutation filtering scheme and constant mutation rate model, the mutation rates were estimated and the possible date of emergence was calculated. Phylogenetic analysis was used to uncover molecular relationships between the isolates.

Results

Close examination of whole genome sequence of 54 SARS-CoV isolates identified before 14^th October 2003, including 22 from patients in Singapore, revealed the mutations engendered during human-to-Vero and Vero-to-human transmission as well as in multiple Vero cell passages in order to refine our analysis of human-to-human transmission. Though co-infection by different quasipecies in individual tissue samples is observed, the in vitro mutation rate of the SARS-CoV in Vero cell passage is negligible. The in vivo mutation rate, however, is consistent with estimates of other RNA viruses at approximately 5.7 × 10^-6 nucleotide substitutions per site per day (0.17 mutations per genome per day), or two mutations per human passage (adjusted R-square = 0.4014). Using the immediate Hotel M contact isolates as roots, we observed that the SARS epidemic has generated four major genetic groups that are geographically associated: two Singapore isolates, one Taiwan isolate, and one North China isolate which appears most closely related to the putative SARS-CoV isolated from a palm civet. Non-synonymous mutations are centered in non-essential ORFs especially in structural and antigenic genes such as the S and M proteins, but these mutations did not distinguish the geographical groupings. However, no non-synonymous mutations were found in the 3CLpro and the polymerase genes.

Conclusions

Our results show that the SARS-CoV is well adapted to growth in culture and did not appear to undergo specific selection in human populations. We further assessed that the putative origin of the SARS epidemic was in late October 2002 which is consistent with a recent estimate using cases from China. The greater sequence divergence in the structural and antigenic proteins and consistent deletions in the 3' – most portion of the viral genome suggest that certain selection pressures are interacting with the functional nature of these validated and putative ORFs.