Background: Respiratory syncytial virus (RSV) is a common cause of upper and lower respiratory tract infections. This study aimed to explore the prevalence of respiratory syncytial virus (RSV) and other respiratory viruses in Bulgaria, characterize the genetic diversity of RSV strains, and perform amino acid sequence analyses of RSV surface and internal proteins.Methods: Clinical and epidemiological data and nasopharyngeal swabs were prospectively collected from patients with acute respiratory infections between October 2020 and May 2023. Real-time PCR for 13 respiratory viruses, whole-genome sequencing, phylogenetic, and amino acid analyses were performed.Results: This study included three epidemic seasons (2020–2021, 2021–2022, and 2022–2023) from week 40 of the previous year to week 20 of the following year. Of the 3,047 patients examined, 1,813 (59.5%) tested positive for at least one viral respiratory pathogen. RSV was the second most detected virus (10.9%) after SARS-CoV-2 (22%). Coinfections between RSV and other respiratory viruses were detected in 68 cases, including 14 with SARS-CoV-2. After two seasons of low circulation, RSV activity increased significantly during the 2022–2023 season. The detection rates of RSV were 3.2, 6.6, and 13.7% in the first, second, and third seasons, respectively. RSV was the most common virus found in children under 5 years old with bronchiolitis (40%) and pneumonia (24.5%). RSV-B drove the 2022–2023 epidemic. Phylogenetic analysis indicated that the sequenced RSV-B strains belonged to the GB5.0.5a and GB5.0.6a genotypes. Amino acid substitutions in the surface and internal proteins, including the F protein antigenic sites were identified compared to the BA prototype strain.Conclusion: This study revealed a strong resurgence of RSV in the autumn of 2022 after the lifting of anti-COVID-19 measures, the leading role of RSV as a causative agent of serious respiratory illnesses in early childhood, and relatively low genetic diversity in circulating RSV strains.

Influenza virological surveillance is an essential tool for early detection of novel genetic variants of epidemiologic and clinical significance. The aim of this study was to determine the antigenic and molecular characteristics of influenza viruses circulating in Bulgaria during the 2015/2016 season. The season was characterized by dominant circulation of A(H1N1)pdm09 viruses, accounting for 66% of detected influenza viruses, followed by B/Victoria-lineage viruses (24%) and A(H3N2) viruses (10%). All sequenced influenza A(H1N1)pdm09, A(H3N2) and B/Victoria-lineage viruses belonged to the 6B.1, 3C.2a and 1A genetic groups, respectively. Amino acid analysis of 57 A(H1N1)pdm09 isolates revealed the presence of 16 changes in hemagglutinin (HA) compared to the vaccine virus, five of which occurred in four antigenic sites, together with 16 changes in neuraminidase (NA) and a number of substitutions in proteins MP, NP, NS and PB2. Despite the many amino acid substitutions, A(H1N1)pdm09 viruses remained antigenically closely related to A/California/7/2009 vaccine virus. Bulgarian A(H3N2) strains (subclade 3C.2a) showed changes at 11 HA positions four of which were located in antigenic sites A and B, together with 6 positions in NA, compared to the subclade 3C.3a vaccine virus. They contained unique HA1 substitutions N171K, S312R and HA2 substitutions I77V and G155E compared to Bulgarian 3C.2a viruses of the previous season. All 20 B/Victoria-lineage viruses sequenced harboured two substitutions in the antigenic 120-loop region of HA, and 5 changes in NA, compared to the B/Brisbane/60/2008 vaccine virus. The results of this study reaffirm the continuous genetic variability of circulating seasonal influenza viruses and the need for continued systematic antigenic and molecular surveillance.

The respiratory syncytial virus (RSV) is a leading cause of acute respiratory illnesses (ARI) in infants and young children. The objectives of this study were to investigate the RSV circulation among children aged <5 years in Bulgaria, to identify the RSV-A and RSV-B genotypes and to perform an amino acid sequence analysis of second hypervariable region (HVR2) of the G gene. During the 2014/15 and 2015/16 winter seasons, nasopharyngeal specimens of 610 children aged <5 years with ARI were tested using Real Time RT-PCR for influenza viruses, RSV, metapneumovirus, parainfluenza viruses, rhinoviruses and adenoviruses. Viral respiratory pathogens were detected in 429 (70%) out of 610 patients examined and RSV was the most frequently identified virus (26%) followed by influenza A(H1N1)pdm09 virus (14%) (p < .05). RSV was the most prevalent pathogen in patients with bronchiolitis (48%) and pneumonia (38%). In the 2014/15 season, RSV-A dominated slightly (53%), while in the next season RSV-B viruses prevailed more strongly (66%). The phylogenetic analysis based on the G gene indicated that all 21 studied RSV-A strains belonged to the ON1 genotype; the vast majority (96%) of the RSV-B strains were classified into BA9 genotype and only one - into BA10 genotype. All Bulgarian RSV-A and RSV-B sequences contained a 72-nt and a 60-nt duplication in the HVR2, respectively. The study showed the leading role of this pathogen as a causative agent of serious respiratory illnesses in early childhood, year-on-year fluctuations in RSV incidence, a shift from RSV-A to RSV-B subgroup dominance and relatively low genetic divergence in the circulating strains.

Influenza viruses have a high potential for genetic changes. The objectives of this study were to analyse influenza virus circulation in Bulgaria during the 2019/2020 season, to perform a phylogenetic and molecular analyses of the haemagglutinin (HA) and neuraminidase (NA) sequences of representative influenza strains, and to identify amino acid substitutions compared to the current vaccine strains. Seasonal influenza viruses A(H3N2), A(H1N1)pdm09 and B/Victoria-lineage were detected using a real-time RT-PCR in 323 (23.3%), 149 (10.7%) and 138 (9.9%) out of 1387 patient samples studied, respectively. The HA genes of A(H3N2) viruses analysed belonged to clades 3C.3a (21 strains) and 3C.2a (5 strains): subclades 3C.2a1b + T131K, 3C.2a1b + T135K-B and 3C.2a1b + T135K-A. The clade 3C.3a and subclade 3C.2a1b viruses carried 5 and 14–17 substitutions in HA, as well as 3 and 9 substitutions in NA, respectively, in comparison with the A/Kansas/14/2017 vaccine virus, including some substitutions in the HA antigenic sites A, B, C and E. All 21 A(H1N1)pdm09 viruses sequenced fell into 6B.1A5A subclade. Amino acid sequence analysis revealed the presence of 7–11 substitutions in HA, compared to the A/Brisbane/02/2018 vaccine virus, three of which occurred in antigenic site Sb, along with 6–9 changes at positions in NA. All 10 B/Victoria-lineage viruses sequenced belonged to clade 1A with a triple deletion in HA1 (genetic group 1A(Δ3)B) and carried 7 and 3 substitutions in HA and NA, respectively, with respect to the B/Colorado/06/2017 vaccine virus. The results of this study confirm the rapid evolution of influenza viruses and the need for continuous antigenic and genetic surveillance.