Molecular evolution of human respiratory syncytial virus attachment glycoprotein (G) gene of new genotype ON1 and ancestor NA1.

We conducted a comprehensive genetic analysis of the C-terminal 3rd hypervariable region of the attachment glycoprotein (G) gene in human respiratory syncytial virus subgroup A (HRSV-A) genotype ON1 (93 strains) and ancestor NA1 (125 strains). Genotype ON1 contains a unique mutation of a 72 nucleotide tandem repeat insertion (corresponding to 24 amino acids) in the hypervariable region. The Bayesian Markov chain Monte Carlo (MCMC) method was used to conduct phylogenetic analysis and a time scale for evolution. We also calculated pairwise distances (p-distances) and estimated the selective pressure. Phylogenetic analysis showed that the analyzed ON1 and NA1 strains formed 4 lineages. A strain belonging to lineage 4 of ON1 showed wide genetic divergence (p-distance, 0.072), which suggests that it might be a candidate new genotype, namely ON2. The emergence of genotype NA1 was estimated to have occurred in 2000 (95% of highest probability density, HPD; 1997-2002) and that of genotype ON1 in 2005 (95% HPD; 2000-2010) based on the time-scaled phylogenetic tree. The evolutionary rate of genotype ON1 was higher than that of ancestral genotype NA1 (6.03×10(-3) vs. 4.61×10(-3) substitutions/site/year, p<0.05). Some positive and many negative selection sites were found in both ON1 and NA1 strains. The results suggested that the new genotype ON1 is rapidly evolving with antigenic changes, leading to epidemics of HRSV infection in various countries.

Hirano E ，Kobayashi M ，Tsukagoshi H ，Yoshida LM ，Kuroda M ，Noda M ，Ishioka T ，Kozawa K ，Ishii H ，Yoshida A ，Oishi K ，Ryo A ，Kimura H ... -  《-》

Molecular evolution of respiratory syncytial virus subgroup A genotype NA1 and ON1 attachment glycoprotein (G) gene in central Vietnam.

We performed molecular evolutionary analyses of the G gene C-terminal 3rd hypervariable region of RSV-A genotypes NA1 and ON1 strains from the paediatric acute respiratory infection patients in central Vietnam during the 2010-2012 study period. Time-scaled phylogenetic analyses were performed using Bayesian Markov Chain Monte Carlo (MCMC) method, and pairwise distances (p-distances) were calculated. Bayesian Skyline Plot (BSP) was constructed to analyze the time-trend relative genetic diversity of central Vietnam RSV-A strains. We also estimated the N-glycosylation sites within G gene hypervariable region. Amino acid substitutions under positive and negative selection pressure were examined using Conservative Single Likelihood Ancestor Counting (SLAC), Fixed Effects Likelihood (FEL), Internal Fixed Effects Likelihood (IFEL) and Mixed Effects Model for Episodic Diversifying Selection (MEME) models. The majority of central Vietnam ON1 strains detected in 2012 were classified into lineage 1 with few positively selected substitutions. As for the Vietnamese NA1 strains, four lineages were circulating during the study period with a few positive selection sites. Shifting patterns of the predominantly circulating NA1 lineage were observed in each year during the investigation period. Median p-distance of central Vietnam NA1 strains was wider (p-distance=0.028) than that of ON1 (p-distance=0.012). The molecular evolutionary rate of central Vietnam ON1 strains was estimated to be 2.55×10-2 (substitutions/site/year) and was faster than NA1 (7.12×10-3 (substitutions/site/year)). Interestingly, the evolutionary rates of both genotypes ON1 and NA1 strains from central Vietnam were faster than the global strains respectively. Furthermore, the shifts of N-glycosylation pattern within the G gene 3rd hypervariable region of Vietnamese NA1 strains were observed in each year. BSP analysis indicated the rapid growth of RSV-A effective population size in early 2012. These results suggested that the molecular evolution of RSV-A G gene detected in central Vietnam was fast with unique evolutionary dynamics.

Yoshihara K ，Le MN ，Nagasawa K ，Tsukagoshi H ，Nguyen HA ，Toizumi M ，Moriuchi H ，Hashizume M ，Ariyoshi K ，Dang DA ，Kimura H ，Yoshida LM ... -  《-》

Rapid replacement of prevailing genotype of human respiratory syncytial virus by genotype ON1 in Beijing, 2012-2014.

Human respiratory syncytial virus (HRSV) is the most common viral pathogen causing lower respiratory infections in infants and young children worldwide. HRSV ON1 genotype in subgroup A with a characteristic of a 72 nucleotide duplication in the second highly variable region of attachment glycoprotein gene, has been reported in some countries since it was first detected in clinical samples collected in Canada in 2010. In this study, 557 HRSV antigen-positive nasopharyngeal aspirates were randomly selected during 2012/2013 to 2013/2014 HRSV seasons in Beijing for subgroup typing and for ON1 genotype screening by using a PCR based method developed for easily identifying genotype ON1 out of strains of subtype A. It was found that subgroup B was dominant in the 2012/2013 season and sudden shift of subgroup dominance from B to A and rapid replacement of previously prevailing NA1 genotype by ON1 genotype occurred in the 2013/2014 season. Reversible amino acid replacement in the G protein gene was found in a new branch of ON1 genotype. The evolutionary rate of the 351 global ON1 sequences was estimated to 7.34 × 10(-3) nucleotide substitutions per site per year (95% highest probability density intervals, HPD, 5.71 × 10(-3) to 9.04 × 10(-3)), with the time of most recent common ancestor dating back to June 2009.

Cui G ，Zhu R ，Deng J ，Zhao L ，Sun Y ，Wang F ，Qian Y ... -  《-》

Molecular evolution of attachment glycoprotein (G) gene in human respiratory syncytial virus detected in Japan 2008-2011.

We investigated the evolution of the C-terminal 3rd hypervariable region of G gene in the prevalent human respiratory syncytial virus (RSV) subgroups A (RSV-A) and B (RSV-B) in Japan in 2008-2011. Phylogenetic analysis and the evolutionary timescale was obtained by the Bayesian Markov Chain Monte Carlo method. All 38 RSV-A strains detected were classified into genotype NA1 and the 17 RSV-B strains detected belonged to genotypes BA and GB2. NA1 subdivided around 1998 in the present phylogenetic tree. Genotype BA subdivided around 1994. The evolutionary rates for RSV-A and RSV-B were estimated at 3.63×10⁻³ and 4.56×10⁻³ substitutions/site/year, respectively. The mean evolutionary rate of RSV-B was significantly faster than that of RSV-A during all seasons. The pairwise distance was relatively short (less than 0.06). In addition, some unique sites under positive selection were found. The results suggested that this region of the RSV strains rapidly evolved with some unique amino acid substitutions due to positive pressure.

Kushibuchi I ，Kobayashi M ，Kusaka T ，Tsukagoshi H ，Ryo A ，Yoshida A ，Ishii H ，Saraya T ，Kurai D ，Yamamoto N ，Kanou K ，Saitoh M ，Noda M ，Kuroda M ，Morita Y ，Kozawa K ，Oishi K ，Tashiro M ，Kimura H ... -  《-》

Genetic analysis of attachment glycoprotein (G) gene in new genotype ON1 of human respiratory syncytial virus detected in Japan.

We studied the evolution of the G gene in the new genotype ON1 of RSV detected from patients with acute respiratory infection in Japan. Phylogenetic analyses and the evolutionary timescale were obtained by the Bayesian MCMC method. We also analyzed p-distance and positive selection sites. A new genotype ON1 emerged around 2001. The evolution rate was rapid (3.57 × 10(-3) substitutions/site per year). The p-distance was short and no positive selection site was found in the present strains. These results suggested that a new genotype ON1 of RSV-A emerged approximately10 years ago and spread to some countries with a high evolution rate.

Tsukagoshi H ，Yokoi H ，Kobayashi M ，Kushibuchi I ，Okamoto-Nakagawa R ，Yoshida A ，Morita Y ，Noda M ，Yamamoto N ，Sugai K ，Oishi K ，Kozawa K ，Kuroda M ，Shirabe K ，Kimura H ... -  《-》