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

Human respiratory syncytial virus diversity and epidemiology among patients hospitalized with severe respiratory illness in South Africa, 2012-2015

BACKGROUND

We aimed to describe the prevalence of human respiratory syncytial virus (HRSV) and evaluate associations between HRSV subgroups and/or genotypes and epidemiologic characteristics and clinical outcomes in patients hospitalized with severe respiratory illness (SRI).

METHODS

Between January 2012 and December 2015, we enrolled patients of all ages admitted to two South African hospitals with SRI in prospective hospital-based syndromic surveillance. We collected respiratory specimens and clinical and epidemiological data. Unconditional random effect multivariable logistic regression was used to assess factors associated with HRSV infection.

RESULTS

HRSV was detected in 11.2% (772/6908) of enrolled patients of which 47.0% (363/772) were under the age of 6 months. There were no differences in clinical outcomes of HRSV subgroup A-infected patients compared with HRSV subgroup B-infected patients but among patients aged <5 years, children with HRSV subgroup A were more likely be coinfected with Streptococcus pneumoniae (23/208, 11.0% vs. 2/90, 2.0%; adjusted odds ratio 5.7). No significant associations of HRSV A genotypes NA1 and ON1 with specific clinical outcomes were observed.

CONCLUSIONS

While HRSV subgroup and genotype dominance shifted between seasons, we showed similar genotype diversity as noted worldwide. We found no association between clinical outcomes and HRSV subgroups or genotypes.

Genetic Diversity and Molecular Epidemiology of Circulating Respiratory Syncytial Virus in Central Taiwan, 2008-2017

In this study, we investigated the molecular evolution and phylodynamics of respiratory syncytial virus (RSV) over 10 consecutive seasons (2008–2017) and the genetic variability of the RSV genotypes ON1 and BA in central Taiwan. The ectodomain region of the G gene was sequenced for genotyping. The nucleotide and deduced amino acid sequences of the second hypervariable region of the G protein in RSV ON1 and BA were analyzed. A total of 132 RSV-A and 81 RSV-B isolates were obtained. Phylogenetic analysis revealed that the NA1, ON1, and BA9 genotypes were responsible for the RSV epidemics in central Taiwan in the study period. For RSV-A, the NA1 genotype predominated during the 2008–2011 seasons. The ON1 genotype was first detected in 2011 and replaced NA1 after 2012. For RSV-B, the BA9 and BA10 genotypes cocirculated from 2008 to 2010, but the BA9 genotype has predominated since 2012. Amino acid sequence alignments revealed the continuous evolution of the G gene in the ectodomain region. The predicted N-glycosylation sites were relatively conserved in the ON1 (site 237 and 318) and BA9 (site 296 and 310) genotype strains. Our results contribute to the understanding and prediction of the temporal evolution of RSV at the local level.

Determination of genetic characterization and circulation pattern of Respiratory Syncytial Virus (RSV) in children with a respiratory infection, Tehran, Iran, during 2018-2019

The RSV-associated disease accounts for a significant health burden particularly in infants and young children who need to be hospitalized. Since continuous surveillance of circulating RSV genotypes is crucial worldwide, this study aimed to investigate the genetic diversity of RSV circulating strains causing upper or lower acute respiratory infection. Our attention was geared towards studying the cases hospitalized or outpatient in children younger than 2 years of age in Iran during 2018/2019. In this study, nasopharyngeal swabs collected from 206 children who presented with respiratory infection symptoms, were admitted to the referral pediatric ward of Bahrami children's hospital in Tehran, Iran. RSV-positive samples were detected via Nested RT-PCR. The glycoprotein gene was sequenced, and virus genotypes were confirmed through phylogenetic analysis by the MEGA X program. A total of 74 (35.92%) samples tested positive for RSV. Among them, sequencing was done in 10 specimens from 2018 (RSV-A: RSV-B=4:6) and 19 specimens from 2019 (RSV-A: RSV-B=16:3). According to phylogenetic analysis, all RSV-A strains were assigned as ON1 genotype and RSV-B strains were assigned as BA9 genotype. A new N-glycosylation site in Iranian BA9 and positive selection in ON1 genotype was observed. Phylogenetic characterization of strains in the current study revealed co-circulation of ON1 and BA9 as the only prevalent genotypes of both RSV-A and -B groups.

Evolutionary analysis of human respiratory syncytial virus collected in Myanmar between 2015 and 2018

We studied genetic variation in the second hypervariable region (HVR) of the G gene of human respiratory syncytial virus (HRSV) from 1701 nasal swab samples collected from outpatients with acute respiratory infections at two general hospitals in the cities Yangon and Pyinmana in Myanmar from 2015 to 2018. HRSV genotypes were characterized using phylogenetic trees constructed using the maximum likelihood method. Time-scale phylogenetic tree analyses were performed using the Bayesian Markov chain Monte Carlo method. In total, 244 (14.3%) samples were HRSV-positive and were classified as HRSV-A (n = 84, 34.4%), HRSV-B (n = 158, 64.8%), and co-detection of HRSV-A/HRSV-B (n = 2, 0.8%). HRSV epidemics occurred seasonally between July (1.9%, 15/785) and August (10.5%, 108/1028), with peak infections in September (35.8%, 149/416) and October (58.2%, 89/153). HRSV infection rate was higher in children ≥1 year of age than in those <1 year of age (70.5% vs. 29.5%). The most common HRSV symptoms in children were cough (80%-90%) and rhinorrhea (70%-100%). The predominant genotypes were ON1for HRSV-A (78%) and BA9 for HRSV-B (64%). Time to the most recent common ancestor was 2014 (95% highest posterior density [HPD], 2012-2015) for HRSV-A ON1 and 2009 (95% HPD, 2004-2012) for HRSV-B BA9. The mean evolutionary rate (substitutions/site/year) for HRSV-B (2.12 × 10^-2, 95% HPD, 8.53 × 10^-3-3.63 × 10^-2) was slightly higher than that for HRSV-A (1.39 × 10^-2, 95% HPD, 6.03 × 10^-3-2.12 × 10^-2). The estimated effective population size (diversity) for HRSV-A increased from 2015 to 2016 and declined in mid-2018, whereas HRSV-B diversity was constant in 2015 and 2016 and increased in mid-2017. In conclusion, the dominant HRSV-A and HRSV-B genotypes in Myanmar were ON1 and BA9, respectively, between 2015 and 2018. HRSV-B evolved slightly faster than HRSV-A and exhibited unique phylogenetic characteristics.

Analysis of phylogenetic diversity and in vitro adherence characteristics of respiratory syncytial virus and Streptococcus pneumoniae clinical isolates obtained during pediatric respiratory co-infections

Respiratory syncytial virus (RSV) and Streptococcus pneumoniae are frequently co-associated during acute respiratory infections, particularly amongst infants and young children. In this study, we aimed to identify strains of RSV and serotypes/sequence types of S. pneumoniae associated with co-infections within a cohort of paediatric patients, and to assess RSV-mediated adhesion of pneumococcal isolates. The RSV glycoprotein sequence was determined for 58 RSV-positive samples and molecular serotyping and MLST was used to analyse 26 pneumococcal isolates. We also compared 23 pneumococcal isolates for their adherence to RSV-infected or mock-infected airway epithelia cells using immunofluorescence microscopy and automated particle counting. The tight association between RSV and S. pneumoniae was also visualized using scanning electron microscopy. This study did not identify any statistically significant trend in the strains of RSV and S. pneumoniae associated with co-infections. Furthermore, almost all isolates (22 of 23) showed significantly increased adherence to RSV-infected cells. The level of adherence did not appear to correlate with pneumococcal strain or sequence type, and isolates obtained from RSV-infected patients displayed a similar level of adherence as those from RSV-negative patients. The absence of particular S. pneumoniae or RSV strains associated with co-infection, together with the near ubiquitous presence of RSV-mediated adhesion throughout the pneumococcal clinical isolates, may indicate that the mechanisms governing the association with RSV are of sufficient importance to be maintained across much of the species.

Molecular characterization of respiratory syncytial virus circulating in Tunisia between 2015 and 2018

Introduction. Respiratory syncytial virus (RSV) is the most frequently identified viral agent in children with lower respiratory tract infection (LRTI). No data are available to date regarding RSV genotypes circulating in Tunisia.Aim. The aim of the present study was to investigate the genetic variability of the glycoprotein G gene in Tunisian RSV strains.Methodology. Nasopharyngeal aspirates were collected from infants hospitalized for LRTI in five Tunisian hospitals. All specimens were screened for RSV by a direct immunofluorescence assay (DIFA). To molecularly characterize Tunisian RSV strains, a phylogenetic analysis was conducted. Randomly selected positive samples were subjected to reverse transcription PCR amplifying the second hyper-variable region (HVR2) of the G gene.Results. Among a total of 1417 samples collected between 2015 and 2018, 394 (27.8 %) were positive for RSV by DIFA. Analysis of 61 randomly selected RSV strains revealed that group A RSV (78.7 %) predominated during the period of study as compared to group B RSV (21.3 %). The phylogenetic analysis showed that two genotypes of RSV-A were co-circulating: the ON1 genotype with a 72-nt duplication in HVR2 of the G gene was predominant (98.0 % of RSV-A strains), while one RSV-A strain clustered with the NA1 genotype (2.0 %). Concerning Tunisian group B RSV strains, all sequences contained a 60-nt insertion in HVR2 and a clustered BA10 genotype.Conclusion. These data suggest that RSV-A genotype ON1 and RSV-B genotype BA10, both with duplications in the G gene, were widely circulating in the Central coastal region of Tunisia between 2015 and 2018.

Epidemiology and Phylogenetic Analysis of Viral Respiratory Infections in Vietnam

Acute respiratory infections (ARIs) impose a major public health burden on fragile healthcare systems of developing Southeast Asian countries such as Vietnam. The epidemiology, genetic diversity and transmission patterns of respiratory viral pathogens that circulate in this region are not well characterized. We used RT-PCR to screen for 14 common respiratory viruses in nasal/throat samples from 4326 ARI patients from 5 sites in Vietnam during 2012-2016. 64% of patients tested positive for viruses; 14% tested positive multiple co-infecting viruses. The most frequently detected viruses were Respiratory syncytial virus (RSV, 23%), Human Rhinovirus (HRV, 13%), Influenza A virus (IAV, 11%) and Human Bocavirus (HBoV, 7%). RSV infections peaked in July to October, were relatively more common in children <1 year and in the northernmost hospital. IAV infections peaked in December to February and were relatively more common in patients >5 years in the central region. Coinfection with IAV or RSV was associated with increased disease severity compared with patients only infected with HBoV or HRV. Over a hundred genomes belonging to 13 families and 24 genera were obtained via metagenomic sequencing, including novel viruses and viruses less commonly associated with ARIs. Phylogenetic and phylogeographic analyses further indicated that neighboring countries were the most likely source of many virus lineages causing ARIs in Vietnam and estimated the period that specific lineages have been circulating. Our study illustrates the value of applying the state-of-the-art virus diagnostic methods (multiplex RT-PCR and metagenomic sequencing) and phylodynamic analyses at a national level to generate an integrated picture of viral ARI epidemiology.

Disease severity of respiratory syncytial virus (RSV) infection correlate to a novel set of five amino acid substitutions in the RSV attachment glycoprotein (G) in China

Human respiratory syncytial virus (RSV) is one of the major viruses of acute respiratory tract disease among infants and young children. We performed molecular epidemiology analysis of RSV among inpatient children in Guangzhou, China. Phylogenetic and Bayesian analysis showed that genotype ON1 was the only subgroup A virus in this study. Interestingly, the majority of Guangzhou ON1 strains formed a well-supported cluster, and these strains shared a novel set of five amino acid substitutions that never illustrated before. Furthermore, the degree of disease severity was assessed using a severity scoring system. The patients carrying the novel RSV A strain were associated with milder respiratory symptoms compared to other RSV A positives. In conclusion, a specific set of five amino acid substitutions was found in China and further analysis showed that disease severity was associated with these alterations. These findings will provide valuable information for the pathogenic mechanism and vaccine development of RSV.

The emergence of subgenotype ON-1 of Human orthopneumovirus type A in Riyadh, Saudi Arabia: A new episode of the virus epidemiological dynamic

Lower respiratory tract infections caused by Human orthopneumovirus are still a threat to the pediatric population worldwide. To date, the molecular epidemiology of the virus in Saudi Arabia has not been adequately charted. In this study, a total of 205 nasopharyngeal aspirate samples were collected from hospitalized children with lower respiratory tract symptoms during the winter seasons of 2014/15 and 2015/16. Human orthopneumovirus was detected in 89 (43.4%) samples, of which 56 (27.3%) were positive for type A and 33 (16.1%) were positive for type B viruses. The fragment that spans the two hypervariable regions (HVR1 and HVR2) of the G gene of Human orthopneumovirus A was amplified and sequenced. Sequence and phylogenetic analyses have revealed a genotype shift from NA1 to ON-1, which was prevalent during the winter seasons of 2007/08 and 2008/09. Based on the intergenotypic p-distance values, ON-1 was reclassified as a subgenotype of the most predominant genotype GA2. Three conserved N-glycosylation sites were observed in the HVR2 of Saudi ON-1 strains. The presence of a 23 amino acid duplicated region in ON-1 strains resulted in a higher number of O-glycosylation sites as compared to other genotypes. The data presented in this report outlined the replacement of NA1 and NA2 subgenotypes in Saudi Arabia with ON-1 within 7 to 8 years. The continuous evolution of Human orthopneumovirus through point mutations and nucleotide duplication may explain its ability to cause recurrent infections.

Global distribution of NA1 genotype of respiratory syncytial virus and its evolutionary dynamics assessed from the past 11 years

Respiratory syncytial virus (RSV) is a potent pathogen having global distribution. The main purpose of this study was to gain an insight into distribution pattern of the NA1 genotype of group A RSV across the globe together with its evolutionary dynamics. We focused on the second hypervariable region of the G protein gene and used the same for Phylogenetic, Bayesian and Network analyses. Eighteen percent of the samples collected from 500 symptomatic pediatric patients with acute respiratory tract infection (ARI) were found to be positive for RSV during 2011-15 from New Delhi, India. Of these, group B RSV was predominant and clustered into two different genotypes (BA and SAB4). Similarly, group A viruses clustered into two genotypes (NA1 and ON1). The data set from the group A viruses included 543 sequences from 23 different countries including 67 strains from India. The local evolutionary dynamics suggested consistent virus population of NA1 genotype in India during 2009 to 2014. The molecular clock analysis suggested that most recent common ancestor of group A and NA1 genotype have emerged in during the years 1953 and 2000, respectively. The global evolutionary rates of group A viruses and NA1 genotype were estimated to be 3.49 × 10^-3 (95% HPD, 2.90-4.17 × 10^-3) and 3.56 × 10^-3 (95% HPD, 2.91 × 10^-3-4.18 × 10^-3) substitution/site/year, respectively. Analysis of the NA1 genotype of group A RSV reported during 11 years i.e. from 2004 to 2014 showed its dominance in 21 different countries across the globe reflecting its evolutionary dynamics. The Network analysis showed highly intricate but an inconsistent pattern of haplotypes of NA1 genotype circulating in the world. Present study seems to be first comprehensive attempt on global distribution and evolution of NA1 genotype augmenting the optimism towards the vaccine development.

Molecular epidemiology of human respiratory syncytial virus among children in Japan during three seasons and hospitalization risk of genotype ON1

We investigated the genetic diversity, the circulation patterns, and risk for hospital admission of human respiratory syncytial virus (HRSV) strains in Japan between 2012 through 2015. During the study period, 744 HRSV-positive cases were identified by rapid diagnostic test. Of these, 572 samples were positive by real-time PCR; 400 (69.9%) were HRSV-A, and 172 (30.1%) were HRSV-B. HRSV-A and -B alternated as the dominant strain in the subsequent seasons. Phylogenetic tree analysis of the second hyper-variable region of the G protein classified the HRSV-A specimens into NA1 (n = 242) and ON1 (n = 114) genotypes and the HRSV-B specimens into BA9 (n = 60), and BA10 (n = 27). The ON1 genotype, containing a 72-nucleotide duplication in the G protein’s second hyper-variable region, was first detected in the 2012–2013 season but it predominated and replaced the older NA1 HRSV-A in the 2014–2015 season, which also coincided with a record number of HRSV cases reported to the National Infectious Disease Surveillance in Japan. The risk of hospitalization was 6.9 times higher for the ON1 genotype compared to NA1. In conclusion, our data showed that the emergence and predominance of the relatively new ON1 genotype in Japan was associated with a record high number of cases and increased risk for hospitalization.

Respiratory syncytial virus genotypes NA1, ON1, and BA9 are prevalent in Thailand, 2012-2015

Respiratory syncytial virus (RSV) causes acute lower respiratory tract infection in infants and young children worldwide. To investigate the RSV burden in Thailand over four consecutive years (January 2012 to December 2015), we screened 3,306 samples obtained from children ≤5 years old with acute respiratory tract infection using semi-nested reverse-transcription polymerase chain reaction (RT-PCR). In all, 8.4% (277/3,306) of the specimens tested positive for RSV, most of which appeared in the rainy months of July to November. We then genotyped RSV by sequencing the G glycoprotein gene and performed phylogenetic analysis to determine the RSV antigenic subgroup. The majority (57.4%, 159/277) of the RSV belonged to subgroup A (RSV-A), of which NA1 genotype was the most common in 2012 while ON1 genotype became prevalent the following year. Among samples tested positive for RSV-B subgroup B (RSV-B) (42.6%, 118/277), most were genotype BA9 (92.6%, 87/94) with some BA10 and BA-C. Predicted amino acid sequence from the partial G region showed highly conserved N-linked glycosylation site at residue N237 among all RSV-A ON1 strains (68/68), and at residues N296 (86/87) and N310 (87/87) among RSV-B BA9 strains. Positive selection of key residues combined with notable sequence variations on the G gene contributed to the continued circulation of this rapidly evolving virus.

Molecular evolution of hepatitis B vaccine escape variants in China, during 2000-2016

Hepatitis B vaccine escape variants are the main threat to hepatitis B virus (HBV) infection in vaccination era worldwide. With 215 genotype B HBV and 313 genotype C HBV vaccine escape variants isolated from China during 2000-2016, we reported that genotype B HBV vaccine escape strains diverged in ∼1997 (95% HPD; 1987-2005), while genotype C HBV vaccine escape strains diverged in ∼1976 (95% HPD; 1955-2003). Additionally, the p-distance of genotype C HBV vaccine escape strains was 0.0291±0.0169, which was significantly higher than that in the genotype B HBV (t=131.02, p<0.05). However, genotype B HBV vaccine escape strains evolved more rapidly than genotype C HBV (2.103×10^-3 vs 1.083×10^-3 substitutions/site/year). Bayesian skyline plot analysis showed that the populations of genotype C HBV vaccine escape strains fluctuated more than those in genotype B HBV. Four sites (A5T/S, L21S, T/A126S and T/N131I/A) and 13 sites (N3S, T5A, G10Q/R/E, L21S, T47K/A/V, L98V/P, I/S126N/V/T, Q129H/R/L, T131P/I/N/A, G145A/R, L175S/F, L213I/S, V224A/G) were found to be under positive selection in genotype B and C HBV vaccine escape strains, respectively. More importantly, N3S, L21S, T47K, L98V, I/S126T and L213I mutations were detected in 1 (2.5%), 1 (2.5%), 1 (2.5%), 3 (7.5%), 1 (2.5%), 1 (2.5%) genotype C HBV infected Chinese younger with neonatal HBV vaccination, respectively. Therefore, our results should be valuable in further understanding the molecular evolution of HBV and providing new ideas for the elimination of HBV infection.

Molecular evolution of the fusion protein (F) gene in human respiratory syncytial virus subgroup B

In this study, we examined the molecular evolution of the fusion protein (F) gene in human respiratory syncytial virus subgroup B (HRSV-B). First, we performed time-scale evolution analyses using the Bayesian Markov chain Monte Carlo (MCMC) method. Next, we performed genetic distance, linear B-cell epitope prediction, N-glycosylation, positive/negative selection site, and Bayesian skyline plot analyses. We also constructed a structural model of the F protein and mapped the amino acid substitutions and the predicted B-cell epitopes. The MCMC-constructed phylogenetic tree indicated that the HRSV F gene diverged from the bovine respiratory syncytial virus gene approximately 580years ago and had a relatively low evolutionary rate (7.14×10^-4substitutions/site/year). Furthermore, a common ancestor of HRSV-A and -B diverged approximately 290years ago, while HRSV-B diverged into three clusters for approximately 60years. The genetic similarity of the present strains was very high. Although a maximum of 11 amino acid substitutions were observed in the structural model of the F protein, only one strain possessed an amino acid substitution located within the palivizumab epitope. Four epitopes were predicted, although these did not correspond to the neutralization sites of the F protein including the palivizumab epitope. In addition, five N-glycosylation sites of the present HRSV-B strains were inferred. No positive selection sites were identified; however, many sites were found to be under negative selection. The effective population size of the gene has remained almost constant. On the basis of these results, it can be concluded that the HRSV-B F gene is highly conserved, as is the F protein of HRSV-A. Moreover, our prediction of B-cell epitopes does not show that the palivizumab reaction site may be recognized as an epitope during naturally occurring infections.