1
|
Jimenez-Vasquez V, Vargas-Herrera N, Bárcena-Flores L, Hurtado V, Padilla-Rojas C, Araujo-Castillo RV. Dispersion of SARS-CoV-2 lineage BA.5.1.25 and its descendants in Peru during two COVID-19 waves in 2022. Genomics Inform 2024; 22:5. [PMID: 38907313 PMCID: PMC11184951 DOI: 10.1186/s44342-024-00006-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Accepted: 03/04/2024] [Indexed: 06/23/2024] Open
Abstract
During the third year of the pandemic in Peru, the persistent transmission of SARS-CoV-2 led to the appearance of more transmissible and immune-evasive Omicron sublineages; in that context, the National Genomic Surveillance of SARS-CoV-2 performed by the Peruvian National Institute of Health detected spike mutations in the circulating Omicron BA.5.1.25 sublineage which was later designated as DJ.1 and increased during the fourth COVID-19 wave, this eventually branched into new sublineages. The introduction, emergence, and timing of the most recent common ancestor (tMRCA) of BA.5.1.25 and its descendants (DJ.1, DJ.1.1, DJ.1.2, and DJ.1.3) were investigated in this paper as well as the time lags between their emergence and identification by the Peruvian National Institute of Health. Our findings show that ongoing genomic surveillance of SARS-CoV-2 is critical for understanding its phylogenetic evolution and the emergence of novel variations.
Collapse
Affiliation(s)
- Victor Jimenez-Vasquez
- Centro Nacional de Salud Pública, Instituto Nacional de Salud, Capac Yupanqui 1400-Jesus Maria, Lima, Peru
| | - Natalia Vargas-Herrera
- Centro Nacional de Salud Pública, Instituto Nacional de Salud, Capac Yupanqui 1400-Jesus Maria, Lima, Peru.
| | - Luis Bárcena-Flores
- Centro Nacional de Salud Pública, Instituto Nacional de Salud, Capac Yupanqui 1400-Jesus Maria, Lima, Peru
| | - Verónica Hurtado
- Centro Nacional de Salud Pública, Instituto Nacional de Salud, Capac Yupanqui 1400-Jesus Maria, Lima, Peru
| | - Carlos Padilla-Rojas
- Centro Nacional de Salud Pública, Instituto Nacional de Salud, Capac Yupanqui 1400-Jesus Maria, Lima, Peru
| | - Roger V Araujo-Castillo
- Centro Nacional de Salud Pública, Instituto Nacional de Salud, Capac Yupanqui 1400-Jesus Maria, Lima, Peru
| |
Collapse
|
2
|
Bukhari AR, Ashraf J, Kanji A, Rahman YA, Trovão NS, Thielen PM, Yameen M, Kanwar S, Khan W, Kabir F, Nisar MI, Merritt B, Hasan R, Spiro D, Rasmussen Z, Aamir UB, Hasan Z. Sequential viral introductions and spread of BA.1 across Pakistan provinces during the Omicron wave. BMC Genomics 2023; 24:432. [PMID: 37532989 PMCID: PMC10399012 DOI: 10.1186/s12864-023-09539-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2023] [Accepted: 07/27/2023] [Indexed: 08/04/2023] Open
Abstract
BACKGROUND COVID-19 waves caused by specific SARS-CoV-2 variants have occurred globally at different times. We focused on Omicron variants to understand the genomic diversity and phylogenetic relatedness of SARS-CoV-2 strains in various regions of Pakistan. METHODS We studied 276,525 COVID-19 cases and 1,031 genomes sequenced from December 2021 to August 2022. Sequences were analyzed and visualized using phylogenetic trees. RESULTS The highest case numbers and deaths were recorded in Sindh and Punjab, the most populous provinces in Pakistan. Omicron variants comprised 93% of all genomes, with BA.2 (32.6%) and BA.5 (38.4%) predominating. The first Omicron wave was associated with the sequential identification of BA.1 in Sindh, then Islamabad Capital Territory, Punjab, Khyber Pakhtunkhwa (KP), Azad Jammu Kashmir (AJK), Gilgit-Baltistan (GB) and Balochistan. Phylogenetic analysis revealed Sindh to be the source of BA.1 and BA.2 introductions into Punjab and Balochistan during early 2022. BA.4 was first introduced in AJK and BA.5 in Punjab. Most recent common ancestor (MRCA) analysis revealed relatedness between the earliest BA.1 genome from Sindh with Balochistan, AJK, Punjab and ICT, and that of first BA.1 from Punjab with strains from KPK and GB. CONCLUSIONS Phylogenetic analysis provides insights into the introduction and transmission dynamics of the Omicron variant in Pakistan, identifying Sindh as a hotspot for viral dissemination. Such data linked with public health efforts can help limit surges of new infections.
Collapse
Affiliation(s)
- Ali Raza Bukhari
- Department of Pathology and Laboratory Medicine, Aga Khan University, Karachi, 74800, Pakistan
| | - Javaria Ashraf
- Department of Pathology and Laboratory Medicine, Aga Khan University, Karachi, 74800, Pakistan
| | - Akbar Kanji
- Department of Pathology and Laboratory Medicine, Aga Khan University, Karachi, 74800, Pakistan
| | - Yusra Abdul Rahman
- Department of Pathology and Laboratory Medicine, Aga Khan University, Karachi, 74800, Pakistan
| | - Nídia S Trovão
- Fogarty International Center, U.S. National Institutes of Health, 16 Center Drive, Bethesda, MD, 20892, USA
| | - Peter M Thielen
- Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, MD, 20723, USA
| | - Maliha Yameen
- Department of Pathology and Laboratory Medicine, Aga Khan University, Karachi, 74800, Pakistan
| | - Samiah Kanwar
- Department of Pediatrics and Child Health, Aga Khan University, Karachi, 74800, Pakistan
| | - Waqasuddin Khan
- Department of Pediatrics and Child Health, Aga Khan University, Karachi, 74800, Pakistan
| | - Furqan Kabir
- Department of Pediatrics and Child Health, Aga Khan University, Karachi, 74800, Pakistan
| | - Muhammad Imran Nisar
- Department of Pediatrics and Child Health, Aga Khan University, Karachi, 74800, Pakistan
| | - Brian Merritt
- Johns Hopkins University Applied Physics Laboratory, 11100 Johns Hopkins Road, Laurel, MD, 20723, USA
| | - Rumina Hasan
- Department of Pathology and Laboratory Medicine, Aga Khan University, Karachi, 74800, Pakistan
| | - David Spiro
- Fogarty International Center, U.S. National Institutes of Health, 16 Center Drive, Bethesda, MD, 20892, USA
| | - Zeba Rasmussen
- Fogarty International Center, U.S. National Institutes of Health, 16 Center Drive, Bethesda, MD, 20892, USA
| | - Uzma Bashir Aamir
- World Health Organization Country Office, Park Road, Chak Shahzad, Islamabad, Pakistan
| | - Zahra Hasan
- Department of Pathology and Laboratory Medicine, Aga Khan University, Karachi, 74800, Pakistan.
| |
Collapse
|
3
|
Zhu M, Zeng Q, Saputro BIL, Chew SP, Chew I, Frendy H, Tan JW, Li L. Tracking the molecular evolution and transmission patterns of SARS-CoV-2 lineage B.1.466.2 in Indonesia based on genomic surveillance data. Virol J 2022; 19:103. [PMID: 35710544 PMCID: PMC9202327 DOI: 10.1186/s12985-022-01830-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2022] [Accepted: 06/02/2022] [Indexed: 12/22/2022] Open
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. Supplementary Information The online version contains supplementary material available at 10.1186/s12985-022-01830-1.
Collapse
Affiliation(s)
- Mingjian Zhu
- State 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 Medicine, Hangzhou, China
| | - Qianli Zeng
- Shanghai Institute of Biological Products, Shanghai, China
| | | | - Sien Ping Chew
- Shanghai Jiao Tong University School of Medicine, Shanghai, China
| | - Ian Chew
- Zhejiang University School of Medicine, Hangzhou, China
| | - Holie Frendy
- Faculty of Medicine and Health Sciences, Krida Wacana Christian University, Jakarta, Indonesia
| | - Joanna Weihui Tan
- Faculty of Arts and Social Sciences, National University of Singapore, Singapore, Singapore
| | - Lanjuan Li
- State 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 Medicine, Hangzhou, China.
| |
Collapse
|
4
|
Valieris R, Drummond RD, Defelicibus A, Dias-Neto E, Rosales RA, Tojal da Silva I. A mixture model for determining SARS-Cov-2 variant composition in pooled samples. Bioinformatics 2022; 38:1809-1815. [PMID: 35104309 DOI: 10.1093/bioinformatics/btac047] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Revised: 12/14/2021] [Accepted: 01/26/2022] [Indexed: 02/03/2023] Open
Abstract
MOTIVATION Despite of the fast development of highly effective vaccines to control the current COVID-19 pandemics, the unequal distribution and availability of these vaccines worldwide and the number of people infected in the world lead to the continuous emergence of Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) variants of concern. Therefore, it is likely that real-time genomic surveillance will be continuously needed as an unceasing monitoring tool, necessary to follow the spread of the disease and the evolution of the virus. In this context, new genomic variants of SARS-CoV-2, including variants refractory to current vaccines, makes genomic surveillance programs tools of utmost importance. Nevertheless, the lack of appropriate analytical tools to quickly and effectively access the viral composition in meta-transcriptomic sequencing data, including environmental surveillance, represent possible challenges that may impact the fast adoption of this approach to mitigate the spread and transmission of viruses. RESULTS We propose a statistical model for the estimation of the relative frequencies of SARS-CoV-2 variants in pooled samples. This model is built by considering a previously defined selection of genomic polymorphisms that characterize SARS-CoV-2 variants. The methods described here support both raw sequencing reads for polymorphisms-based markers calling and predefined markers in the variant call format. Results obtained using simulated data show that our method is quite effective in recovering the correct variant proportions. Further, results obtained by considering longitudinal data from wastewater samples of two locations in Switzerland agree well with those describing the epidemiological evolution of COVID-19 variants in clinical samples of these locations. Our results show that the described method can be a valuable tool for tracking the proportions of SARS-CoV-2 variants in complex mixtures such as waste water and environmental samples. AVAILABILITY AND IMPLEMENTATION http://github.com/rvalieris/LCS. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Renan Valieris
- Laboratory of Computational Biology and Bioinformatics, CIPE/A.C. Camargo Cancer Center, São Paulo 01508-010, Brazil
| | - Rodrigo D Drummond
- Laboratory of Computational Biology and Bioinformatics, CIPE/A.C. Camargo Cancer Center, São Paulo 01508-010, Brazil
| | - Alexandre Defelicibus
- Laboratory of Computational Biology and Bioinformatics, CIPE/A.C. Camargo Cancer Center, São Paulo 01508-010, Brazil
| | - Emmanuel Dias-Neto
- Laboratory of Medical Genomics, CIPE/A.C. Camargo Cancer Center, São Paulo 01508-010, Brazil
| | - Rafael A Rosales
- Departamento de Computação e Matemática, Universidade de São Paulo, Ribeirão Preto, São Paulo 14040-901, Brazil
| | - Israel Tojal da Silva
- Laboratory of Computational Biology and Bioinformatics, CIPE/A.C. Camargo Cancer Center, São Paulo 01508-010, Brazil
| |
Collapse
|
5
|
Limaye S, Kasibhatla SM, Ramtirthkar M, Kinikar M, Kale MM, Kulkarni-Kale U. Circulation and Evolution of SARS-CoV-2 in India: Let the Data Speak. Viruses 2021; 13:2238. [PMID: 34835044 PMCID: PMC8619538 DOI: 10.3390/v13112238] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2021] [Revised: 10/28/2021] [Accepted: 10/28/2021] [Indexed: 01/08/2023] Open
Abstract
The COVID-19 pandemic is a global challenge that impacted 200+ countries. India ranks in the second and third positions in terms of number of reported cases and deaths. Being a populous country with densely packed cities, SARS-CoV-2 spread exponentially. India sequenced ≈0.14% isolates from confirmed cases for pandemic surveillance and contributed ≈1.58% of complete genomes sequenced globally. This study was designed to map the circulating lineage diversity and to understand the evolution of SARS-CoV-2 in India using comparative genomics and population genetics approaches. Despite varied sequencing coverage across Indian States and Union Territories, isolates belonging to variants of concern (VoC) and variants of interest (VoI) circulated, persisted, and diversified during the first seventeen months of the pandemic. Delta and Kappa lineages emerged in India and spread globally. The phylogenetic tree shows lineage-wise monophyletic clusters of VoCs/VoIs and diversified tree topologies for non-VoC/VoI lineages designated as 'Others' in this study. Evolutionary dynamics analyses substantiate a lack of spatio-temporal clustering, which is indicative of multiple global and local introductions. Sites under positive selection and significant variations in spike protein corroborate with the constellation of mutations to be monitored for VoC/VoI as well as substitutions that are characteristic of functions with implications in virus-host interactions, differential glycosylation, immune evasion, and escape from neutralization.
Collapse
Affiliation(s)
- Sanket Limaye
- Bioinformatics Centre, Savitribai Phule Pune University (Formerly University of Pune), Pune 411007, India; (S.L.); (S.M.K.); (M.K.)
| | - Sunitha M. Kasibhatla
- Bioinformatics Centre, Savitribai Phule Pune University (Formerly University of Pune), Pune 411007, India; (S.L.); (S.M.K.); (M.K.)
- HPC-Medical & Bioinformatics Applications Group, Centre for Development of Advanced Computing, Pune 411008, India
| | - Mukund Ramtirthkar
- Department of Statistics, Savitribai Phule Pune University (Formerly University of Pune), Pune 411007, India; (M.R.); (M.M.K.)
| | - Meenal Kinikar
- Bioinformatics Centre, Savitribai Phule Pune University (Formerly University of Pune), Pune 411007, India; (S.L.); (S.M.K.); (M.K.)
| | - Mohan M. Kale
- Department of Statistics, Savitribai Phule Pune University (Formerly University of Pune), Pune 411007, India; (M.R.); (M.M.K.)
| | - Urmila Kulkarni-Kale
- Bioinformatics Centre, Savitribai Phule Pune University (Formerly University of Pune), Pune 411007, India; (S.L.); (S.M.K.); (M.K.)
| |
Collapse
|
6
|
Zhu M, Kleepbua J, Guan Z, Chew SP, Tan JW, Shen J, Latthitham N, Hu J, Law JX, Li L. Early Spatiotemporal Patterns and Population Characteristics of the COVID-19 Pandemic in Southeast Asia. Healthcare (Basel) 2021; 9:1220. [PMID: 34574997 PMCID: PMC8466219 DOI: 10.3390/healthcare9091220] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2021] [Revised: 09/02/2021] [Accepted: 09/10/2021] [Indexed: 12/28/2022] Open
Abstract
This observational study aims to investigate the early disease patterns of coronavirus disease 2019 (COVID-19) in Southeast Asia, consequently providing historical experience for further interventions. Data were extracted from official websites of the WHO and health authorities of relevant countries. A total of 1346 confirmed cases of COVID-19, with 217 recoveries and 18 deaths, were reported in Southeast Asia as of 16 March 2020. The basic reproductive number (R0) of COVID-19 in the region was estimated as 2.51 (95% CI:2.31 to 2.73), and there were significant geographical variations at the subregional level. Early transmission dynamics were examined with an exponential regression model: y = 0.30e0.13x (p < 0.01, R2 = 0.96), which could help predict short-term incidence. Country-level disease burden was positively correlated with Human Development Index (r = 0.86, p < 0.01). A potential early shift in spatial diffusion patterns and a spatiotemporal cluster occurring in Malaysia and Singapore were detected. Demographic analyses of 925 confirmed cases indicated a median age of 44 years and a sex ratio (male/female) of 1.25. Age may play a significant role in both susceptibilities and outcomes. The COVID-19 situation in Southeast Asia is challenging and unevenly geographically distributed. Hence, enhanced real-time surveillance and more efficient resource allocation are urgently needed.
Collapse
Affiliation(s)
- Mingjian Zhu
- State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, National Clinical Research Center for Infectious Diseases, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou 310003, China; (M.Z.); (Z.G.); (J.S.)
| | - Jirapat Kleepbua
- Thammasat University Hospital, Pathum Thani 12120, Thailand; (J.K.); (N.L.)
| | - Zhou Guan
- State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, National Clinical Research Center for Infectious Diseases, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou 310003, China; (M.Z.); (Z.G.); (J.S.)
| | - Sien Ping Chew
- School of Medicine, Shanghai Jiao Tong University, Shanghai 200025, China;
| | - Joanna Weihui Tan
- Faculty of Arts and Social Sciences, National University of Singapore, Singapore 117570, Singapore;
| | - Jian Shen
- State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, National Clinical Research Center for Infectious Diseases, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou 310003, China; (M.Z.); (Z.G.); (J.S.)
| | | | - Jianxiong Hu
- Guangdong Provincial Institute of Public Health, Guangdong Provincial Center for Disease Control and Prevention, Guangzhou 511430, China;
| | - Jia Xian Law
- Tuanku Ja’afar Hospital, Seremban 70300, Malaysia;
| | - Lanjuan Li
- State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, National Clinical Research Center for Infectious Diseases, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou 310003, China; (M.Z.); (Z.G.); (J.S.)
| |
Collapse
|