Abstract:Background: The ongoing coronavirus disease 2019 (COVID-19) pandemic has posed an unprecedented challenge to public health in Southeast Asia, a tropical region with limited resources. This study aimed to investigate the evolutionary dynamics and spatiotemporal patterns of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the region. Materials and Methods: A total of 1491 complete SARS-CoV-2 genome sequences from 10 Southeast Asian countries were downloaded from the Global Initiative on Sharing Avian Influenza Data (GISAID) database on November 17, 2020. The evolutionary relationships were assessed using maximum likelihood (ML) and time-scaled Bayesian phylogenetic analyses, and the phylogenetic clustering was tested using principal component analysis (PCA). The spatial patterns of SARS-CoV-2 spread within Southeast Asia were inferred using the Bayesian stochastic search variable selection (BSSVS) model. The effective population size (Ne) trajectory was inferred using the Bayesian Skygrid model. Results: Four major clades (including one potentially endemic) were identified based on the maximum clade credibility (MCC) tree. Similar clustering was yielded by PCA; the first three PCs explained 46.9% of the total genomic variations among the samples. The time to the most recent common ancestor (tMRCA) and the evolutionary rate of SARS-CoV-2 circulating in Southeast Asia were estimated to be November 28, 2019 (September 7, 2019 to January 4, 2020) and 1.446 x 10(-3) (1.292 x 10(-3) to 1.613 x 10(-3)) substitutions per site per year, respectively. Singapore and Thailand were the two most probable root positions, with posterior probabilities of 0.549 and 0.413, respectively. There were high-support transmission links (Bayes factors exceeding 1,000) in Singapore, Malaysia, and Indonesia; Malaysia involved the highest number (7) of inferred transmission links within the region. A twice-accelerated viral population expansion, followed by a temporary setback, was inferred during the early stages of the pandemic in Southeast Asia. Conclusions: With available genomic data, we illustrate the phylogeography and phylodynamics of SARS-CoV-2 circulating in Southeast Asia. Continuous genomic surveillance and enhanced strategic collaboration should be listed as priorities to curb the pandemic, especially for regional communities dominated by developing countries.

Comprehensive detection and dissection of interlineage recombination events in the SARS-CoV-2 pandemic

VirusRecom: an information-theory-based method for recombination detection of viral lineages and its application on SARS-CoV-2

Unraveling the genetic evolution of SARS-CoV-2 Recombinants using mutational dynamics across the different lineages

Systematic detection of co-infection and intra-host recombination in more than 2 million global SARS-CoV-2 samples

Data-driven recombination detection in viral genomes

The nuts and bolts of recombination in the generation of SARS-CoV-2 variants; from XA to XBB

Tracking SARS-CoV-2 Omicron diverse spike gene mutations identifies multiple inter-variant recombination events

New framework for recombination and adaptive evolution analysis with application to the novel coronavirus SARS-CoV-2

Genomic Feature Analysis of Betacoronavirus Provides Insights Into SARS and COVID-19 Pandemics

Identification of Various Recombinants in a Patient Coinfected With the Different SARS‐CoV‐2 Variants

Recombinant SARS-CoV-2 Delta/Omicron BA.5 emerging in an immunocompromised long-term infected COVID-19 patient

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

Recombination-aware phylogenetic analysis sheds light on the evolutionary origin of SARS-CoV-2

Evolutionary analysis and lineage designation of SARS-CoV-2 genomes

On the Origin and Continuing Evolution of SARS-CoV-2

Analysis of Recombinant Characteristics Based on 949 PRRSV-2 Genomic Sequences Obtained from 1991 to 2021 Shows That Viral Multiplication Ability Contributes to Dominant Recombination

Molecular Phylogenesis and Spatiotemporal Spread of SARS-CoV-2 in Southeast Asia

The dominant strain of SARS-CoV-2 is a mosaicism

Classification of SARS-CoV-2 sequences as recombinants via a pre-trained CNN and identification of a mathematical signature relative to recombinant feature at Spike, via interpretability

Mining genomic repositories to further our knowledge of the extent of SARS-CoV-2 co-infections