Phylogenetic Patterns of Human Coxsackievirus B5 Arise from Population Dynamics between Two Genogroups and Reveal Evolutionary Factors of Molecular Adaptation and Transmission

Abstract

The aim of this study was to gain insights into the tempo and mode of the evolutionary processes that sustain genetic diversity in coxsackievirus B5 (CVB5) and into the interplay with virus transmission. We estimated phylodynamic patterns with a large sample of virus strains collected in Europe by Bayesian statistical methods, reconstructed the ancestral states of genealogical nodes, and tested for selection. The genealogies estimated with the structural one-dimensional gene encoding the VP1 protein and nonstructural 3CD locus allowed the precise description of lineages over time and cocirculating virus populations within the two CVB5 clades, genogroups A and B. Strong negative selection shaped the evolution of both loci, but compelling phylogenetic data suggested that immune selection pressure resulted in the emergence of the two genogroups with opposed evolutionary pathways. The genogroups also differed in the temporal occurrence of the amino acid changes. The virus strains of genogroup A were characterized by sequential acquisition of nonsynonymous changes in residues exposed at the virus 5-fold axis. The genogroup B viruses were marked by selection of three changes in a different domain (VP1 C terminus) during its early emergence. These external changes resulted in a selective sweep, which was followed by an evolutionary stasis that is still ongoing after 50 years. The inferred population history of CVB5 showed an alternation of the prevailing genogroup during meningitis epidemics across Europe and is interpreted to be a consequence of partial cross-immunity.