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...

Full description

Bibliographic Details
Main Authors: Mingjian Zhu, Qianli Zeng, Bryanna Infinita Laviashna Saputro, Sien Ping Chew, Ian Chew, Holie Frendy, Joanna Weihui Tan, Lanjuan Li
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