Dataset Information

SNPs or RNA modifications? Concerns on mutation-based evolutionary studies of SARS-CoV-2.

ABSTRACT: SARS-CoV-2 is still rampaging throughout the world while the many evolutionary studies on it are simultaneously springing up. Researchers have simply utilized the public RNA-seq data to find out the so-called SNPs in the virus genome. The evolutionary analyses were largely based on these mutations. Here, we claim that we reliably detected A-to-G RNA modifications in the RNA-seq data of SARS-CoV-2 with high signal to noise ratios, presumably caused by the host's deamination enzymes. Intriguingly, since SARS-CoV-2 is an RNA virus, it is technically impossible to distinguish SNPs and RNA modifications from the RNA-seq data alone without solid evidence, making it difficult to tell the evolutionary patterns behind the mutation spectrum. Researchers should clarify their biological significance before they automatically regard the mutations as SNPs or RNA modifications. This is not a problem for DNA organisms but should be seriously considered when we are investigating the RNA viruses.

SUBMITTER: Li Y

PROVIDER: S-EPMC7454988 | biostudies-literature | 2020

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

SNPs or RNA modifications? Concerns on mutation-based evolutionary studies of SARS-CoV-2.

Li Yue Y Yang Xinai X Wang Na N Wang Haiyan H Yin Bin B Yang Xiaoping X Jiang Wenqing W

PloS one 20200828 8

SARS-CoV-2 is still rampaging throughout the world while the many evolutionary studies on it are simultaneously springing up. Researchers have simply utilized the public RNA-seq data to find out the so-called SNPs in the virus genome. The evolutionary analyses were largely based on these mutations. Here, we claim that we reliably detected A-to-G RNA modifications in the RNA-seq data of SARS-CoV-2 with high signal to noise ratios, presumably caused by the host's deamination enzymes. Intriguingly, ...[more]

PMID: 32857808

Similar Datasets

Project description:Since the emergence of SARS-CoV-2 in Wuhan, China more than a year ago, it has spread across the world in a very short span of time. Although, different forms of vaccines are being rolled out for vaccination programs around the globe, the mutation of the virus is still a cause of concern among the research communities. Hence, it is important to study the constantly evolving virus and its strains in order to provide a much more stable form of cure. This fact motivated us to conduct this research where we have initially carried out multiple sequence alignment of 15359 and 3033 global dataset without Indian and the dataset of exclusive Indian SARS-CoV-2 genomes respectively, using MAFFT. Subsequently, phylogenetic analyses are performed using Nextstrain to identify virus clades. Consequently, the virus strains are found to be distributed among 5 major clades or clusters viz. 19A, 19B, 20A, 20B and 20C. Thereafter, mutation points as SNPs are identified in each clade. Henceforth, from each clade top 10 signature SNPs are identified based on their frequency i.e. number of occurrences in the virus genome. As a result, 50 such signature SNPs are individually identified for global dataset without Indian and dataset of exclusive Indian SARS-CoV-2 genomes respectively. Out of each 50 signature SNPs, 39 and 41 unique SNPs are identified among which 25 non-synonymous signature SNPs (out of 39) resulted in 30 amino acid changes in protein while 27 changes in amino acid are identified from 22 non-synonymous signature SNPs (out of 41). These 30 and 27 amino acid changes for the non-synonymous signature SNPs are visualised in their respective protein structure as well. Finally, in order to judge the characteristics of the identified clades, the non-synonymous signature SNPs are considered to evaluate the changes in proteins as biological functions with the sequences using PROVEAN and PolyPhen-2 while I-Mutant 2.0 is used to evaluate their structural stability. As a consequence, for global dataset without Indian sequences, G251V in ORF3a in clade 19A, F308Y and G196V in NSP4 and ORF3a in 19B are the unique amino acid changes which are responsible for defining each clade as they are all deleterious and unstable. Such changes which are common for both global dataset without Indian and dataset of exclusive Indian sequences are R203M in Nucleocapsid for 20B, T85I and Q57H in NSP2 and ORF3a respectively for 20C while for exclusive Indian sequences such unique changes are A97V in RdRp, G339S and G339C in NSP2 in 19A and Q57H in ORF3a in 20A.

Project description:BackgroundSARS-CoV-2 is a RNA coronavirus responsible for the pandemic of the Severe Acute Respiratory Syndrome (COVID-19). RNA viruses are characterized by a high mutation rate, up to a million times higher than that of their hosts. Virus mutagenic capability depends upon several factors, including the fidelity of viral enzymes that replicate nucleic acids, as SARS-CoV-2 RNA dependent RNA polymerase (RdRp). Mutation rate drives viral evolution and genome variability, thereby enabling viruses to escape host immunity and to develop drug resistance.MethodsWe analyzed 220 genomic sequences from the GISAID database derived from patients infected by SARS-CoV-2 worldwide from December 2019 to mid-March 2020. SARS-CoV-2 reference genome was obtained from the GenBank database. Genomes alignment was performed using Clustal Omega. Mann-Whitney and Fisher-Exact tests were used to assess statistical significance.ResultsWe characterized 8 novel recurrent mutations of SARS-CoV-2, located at positions 1397, 2891, 14408, 17746, 17857, 18060, 23403 and 28881. Mutations in 2891, 3036, 14408, 23403 and 28881 positions are predominantly observed in Europe, whereas those located at positions 17746, 17857 and 18060 are exclusively present in North America. We noticed for the first time a silent mutation in RdRp gene in England (UK) on February 9th, 2020 while a different mutation in RdRp changing its amino acid composition emerged on February 20th, 2020 in Italy (Lombardy). Viruses with RdRp mutation have a median of 3 point mutations [range: 2-5], otherwise they have a median of 1 mutation [range: 0-3] (p value < 0.001).ConclusionsThese findings suggest that the virus is evolving and European, North American and Asian strains might coexist, each of them characterized by a different mutation pattern. The contribution of the mutated RdRp to this phenomenon needs to be investigated. To date, several drugs targeting RdRp enzymes are being employed for SARS-CoV-2 infection treatment. Some of them have a predicted binding moiety in a SARS-CoV-2 RdRp hydrophobic cleft, which is adjacent to the 14408 mutation we identified. Consequently, it is important to study and characterize SARS-CoV-2 RdRp mutation in order to assess possible drug-resistance viral phenotypes. It is also important to recognize whether the presence of some mutations might correlate with different SARS-CoV-2 mortality rates.

Dataset Information

SNPs or RNA modifications? Concerns on mutation-based evolutionary studies of SARS-CoV-2.

Publications

SNPs or RNA modifications? Concerns on mutation-based evolutionary studies of SARS-CoV-2.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets