Persistence and continuous evolution of the human respiratory syncytial virus in northern Taiwan for two decades

Evolutionary trends of respiratory syncytial viruses: Insights from large-scale surveillance and molecular dynamics of G glycoprotein

Human respiratory syncytial virus (RSV) is an underlying cause of lower respiratory illnesses in children, elderly and immunocompromised adults. RSV contains multiple structural and non-structural proteins with two major glycoproteins that control the initial phase of infection, fusion glycoprotein and the attachment (G) glycoprotein. G protein attaches to the ciliated cells of airways initiating the infection. The hypervariable G protein plays a vital role in evolution of RSV strains. We employed multiple bioinformatics tools on systematically accessed large-scale data to evaluate mutations, evolutionary history, and phylodynamics of RSV. Mutational analysis of central conserved region (CCR) on G protein-coding sequences between 163 and 189 positions revealed frequent mutations at site 178 in human RSV (hRSV) A while arginine to glutamine substitutions at site 180 positions in hRSV B, remained prevalent from 2009 to 2014. Phylogenetic analysis indicates multiple signature mutations within G protein responsible for diversification of clades. The USA and China have highest number of surveillance records, followed by Kenya. Markov Chain Monte Carlo Bayesian skyline plot revealed that RSV A evolved steadily from 1990 to 2000, and rapidly between 2003 and 2005. Evolution of RSV B continued from 2003 to 2022, with a high evolution stage from 2016 to 2020. Throughout evolution, cysteine residues maintained their strict conserved states while CCR has an entropy value of 0.0039(±0.0005). This study concludes the notion that RSV G glycoprotein is continuously evolving while the CCR region of G protein maintains its conserved state providing an opportunity for CCR-specific monoclonal antibodys (mAbs) and inhibitors as potential candidates for immunoprophylaxis.

Comparative analysis of human respiratory syncytial virus evolutionary patterns during the COVID-19 pandemic and pre-pandemic periods

The COVID-19 pandemic has resulted in the implementation of strict mitigation measures that have impacted the transmission dynamics of human respiratory syncytial virus (HRSV). The measures also have the potential to influence the evolutionary patterns of the virus. In this study, we conducted a comprehensive analysis comparing genomic variations and evolving characteristics of its neutralizing antigens, specifically F and G proteins, before and during the COVID-19 pandemic. Our findings showed that both HRSV A and B exhibited an overall chronological evolutionary pattern. For the sequences obtained during the pandemic period (2019-2022), we observed that the HRSV A distributed in A23 genotype, but formed into three subclusters; whereas the HRSV B sequences were relatively concentrated within genotype B6. Additionally, multiple positively selected sites were detected on F and G proteins but none were located at neutralizing antigenic sites of the F protein. Notably, amino acids within antigenic site III, IV, and V of F protein remained strictly conserved, while some substitutions occurred over time on antigenic site Ø, I, II and VIII; substitution S389P on antigenic site I of HRSV B occurred during the pandemic period with nearly 50% frequency. However, further analysis revealed no substitutions have altered the structural conformations of the antigenic sites, the vial antigenicity has not been changed. We inferred that the intensive public health interventions during the COVID-19 pandemic did not affect the evolutionary mode of HRSV.

Unprecedented outbreak of respiratory syncytial virus in Taiwan associated with ON1 variant emergence between 2010 and 2020

An outbreak of respiratory syncytial virus (RSV) has been observed in Taiwan since August 2020. We reviewed a central laboratory-based surveillance network established over 20 years by Taiwan Centres for Disease Control for respiratory viral pathogens between 2010 and 2020.A retrospective study of children <5 years old hospitalized with RSV infection at Chang Gung Memorial Hospital between 2018 and 2020 was conducted, and samples positive for RSV-A were sequenced. Clinical data were obtained and stratified by genotype and year.Data from 2020 showed an approximately 4-fold surge in RSV cases compared to 2010 in Taiwan, surpassing previous years during which ON1 was prevalent. Phylogenetic analysis of G protein showed that novel ON1 variants were clustered separately from those of 2018 and 2019 seasons and ON1 reference strains. The variant G protein carried six amino acid changes that emerged gradually in 2019; high consistency was observed in 2020. A unique substitution, E257K, was observed in 2020 exclusively. The F protein of the variant carried T12I and H514N substitutions, which weren't at antigenic sites. In terms of multivariate analysis, age (OR: 0.97; 95% CI: 0.94-0.99; p  = 0.02) and 2020 ON1 variant (OR:2.52; 95% CI:1.13-5.63; p = 0.025) were independently associated with oxygen saturation <94% during hospitalization.The 2020 ON1 variant didn't show higher replication or virulence compared with those in 2018 in our study. The unprecedented 2020 RSV epidemic may attribute to antigenic changes and lack of interferon-stimulated immunity induced by seasonal circulating virus under non-pharmaceutical intervention.

Outbreak of respiratory syncytial virus subtype ON1 among children during COVID-19 pandemic in Southern Taiwan

BACKGROUND

The regional respiratory syncytial virus (RSV) outbreak in southern Taiwan in late 2020 followed the surge of RSV cases in the national surveillance data and displayed distinct clinical features. This study investigated RSV epidemiology in the most recent five years and compared the clinical manifestations of this outbreak with non-outbreak period.

METHODS

Medical records of RSV-infected children at the National Cheng Kung University Hospital from January 2016 to December 2020 were retrospectively retrieved from hospital-based electronic medical database. Cases of RSV infection were identified by RSV antigen positive and/or RSV isolated from respiratory specimens. The demographic, clinical presentations, and laboratory data were recorded. The RSV isolates in 2020 was sequenced for phylogenetic analysis.

RESULTS

Overall, 442 RSV-infected cases were retrieved and 42.1% (186 cases) clustered in late 2020. The 2020 outbreak started in September, peaked in November, and lasted for 3 months. 2020 RSV-infected children were older (2.3 ± 2.2 years vs. 1.0 ± 1.0 years), more likely to be diagnosed with bronchopneumonia (57.5% vs. 31.6%), but also had a lower hospitalization rate, shorter hospital stay, less oxygen use, and less respiratory distress than those in 2016-2019 (all p value < 0.05). The RSV isolates in 2020 belonged to RSV-A subtype ON1 but were phylogenetically distinct from the ON1 strains prevalent in Taiwan previously.

CONCLUSION

The 2020 RSV outbreak was led by the novel RSV-A subtype ON1 variant with clinical manifestations distinct from previous years. Continuous surveillance of new emerging variants of respiratory viruses in the post-pandemic era is warranted.

Global distribution of respiratory syncytial virus A and B infections: a systematic review

Respiratory syncytial virus (RSV) is the main cause of severe respiratory infections in young children. The need for global epidemiologic data regarding RSV has been increasingly recognized. RSV A infections are reported more frequently than RSV B. Nonetheless, the temporal distribution of infections caused by both RSV groups has not been investigated globally. A systematic review was carried out regarding published studies on RSV A and B epidemiology, as well as RSV G gene ectodomain sequence data available at GenBank. A total of 76,668 [45,990 (60%) RSV A and 30,678 (40%) RSV B] positive samples from 83 countries were identified and included in the analysis. Genotype assignment was obtained in 5,340 RSV A and 2,518 RSV B sequences. Two patterns of RSV circulation were observed: continuous seasons with RSV A predominance and alternate predominance of RSV A and B. These patterns were observed in all regions, but the predominant RSV group seldom coincided in all continents during a given year or season. The most frequently identified RSV A genotype was NA1 (including ON1 viruses) (76.30%), and the most frequently identified RSV B genotype was BA (70.65%). Multiple genotypes circulated simultaneously throughout the evolutionary history of RSV, but genotype diversity decreased after the year 2000. The classification of RSV group and genotype is important for the development of vaccines, as well as to understand viral dynamics. This study displays the global and continental RSV circulation patterns from the first report of human RSV infection until the end of 2020.

IRIS: Infection with RespIratory Syncytial Virus in infants-a prospective observational cohort study

BACKGROUND

Respiratory syncytial virus (RSV) is the most common cause of acute lower respiratory tract infection in infants. Globally, RSV is responsible for approximately 3.2 million hospital admissions and about 60,000 in-hospital deaths per year.

METHODS

Infection with RespIratory Syncytial Virus (IRIS) is an observational, multi-centre study enrolling infants with severe RSV infection and healthy controls. Inclusion criteria are age between 0 and 36 months and hospitalisation due to RSV infection at three German sites. Exclusion criteria are premature birth, congenital or acquired bronchopulmonary or cardiac diseases, and immunodeficiency. Healthy control probands are enrolled via recruitment of patients undergoing routine surgical procedures. Blood and respiratory specimens are collected upon admission, and RSV and other pathogens are analysed by multiplex polymerase chain reaction. Different biomaterials, including plasma, nasal lining fluid, blood cells, DNA, and RNA specimens, are sampled in a dedicated biobank. Detailed information on demographic characteristics and medical history is recorded, and comprehensive clinical data, including vital signs, medication, and interventions.

DISCUSSION

The IRIS study aims to discover host and viral factors controlling RSV disease courses in infants. The approach including multi-omics characterisation in clinically well-characterized children with RSV bronchiolitis seeks to improve our understanding of the immune response against this virus. It may disclose novel diagnostic and treatment approaches for respiratory infections in infants.

TRIAL REGISTRATION

ClinicalTrials.gov, NCT04925310. Registered 01 October 2021-Retrospectively registered. https://clinicaltrials.gov/ct2/show/NCT04925310?cond=NCT04925310&draw=2&rank=1.

Spatial-temporal distribution and sequence diversity of group a human respiratory syncytial viruses in Kenya preceding the emergence of ON1 genotype

Background

Human respiratory syncytial virus (HRSV) is a major cause of severe viral acute respiratory illness and contributes significantly to severe pneumonia cases in Africa. Little is known about its spatial–temporal distribution as defined by its genetic diversity.

Methods

A retrospective study conducted utilizing archived nasopharyngeal specimens from patients attending outpatient clinics in hospitals located in five demographically and climatically distinct regions of Kenya; Coast, Western, Highlands, Eastern and Nairobi. The viral total RNA was extracted and tested using multiplex real time RT‐PCR (reverse transcriptase polymerase chain reaction). A segment of the G‐gene was amplified using one‐step RT‐PCR and sequenced by Sanger di‐deoxy method. Bayesian analysis of phylogeny was utilized and subsequently median joining methods for haplotype network reconstruction.

Results

Three genotypes of HRSVA were detected; GA5 (14.0%), GA2 (33.1%), and NA1 (52.9%). HRSVA prevalence varied by location from 33% to 13.2% in the Highlands and the Eastern regions respectively. The mean nucleotide diversity (Pi[π]) varied by genotype: highest of 0.018 for GA5 and lowest of 0.005 for NA1. A total of 58 haplotypes were identified (GA5 10; GA2 20; NA1 28). These haplotypes were introduced into the population locally by single haplotypes and additional subsidiary seeds amongst the GA2 and the NA1 haplotypes.

Conclusions

HRSVA was found across all the regions throughout the study period and comprised three genotypes; GA5, GA2, and NA1 genotypes. The genotypes were disproportionately distributed across the regions with GA5 gradually increasing toward the Western zones and decreasing toward the Eastern zones of the country.

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.

Delayed respiratory syncytial virus outbreak in 2020 in Taiwan was correlated with two novel RSV-A genotype ON1 variants

Background

Human respiratory syncytial virus (RSV) is a leading pathogen of acute respiratory tract disease among infants and young children. Compared with previous seasons, RSV outbreaks in Taiwan during the 2020–2021 season were delayed because of COVID‐19 mitigation measures. We conducted this study to determine the association of viral factors with clinical characteristics of preschool children with RSV infection.

Methods

We performed a molecular epidemiology analysis of RSV among inpatient preschool children in Taiwan. In 80 nasopharyngeal samples positive for RSV, we sequenced and analyzed viral genotypes according to patient data. Patients' clinical data were obtained from medical files, and their clinical profiles were compared with those of RSV cases recorded during the 2014–2017 seasons.

Results

Phylogenetic analysis revealed that among the RSV‐positive samples, all RSV strains identified during the 2020–2021 season belonged to the ON1 genotype. Most of the Taiwan ON1 strains were categorized into two well‐supported clusters with distinct G protein amino acid substitution patterns that had never been demonstrated previously. Furthermore, the proportion of cases among children aged >24 months increased (P < 0.001). Compared with patients infected during the 2014–2017 seasons, patients infected during the 2020–2021 season were hospitalized for shorter days from hospital admission to dereference (P = 0.004) and had a greater need for oxygen supplements (P = 0.021) and systemic steroid therapy (P = 0.026).

Conclusion

The delayed 2020–2021 RSV outbreak in Taiwan was caused by two novel RSV ON1.1 variants. How the change in RSV epidemiology affects future RSV outbreaks warrants exploration.

Evolutionary dynamics of group A and B respiratory syncytial virus in China, 2009-2018

Respiratory syncytial virus (RSV) is a major cause of acute respiratory tract infections in children and is a public health threat globally. To investigate the spatiotemporal dynamics of RSV evolution, we performed systematic phylogenetic analysis using all available sequences from the GenBank database, together with sequences from Shanghai, China. Both RSV-A and RSV-B appear to have originated in North America, with an inferred origin time of 1954.0 (1938.7-1967.6) and 1969.7 (1962.6-1975.5), respectively. BA-like strains of RSV-B, with a 60-nt insertion, and the ON1 strain of RSV-A, with a 72-nt insertion, emerged in 1997.6 (1996.2-1998.6) and 2010.1 (2009.1-2010.3), respectively. Since their origin, both genotypes have gradually replaced the former circulating genotypes to become the dominant strain. The population dynamic of RSV-A showed a seasonal epidemic pattern with obvious expansion in the periods of 2006-2007, 2010-2011, 2011-2012, and 2013-2014. Thirty fixed amino acid substitutions were identified during the divergence of NA4 from GA1 genotypes of RSV-A, and 13 were found during the divergence of SAB4 from GB1 of RSV-B. Importantly, ongoing evolution has occurred since the emergence of ON1, including four amino acid substitutions (I208L, E232G, T253K, and P314L). RSV-A genotypes GA5, NA4, NA1, and ON1 and RSV-B genotypes CB1, SAB4, BA-C, BA10, BA7, and BA9 were co-circulating in China from 2005 to 2015. In particular, RSV-A genotype ON1 was first detected in China in 2011, and it completely replaced GA2 to become the predominant strain after 2016. These data provide important insights into the evolution and epidemiology of RSV.

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.