Molecular evolution of the fusion protein gene in human respiratory syncytial virus subgroup A

We studied the molecular evolution of the fusion protein (F) gene in the human respiratory syncytial virus subgroup A (HRSV-A). We performed time-scaled phylogenetic analyses using the Bayesian Markov chain Monte Carlo (MCMC) method. We also conducted genetic distance (p-distance), positive/negative selection, and Bayesian skyline plot analyses. Furthermore, we mapped the amino acid substitutions of the protein. The MCMC-constructed tree indicated that the HRSV F gene diverged from the bovine RSV (BRSV) gene approximately 550years ago and had a relatively low substitution rate (7.59×10(-4) substitutions/site/year). Moreover, a common ancestor of HRSV-A and -B diverged approximately 280years ago, which has since formed four distinct clusters. The present HRSV-A strains were assigned six genotypes based on F gene sequences and attachment glycoprotein gene sequences. The present strains exhibited high F gene sequence similarity values and low genetic divergence. No positive selection sites were identified; however, 50 negative selection sites were identified. F protein amino acid substitutions at 17 sites were distributed in the F protein. The effective population size of the gene has remained relatively constant, but the population size of the prevalent genotype (GA2) has increased in the last 10years. These results suggest that the HRSV-AF gene has evolved independently and formed some genotypes.