The footprint of genome architecture in the largest genome expansion in RNA viruses

PLoS Pathog. 2013;9(7):e1003500. doi: 10.1371/journal.ppat.1003500. Epub 2013 Jul 18.

Abstract

The small size of RNA virus genomes (2-to-32 kb) has been attributed to high mutation rates during replication, which is thought to lack proof-reading. This paradigm is being revisited owing to the discovery of a 3'-to-5' exoribonuclease (ExoN) in nidoviruses, a monophyletic group of positive-stranded RNA viruses with a conserved genome architecture. ExoN, a homolog of canonical DNA proof-reading enzymes, is exclusively encoded by nidoviruses with genomes larger than 20 kb. All other known non-segmented RNA viruses have smaller genomes. Here we use evolutionary analyses to show that the two- to three-fold expansion of the nidovirus genome was accompanied by a large number of replacements in conserved proteins at a scale comparable to that in the Tree of Life. To unravel common evolutionary patterns in such genetically diverse viruses, we established the relation between genomic regions in nidoviruses in a sequence alignment-free manner. We exploited the conservation of the genome architecture to partition each genome into five non-overlapping regions: 5' untranslated region (UTR), open reading frame (ORF) 1a, ORF1b, 3'ORFs (encompassing the 3'-proximal ORFs), and 3' UTR. Each region was analyzed for its contribution to genome size change under different models. The non-linear model statistically outperformed the linear one and captured >92% of data variation. Accordingly, nidovirus genomes were concluded to have reached different points on an expansion trajectory dominated by consecutive increases of ORF1b, ORF1a, and 3'ORFs. Our findings indicate a unidirectional hierarchical relation between these genome regions, which are distinguished by their expression mechanism. In contrast, these regions cooperate bi-directionally on a functional level in the virus life cycle, in which they predominantly control genome replication, genome expression, and virus dissemination, respectively. Collectively, our findings suggest that genome architecture and the associated region-specific division of labor leave a footprint on genome expansion and may limit RNA genome size.

Publication types

  • Comparative Study
  • Research Support, Non-U.S. Gov't

MeSH terms

  • 3' Untranslated Regions
  • 5' Untranslated Regions
  • Base Sequence
  • Conserved Sequence
  • Databases, Nucleic Acid
  • Exoribonucleases / chemistry
  • Exoribonucleases / genetics
  • Exoribonucleases / metabolism
  • Genome Size
  • Genome, Viral*
  • Models, Biological*
  • Nidovirales / enzymology
  • Nidovirales / metabolism
  • Open Reading Frames
  • Phylogeny*
  • RNA Viruses / enzymology
  • RNA Viruses / genetics*
  • Viral Proteins / chemistry
  • Viral Proteins / genetics
  • Viral Proteins / metabolism

Substances

  • 3' Untranslated Regions
  • 5' Untranslated Regions
  • Viral Proteins
  • Exoribonucleases

Grants and funding

This research has received funding from the Program of Japan Initiative for Global Research Network on Infectious Diseases (J-GRID), MEXT, Japan, the European Union Seventh Framework Programme (FP7/2007-2013) under the program SILVER (grant agreement no. 260644), The Netherlands Bioinformatics Centre (BioRange SP3.2.2), the Collaborative Agreement in Bioinformatics between Leiden University Medical Center and Moscow State University (MoBiLe), and Leiden University Fund. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.