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.

Molecular interactions between respiratory syncytial virus (RSV) fusion protein (F protein) and the cellular receptor Toll-like receptor 4 (TLR4) and myeloid differentiation factor-2 (MD-2) protein complex are unknown. Thus, to reveal the detailed molecular interactions between them, in silico analyses were performed using various bioinformatics techniques. The present simulation data showed that the neutralizing antibody (NT-Ab) binding sites in both prefusion and postfusion proteins at sites II and IV were involved in the interactions between them and the TLR4 molecule. Moreover, the binding affinity between postfusion proteins and the TLR4/MD-2 complex was higher than that between prefusion proteins and the TLR4/MD-2 complex. This increased binding affinity due to conformational changes in the F protein may be able to form syncytium in RSV-infected cells. These results may contribute to better understand the infectivity and pathogenicity (syncytium formation) of RSV.

To understand the evolution of GII.P6-GII.6 and GII.P7-GII.6 strains, the prevalent human norovirus genotypes, we analysed both the RdRp region and VP1 gene in globally collected strains using authentic bioinformatics technologies. A common ancestor of the P6- and P7-type RdRp region emerged approximately 50 years ago and a common ancestor of the P6- and P7-type VP1 gene emerged approximately 110 years ago. Subsequently, the RdRp region and VP1 gene evolved. Moreover, the evolutionary rates were significantly faster for the P6-type RdRp region and VP1 gene than for the P7-type RdRp region and VP1 genes. Large genetic divergence was observed in the P7-type RdRp region and VP1 gene compared with the P6-type RdRp region and VP1 gene. The phylodynamics of the RdRp region and VP1 gene fluctuated after the year 2000. Positive selection sites in VP1 proteins were located in the antigenicity-related protruding 2 domain, and these sites overlapped with conformational epitopes. These results suggest that the GII.6 VP1 gene and VP1 proteins evolved uniquely due to recombination between the P6- and P7-type RdRp regions in the HuNoV GII.P6-GII.6 and GII.P7-GII.6 virus strains.

To understand the evolution of the human norovirus GII.P6-GII.6 and GII.P7-GII.6 strains, we analyzed both the RdRp region and VP1 gene in globally collected strains using authentic bioinformatics technologies. A common ancestor of the P6- and P7-type RdRp region emerged approximately 50 years ago and a common ancestor of the P6- and P7-type VP1 gene emerged approximately 110 years ago. Subsequently, the RdRp region and VP1 gene evolved. Moreover, the evolutionary rates were significantly faster for the P6-type RdRp region and VP1 gene than the P7-type RdRp region and VP1 genes. Large genetic divergence was observed in the P7-type RdRp region and VP1 gene compared with the P6-type RdRp region and VP1 gene. The phylodynamics of the RdRp region and VP1 gene fluctuated after the year 2000. Positive selection sites in VP1 proteins were located in the antigenicity-related protruding 2 domain, and these sites overlapped with conformational epitopes. These results suggest that the GII.6 VP1 gene and VP1 proteins evolved uniquely due to recombination between the P6- and P7-type RdRp regions in the HuNoV GII.P6-GII.6 and GII.P7-GII.6 virus strains.