Abstract

Proteins that are exposed on the surface of a virus are frequently subject to strong selection to escape from neutralizing antibodies. To investigate whether surface-exposed (SE) and non-exposed (NE) proteins encoded by RNA viruses exhibit different patterns of evolution under selection, we analyzed 244 protein-coding genes from 28 species of RNA viruses representing 15 taxonomic families. First, we show that gene-wide rates of non-synonymous (dN) and synonymous (dS) substitutions do not differentiate between SE and NE proteins. To incorporate variation in substitution rates among codon sites, we inferred the posterior distribution over a fixed grid of dN and dS rates for each alignment. This ‘evolutionary fingerprint’ provides a common framework for comparing the selection profiles of non-homologous genes. Next, we computed the Wasserstein distance for every pair of fingerprints, which is analogous to amount of work required to reshape one distribution to another. After compensating for differences in genetic variation among alignments, we found a small but significant difference between the fingerprints of SE and NE proteins (PERMANOVA, P = 0.03). However, we observed larger and more significant effects of whether the virus is enveloped (P < 10-5) and the interaction between these factors (P=6.9×10-4). The latter effects were driven by high levels of purifying selection in capsid proteins of Picornaviruses. Furthermore, greater amounts of variation in fingerprints were explained by significant differences among virus families and modes of transmission (P < 10-5). These results imply the pattern of selection on a virus protein is shaped more by characteristics of the virus than the protein itself.