Investigation of the Molecular Epidemiology and Evolution of Circulating Severe Acute Respiratory Syndrome Coronavirus 2 in Thailand from 2020 to 2022 via Next-Generation Sequencing.

Coronavirus disease 2019 (COVID-19) is an infectious condition caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), which surfaced in Thailand in early 2020. The current study investigated the SARS-CoV-2 lineages circulating in Thailand and their evolutionary history. Complete genome sequencing of 210 SARS-CoV-2 samples collected from collaborating hospitals and the Institute of Urban Disease Control and Prevention over two years, from December 2020 to July 2022, was performed using next-generation sequencing technology. Multiple lineage introductions were observed before the emergence of the B.1.1.529 omicron variant, including B.1.36.16, B.1.351, B.1.1, B.1.1.7, B.1.524, AY.30, and B.1.617.2. The B.1.1.529 omicron variant was subsequently detected between January 2022 and June 2022. The evolutionary rate for the spike gene of SARS-CoV-2 was estimated to be between 0.87 and 1.71 × 10-3 substitutions per site per year. There was a substantial prevalence of the predominant mutations C25672T (L94F), C25961T (T190I), and G26167T (V259L) in the ORF3a gene during the Thailand outbreaks. Complete genome sequencing can enhance the prediction of future variant changes in viral genomes, which is crucial to ensuring that vaccine strains are protective against worldwide outbreaks.

Molecular characterisation and tracking of severe acute respiratory syndrome coronavirus 2 in Thailand, 2020-2022.

The global COVID-19 pandemic, caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was first detected in China in December 2019. To date, there have been approximately 3.4 million reported cases of COVID-19 and over 24,000 deaths in Thailand. In this study, we investigated the molecular characteristics and evolution of SARS-CoV-2 in Thailand from 2020 to 2022. Two hundred sixty-eight SARS-CoV-2 isolates, collected mostly in Bangkok from COVID-19 patients, were characterised by partial genome sequencing. Moreover, the viruses in 5,627 positive SARS-CoV-2 samples were identified as viral variants - B.1.1.7 (Alpha), B.1.617.2 (Delta), B.1.1.529 (Omicron/BA.1), or B.1.1.529 (Omicron/BA.2) - by multiplex real-time reverse transcription polymerase chain reaction (RT-PCR) assays. The results revealed that B.1.36.16 caused the predominant outbreak in the second wave (December 2020-January 2021), B.1.1.7 (Alpha) in the third wave (April-June 2021), B.1.617.2 (Delta) in the fourth wave (July-December 2021), and B.1.1.529 (Omicron) in the fifth wave (January-March 2022). The evolutionary rate of the viral genome was 2.60 × 10-3 (95% highest posterior density [HPD], 1.72 × 10-3 to 3.62 × 10-3) nucleotide substitutions per site per year. Continued molecular surveillance of SARS-CoV-2 is crucial for monitoring emerging variants with the potential to cause new COVID-19 outbreaks.

Evolutionary and Genetic Recombination Analyses of Coxsackievirus A6 Variants Associated with Hand, Foot, and Mouth Disease Outbreaks in Thailand between 2019 and 2022.

Coxsackievirus (CV)-A6 infections cause hand, foot, and mouth disease (HFMD) in children and adults. Despite the serious public health threat presented by CV-A6 infections, our understanding of the mechanisms by which new CV-A6 strains emerge remains limited. This study investigated the molecular epidemiological trends, evolutionary dynamics, and recombination characteristics of CV-A6-associated HFMD in Thailand between 2019 and 2022. In the HFMD patient samples collected during the 4-year study period, we identified enterovirus (EV) RNA in 368 samples (48.7%), of which CV-A6 (23.7%) was the predominant genotype, followed by CV-A4 (6%), EV-A71 (3.7%), and CV-A16 (3.4%). According to the partial viral protein (VP) 1 sequences, all these CV-A6 strains belonged to the D3 clade. Based on the viral-RNA-dependent RNA polymerase (RdRp) gene, four recombinant forms (RFs), RF-A (147, 84.5%), RF-N (11, 6.3%), RF-H (1, 0.6%), and newly RF-Y (15, 8.6%), were identified throughout the study period. Results from the similarity plot and bootscan analyses revealed that the 3D polymerase (3Dpol) region of the D3/RF-Y subclade consists of sequences highly similar to CV-A10. We envisage that the epidemiological and evolutionarily insights presented in this manuscript will contribute to the development of vaccines to prevent the spread of CV-A6 infection.