MBE Advance Access published online on January 30, 2009
Molecular Biology and Evolution, doi:10.1093/molbev/msp009
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Article |
Phylogenetic Relationships and Molecular Adaptation Dynamics of Human Rhinoviruses
a Department of Microbiology and Molecular Biology, Brigham Young University, Provo, UT 84602, USA
b Department of Biology, Brigham Young University, Provo, UT 84602, USA
* Corresponding author. Tel.: +1 801-422-1733; fax: +1 801-422-0519. E-mail address: nlewisrogers{at}yahoo.com (N. Lewis-Rogers).
Received for publication August 8, 2008. Revision received December 5, 2008. Accepted for publication December 10, 2008.
Human rhinoviruses (HRVs) are responsible for nearly 50% of all common cold infections. Ordinarily, HRV infections are mild and self-limiting; nonetheless, every year they result in significant loss of economic productivity and substantial inappropriate antibiotic use. Development of effective vaccine and antiviral prophylaxis against HRV has been hampered by the extensive antigenic diversity present among the nearly 100 serotypes. To gain new insights into the evolutionary processes that create the genetic diversity present among HRVs, we tested for recombination and selection for individual genes and the coding genome for 45 HRV serotypes using estimated phylogenetic relationships. Although the structural capsid genes and nonstructural genes recovered incongruent tree topologies, no recombination was detected using substitution methods. Therefore, the coding genome was determined appropriate for phylogenetic tests. Results of the Shimodaira-Hasegawa (SH) test support the hypothesis that the capsid genes recover a different evolutionary history than the nonstructural genes. Our best phylogenetic estimate based on the coding genome suggests that HRV-B is more closely related to enterovirus than to HRV-A; however, several alternative phylogenetic hypotheses were not rejected by the SH test. Positive selection was examined by using two different approaches; dN/dS rate ratio and the physico-chemical phenotypes for 31 amino acid properties. Analyses using dN/dS failed to detect positive selection. However, protein phenotypic expression appears to be a more sensitive approach. There was extensive stabilizing and destabilizing-positive selection in HRV-A major and HRV-B serotypes for all proteins, except in 3A in HRV-B, which overlapped with functional, structural, and to a greater extent in uncharacterized genomic regions. In contrast, the evolution of HRV-A minor serotypes appears to be driven primarily by destabilizing selection. Our results demonstrate that HRV-A major, HRV-A minor and HRV-B serotypes have not been similarly influenced by purifying selection.
Key Words: human rhinovirus • enterovirus • selection • adaptive evolution • recombination