Genome-wide analysis of SARS-CoV-2 virus strains circulating worldwide implicates heterogeneity

Sci Rep. 2020 Aug 19;10(1):14004. doi: 10.1038/s41598-020-70812-6.

Abstract

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), a novel evolutionary divergent RNA virus, is responsible for the present devastating COVID-19 pandemic. To explore the genomic signatures, we comprehensively analyzed 2,492 complete and/or near-complete genome sequences of SARS-CoV-2 strains reported from across the globe to the GISAID database up to 30 March 2020. Genome-wide annotations revealed 1,516 nucleotide-level variations at different positions throughout the entire genome of SARS-CoV-2. Moreover, nucleotide (nt) deletion analysis found twelve deletion sites throughout the genome other than previously reported deletions at coding sequence of the ORF8 (open reading frame), spike, and ORF7a proteins, specifically in polyprotein ORF1ab (n = 9), ORF10 (n = 1), and 3´-UTR (n = 2). Evidence from the systematic gene-level mutational and protein profile analyses revealed a large number of amino acid (aa) substitutions (n = 744), demonstrating the viral proteins heterogeneous. Notably, residues of receptor-binding domain (RBD) showing crucial interactions with angiotensin-converting enzyme 2 (ACE2) and cross-reacting neutralizing antibody were found to be conserved among the analyzed virus strains, except for replacement of lysine with arginine at 378th position of the cryptic epitope of a Shanghai isolate, hCoV-19/Shanghai/SH0007/2020 (EPI_ISL_416320). Furthermore, our results of the preliminary epidemiological data on SARS-CoV-2 infections revealed that frequency of aa mutations were relatively higher in the SARS-CoV-2 genome sequences of Europe (43.07%) followed by Asia (38.09%), and North America (29.64%) while case fatality rates remained higher in the European temperate countries, such as Italy, Spain, Netherlands, France, England and Belgium. Thus, the present method of genome annotation employed at this early pandemic stage could be a promising tool for monitoring and tracking the continuously evolving pandemic situation, the associated genetic variants, and their implications for the development of effective control and prophylaxis strategies.

MeSH terms

  • Amino Acid Sequence / genetics
  • Angiotensin-Converting Enzyme 2
  • Antibodies, Neutralizing / immunology
  • Base Pair Mismatch
  • Base Sequence / genetics
  • Betacoronavirus / classification*
  • Betacoronavirus / genetics*
  • COVID-19
  • Climate
  • Coronavirus Infections / epidemiology*
  • Coronavirus Infections / virology
  • Genetic Heterogeneity*
  • Genome, Viral / genetics*
  • Genome-Wide Association Study / methods*
  • Global Health*
  • Humans
  • Open Reading Frames / genetics
  • Pandemics
  • Peptidyl-Dipeptidase A / metabolism
  • Phylogeny
  • Pneumonia, Viral / epidemiology*
  • Pneumonia, Viral / virology
  • Protein Domains / genetics
  • Protein Domains / immunology
  • SARS-CoV-2
  • Sequence Deletion
  • Spike Glycoprotein, Coronavirus / chemistry
  • Spike Glycoprotein, Coronavirus / genetics
  • Spike Glycoprotein, Coronavirus / metabolism

Substances

  • Antibodies, Neutralizing
  • Spike Glycoprotein, Coronavirus
  • spike protein, SARS-CoV-2
  • Peptidyl-Dipeptidase A
  • ACE2 protein, human
  • Angiotensin-Converting Enzyme 2