Genomic evolution of human respiratory syncytial virus during a decade (2013–2023): bridging the path to monoclonal antibody surveillance

Abstract

This study investigated the prevalence, genetic diversity, and evolution of human respiratory syncytial virus (HRSV) in Barcelona from 2013 to 2023. Respiratory specimens from patients with RTI suspicion at Hospital Universitari Vall d'Hebron were collected from October 2013 to May 2023 for laboratory-confirmation of respiratory viruses. Next-generation sequencing was performed in randomly-selected samples with Illumina technology. Phylogenetic analyses of whole genome sequences were performed with BEAST v1.10.4. Signals of selection and evolutionary pressures were inferred by population dynamics and evolutionary analyses. Mutations in major surface proteins were genetic and structurally characterised, emphasizing those within antigenic epitopes. Analyzing 139,625 samples, 5.3% were HRSV-positive (3008 HRSV-A, 3882 HRSV-B, 56 HRSV-A and -B, and 495 unsubtyped HRSV), with a higher prevalence observed in the paediatric population. Pandemic-related shifts in seasonal patterns returned to normal in 2022-2023. A total of 198 whole-genome sequences were obtained for HRSV-A (6.6% of the HRSV-A positive samples) belonging to GA2.3.5 lineage. For HRSV-B, 167 samples were sequenced (4.3% of the HRSV-B positive samples), belonging to GB5.0.2, GB5.0.4a and GB5.0.5a. HRSV-B exhibited a higher evolution rate. Post-SARS-CoV-2 pandemic, both subtypes showed increased evolutionary rates and decreased effective population size initially, followed by a sharp increase. Analyses indicated negative selective pressure on HRSV. Mutations in antigenic epitopes, including S276N and M274I in palivizumab-targeted site II, and I206M, Q209R, and S211N in nirsevimab-targeted site Ø, were identified. Particularly in the context of the large-scale use in 2023-2024 season of nirsevimab, continuous epidemiological and genomic surveillance is crucial.