Identification of group B respiratory syncytial viruses that lack the 60-nucleotide duplication after six consecutive epidemics of total BA dominance at coastal Kenya

Epidemiology of Respiratory Syncytial Virus Across Five Influenza Seasons Among Adults and Children One Year of Age and Older-Washington State, 2011/2012-2015/2016

BACKGROUND

Vaccines and novel prophylactics against respiratory syncytial virus (RSV) are in development. To provide a baseline for evaluating these interventions, we characterized the incidence and molecular epidemiology of RSV in persons aged ≥1 year.

METHODS

We identified patients with medically attended acute respiratory illness (MAARI) from the 2011/2012 through 2015/2016 influenza seasons among members of Kaiser Permanente Washington. We estimated the cumulative incidence of MAARI for laboratory-confirmed RSV or influenza infection.

RESULTS

Annual cohorts ranged from 82 266 to 162 633 individuals, 14% of whom were children aged 1 to 17 years. Cumulative incidence of RSV each season ranged from 14 per 1000 population (95% confidence interval [CI], 12-16) to 22 per 1000 (95% CI, 19-25). Incidence of RSV was greater than influenza in children aged 12-23 months and 2-4 years; incidence of influenza was greater in other age groups. Respiratory syncytial virus subtype A dominated in 2011/2012, 2012/2013, and 2015/2016, with ON1 being the most common genotype. Respiratory syncytial virus subtype B dominated in 2013/2014 and 2014/2015, primarily of the BA genotype.

CONCLUSIONS

The burden of RSV is comparable to that of influenza across the life course. These results provide a baseline for evaluating the impact of new RSV interventions on the epidemiology of RSV.

Results from the WHO external quality assessment for the respiratory syncytial virus pilot, 2016-17

Background

External quality assessments (EQAs) for the molecular detection of respiratory syncytial virus (RSV) are necessary to ensure the provision of reliable and accurate results. One of the objectives of the pilot of the World Health Organization (WHO) Global RSV Surveillance, 2016‐2017, was to evaluate and standardize RSV molecular tests used by participating countries. This paper describes the first WHO RSV EQA for the molecular detection of RSV.

Methods

The WHO implemented the pilot of Global RSV Surveillance based on the WHO Global Influenza Surveillance and Response System (GISRS) from 2016 to 2018 in 14 countries. To ensure standardization of tests, 13 participating laboratories were required to complete a 12 panel RSV EQA prepared and distributed by the Centers for Disease Control and Prevention (CDC), USA. The 14th laboratory joined the pilot late and participated in a separate EQA. Laboratories evaluated a RSV rRT‐PCR assay developed by CDC and compared where applicable, other Laboratory Developed Tests (LDTs) or commercial assays already in use at their laboratories.

Results

Laboratories performed well using the CDC RSV rRT‐PCR in comparison with LDTs and commercial assays. Using the CDC assay, 11 of 13 laboratories reported correct results. Two laboratories each reported one false‐positive finding. Of the laboratories using LDTs or commercial assays, results as assessed by Ct values were 100% correct for 1/5 (20%). With corrective actions, all laboratories achieved satisfactory outputs.

Conclusions

These findings indicate that reliable results can be expected from this pilot. Continued participation in EQAs for the molecular detection of RSV is recommended.

Optimization of parameters for expression and purification of G glycoprotein ectodomain of respiratory syncytial virus

Aim: G glycoprotein ectodomain (Ge) of BA genotype of group B respiratory syncytial virus was expressed and purified to achieve maximum yield of the protein. Materials & methods: We optimized different parameters like strains, temperature, inducer concentration and post induction time period for efficient protein expression in Escherichia coli. The protein was purified using affinity chromatography and confirmed by western blotting. Results: It was concluded that a 5-h induction with 0.75 mM isopropyl β-D-1-thiogalactopyranoside at 37°C in BL21(DE3) cells was the most favorable condition for maximal protein expression. The far-UV circular dichroism spectroscopy suggested that it is an α-helical protein. Conclusion: The purified Ge protein can be characterized by antigenic and biophysical methods in future studies, which will probably assist in vaccine development.

Genomic analysis of respiratory syncytial virus infections in households and utility in inferring who infects the infant

Infants (under 1-year-old) are at most risk of life threatening respiratory syncytial virus (RSV) disease. RSV epidemiological data alone has been insufficient in defining who acquires infection from whom (WAIFW) within households. We investigated RSV genomic variation within and between infected individuals and assessed its potential utility in tracking transmission in households. Over an entire single RSV season in coastal Kenya, nasal swabs were collected from members of 20 households every 3-4 days regardless of symptom status and screened for RSV nucleic acid. Next generation sequencing was used to generate >90% RSV full-length genomes for 51.1% of positive samples (191/374). Single nucleotide polymorphisms (SNPs) observed during household infection outbreaks ranged from 0-21 (median: 3) while SNPs observed during single-host infection episodes ranged from 0-17 (median: 1). Using the viral genomic data alone there was insufficient resolution to fully reconstruct within-household transmission chains. For households with clear index cases, the most likely source of infant infection was via a toddler (aged 1 to <3 years-old) or school-aged (aged 6 to <12 years-old) co-occupant. However, for best resolution of WAIFW within households, we suggest an integrated analysis of RSV genomic and epidemiological data.

Model-based estimates of transmission of respiratory syncytial virus within households

Introduction

Respiratory syncytial virus (RSV) causes a significant respiratory disease burden in the under 5 population. The transmission pathway to young children is not fully quantified in low-income settings, and this information is required to design interventions.

Methods

We used an individual level transmission model to infer transmission parameters using data collected from 493 individuals distributed across 47 households over a period of 6 months spanning the 2009/2010 RSV season. A total of 208 episodes of RSV were observed from 179 individuals. We model competing transmission risk from within household exposure and community exposure while making a distinction between RSV groups A and B.

Results

We find that 32–53% of all RSV transmissions are between members of the same household; the rate of pair-wise transmission is 58% (95% CrI: 30–74%) lower in larger households (≥8 occupants) than smaller households; symptomatic individuals are 2–7 times more infectious than asymptomatic individuals i.e. 2.48 (95% CrI: 1.22–5.57) among symptomatic individuals with low viral load and 6.7(95% CrI: 2.56–16) among symptomatic individuals with high viral load; previous infection reduces susceptibility to re-infection within the same epidemic by 47% (95% CrI: 17%–68%) for homologous RSV group and 39% (95%CrI: -8%-69%) for heterologous group; RSV B is more frequently introduced into the household, and RSV A is more rapidly transmitted once in the household.

Discussion

Our analysis presents the first transmission modelling of cohort data for RSV and we find that it is important to consider the household social structuring and household size when modelling transmission. The increased infectiousness of symptomatic individuals implies that a vaccine against RSV related disease would also have an impact on infection transmission. Together, the weak cross immunity between RSV groups and the possibility of different transmission niches could form part of the explanation for the group co-existence.

Emergence of BA9 genotype of human respiratory syncytial virus subgroup B in China from 2006 to 2014

A study was conducted to investigate the circulation of HRSV subgroup B (HRSVB) in China in recent years. HRSVB sequences from 365 samples collected in 1991, 2004 and 2008-2014 in China, together with 332 Chinese HRSVB sequences obtained from GenBank were analyzed to determine the geographic and yearly distribution of HRSVB. Phylogenetic analysis revealed these HRSVB sequences clustered into 4 genotypes with different frequencies: BA (83%), CB1 (11%), SAB (3.0%) and GB3 (0.7%). Between 2005 and 2013, there was a co-circulation of BA and non-BA genotypes in China. Genotypes BA9 and BA10 were two of the main BA genotypes detected in this study. Genotype BA9 was first detected in China in 2006 and became the predominant HRSVB genotype circulating in China from 2008 to 2014. Three different lineages were detected for both genotypes BA9 and BA10. Time to the most recent common ancestor for genotypes BA9 and BA10 was estimated for years 1997 and 1996, respectively. Results of this study not only contribute to the understanding of the circulation pattern, but also the phylogenetic pattern and evolution of HRSVB in China from 1991 to 2014.

Biophysical characterization of G protein ectodomain of group B human respiratory syncytial virus from E. coli

Human respiratory syncytial virus (hRSV) is an important pathogen of acute respiratory tract infection. The G protein of hRSV is a transmembrane glycoprotein that is a neutralizing antigen and is thus a vaccine candidate. In this study, synthetic codon optimized ectodomain G protein [G(ΔTM)] of BA genotype of group B hRSV was cloned, expressed, and characterized using biophysical techniques. The molar absorption coefficient and mean residue ellipticity at 222 nm ([θ]222) of G (ΔTM) was found to be 7950 M(-1) cm(-1) and -19701.7 deg cm(2) dmol(-1) respectively. It was concluded that G(ΔTM) mainly consist of α-helix (74.9%) with some amount of β-sheet (4%). The protein was stable up to 85°C without any transition curve. However, heat-induced denaturation of G(ΔTM) resulted in total loss of β-sheet whereas not much change was observed in the α-helix part of the secondary structure. It was concluded that G(ΔTM) is an α-helical protein and it is highly stable at high temperature, but could be easily denatured using high concentrations of GdmCl/urea or acidic condition. This is the first investigation of cloning, expression, and characterization of G(ΔTM) of BA viruses from India. Structural characterization of G protein will assist in drug designing and vaccine development for hRSV.

Molecular detection and characterization of respiratory syncytial virus B genotypes circulating in Pakistani children

Respiratory syncytial virus (RSV) is the major cause of acute lower respiratory tract infections in young children, but very little is known about its epidemiology and circulating genotypes in Pakistan. This study analyzed the epidemiological and molecular characteristics of RSV B genotypes in Pakistani children below 5years with acute respiratory tract infections (ARIs) during three consecutive winter seasons from 2010 to 2013. A total of 1941 samples were analyzed for RSV infection by real time PCR and 24% (472/1941) samples were found positive out of which 22.3% (105/472) were sub-typed as RSV-B. The frequency of outpatient cases was higher (62.5%; 295/472) as compared to hospitalized patients (37.5%; 177/472). Patient ages ranged from 2month to 5years with a mean age of 1.48±1.2 (years) and a median age of 1year. Children below one year made up the highest percentage of enrolled subjects and male to female ratio of RSVB positive cases was nearly equivalent (1:1.1). The most common clinical symptoms were cough (96%), fever (80%) and sore throat (50%). All Pak RSVB strains ascribed to the BA genotype showing 91.9-97.1% and 86.2-95.3% homology at the nucleotide and amino acid levels respectively in comparison to BA prototype strain. On phylogenetic analysis, three genotypes of Pakistan RSV B viruses were observed; BA-9 and BA-10 which have been reported previously from other regions, and a third novel genotype assigned as BA-13 which formed a distinct cluster with protein length of 319 AA and showed 9-11 unique AA substitutions. All the RSV B isolates had two potential N-glycosylation sites in HVR2 of G protein and with heavy O-glycosylation of serine and threonine residues (G scores of 0.5-0.7). This study highlights the diversity of RSVB viruses and the significance of RSV as a dominant viral etiologic agent of pediatric ARI. It also emphasizes the need for continued molecular surveillance for early detection of prevalent and newly emerging genotypes to understand epidemiology of RSV infections in various regions of Pakistan.

Co-Circulation of 72bp Duplication Group A and 60bp Duplication Group B Respiratory Syncytial Virus (RSV) Strains in Riyadh, Saudi Arabia during 2014

Respiratory syncytial virus (RSV) is an important viral pathogen of acute respiratory tract infection (ARI). Limited data are available on molecular epidemiology of RSV from Saudi Arabia. A total of 130 nasopharyngeal aspirates were collected from children less than 5 years of age with ARI symptoms attending the Emergency Department at King Khalid University Hospital and King Fahad Medical City, Riyadh, Saudi Arabia between October and December, 2014. RSV was identified in the 26% of the hospitalized children by reverse transcriptase PCR. Group A RSV (77%) predominated during the study as compared to group B RSV (23%). The phylogenetic analysis of 28 study strains clustered group A RSV in NA1 and ON1 genotypes and group B viruses in BA (BA9) genotype. Interestingly, 26% of the positive samples clustered in genotypes with duplication in the G protein gene (ON1 for group A and BA for group B). Both the genotypes showed enhanced O-linked glycosylation in the duplicated region, with 10 and 2 additional sites in ON1 and BA respectively. Selection pressure analysis revealed purifying selection in both the ON1 and BA genotypes. One codon each in the ON1 (position 274) and BA genotypes (position 219) were positively selected and had high entropy values indicating variations at these amino acid positions. This is the first report describing the presence of ON1 genotype and the first report on co-circulation of two different genotypes of RSV with duplication in the G protein gene from Saudi Arabia. The clinical implications of the simultaneous occurrence of genotypes with duplication in G protein gene in a given population especially in the concurrent infections should be investigated in future. Further, the ongoing surveillance of RSV in this region will reveal the evolutionary trajectory of these two genotypes with duplication in G protein gene from largest country in the Middle East.

Molecular Characterization of Human Respiratory Syncytial Virus in the Philippines, 2012-2013

Human respiratory syncytial virus (HRSV) is a major cause of acute lower respiratory tract infections in infants and children worldwide. We performed molecular analysis of HRSV among infants and children with clinical diagnosis of severe pneumonia in four study sites in the Philippines, including Biliran, Leyte, Palawan, and Metro Manila from June 2012 to July 2013. Nasopharyngeal swabs were collected and screened for HRSV using real-time polymerase chain reaction (PCR). Positive samples were tested by conventional PCR and sequenced for the second hypervariable region (2^nd HVR) of the G gene. Among a total of 1,505 samples, 423 samples were positive for HRSV (28.1%), of which 305 (72.1%) and 118 (27.9%) were identified as HRSV-A and HRSV-B, respectively. Two genotypes of HRSV-A, NA1 and ON1, were identified during the study period. The novel ON1 genotype with a 72-nucleotide duplication in 2^nd HVR of the G gene increased rapidly and finally became the predominant genotype in 2013 with an evolutionary rate higher than the NA1 genotype. Moreover, in the ON1 genotype, we found positive selection at amino acid position 274 (p<0.05) and massive O- and N-glycosylation in the 2^nd HVR of the G gene. Among HRSV-B, BA9 was the predominant genotype circulating in the Philippines. However, two sporadic cases of GB2 genotype were found, which might share a common ancestor with other Asian strains. These findings suggest that HRSV is an important cause of severe acute respiratory infection among children in the Philippines and revealed the emergence and subsequent predominance of the ON1 genotype and the sporadic detection of the GB2 genotype. Both genotypes were detected for the first time in the Philippines.

Retrospective phylogenetic analysis of circulating BA genotype of human respiratory syncytial virus with 60 bp duplication from New Delhi, India during 2007-2010

Human respiratory syncytial virus (hRSV) is the most common viral pathogen of acute lower respiratory tract infection in infants and young children. The G protein of hRSV is the trans-membrane glycoprotein that is involved in the attachment of virion with the host cell. The nasopharyngeal aspirates were subjected to RT-PCR for the second hypervariable region of the G protein gene in the present investigation. Sequencing and phylogenetic analysis revealed that all the study strains clustered within the BA genotype. The study sequences further clustered in BA-9, BA-7, BA-10 and BA-12 subgroups within the BA genotype. The G proteins of the study sequences were predicted to encode 312 and 319 amino acids. Three different N-linked glycosylation sites were observed in the deduced 93-100 amino acid region. There were 40-43 serine and threonine residues that are the potential O-linked glycosylation sites. The non-synonymous/synonymous (dN/dS) ratio was less than one indicating negative selection pressure for amino acid change in the analyzed region of the G protein. The present investigation provides information on circulating strains of BA genotype from New Delhi, India. Further elaborate investigations of the BA viruses from different regions of the world will establish the basis of the rapid global spread and evolutionary pattern of this expanding genotype.

Conservation of G-Protein Epitopes in Respiratory Syncytial Virus (Group A) Despite Broad Genetic Diversity: Is Antibody Selection Involved in Virus Evolution?

Worldwide G-glycoprotein phylogeny of human respiratory syncytial virus (hRSV) group A sequences revealed diversification in major clades and genotypes over more than 50 years of recorded history. Multiple genotypes cocirculated during prolonged periods of time, but recent dominance of the GA2 genotype was noticed in several studies, and it is highlighted here with sequences from viruses circulating recently in Spain and Panama. Reactivity of group A viruses with monoclonal antibodies (MAbs) that recognize strain-variable epitopes of the G glycoprotein failed to correlate genotype diversification with antibody reactivity. Additionally, no clear correlation was found between changes in strain-variable epitopes and predicted sites of positive selection, despite both traits being associated with the C-terminal third of the G glycoprotein. Hence, our data do not lend support to the proposed antibody-driven selection of variants as a major determinant of hRSV evolution. Other alternative mechanisms are considered to account for the high degree of hRSV G-protein variability.

IMPORTANCE

An unusual characteristic of the G glycoprotein of human respiratory syncytial virus (hRSV) is the accumulation of nonsynonymous (N) changes at higher rates than synonymous (S) changes, reaching dN/dS values at certain sites predictive of positive selection. Since these sites cluster preferentially in the C-terminal third of the G protein, like certain epitopes recognized by murine antibodies, it was proposed that immune (antibody) selection might be driving the apparent positive selection, analogous to the antigenic drift observed in the influenza virus hemagglutinin (HA). However, careful antigenic and genetic comparison of the G glycoprotein does not provide evidence of antigenic drift in the G molecule, in agreement with recently published data which did not indicate antigenic drift in the G protein with human sera. Alternative explanations to the immune-driven selection hypothesis are offered to account for the high level of G-protein genetic diversity highlighted in this study.

Local evolutionary patterns of human respiratory syncytial virus derived from whole-genome sequencing

Human respiratory syncytial virus (RSV) is associated with severe childhood respiratory infections. A clear description of local RSV molecular epidemiology, evolution, and transmission requires detailed sequence data and can inform new strategies for virus control and vaccine development. We have generated 27 complete or nearly complete genomes of RSV from hospitalized children attending a rural coastal district hospital in Kilifi, Kenya, over a 10-year period using a novel full-genome deep-sequencing process. Phylogenetic analysis of the new genomes demonstrated the existence and cocirculation of multiple genotypes in both RSV A and B groups in Kilifi. Comparison of local versus global strains demonstrated that most RSV A variants observed locally in Kilifi were also seen in other parts of the world, while the Kilifi RSV B genomes encoded a high degree of variation that was not observed in other parts of the world. The nucleotide substitution rates for the individual open reading frames (ORFs) were highest in the regions encoding the attachment (G) glycoprotein and the NS2 protein. The analysis of RSV full genomes, compared to subgenomic regions, provided more precise estimates of the RSV sequence changes and revealed important patterns of RSV genomic variation and global movement. The novel sequencing method and the new RSV genomic sequences reported here expand our knowledge base for large-scale RSV epidemiological and transmission studies.

IMPORTANCE The new RSV genomic sequences and the novel sequencing method reported here provide important data for understanding RSV transmission and vaccine development. Given the complex interplay between RSV A and RSV B infections, the existence of local RSV B evolution is an important factor in vaccine deployment.

Rapid spread and diversification of respiratory syncytial virus genotype ON1, Kenya

Surveillance of this new genotype helps clarify the mechanisms of rapid emergence of respiratory viruses.

Respiratory syncytial virus genotype ON1, which is characterized by a 72-nt duplication in the attachment protein gene, has been detected in >10 countries since first identified in Ontario, Canada, in 2010. We describe 2 waves of genotype ON1 infections among children admitted to a rural hospital in Kenya during 2012. Phylogenetic analysis of attachment protein gene sequences showed multiple introductions of genotype ON1; variants distinct from the original Canadian viruses predominated in both infection waves. The genotype ON1 dominated over the other group A genotypes during the second wave, and some first wave ON1 variants reappeared in the second wave. An analysis of global genotype ON1 sequences determined that this genotype has become considerably diversified and has acquired signature coding mutations within immunogenic regions, and its most recent common ancestor dates to ≈2008–2009. Surveillance of genotype ON1 contributes to an understanding of the mechanisms of rapid emergence of respiratory viruses.

Molecular characterization of human respiratory syncytial virus subtype B: a novel genotype of subtype B circulating in China

Human respiratory syncytial virus (HRSV) is major pathogen of lower respiratory tract infections in infants and young children worldwide. There have been many studies regarding HRSV subgroup A (HRSV-A) G protein genetic variability but little information about HRSV subtype B (HRSV-B) G protein genetic diversity and molecular evolution in China. Thus, a survey of the molecular epidemiology and evolution of the G protein in China is of high importance. In this study, the circulation and genetic diversity of HRSV in Chongqing, Southwestern China, from June 2009 to May 2013, were investigated. A total of 3,167 nasopharyngeal aspirates were obtained in this study, and it was found that HRSV-B predominated in the 2009-2010 and 2012-2013 epidemic seasons. This study identified the genetic variability of the glycoprotein G gene among 102 HRSV-B strains isolated by cell culture from Chongqing nasopharyngeal aspirates, and 68 Chinese HRSV-B sequences were deposited in GenBank. Genotyping and phylogenetic analysis revealed that the HRSV-B strains were clustered into three genotypes: BA (n = 111, 65.29%), GB3 (n = 5, 2.94%), and a new GB genotype (n = 54, 31.77%) named GB5. The GB5 strains varied from other genotypes in the central conserved region and N-glycosylation sites. The estimated evolutionary rate of Chinese HRSV-B was 2.01 × 10(-3) nucleotide substitutions/site/year, which is similar to the reports from Belgium and the Netherlands with 1.95 × 10(-3) and 2.78 × 10(-3) nucleotide substitutions/site/year, respectively. This study provides data on the circulating pattern and molecular characterization of HRSV-B genotypes in China during four consecutive years and may contribute to HRSV vaccine development.

Examining strain diversity and phylogeography in relation to an unusual epidemic pattern of respiratory syncytial virus (RSV) in a long-term refugee camp in Kenya

Background

A recent longitudinal study in the Dadaab refugee camp near the Kenya-Somalia border identified unusual biannual respiratory syncytial virus (RSV) epidemics. We characterized the genetic variability of the associated RSV strains to determine if viral diversity contributed to this unusual epidemic pattern.

Methods

For 336 RSV positive specimens identified from 2007 through 2011 through facility-based surveillance of respiratory illnesses in the camp, 324 (96.4%) were sub-typed by PCR methods, into 201 (62.0%) group A, 118 (36.4%) group B and 5 (1.5%) group A-B co-infections. Partial sequencing of the G gene (coding for the attachment protein) was completed for 290 (89.5%) specimens. These specimens were phylogenetically analyzed together with 1154 contemporaneous strains from 22 countries.

Results

Of the 6 epidemic peaks recorded in the camp over the period, the first and last were predominantly made up of group B strains, while the 4 in between were largely composed of group A strains in a consecutive series of minor followed by major epidemics. The Dadaab group A strains belonged to either genotype GA2 (180, 98.9%) or GA5 (2, < 1%) while all group B strains (108, 100%) belonged to BA genotype. In sequential epidemics, strains within these genotypes appeared to be of two types: those continuing from the preceding epidemics and those newly introduced. Genotype diversity was similar in minor and major epidemics.

Conclusion

RSV strain diversity in Dadaab was similar to contemporaneous diversity worldwide, suggested both between-epidemic persistence and new introductions, and was unrelated to the unusual epidemic pattern.

Identification of group B respiratory syncytial viruses that lack the 60-nucleotide duplication after six consecutive epidemics of total BA dominance at coastal Kenya

Respiratory syncytial virus BA genotype has reportedly replaced other group B genotypes worldwide. We report the observation of three group B viruses, all identical in G sequence but lacking the BA duplication, at a coastal district hospital in Kenya in early 2012. This follows a period of six consecutive respiratory syncytial virus (RSV) epidemics with 100% BA dominance among group B isolates. The new strains appear only distantly related to BA variants and to previously circulating SAB1 viruses last seen in the district in 2005, suggesting that they were circulating elsewhere undetected. These results are of relevance to an understanding of RSV persistence.