Dataset Information

Bioinformatic analysis of the whole genome sequences of SARS-CoV-2 from Indonesia.

ABSTRACT:

Background and objectives

In first May 2020, Indonesia has been successfully submitted the first three full-length sequence of SARS-CoV-2 to GISAID database. Until September 10^th, 2020, Indonesia had submitted 54 WGS. In this study, we have analyzed and annotated SARS-CoV-2 mutations in spike protein and main proteases.

Materials and methods

The Whole Genome Sequence (WGS) of Indonesia were obtained from GISAID data base. The 54 data were taken from March to September 10^th, 2020. The sequences corresponded to Spike Protein (SP), 3-chymotrypsin like protease (3CLpro), and papain like protease (PLpro) were selected. The Wuhan genome was used as reference.

Results

In total WGS from Indonesia, we found 5 major clades, which dominated as G clade, where the mutation of D614G was found. This D614G was identified as much as 59%, which mostly reported in late samples submitted. Beside D614G mutation, we report three unique mutations: A352S, S477I, and Q677H. Besides, some mutations were also detected in two domains that were expected to be conserved region, the main viral proteases: PLpro (P77L and V205I), 3CLpro (M49I and L50F).

Conclusion

The analysis of SARS-CoV-2 from WGS Indonesia showed a high genetic variation. The diversity in SARS-CoV-2 may epidemiologically enhance virulence and transmission of this virus. The prevalence of D614G over the time in different locations, indicating that changes in this mutation may related to host infection and the viral transmission. However, some mutations that have been reported in this study were not eligible for the most stable conformation.

SUBMITTER: Ulfah M

PROVIDER: S-EPMC8408024 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Background & objectivesSeveral phylogenetic classification systems have been devised to trace the viral lineages of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, inconsistency in the nomenclature limits uniformity in its epidemiological understanding. This study provides an integration of existing classifications and describes evolutionary trends of the SARS-CoV-2 strains circulating in India.MethodsThe whole genomes of 330 SARS-CoV-2 samples were sequenced using next-generation sequencing (NGS). Phylogenetic and sequence analysis of a total of 3014 Indian SARS-CoV-2 sequences from 20 different States/Union Territories (January to September 2020) from the Global Initiative on Sharing All Influenza Data (GISAID) database was performed to observe the clustering of Nextstrain and Phylogenetic Assignment of Named Global Outbreak LINeages (Pangolin) lineages with the GISAID clades. The identification of mutational sites under selection pressure was performed using Mixed Effects Model of Evolution and Single-Likelihood Ancestor Counting methods available in the Datamonkey server.ResultsTemporal data of the Indian SARS-CoV-2 genomes revealed that except for Uttarakhand, West Bengal and Haryana that showed the circulation of GISAID clade O even after July 2020, the rest of the States showed a complete switch to GR/GH clades. Pangolin lineages B.1.1.8 and B.1.113 identified within GR and GH clades, respectively, were noted to be indigenous evolutions. Sites identified to be under positive selection pressure within these clades were found to occur majorly in the non-structural proteins coded by ORF1a and ORF1b.Interpretation & conclusionsThis study interpreted the geographical and temporal dominance of SARS-CoV-2 strains in India over a period of nine months based on the GISAID classification. An integration of the GISAID, Nextstrain and Pangolin classifications is also provided. The emergence of new lineages B.1.1.8 and B.1.113 was indicative of host-specific evolution of the SARS-CoV-2 strains in India. The hotspot mutations such as those driven by positive selection need to be further characterized.