We performed evolution, phylodynamics, and reinfection-related antigenicity analyses of respiratory syncytial virus subgroup A (RSV-A) fusion (F) gene in globally collected strains (1465 strains) using authentic bioinformatics methods. The time-scaled evolutionary tree using the Bayesian Markov chain Monte Carlo method estimated that a common ancestor of the RSV-A, RSV-B, and bovine-RSV diverged at around 450 years ago, and RSV-A and RSV-B diverged around 250 years ago. Finally, the RSV-A F gene formed eight genotypes (GA1-GA7 and NA1) over the last 80 years. Phylodynamics of RSV-A F gene, including all genotype strains, increased twice in the 1990s and 2010s, while patterns of each RSV-A genotype were different. Phylogenetic distance analysis suggested that the genetic distances of the strains were relatively short (less than 0.05). No positive selection sites were estimated, while many negative selection sites were found. Moreover, the F protein 3D structure mapping and conformational epitope analysis implied that the conformational epitopes did not correspond to the neutralizing antibody binding sites of the F protein. These results suggested that the RSV-A F gene is relatively conserved, and mismatches between conformational epitopes and neutralizing antibody binding sites of the F protein are responsible for the virus reinfection.
