Abstract:Background: The coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome-related coronavirus-2 (SARS-CoV-2) is pandemic. However, the origins and global transmission pattern of SARS-CoV-2 remain largely unknown. We aimed to characterize the origination and transmission of SARS-CoV-2 based on evolutionary dynamics. Methods: Using the full-length sequences of SARS-CoV-2 with intact geographic, demographic, and temporal information worldwide from the GISAID database during 26 December 2019 and 30 November 2020, we constructed the transmission tree to depict the evolutionary process by the R package "outbreaker". The affinity of the mutated receptor-binding region of the spike protein to angiotensin-converting enzyme 2 (ACE2) was predicted using mCSM-PPI2 software. Viral infectivity and antigenicity were tested in ACE2-transfected HEK293T cells by pseudovirus transfection and neutralizing antibody test. Results: From 26 December 2019 to 8 March 2020, early stage of the COVID-19 pandemic, SARS-CoV-2 strains identified worldwide were mainly composed of three clusters: the Europe-based cluster including two USA-based sub-clusters; the Asia-based cluster including isolates in China, Japan, the USA, Singapore, Australia, Malaysia, and Italy; and the USA-based cluster. The SARS-CoV-2 strains identified in the USA formed four independent clades while those identified in China formed one clade. After 8 March 2020, the clusters of SARS-CoV-2 strains tended to be independent and became "pure" in each of the major countries. Twenty-two of 60 mutations in the receptor-binding domain of the spike protein were predicted to increase the binding affinity of SARS-CoV-2 to ACE2. Of all predicted mutants, the number of E484K was the largest one with 86 585 sequences, followed by S477N with 55 442 sequences worldwide. In more than ten countries, the frequencies of the isolates with E484K and S477N increased significantly. V367F and N354D mutations increased the infectivity of SARS-CoV-2 pseudoviruses (P < 0.001). SARS-CoV-2 with V367F was more sensitive to the S1-targeting neutralizing antibody than the wild-type counterpart (P < 0.001). Conclusions: SARS-CoV-2 strains might have originated in several countries simultaneously under certain evolutionary pressure. Travel restrictions might cause location-specific SARS-CoV-2 clustering. The SARS-CoV-2 evolution appears to facilitate its transmission via altering the affinity to ACE2 or immune evasion.

Temporal epistasis inference from more than 3,500,000 SARS-CoV-2 Genomic Sequences

Global analysis of more than 50,000 SARS-Cov-2 genomes reveals epistasis between 8 viral genes

Real-time identification of epistatic interactions in SARS-CoV-2 from large genome collections

Genetic drift of human coronavirus OC43 spike gene during adaptive evolution

DCA for genome-wide epistasis analysis: the statistical genetics perspective

Efficient epistasis inference via higher-order covariance matrix factorization

Molecular Evolution and Multilocus Sequence Typing of 145 Strains of Sars-Cov

Genomic epidemiology reveals early transmission of SARS-CoV-2 and mutational dynamics in Nanning, China

Community structure and temporal dynamics of SARS-CoV-2 epistatic network allow for early detection of emerging variants with altered phenotypes

Genome-wide association study between SARS-CoV-2 single nucleotide polymorphisms and virus copies during infections

Epi-Clock: A sensitive platform to help understand pathogenic disease outbreaks and facilitate the response to future outbreaks of concern

The effects of amino acid substitution of spike protein and genomic recombination on the evolution of SARS-CoV-2

Evaluating the effect of SARS-CoV-2 spike mutations with a linear doubly robust learner

Characterization of SARS-CoV-2 worldwide transmission based on evolutionary dynamics and specific viral mutations in the spike protein

Using big sequencing data to identify chronic SARS-Coronavirus-2 infections

Phylogenetic signatures reveal multilevel selection and fitness costs in SARS-CoV-2

Topological data analysis identifies emerging adaptive mutations in SARS-CoV-2

Prediction of Prospective Mutational Landscape of SARS-CoV-2 Spike ssRNA and Evolutionary Basis of Its Host Interaction

Intra-host variation and evolutionary dynamics of SARS-CoV-2 populations in COVID-19 patients

Two fitness inference schemes compared using allele frequencies from 1,068,391 sequences sampled in the UK during the COVID-19 pandemic

Structural and molecular basis of the epistasis effect in enhanced affinity between SARS-CoV-2 KP.3 and ACE2