The pandemic caused by the spread of SARS-CoV-2 has led to considerable interest in its evolutionary origin and genome structure. Here, we analyzed mutation patterns in 34 human SARS-CoV-2 isolates and a closely related RaTG13 isolated from Rhinolophus affinis (a horseshoe bat). We also evaluated the CpG dinucleotide contents in SARS-CoV-2 and other human and animal coronavirus genomes. Out of 1136 single nucleotide variations (~4% divergence) between human SARS-CoV-2 and bat RaTG13, 682 (60%) can be attributed to C>U and U>C substitutions, far exceeding other types of substitutions. An accumulation of C>U mutations was also observed in SARS-CoV2 variants that arose within the human population. Globally, the C>U substitutions increased the frequency of codons for hydrophobic amino acids in SARS-CoV-2 peptides, while U>C substitutions decreased it. In contrast to most other coronaviruses, both SARS-CoV-2 and RaTG13 exhibited CpG depletion in their genomes. The data suggest that C-to-U conversion mediated by C deamination played a significant role in the evolution of the SARS-CoV-2 coronavirus. We hypothesize that the high frequency C>U transitions reflect virus adaptation processes in their hosts, and that SARS-CoV-2 could have been evolving for a relatively long period in humans following the transfer from animals before spreading worldwide.

Laboratory of Molecular Epigenetics Institute of Biophysics Academy of Sciences of the Czech Republic Královopolská 135 61265 Brno Czech Republic

Zobrazit více v PubMed

Zhou P., Yang X.L., Wang X.G., Hu B., Zhang L., Zhang W., Si H.R., Zhu Y., Li B., Huang C.L., et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020;579:270–273. doi: 10.1038/s41586-020-2012-7. PubMed DOI PMC

Kandeel M., Ibrahim A., Fayez M., Al-Nazawi M. From SARS and MERS CoVs to SARS-CoV-2: Moving toward more biased codon usage in viral structural and nonstructural genes. J. Med. Virol. 2020 doi: 10.1002/jmv.25754. PubMed DOI PMC

Wan Y., Shang J., Graham R., Baric R.S., Li F. Receptor recognition by the novel coronavirus from Wuhan: An analysis based on decade-long structural studies of SARS coronavirus. J. Virol. 2020;94 doi: 10.1128/JVI.00127-20. PubMed DOI PMC

Wrapp D., Wang N., Corbett K.S., Goldsmith J.A., Hsieh C.L., Abiona O., Graham B.S., McLellan J.S. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science. 2020;367:1260–1263. doi: 10.1126/science.abb2507. PubMed DOI PMC

Zhang T., Wu Q., Zhang Z. Probable pangolin origin of SARS-CoV-2 associated with the COVID-19 outbreak. Curr. Biol. 2020 doi: 10.1016/j.cub.2020.03.022. PubMed DOI PMC

Andersen K.G., Rambaut A., Lipkin W.I., Holmes E.C., Garry R.F. The proximal origin of SARS-CoV-2. Nat. Med. 2020;26:450–452. doi: 10.1038/s41591-020-0820-9. PubMed DOI PMC

Frederico L.A., Kunkel T.A., Shaw B.R. Cytosine deamination in mismatched base pairs. Biochemistry. 1993;32:6523–6530. doi: 10.1021/bi00077a005. PubMed DOI

Poole A., Penny D., Sjoberg B.M. Confounded cytosine! Tinkering and the evolution of DNA. Nat. Rev. Mol. Cell Biol. 2001;2:147–151. doi: 10.1038/35052091. PubMed DOI

Jabbari K., Bernardi G. Cytosine methylation and CpG, TpG (CpA) and TpA frequencies. Gene. 2004;333:143–149. doi: 10.1016/j.gene.2004.02.043. PubMed DOI

Stevens M., Cheng J.B., Li D., Xie M., Hong C., Maire C.L., Ligon K.L., Hirst M., Marra M.A., Costello J.F., et al. Estimating absolute methylation levels at single-CpG resolution from methylation enrichment and restriction enzyme sequencing methods. Genome Res. 2013;23:1541–1553. doi: 10.1101/gr.152231.112. PubMed DOI PMC

Bird A. DNA methylation patterns and epigenetic memory. Genes Dev. 2002;16:6–21. doi: 10.1101/gad.947102. PubMed DOI

Kypr J., Mrazek J., Reich J. Nucleotide composition bias and CpG dinucleotide content in the genomes of HIV and HTLV 1/2. Biochim. Biophys. Acta. 1989;1009:280–282. doi: 10.1016/0167-4781(89)90114-0. PubMed DOI

Cheng X., Virk N., Chen W., Ji S., Sun Y., Wu X. CpG usage in RNA viruses: Data and hypotheses. PLoS ONE. 2013;8:e74109. doi: 10.1371/journal.pone.0074109. PubMed DOI PMC

Upadhyay M., Samal J., Kandpal M., Vasaikar S., Biswas B., Gomes J., Vivekanandan P. CpG dinucleotide frequencies reveal the role of host methylation capabilities in parvovirus evolution. J. Virol. 2013;87:13816–13824. doi: 10.1128/JVI.02515-13. PubMed DOI PMC

Alinejad-Rokny H., Anwar F., Waters S.A., Davenport M.P., Ebrahimi D. Source of CpG depletion in the HIV-1 genome. Mol. Biol. Evol. 2016;33:3205–3212. doi: 10.1093/molbev/msw205. PubMed DOI

Gojobori T., Li W.H., Graur D. Patterns of nucleotide substitution in pseudogenes and functional genes. J. Mol. Evol. 1982;18:360–369. doi: 10.1007/BF01733904. PubMed DOI

Ochman H. Neutral mutations and neutral substitutions in bacterial genomes. Mol. Biol. Evol. 2003;20:2091–2096. doi: 10.1093/molbev/msg229. PubMed DOI

Wang W.C., Ma L., Becher H., Garcia S., Kovarikova A., Leitch I.J., Leitch A.R., Kovarik A. Astonishing 35S rDNA diversity in the gymnosperm species Cycas revoluta Thunb. Chromosoma. 2016;125:683–699. doi: 10.1007/s00412-015-0556-3. PubMed DOI PMC

Bulmer M. Neighboring base effects on substitution rates in pseudogenes. Mol. Biol. Evol. 1986;3:322–329. PubMed

Keller I., Bensasson D., Nichols R.A. Transition-transversion bias is not universal: A counter example from grasshopper pseudogenes. PLoS Genet. 2007;3:e22. doi: 10.1371/journal.pgen.0030022. PubMed DOI PMC

Duchene S., Ho S.Y., Holmes E.C. Declining transition/transversion ratios through time reveal limitations to the accuracy of nucleotide substitution models. BMC Evol. Biol. 2015;15:36. doi: 10.1186/s12862-015-0312-6. PubMed DOI PMC

Afgan E., Baker D., Batut B., van den Beek M., Bouvier D., Cech M., Chilton J., Clements D., Coraor N., Grüning B.A., et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucl. Acids Res. 2018;46:W537–W544. doi: 10.1093/nar/gky379. PubMed DOI PMC

Rozas J., Rozas R. DnaSP version 2.0: A novel software package for extensive molecular population genetics analysis. Comput. Appl. Biosci. 1997;13:307–311. PubMed

Pettersen E.F., Goddard T.D., Huang C.C., Couch G.S., Greenblatt D.M., Meng E.C., Ferrin T.E. UCSF Chimera–A visualization system for exploratory research and analysis. J. Comput. Chem. 2004;25:1605–1612. doi: 10.1002/jcc.20084. PubMed DOI

Kyte J., Doolittle R.F. A simple method for displaying the hydropathic character of a protein. J. Mol. Biol. 1982;157:105–132. doi: 10.1016/0022-2836(82)90515-0. PubMed DOI

Gouy M., Guindon S., Gascuel O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol. Biol. Evol. 2010;27:221–224. doi: 10.1093/molbev/msp259. PubMed DOI

BoxPlotR: A Web-Tool for Generation of Box Plots. [(accessed on 2 May 2020)]; Available online: http://shiny.chemgrid.org/boxplotr/

Tang X., Wu C., Li X., Song Y., Yao X., Wu X., Duan Y., Zhang H., Wang Y., Qian Z., et al. On the origin and continuing evolution of SARS-CoV-2. Natl. Sci. Rev. 2020;7:1012–1023. doi: 10.1093/nsr/nwaa036. PubMed DOI PMC

Anisimova M., Gascuel O. Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst. Biol. 2006;55:539–552. doi: 10.1080/10635150600755453. PubMed DOI

Nei M., Gojobori T. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol. Biol. Evol. 1986;3:418–426. PubMed

Duffy S. Why are RNA virus mutation rates so damn high? PLoS Biol. 2018;16 doi: 10.1371/journal.pbio.3000003. PubMed DOI PMC

Pachetti M., Marini B., Benedetti F., Giudici F., Mauro E., Storici P., Masciovecchio C., Angeletti S., Ciccozzi M., Gallo R.C., et al. Emerging SARS-CoV-2 mutation hot spots include a novel RNA-dependent-RNA polymerase variant. J. Transl. Med. 2020;18:179. doi: 10.1186/s12967-020-02344-6. PubMed DOI PMC

Minskaia E., Hertzig T., Gorbalenya A.E., Campanacci V., Cambillau C., Canard B., Ziebuhr J. Discovery of an RNA virus 3′->5′ exoribonuclease that is critically involved in coronavirus RNA synthesis. Proc. Natl. Acad. Sci. USA. 2006;103:5108–5113. doi: 10.1073/pnas.0508200103. PubMed DOI PMC

Sethna P.B., Hofmann M.A., Brian D.A. Minus-strand copies of replicating coronavirus mRNAs contain antileaders. J. Virol. 1991;65:320–325. doi: 10.1128/JVI.65.1.320-325.1991. PubMed DOI PMC

Zhao Z., Li H., Wu X., Zhong Y., Zhang K., Zhang Y.P., Boerwinkle E., Fu Y.X. Moderate mutation rate in the SARS coronavirus genome and its implications. BMC Evol. Biol. 2004;4:21. doi: 10.1186/1471-2148-4-21. PubMed DOI PMC

Roberts S.A., Lawrence M.S., Klimczak L.J., Grimm S.A., Fargo D., Stojanov P., Kiezun A., Kryukov G.V., Carter S.L., Saksena G., et al. An APOBEC cytidine deaminase mutagenesis pattern is widespread in human cancers. Nat. Genet. 2013;45:970–976. doi: 10.1038/ng.2702. PubMed DOI PMC

Opi S., Takeuchi H., Kao S., Khan M.A., Miyagi E., Goila-Gaur R., Iwatani Y., Levin J.G., Strebel K. Monomeric APOBEC3G is catalytically active and has antiviral activity. J. Virol. 2006;80:4673–4682. doi: 10.1128/JVI.80.10.4673-4682.2006. PubMed DOI PMC

Upadhyay M., Vivekanandan P. Depletion of CpG Dinucleotides in papillomaviruses and polyomaviruses: A role for divergent evolutionary pressures. PLoS ONE. 2015;10:e0142368. doi: 10.1371/journal.pone.0142368. PubMed DOI PMC

Klinman D.M., Yamshchikov G., Ishigatsubo Y. Contribution of CpG motifs to the immunogenicity of DNA vaccines. J. Immunol. 1997;158:3635–3639. PubMed

Zhang J. Rates of conservative and radical nonsynonymous nucleotide substitutions in mammalian nuclear genes. J. Mol. Evol. 2000;50:56–68. doi: 10.1007/s002399910007. PubMed DOI

Cytosine Deamination and Evolution. [(accessed on 9 January 2020)]; Available online: https://designmatrix.wordpress.com/2009/02/22/cytosine-deamination-and-evolution/

G-quadruplexes in the evolution of hepatitis B virus

Letter to the Editor: Significant mutation enrichment in inverted repeat sites of new SARS-CoV-2 strains

In silico analysis of RNA-dependent RNA polymerase of the SARS-CoV-2 and therapeutic potential of existing antiviral drugs

Mutational Asymmetries in the SARS-CoV-2 Genome May Lead to Increased Hydrophobicity of Virus Proteins

SARS-CoV-2 hot-spot mutations are significantly enriched within inverted repeats and CpG island loci