Studies on the etiology of respiratory infections among children in Qatar and surrounding countries are limited. To describe the prevalence and seasonality of RSV, influenza, and other respiratory pathogens among children in Qatar. We retrospectively collected and analyzed data of 33,404 children (<15 years) presented with influenza-like illness from 2012 to 2017. At least one respiratory pathogen was detected in 26,138 (78%) of patients. Together, human rhinoviruses (HRV), respiratory syncytial virus (RSV), and influenza viruses comprised nearly two-thirds of all cases, affecting 24%, 19.7%, and 18.5%, respectively. A prevalence of 5-10% was recorded for adenovirus, parainfluenza viruses (PIVs), human bocavirus (HboV), and human coronaviruses (HCoVs). Human metapneumovirus (HMPV), enteroviruses, M. pneumonia, and parechovirus had prevalences below 5%. While RSV, influenza, and HMPV exhibited strong seasonal activity in the winter, HRV was active during low RSV and influenza circulation. The burden of RSV exceeds that of influenza among young age groups, whereas influenza correlated positively with age. Further, HRV, adenovirus, influenza, and RSV infection rates varied significantly between male and females. This comprehensive multi-year study provides insights into the etiology of ILI among children in Qatar, which represents the Gulf region. Our results reinforce the significance of active surveillance of respiratory pathogens to improve infection prevention and control strategies, particularly among children.

Influenza surveillance is necessary for detection of emerging variants of epidemiologic and clinical significance. This study describes the epidemiology of influenza types A and B, and molecular characteristics of surface glycoproteins (hemagglutinin [HA] and neuraminidase [NA]) of influenza A subtypes: pH1N1 and H3N2 circulated in Arabian Gulf, Levant and North Africa regions during 2009–2017. Analysis of phylogenetics and evolution of HA and NA genes was done using full HA and NA sequences (n = 1229) downloaded from Influenza Research Database (IRD). In total, 130,354 influenza positive cases were reported to WHO during study period. Of these, 50.8% were pH1N1 positive, 15.9% were H3N2 positives and 17.2% were influenza B positive. With few exceptions, all three regions were showing the typical seasonal influenza peak similar to that reported in Northern hemisphere (December–March). However, influenza activity started earlier (October) in both Gulf and North Africa while commenced later during November in Levant countries. The molecular analysis of the HA genes (influenza A subtypes) revealed similar mutations to those reported worldwide. Generally, amino acid substitutions were most frequently found in head domain in H1N1 pandemic viruses, while localized mainly in the stem region in H3N2 viruses. Expectedly, seasons with high pH1N1 influenza activity was associated with a relatively higher number of substitutions in the head domain of the HA in pH1N1 subtype. Furthermore, nucleotide variations were lower at the antigenic sites of pH1N1 viruses compared to H3N2 viruses, which experienced higher variability at the antigenic sites, reflecting the increased immunological pressure because of longer circulation and continuous vaccine changes. Analysis of NA gene of pH1N1 viruses revealed sporadic detections of oseltamivir-resistance mutation, H275Y, in 4% of reported sequences, however, none of NAI resistance mutations were found in the NA of H3N2 viruses. Molecular characterization of H1N1 and H3N2 viruses over 9 years revealed significant differences with regard to position and function of characterized substitutions. While pH1N1 virus substitutions were mainly found in HA head domain, H3N2 virus substitutions were mostly found in HA stem domain. Additionally, more fixed substitutions were encountered in H3N2 virus compared to larger number of non-fixed substitutions in pH1N1.

The diversity of RNA viruses dictates their evolution in a particular host, community or environment. Here, we reported within- and between-host pH1N1virus diversity at consensus and sub-consensus levels over a three-year period (2015–2017) and its implications on disease severity. A total of 90 nasal samples positive for the pH1N1 virus were deep-sequenced and analyzed to detect low-frequency variants (LFVs) and haplotypes. Parallel evolution of LFVs was seen in the hemagglutinin (HA) gene across three scales: among patients (33%), across years (22%), and at global scale. Remarkably, investigating the emergence of LFVs at the consensus level demonstrated that within-host virus evolution recapitulates evolutionary dynamics seen at the global scale. Analysis of virus diversity at the HA haplotype level revealed the clustering of low-frequency haplotypes from early 2015 with dominant strains of 2016, indicating rapid haplotype evolution. Haplotype sharing was also noticed in all years, strongly suggesting haplotype transmission among patients infected during a specific influenza season. Finally, more than half of patients with severe symptoms harbored a larger number of haplotypes, mostly in patients under the age of five. Therefore, patient age, haplotype diversity, and the presence of certain LFVs should be considered when interpreting illness severity. In addition to its importance in understanding virus evolution, sub-consensus virus diversity together with whole genome sequencing is essential to explain variabilities in clinical outcomes that cannot be explained by either analysis alone.