Dataset Information

Comparative Genomic Analyses Reveal a Specific Mutation Pattern Between Human Coronavirus SARS-CoV-2 and Bat-CoV RaTG13

ABSTRACT:

Background

The outbreak of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Wuhan, China, rapidly grew into a global pandemic. How SARS-CoV-2 evolved remains unclear.

Methods

We performed a comprehensive analysis using the available genomes of SARS-CoV-2 and its closely related coronaviruses.

Results

The ratio of nucleotide substitutions to amino acid substitutions of the spike gene (9.07) between SARS-CoV-2 WIV04 and Bat-CoV RaTG13 was markedly higher than that between other coronaviruses (range, 1.29–4.81); the ratio of non-synonymous to synonymous substitution rates (dN/dS) between SARS-CoV-2 WIV04 and Bat-CoV RaTG13 was the lowest among all the performed comparisons, suggesting evolution under stringent selective pressure. Notably, the relative proportion of the T:C transition was markedly higher between SARS-CoV-2 WIV04 and Bat-CoV RaTG13 than between other compared coronaviruses. Codon usage is similar across these coronaviruses and is unlikely to explain the increased number of synonymous mutations. Moreover, some sites of the spike protein might be subjected to positive selection.

Conclusions

Our results showed an increased proportion of synonymous substitutions and the T:C transition between SARS-CoV-2 and RaTG13. Further investigation of the mutation pattern mechanism would contribute to understanding viral pathogenicity and its adaptation to hosts.

SUBMITTER: Lv L

PROVIDER: S-EPMC7793997 | biostudies-literature | 2020 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The evolutional process of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) development remains inconclusive. This study compared the genome sequences of severe acute respiratory syndrome coronavirus (SARS-CoV), bat coronavirus RaTG13, and SARS-CoV-2. In total, the genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. A total of 3.6% (1,068 bases) of the SARS-CoV-2 genome was derived from insertion and/or deletion (indel) mutations, and 18.6% (5,548 bases) was from point mutations from the genome of SARS-CoV. At least 35 indel sites were confirmed in the genome of SARS-CoV-2, in which 17 were with ≥10 consecutive bases long. Ten of these relatively long indels were located in the spike (S) gene, five in nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, and one in ORF8 and noncoding region. Seventeen (48.6%) of the 35 indels were based on insertion-and-deletion mutations with exchanged gene sequences of 7-325 consecutive bases. Almost the complete ORF8 gene was replaced by a single 325 consecutive base-long indel. The distribution of these indels was roughly in accordance with the distribution of the rate of point mutation rate around the indels. The genome sequence of SARS-CoV-2 was 96.0% identical to that of RaTG13. There was no long insertion-and-deletion mutation between the genomes of RaTG13 and SARS-CoV-2. The findings of the uneven distribution of multiple indels and the presence of multiple long insertion-and-deletion mutations with exchanged consecutive base sequences in the viral genome may provide insights into SARS-CoV-2 development. IMPORTANCE The developmental mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains inconclusive. This study compared the base sequence one-by-one between severe acute respiratory syndrome coronavirus (SARS-CoV) or bat coronavirus RaTG13 and SARS-CoV-2. The genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. Seventeen of the 35 sites with insertion and/or deletion mutations between SARS-CoV-2 and SARS-CoV were based on insertion-and-deletion mutations with the replacement of 7-325 consecutive bases. Most of these long insertion-and-deletion sites were concentrated in the nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, S1 domain of the spike protein, and ORF8 genes. Such long insertion-and-deletion mutations were not observed between the genomes of RaTG13 and SARS-CoV-2. The presence of multiple long insertion-and-deletion mutations in the genome of SARS-CoV-2 and their uneven distributions may provide further insights into the development of the virus.

Project description:Coronaviruses (CoVs), including severe acute respiratory syndrome CoV (SARS-CoV), Middle East respiratory syndrome CoV (MERS-CoV), and SARS-CoV-2, produce double-stranded RNA (dsRNA) that activates antiviral pathways such as PKR and OAS/RNase L. To successfully replicate in hosts, viruses must evade such antiviral pathways. Currently, the mechanism of how SARS-CoV-2 antagonizes dsRNA-activated antiviral pathways is unknown. In this study, we demonstrate that the SARS-CoV-2 nucleocapsid (N) protein, the most abundant viral structural protein, is capable of binding to dsRNA and phosphorylated PKR, inhibiting both the PKR and OAS/RNase L pathways. The N protein of the bat coronavirus (bat-CoV) RaTG13, the closest relative of SARS-CoV-2, has a similar ability to inhibit the human PKR and RNase L antiviral pathways. Via mutagenic analysis, we found that the C-terminal domain (CTD) of the N protein is sufficient for binding dsRNA and inhibiting RNase L activity. Interestingly, while the CTD is also sufficient for binding phosphorylated PKR, the inhibition of PKR antiviral activity requires not only the CTD but also the central linker region (LKR). Thus, our findings demonstrate that the SARS-CoV-2 N protein is capable of antagonizing the two critical antiviral pathways activated by viral dsRNA and that its inhibition of PKR activities requires more than dsRNA binding mediated by the CTD. IMPORTANCE The high transmissibility of SARS-CoV-2 is an important viral factor defining the coronavirus disease 2019 (COVID-19) pandemic. To transmit efficiently, SARS-CoV-2 must be capable of disarming the innate immune response of its host efficiently. Here, we describe that the nucleocapsid protein of SARS-CoV-2 is capable of inhibiting two critical innate antiviral pathways, PKR and OAS/RNase L. Moreover, the counterpart of the closest animal coronavirus relative of SARS-CoV-2, bat-CoV RaTG13, can also inhibit human PKR and OAS/RNase L antiviral activities. Thus, the importance of our discovery for understanding the COVID-19 pandemic is 2-fold. First, the ability of SARS-CoV-2 N to inhibit innate antiviral activity is likely a factor contributing to the transmissibility and pathogenicity of the virus. Second, the bat relative of SARS-CoV-2 has the capacity to inhibit human innate immunity, which thus likely contributed to the establishment of infection in humans. The findings described in this study are valuable for developing novel antivirals and vaccines.

Project description:BackgroundClostridium botulinum is a diverse group of bacteria characterized by the production of botulinum neurotoxin. Botulinum neurotoxins are classified into serotypes (BoNT/A-G), which are produced by six species/Groups of Clostridia, but the genetic background of the bacteria remains poorly understood. The purpose of this study was to use comparative genomics to provide insights into the genetic diversity and evolutionary history of bacteria that produce the potent botulinum neurotoxin.ResultsComparative genomic analyses of over 170 Clostridia genomes, including our draft genome assemblies for 59 newly sequenced Clostridia strains from six continents and publicly available genomic data, provided in-depth insights into the diversity and distribution of BoNT-producing bacteria. These newly sequenced strains included Group I and II strains that express BoNT/A,/B,/E, or/F as well as bivalent strains. BoNT-producing Clostridia and closely related Clostridia species were delineated with a variety of methods including 16S rRNA gene, concatenated marker genes, core genome and concatenated multi-locus sequencing typing (MLST) gene phylogenies that related whole genome sequenced strains to publicly available strains and sequence types. These analyses illustrated the phylogenetic diversity in each Group and the diversity of genomic backgrounds that express the same toxin type or subtype. Comparisons of the botulinum neurotoxin genes did not identify novel toxin types or variants.ConclusionsThis study represents one of the most comprehensive analyses of whole genome sequence data for Group I and II BoNT-producing strains. Read data and draft genome assemblies generated for 59 isolates will be a resource to the research community. Core genome phylogenies proved to be a powerful tool for differentiating BoNT-producing strains and can provide a framework for the study of these bacteria. Comparative genomic analyses of Clostridia species illustrate the diversity of botulinum-neurotoxin-producing strains and the plasticity of the genomic backgrounds in which bont genes are found.