Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data
Abstract Background As a new epi-center of COVID-19 in Asia and a densely populated developing country, Indonesia is facing unprecedented challenges in public health. SARS-CoV-2 lineage B.1.466.2 was reported to be an indigenous dominant strain in Indonesia (once second only to the Delta variant). H...
Main Authors: | , , , , , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
BMC
2022-06-01
|
Series: | Virology Journal |
Subjects: | |
Online Access: | https://doi.org/10.1186/s12985-022-01830-1 |
_version_ | 1811344954745683968 |
---|---|
author | Mingjian Zhu Qianli Zeng Bryanna Infinita Laviashna Saputro Sien Ping Chew Ian Chew Holie Frendy Joanna Weihui Tan Lanjuan Li |
author_facet | Mingjian Zhu Qianli Zeng Bryanna Infinita Laviashna Saputro Sien Ping Chew Ian Chew Holie Frendy Joanna Weihui Tan Lanjuan Li |
author_sort | Mingjian Zhu |
collection | DOAJ |
description | Abstract Background As a new epi-center of COVID-19 in Asia and a densely populated developing country, Indonesia is facing unprecedented challenges in public health. SARS-CoV-2 lineage B.1.466.2 was reported to be an indigenous dominant strain in Indonesia (once second only to the Delta variant). However, it remains unclear how this variant evolved and spread within such an archipelagic nation. Methods For statistical description, the spatiotemporal distributions of the B.1.466.2 variant were plotted using the publicly accessible metadata in GISAID. A total of 1302 complete genome sequences of Indonesian B.1.466.2 strains with high coverage were downloaded from the GISAID’s EpiCoV database on 28 August 2021. To determine the molecular evolutionary characteristics, we performed a time-scaled phylogenetic analysis using the maximum likelihood algorithm and called the single nucleotide variants taking the Wuhan-Hu-1 sequence as reference. To investigate the spatiotemporal transmission patterns, we estimated two dynamic parameters (effective population size and effective reproduction number) and reconstructed the phylogeography among different islands. Results As of the end of August 2021, nearly 85% of the global SARS-CoV-2 lineage B.1.466.2 sequences (including the first one) were obtained from Indonesia. This variant was estimated to account for over 50% of Indonesia’s daily infections during the period of March–May 2021. The time-scaled phylogeny suggested that SARS-CoV-2 lineage B.1.466.2 circulating in Indonesia might have originated from Java Island in mid-June 2020 and had evolved into two disproportional and distinct sub-lineages. High-frequency non-synonymous mutations were mostly found in the spike and NSP3; the S-D614G/N439K/P681R co-mutations were identified in its larger sub-lineage. The demographic history was inferred to have experienced four phases, with an exponential growth from October 2020 to February 2021. The effective reproduction number was estimated to have reached its peak (11.18) in late December 2020 and dropped to be less than one after early May 2021. The relevant phylogeography showed that Java and Sumatra might successively act as epi-centers and form a stable transmission loop. Additionally, several long-distance transmission links across seas were revealed. Conclusions SARS-CoV-2 variants circulating in the tropical archipelago may follow unique patterns of evolution and transmission. Continuous, extensive and targeted genomic surveillance is essential. |
first_indexed | 2024-04-13T19:55:46Z |
format | Article |
id | doaj.art-ada4918767ec4babb5503e12174fac54 |
institution | Directory Open Access Journal |
issn | 1743-422X |
language | English |
last_indexed | 2024-04-13T19:55:46Z |
publishDate | 2022-06-01 |
publisher | BMC |
record_format | Article |
series | Virology Journal |
spelling | doaj.art-ada4918767ec4babb5503e12174fac542022-12-22T02:32:21ZengBMCVirology Journal1743-422X2022-06-0119111310.1186/s12985-022-01830-1Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance dataMingjian Zhu0Qianli Zeng1Bryanna Infinita Laviashna Saputro2Sien Ping Chew3Ian Chew4Holie Frendy5Joanna Weihui Tan6Lanjuan Li7State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, National Clinical Research Center for Infectious Diseases, National Medical Center for Infectious Diseases, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Zhejiang University School of MedicineShanghai Institute of Biological ProductsFaculty of Medicine, University of IndonesiaShanghai Jiao Tong University School of MedicineZhejiang University School of MedicineFaculty of Medicine and Health Sciences, Krida Wacana Christian UniversityFaculty of Arts and Social Sciences, National University of SingaporeState Key Laboratory for Diagnosis and Treatment of Infectious Diseases, National Clinical Research Center for Infectious Diseases, National Medical Center for Infectious Diseases, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Zhejiang University School of MedicineAbstract Background As a new epi-center of COVID-19 in Asia and a densely populated developing country, Indonesia is facing unprecedented challenges in public health. SARS-CoV-2 lineage B.1.466.2 was reported to be an indigenous dominant strain in Indonesia (once second only to the Delta variant). However, it remains unclear how this variant evolved and spread within such an archipelagic nation. Methods For statistical description, the spatiotemporal distributions of the B.1.466.2 variant were plotted using the publicly accessible metadata in GISAID. A total of 1302 complete genome sequences of Indonesian B.1.466.2 strains with high coverage were downloaded from the GISAID’s EpiCoV database on 28 August 2021. To determine the molecular evolutionary characteristics, we performed a time-scaled phylogenetic analysis using the maximum likelihood algorithm and called the single nucleotide variants taking the Wuhan-Hu-1 sequence as reference. To investigate the spatiotemporal transmission patterns, we estimated two dynamic parameters (effective population size and effective reproduction number) and reconstructed the phylogeography among different islands. Results As of the end of August 2021, nearly 85% of the global SARS-CoV-2 lineage B.1.466.2 sequences (including the first one) were obtained from Indonesia. This variant was estimated to account for over 50% of Indonesia’s daily infections during the period of March–May 2021. The time-scaled phylogeny suggested that SARS-CoV-2 lineage B.1.466.2 circulating in Indonesia might have originated from Java Island in mid-June 2020 and had evolved into two disproportional and distinct sub-lineages. High-frequency non-synonymous mutations were mostly found in the spike and NSP3; the S-D614G/N439K/P681R co-mutations were identified in its larger sub-lineage. The demographic history was inferred to have experienced four phases, with an exponential growth from October 2020 to February 2021. The effective reproduction number was estimated to have reached its peak (11.18) in late December 2020 and dropped to be less than one after early May 2021. The relevant phylogeography showed that Java and Sumatra might successively act as epi-centers and form a stable transmission loop. Additionally, several long-distance transmission links across seas were revealed. Conclusions SARS-CoV-2 variants circulating in the tropical archipelago may follow unique patterns of evolution and transmission. Continuous, extensive and targeted genomic surveillance is essential.https://doi.org/10.1186/s12985-022-01830-1Severe acute respiratory syndrome coronavirus 2B.1.466.2 variantGenomic surveillancePrevalencePhylogenesisSingle nucleotide variant |
spellingShingle | Mingjian Zhu Qianli Zeng Bryanna Infinita Laviashna Saputro Sien Ping Chew Ian Chew Holie Frendy Joanna Weihui Tan Lanjuan Li Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data Virology Journal Severe acute respiratory syndrome coronavirus 2 B.1.466.2 variant Genomic surveillance Prevalence Phylogenesis Single nucleotide variant |
title | Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data |
title_full | Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data |
title_fullStr | Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data |
title_full_unstemmed | Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data |
title_short | Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data |
title_sort | tracking the molecular evolution and transmission patterns of sars cov 2 lineage b 1 466 2 in indonesia based on genomic surveillance data |
topic | Severe acute respiratory syndrome coronavirus 2 B.1.466.2 variant Genomic surveillance Prevalence Phylogenesis Single nucleotide variant |
url | https://doi.org/10.1186/s12985-022-01830-1 |
work_keys_str_mv | AT mingjianzhu trackingthemolecularevolutionandtransmissionpatternsofsarscov2lineageb14662inindonesiabasedongenomicsurveillancedata AT qianlizeng trackingthemolecularevolutionandtransmissionpatternsofsarscov2lineageb14662inindonesiabasedongenomicsurveillancedata AT bryannainfinitalaviashnasaputro trackingthemolecularevolutionandtransmissionpatternsofsarscov2lineageb14662inindonesiabasedongenomicsurveillancedata AT sienpingchew trackingthemolecularevolutionandtransmissionpatternsofsarscov2lineageb14662inindonesiabasedongenomicsurveillancedata AT ianchew trackingthemolecularevolutionandtransmissionpatternsofsarscov2lineageb14662inindonesiabasedongenomicsurveillancedata AT holiefrendy trackingthemolecularevolutionandtransmissionpatternsofsarscov2lineageb14662inindonesiabasedongenomicsurveillancedata AT joannaweihuitan trackingthemolecularevolutionandtransmissionpatternsofsarscov2lineageb14662inindonesiabasedongenomicsurveillancedata AT lanjuanli trackingthemolecularevolutionandtransmissionpatternsofsarscov2lineageb14662inindonesiabasedongenomicsurveillancedata |