Bayesian phylodynamic inferences on the temporal evolution and global transmission of SARS-CoV-2
Abstract Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread widely from China to the world. Although the viral genome has been well characterized, the evolutionary origin and global transmission dynamics of SARS-CoV-2 remain poorly investigated. To address this, we retrieved 313 SARS-CoV-2 genomes from the GISAID database (https://www.gisaid.org), from which 99 genomes generated from original clinical specimens with exact collection dates from 16 countries were selected and enrolled for Bayesian phylodynamic analysis. Here we show that the time to the Most Recent Common Ancestor (tMRCA) of SARS-CoV-2 is Dec 11, 2019 (95%HPD, Nov 21 - Dec 24). Two clades of global circulating strains of SARS-CoV-2 were suggested by Bayesian Maximum Clade Credibility (MCC) tree. The USA circulating strains of SARS-CoV-2 seemed to be from both of the two clades, the UK and Australia circulating strains were from Clade 1, the circulating strains in Singapore, Japan, Germany, France, and Italy were from Clade 2. Although we have not found any obvious bottle-neck-effect from the Bayesian Skyline Plot of the viral population dynamics reconstruction, a sharp reduction of the lower 95% HPD of the relative genetic diversity was observed from Feb 5, 2020, suggesting a possible initiation of a bottle-neck-effect. Thirteen (6 synonymous and 7 non-synonymous) mutations in the viral genome were observed, including two clade-specific mutations (C8782T and T1844C in Clade 1 rather than Clade 2) and eleven sub-clade specific mutations. All of the observed mutations occurred in the USA circulating strains, except one mutation T18488C only occurred in the UK circulating strains. A non-synonymous mutation in the 3’-UTR was also observed, suggesting an altered RNA replication capacity of SARS-CoV-2. We thus came to the conclusion that continuous evolution occurred in almost all regions of the SARS-CoV-2 genome and potentially in a country-specific manner. Further efforts on monitoring the genomic mutations of SARS-CoV-2 from different countries are recommended.