Surveillance of respiratory viruses in the outpatient setting in rural coastal Kenya: baseline epidemiological observations

Local patterns of spread of influenza A H3N2 virus in coastal Kenya over a 1-year period revealed through virus sequence data

The patterns of spread of influenza A viruses in local populations in tropical and sub-tropical regions are unclear due to sparsity of representative spatiotemporal sequence data. We sequenced and analyzed 58 influenza A(H3N2) virus genomes sampled between December 2015 and December 2016 from nine health facilities within the Kilifi Health and Demographic Surveillance System (KHDSS), a predominantly rural region, covering approximately 891 km2 along the Kenyan coastline. The genomes were compared with 1571 contemporaneous global sequences from 75 countries. We observed at least five independent introductions of A(H3N2) viruses into the region during the one-year period, with the importations originating from Africa, Europe, and North America. We also inferred 23 virus location transition events between the nine facilities included in the study. International virus imports into the study area were captured at the facilities of Chasimba, Matsangoni, Mtondia, and Mavueni, while all four exports from the region were captured from the Chasimba facility, all occurring to Africa destinations. A strong spatial clustering of virus strains at all locations was observed associated with local evolution. Our study shows that influenza A(H3N2) virus epidemics in local populations appear to be characterized by limited introductions followed by significant local spread and evolution. Knowledge of the viral lineages that circulate within specific populations in understudied tropical and subtropical regions is required to understand the full diversity and global ecology of influenza viruses and to inform vaccination strategies within these populations.

Global Epidemiology and Seasonality of Human Seasonal Coronaviruses: A Systematic Review

Background

We characterized the global epidemiology and seasonality of human coronaviruses (HCoVs) OC43, NL63, 229E, and HKU1.

Methods

In this systematic review, we searched MEDLINE, EMBASE, Web of Science, SCOPUS, CINAHL, and backward citations for studies published until 1 September 2023. We included studies with ≥12 months of consecutive data and tested for ≥1 HCoV species. Case reports, review articles, animal studies, studies focusing on SARS-CoV-1, SARS-CoV-2, and/or Middle East respiratory syndrome, and those including <100 cases were excluded. Study quality and risk of bias were assessed using Joanna Briggs Institute Critical Appraisal Checklist tools. We reported the prevalence of all HCoVs and individual species. Seasonality was reported for studies that included ≥100 HCoVs annually. This study is registered with PROSPERO, CRD42022330902.

Results

A total of 201 studies (1 819 320 samples) from 68 countries were included. A high proportion were from China (19.4%; n = 39), whereas the Southern Hemisphere was underrepresented. Most were case series (77.1%, n = 155) with samples from secondary care (74.1%, n = 149). Seventeen (8.5%) studies included asymptomatic controls, whereas 76 (37.8%) reported results for all 4 HCoV species. Overall, OC43 was the most prevalent HCoV. Median test positivity of OC43 and NL63 was higher in children, and 229E and HKU1 in adults. Among 18 studies that described seasonality (17 from the Northern Hemisphere), circulation of all HCoVs mostly peaked during cold months.

Conclusions

In our comprehensive review, few studies reported the prevalence of individual HCoVs or seasonality. Further research on the burden and circulation of HCoVs is needed, particularly from Africa, South Asia, and Central/South America.

Rhinovirus dynamics across different social structures

Rhinoviruses (RV), common human respiratory viruses, exhibit significant antigenic diversity, yet their dynamics across distinct social structures remain poorly understood. Our study delves into RV dynamics within Kenya by analysing VP4/2 sequences across four different social structures: households, a public primary school, outpatient clinics in the Kilifi Health and Demographics Surveillance System (HDSS), and countrywide hospital admissions and outpatients. The study revealed the greatest diversity of RV infections at the countrywide level (114 types), followed by the Kilifi HDSS (78 types), the school (47 types), and households (40 types), cumulatively representing >90% of all known RV types. Notably, RV diversity correlated directly with the size of the population under observation, and several RV type variants occasionally fuelled RV infection waves. Our findings highlight the critical role of social structures in shaping RV dynamics, information that can be leveraged to enhance public health strategies. Future research should incorporate whole-genome analysis to understand fine-scale evolution across various social structures.

Prevalence and Molecular Epidemiology of Human Coronaviruses in Africa Prior to the SARS-CoV-2 Outbreak: A Systematic Review

Coronaviruses, re-emerging in human populations, cause mild or severe acute respiratory diseases, and occasionally epidemics. This study systematically reviewed human coronavirus (HCoVs) infections in Africa prior to the SARS-CoV-2 outbreak. Forty studies on the prevalence or molecular epidemiology of HCoVs were available from 13/54 African countries (24%). The first published data on HCoV was from South Africa in 2008. Eight studies (20%) reported on HCoV molecular epidemiology. Endemic HCoV prevalence ranged from 0.0% to 18.2%. The prevalence of zoonotic MERS-CoV ranged from 0.0% to 83.5%. Two studies investigated SARS-CoV infection, for which a prevalence of 0.0% was reported. There was heterogeneity in the type of tests used in determining HCoV prevalence. Two studies reported that risk factors for HCoV include exposure to infected animals or humans. The quantity of virologic investigations on HCoV on the African continent was scant, and Africa was not prepared for SARS-CoV-2.

The burden of RSV-associated illness in children aged < 5 years, South Africa, 2011 to 2016

BACKGROUND

Vaccines and monoclonal antibodies to protect the very young infant against the respiratory syncytial virus (RSV)-associated illness are effective for limited time periods. We aimed to estimate age-specific burden to guide implementation strategies and cost-effectiveness analyses.

METHODS

We combined case-based surveillance and ecological data to generate a national estimate of the burden of RSV-associated acute respiratory illness (ARI) and severe acute respiratory illness (SARI) in South African children aged < 5 years (2011-2016), including adjustment for attributable fraction. We estimated the RSV burden by month of life in the < 1-year age group, by 3-month intervals until 2 years, and then 12 monthly intervals to < 5 years for medically and non-medically attended illness.

RESULTS

We estimated a mean annual total (medically and non-medically attended) of 264,112 (95% confidence interval (CI) 134,357-437,187) cases of RSV-associated ARI and 96,220 (95% CI 66,470-132,844) cases of RSV-associated SARI (4.7% and 1.7% of the population aged < 5 years, respectively). RSV-associated ARI incidence was highest in 2-month-old infants (18,361/100,000 population, 95% CI 9336-28,466). The highest incidence of RSV-associated SARI was in the < 1-month age group 14,674/100,000 (95% CI 11,175-19,645). RSV-associated deaths were highest in the first and second month of life (110.8 (95% CI 74.8-144.5) and 111.3 (86.0-135.8), respectively).

CONCLUSIONS

Due to the high burden of RSV-associated illness, specifically SARI cases in young infants, maternal vaccination and monoclonal antibody products delivered at birth could prevent significant RSV-associated disease burden.

Estimating the cost-effectiveness of maternal vaccination and monoclonal antibodies for respiratory syncytial virus in Kenya and South Africa

BACKGROUND

Respiratory syncytial virus (RSV) causes a substantial burden of acute lower respiratory infection in children under 5 years, particularly in low- and middle-income countries (LMICs). Maternal vaccine (MV) and next-generation monoclonal antibody (mAb) candidates have been shown to reduce RSV disease in infants in phase 3 clinical trials. The cost-effectiveness of these biologics has been estimated using disease burden data from global meta-analyses, but these are sensitive to the detailed age breakdown of paediatric RSV disease, for which there have previously been limited data.

METHODS

We use original hospital-based incidence data from South Africa (ZAF) and Kenya (KEN) collected between 2010 and 2018 of RSV-associated acute respiratory infection (ARI), influenza-like illness (ILI), and severe acute respiratory infection (SARI) as well as deaths with monthly age-stratification, supplemented with data on healthcare-seeking behaviour and costs to the healthcare system and households. We estimated the incremental cost per DALY averted (incremental cost-effectiveness ratio or ICER) of public health interventions by MV or mAb for a plausible range of prices (5-50 USD for MV, 10-125 USD for mAb), using an adjusted version of a previously published health economic model of RSV immunisation.

RESULTS

Our data show higher disease incidence for infants younger than 6 months of age in the case of Kenya and South Africa than suggested by earlier projections from community incidence-based meta-analyses of LMIC data. Since MV and mAb provide protection for these youngest age groups, this leads to a substantially larger reduction of disease burden and, therefore, more favourable cost-effectiveness of both interventions in both countries. Using the latest efficacy data and inferred coverage levels based on antenatal care (ANC-3) coverage (KEN: 61.7%, ZAF: 75.2%), our median estimate of the reduction in RSV-associated deaths in children under 5 years in Kenya is 10.5% (95% CI: 7.9, 13.3) for MV and 13.5% (10.7, 16.4) for mAb, while in South Africa, it is 27.4% (21.6, 32.3) and 37.9% (32.3, 43.0), respectively. Starting from a dose price of 5 USD, in Kenya, net cost (for the healthcare system) per (undiscounted) DALY averted for MV is 179 (126, 267) USD, rising to 1512 (1166, 2070) USD at 30 USD per dose; for mAb, it is 684 (543, 895) USD at 20 USD per dose and 1496 (1203, 1934) USD at 40 USD per dose. In South Africa, a MV at 5 USD per dose would be net cost-saving for the healthcare system and net cost per DALY averted is still below the ZAF's GDP per capita at 40 USD dose price (median: 2350, 95% CI: 1720, 3346). For mAb in ZAF, net cost per DALY averted is 247 (46, 510) USD at 20 USD per dose, rising to 2028 (1565, 2638) USD at 50 USD per dose and to 6481 (5364, 7959) USD at 125 USD per dose.

CONCLUSIONS

Incorporation of new data indicating the disease burden is highly concentrated in the first 6 months of life in two African settings suggests that interventions against RSV disease may be more cost-effective than previously estimated.

Surveillance of respiratory viruses at health facilities from across Kenya, 2014

Background: Acute respiratory illnesses (ARI) are a major cause of morbidity and mortality globally.  With (re)emergence of novel viruses and increased access to childhood bacterial vaccines, viruses have assumed greater importance in the aetiology of ARI. There are now promising candidate vaccines against some of the most common endemic respiratory viruses. Optimal delivery strategies for these vaccines, and the need for interventions against other respiratory viruses, requires geographically diverse data capturing temporal variations in virus circulation. Methods: We leveraged three health facility-based respiratory illness surveillance platforms operating in 11 sites across Kenya. Nasopharyngeal (NP) and/or oropharyngeal (OP) specimens, patient demographic, and clinical characteristics were collected in 2014 from individuals of various ages presenting with respiratory symptoms at the surveillance facilities. Real time multiplex polymerase chain reaction was used to detect rhinoviruses, respiratory syncytial virus (RSV), influenza virus, human coronaviruses (hCoV), and adenoviruses. Results: From 11 sites, 5451 NP/OP specimens were collected and tested from patients. Of these, 40.2% were positive for at least one of the targeted respiratory viruses. The most frequently detected were rhinoviruses (17.0%) and RSV A/B (10.5%), followed by influenza A (6.2%), adenovirus (6.0%) and hCoV (4.2%). RSV was most prevalent among infants aged <12 months old (18.9%), adenovirus among children aged 12–23 months old (11.0%), influenza A among children aged 24–59 months (9.3%), and rhinovirus across all age groups (range, 12.7–19.0%). The overall percent virus positivity varied by surveillance site, health facility type and case definition used in surveillance. Conclusions:  We identify rhinoviruses, RSV, and influenza A as the most prevalent respiratory viruses. Higher RSV positivity in inpatient settings compared to outpatient clinics strengthen the case for RSV vaccination. To inform the design and delivery of public health interventions, long-term surveillance is required to establish regional heterogeneities in respiratory virus circulation and seasonality.

Surveillance of respiratory viruses at health facilities from across Kenya, 2014

Background: Acute respiratory illnesses (ARI) are a major cause of morbidity and mortality globally.  With (re)emergence of novel viruses and increased access to childhood bacterial vaccines, viruses have assumed greater importance in the aetiology of ARI. There are now promising candidate vaccines against some of the most common endemic respiratory viruses. Optimal delivery strategies for these vaccines, and the need for interventions against other respiratory viruses, requires geographically diverse data capturing temporal variations in virus circulation. Methods: We leveraged three health facility-based respiratory illness surveillance platforms operating in 11 sites across Kenya. Nasopharyngeal (NP) and/or oropharyngeal (OP) specimens, patient demographic, and clinical characteristics were collected in 2014 from individuals of various ages presenting with respiratory symptoms at the surveillance facilities. Real time multiplex polymerase chain reaction was used to detect rhinoviruses, respiratory syncytial virus (RSV), influenza virus, human coronaviruses (hCoV), and adenoviruses. Results: From 11 sites, 5451 NP/OP specimens were collected and tested from patients. Of these, 40.2% were positive for at least one of the targeted respiratory viruses. The most frequently detected were rhinoviruses (17.0%) and RSV A/B (10.5%), followed by influenza A (6.2%), adenovirus (6.0%) and hCoV (4.2%). RSV was most prevalent among infants aged <12 months old (18.9%), adenovirus among children aged 12–23 months old (11.0%), influenza A among children aged 24–59 months (9.3%), and rhinovirus across all age groups (range, 12.7–19.0%). The overall percent virus positivity varied by surveillance site, health facility type and case definition used in surveillance. Conclusions:  We identify rhinoviruses, RSV, and influenza A as the most prevalent respiratory viruses. Higher RSV positivity in inpatient settings compared to outpatient clinics strengthen the case for RSV vaccination. To inform the design and delivery of public health interventions, long-term surveillance is required to establish regional heterogeneities in respiratory virus circulation and seasonality.

The national burden of influenza-like illness and severe respiratory illness overall and associated with nine respiratory viruses in South Africa, 2013-2015

BACKGROUND

Estimates of the disease burden associated with different respiratory viruses are severely limited in low- and middle-income countries, especially in Africa.

METHODS

We estimated age-specific numbers and rates of medically and non-medically attended influenza-like illness (ILI) and severe respiratory illness (SRI) that were associated with influenza, respiratory syncytial virus (RSV), rhinovirus, human metapneumovirus, adenovirus, enterovirus and parainfluenza virus types 1-3 after adjusting for the attributable fraction (AF) of virus detection to illness in South Africa during 2013-2015. The base rates were estimated from five surveillance sites and extrapolated nationally.

RESULTS

The mean annual rates per 100,000 population were 51,383 and 4196 for ILI and SRI, respectively. Of these, 26% (for ILI) and 46% (for SRI) were medically attended. Among outpatients with ILI, rhinovirus had the highest AF-adjusted rate (7221), followed by influenza (6443) and adenovirus (1364); whereas, among inpatients with SRI, rhinovirus had the highest AF-adjusted rate (400), followed by RSV (247) and influenza (130). Rhinovirus (9424) and RSV (2026) had the highest AF-adjusted rates among children aged <5 years with ILI or SRI, respectively, whereas rhinovirus (757) and influenza (306) had the highest AF-adjusted rates among individuals aged ≥65 years with ILI or SRI, respectively.

CONCLUSIONS

There was a substantial burden of ILI and SRI in South Africa during 2013-2015. Rhinovirus and influenza had a prominent disease burden among patients with ILI. RSV and influenza were the most prominent causes of SRI in children and the elderly, respectively.

Identification of missed viruses by metagenomic sequencing of clinical respiratory samples from Kenya

Pneumonia remains a major cause of mortality and morbidity. Most molecular diagnoses of viruses rely on polymerase chain reaction (PCR) assays that however can fail due to primer mismatch. We investigated the performance of routine virus diagnostics in Kilifi, Kenya, using random-primed viral next generation sequencing (viral NGS) on respiratory samples which tested negative for the common viral respiratory pathogens by a local standard diagnostic panel. Among 95 hospitalised pneumonia patients and 95 household-cohort individuals, analysis of viral NGS identified at least one respiratory-associated virus in 35 (37%) and 23 (24%) samples, respectively. The majority (66%; 42/64) belonged to the Picornaviridae family. The NGS data analysis identified a number of viruses that were missed by the diagnostic panel (rhinovirus, human metapneumovirus, respiratory syncytial virus and parainfluenza virus), and these failures could be attributed to PCR primer/probe binding site mismatches. Unexpected viruses identified included parvovirus B19, enterovirus D68, coxsackievirus A16 and A24 and rubella virus. The regular application of such viral NGS could help evaluate assay performance, identify molecular causes of missed diagnoses and reveal gaps in the respiratory virus set used for local screening assays. The results can provide actionable information to improve the local pneumonia diagnostics and reveal locally important viral pathogens.

Whole genome sequencing of two human rhinovirus A types (A101 and A15) detected in Kenya, 2016-2018

Background: Virus genome sequencing is increasingly utilized in epidemiological surveillance. Genomic data allows comprehensive evaluation of underlying viral diversity and epidemiology to inform control. For human rhinovirus (HRV), genomic amplification and sequencing is challenging due to numerous types, high genetic diversity and inadequate reference sequences. Methods: We developed a tiled amplicon type-specific protocol for genome amplification and sequencing on the Illumina MiSeq platform of two HRV types, A15 and A101. We then assessed added value in analyzing whole genomes relative to the VP4/2 region only in the investigation of HRV molecular epidemiology within the community in Kilifi, coastal Kenya. Results: We processed 73 nasopharyngeal swabs collected between 2016-2018, and 48 yielded at least 70% HRV genome coverage. These included all A101 samples (n=10) and 38 (60.3%) A15 samples.  Phylogenetic analysis revealed that the Kilifi A101 sequences interspersed with global A101 genomes available in GenBank collected between 1999-2016. On the other hand, our A15 sequences formed a monophyletic group separate from the global genomes collected in 2008 and 2019. An improved phylogenetic resolution was observed with the genome phylogenies compared to the VP4/2 phylogenies. Conclusions: We present a type-specific full genome sequencing approach for obtaining HRV genomic data and characterizing infections.

Whole genome sequencing of two human rhinovirus A types (A101 and A15) detected in Kenya, 2016-2018

Background: Virus genome sequencing is increasingly utilized in epidemiological surveillance. Genomic data allows comprehensive evaluation of underlying viral diversity and epidemiology to inform control. For human rhinovirus (HRV), genomic amplification and sequencing is challenging due to numerous types, high genetic diversity and inadequate reference sequences. Methods: We developed a tiled amplicon type-specific protocol for genome amplification and sequencing on the Illumina MiSeq platform of two HRV types, A15 and A101. We then assessed added value in analyzing whole genomes relative to the VP4/2 region only in the investigation of HRV molecular epidemiology within the community in Kilifi, coastal Kenya. Results: We processed 73 samples collected between 2016-2018, and 48 yielded at least 70% HRV genome coverage. These included all A101 samples (n=10) and 38 (60.3%) A15 samples.  Phylogenetic analysis revealed that the Kilifi A101 sequences interspersed with global A101 genomes available in GenBank collected between 1999-2016. On the other hand, our A15 sequences formed a monophyletic group separate from the global genomes collected in 2008 and 2019. Improved phylogenetic resolution was observed with the genome phylogenies compared to the VP4/2 phylogenies. Conclusions: We present a type-specific full genome sequencing approach for obtaining HRV genomic data and characterizing infections.

Prevalence of human respiratory syncytial virus, parainfluenza and adenoviruses in East Africa Community partner states of Kenya, Tanzania, and Uganda: A systematic review and meta-analysis (2007-2020)

Background

Viruses are responsible for a large proportion of acute respiratory tract infections (ARTIs). Human influenza, parainfluenza, respiratory-syncytial-virus, and adenoviruses are among the leading cause of ARTIs. Epidemiological evidence of those respiratory viruses is limited in the East Africa Community (EAC) region. This review sought to identify the prevalence of respiratory syncytial virus, parainfluenza, and adenoviruses among cases of ARTI in the EAC from 2007 to 2020.

Methods

A literature search was conducted in Medline, Global Index Medicus, and the grey literature from public health institutions and programs in the EAC. Two independent reviewers performed data extraction. We used a random effects model to pool the prevalence estimate across studies. We assessed heterogeneity with the I² statistic, and Cochran’s Q test, and further we did subgroup analysis. This review was registered with PROSPERO under registration number CRD42018110186.

Results

A total of 12 studies met the eligibility criteria for the studies documented from 2007 to 2020. The overall pooled prevalence of adenoviruses was 13% (95% confidence interval [CI]: 6–21, N = 28829), respiratory syncytial virus 11% (95% CI: 7–15, N = 22627), and parainfluenza was 9% (95% CI: 7–11, N = 28363). Pooled prevalence of reported ARTIs, all ages, and locality varied in the included studies. Studies among participants with severe acute respiratory disease had a higher pooled prevalence of all the three viruses. Considerable heterogeneity was noted overall and in subgroup analysis.

Conclusion

Our findings indicate that human adenoviruses, respiratory syncytial virus and parainfluenza virus are prevalent in Kenya, Tanzania, and Uganda. These three respiratory viruses contribute substantially to ARTIs in the EAC, particularly among those with severe disease and those aged five and above.

Qingfei oral liquid inhibited autophagy to alleviate inflammation via mTOR signaling pathway in RSV-infected asthmatic mice

Qingfei oral liquid (QF) is a traditional Chinese medicine that has been used to treat patients with viral pneumonia and asthma for decades. Our previous study revealed that QF prevents airway inflammation and reduces airway hyperresponsiveness (AHR) in respiratory syncytial virus (RSV)-infected asthmatic mice. RSV infection can exacerbate asthma in pediatric patients and induce autophagy, which leads to the promotion of inflammatory cytokine production in the pathology of this disease. The effect of QF on regulating autophagy in RSV-infected asthma patients has not been fully elucidated. In this study, we identified compounds of QF by HPLC-DAD-Q-TOF-MS/MS. The RSV infected OVA challenged mice, we evaluated the RSV-infected asthma model. We found that treatment with QF alleviated airway inflammation and mitigated airway AHR in RSV-infected asthmatic mice. In addition, we found that QF inhibited autophagosome formation and the expression of LC3 protein by using electron and laser confocal microscopy, respectively, to assess RSV-infected asthmatic mice lung tissues. Furthermore, QF was found to reduce the quantity of autophagy and its related proteins LC3B (light chain 3B), Beclin-1, p62 and Atg5 (autophagy-related gene 5) and downstream inflammatory cytokines TNF-α, IL-4, IL-6, and IL-13 via an action in mTOR-dependent signaling in vivo and in vitro. These findings suggest that QF can alleviate the inflammation caused by RSV infection in asthmatic mice, and its mechanism may be involved in the regulation of autophagy via the mTOR signaling pathway.

Evolution of respiratory syncytial virus genotype BA in Kilifi, Kenya, 15 years on

Respiratory syncytial virus (RSV) is recognised as a leading cause of severe acute respiratory disease and deaths among infants and vulnerable adults. Clinical RSV isolates can be divided into several known genotypes. RSV genotype BA, characterised by a 60-nucleotide duplication in the G glycoprotein gene, emerged in 1999 and quickly disseminated globally replacing other RSV group B genotypes. Continual molecular epidemiology is critical to understand the evolutionary processes maintaining the success of the BA viruses. We analysed 735 G gene sequences from samples collected from paediatric patients in Kilifi, Kenya, between 2003 and 2017. The virus population comprised of several genetically distinct variants (n = 56) co-circulating within and between epidemics. In addition, there was consistent seasonal fluctuations in relative genetic diversity. Amino acid changes increasingly accumulated over the surveillance period including two residues (N178S and Q180R) that mapped to monoclonal antibody 2D10 epitopes, as well as addition of putative N-glycosylation sequons. Further, switching and toggling of amino acids within and between epidemics was observed. On a global phylogeny, the BA viruses from different countries form geographically isolated clusters suggesting substantial localized variants. This study offers insights into longitudinal population dynamics of a globally endemic RSV genotype within a discrete location.

Implications of gestational age at antenatal care attendance on the successful implementation of a maternal respiratory syncytial virus (RSV) vaccine program in coastal Kenya

BACKGROUND

Maternal immunisation to boost respiratory syncytial virus (RSV) specific antibodies in pregnant women is a strategy to enhance infant protection. The timing of maternal vaccination during pregnancy may be critical for its effectiveness. However, Kenya has no documented published data on gestational age distribution of pregnant women attending antenatal care (ANC), or the proportion of women attending ANC during the proposed window period for vaccination, to inform appropriate timing for delivery or estimate potential uptake of this vaccine.

METHODS

A cross-sectional survey was conducted within the Kilifi Health and Demographic Surveillance System (KHDSS), coastal Kenya. A simple random sample of 1000 women who had registered pregnant in 2017 to 2018 and with a birth outcome by the time of data collection was taken. The selected women were followed at their homes, and individually written informed consent was obtained. Records of their antenatal attendance during pregnancy were abstracted from their ANC booklet. The proportion of all pregnant women from KHDSS (55%) who attended for one or more ANC in 2018 was used to estimate vaccine coverage.

RESULTS

Of the 1000 women selected, 935 were traced with 607/935 (64.9%) available for interview, among whom 470/607 (77.4%) had antenatal care booklets. The median maternal age during pregnancy was 28.6 years. The median (interquartile range) gestational age in weeks at the first to fifth ANC attendance was 26 (21-28), 29 (26-32), 32 (28-34), 34 (32-36) and 36 (34-38), respectively. The proportion of women attending for ANC during a gestational age window for vaccination of 28-32 weeks (recommended), 26-33 weeks and 24-36 weeks was 76.6% (360/470), 84.5% (397/470) and 96.2% (452/470), respectively. Estimated vaccine coverage was 42.1, 46.5 and 52.9% within the narrow, wide and wider gestational age windows, respectively.

CONCLUSIONS

In a random sample of pregnant women from Kilifi HDSS, Coastal Kenya with card-confirmed ANC clinic attendance, 76.6% would be reached for maternal RSV vaccination within the gestational age window of 28-32 weeks. Widening the vaccination window (26-33 weeks) or (24-36 weeks) would not dramatically increase vaccine coverage and would require consideration of antibody kinetics data that could affect vaccine efficacy.

Surveillance of respiratory viruses among children attending a primary school in rural coastal Kenya

Background: Respiratory viruses are primary agents of respiratory tract diseases. Knowledge on the types and frequency of respiratory viruses affecting school-children is important in determining the role of schools in transmission in the community and identifying targets for interventions. Methods: We conducted a one-year (term-time) surveillance of respiratory viruses in a rural primary school in Kilifi County, coastal Kenya between May 2017 and April 2018. A sample of 60 students with symptoms of ARI were targeted for nasopharyngeal swab (NPS) collection weekly.  Swabs were screened for 15 respiratory virus targets using real time PCR diagnostics. Data from respiratory virus surveillance at the local primary healthcare facility was used for comparison. Results: Overall, 469 students aged 2-19 years were followed up for 220 days. A total of 1726 samples were collected from 325 symptomatic students; median age of 7 years (IQR 5-11). At least one virus target was detected in 384 (22%) of the samples with a frequency of 288 (16.7%) for rhinovirus, 47 (2.7%) parainfluenza virus, 35 (2.0%) coronavirus, 15 (0.9%) adenovirus, 11 (0.6%) respiratory syncytial virus (RSV) and 5 (0.3%) influenza virus.  The proportion of virus positive samples was higher among lower grades compared to upper grades (25.9% vs 17.5% respectively; χ 2 = 17.2, P -value <0.001). Individual virus target frequencies did not differ by age, sex, grade, school term or class size. Rhinovirus was predominant in both the school and outpatient setting. Conclusion: Multiple respiratory viruses circulated in this rural school population.  Rhinovirus was dominant in both the school and outpatient setting and RSV was of notably low frequency in the school. The role of school children in transmitting viruses to the household setting is still unclear and further studies linking molecular data to contact patterns between the school children and their households are required.

Surveillance of endemic human coronaviruses (HCoV-NL63, OC43 and 229E) associated with childhood pneumonia in Kilifi, Kenya

Introduction: Human coronaviruses (HCoVs) circulate endemically in human populations, often with seasonal variation. We describe the long-term patterns of paediatric disease associated with three of these viruses, HCoV-NL63, OC43 and 229E, in coastal Kenya. Methods: Continuous surveillance of pneumonia admissions was conducted at the Kilifi county hospital (KCH) located in the northern coastal region of Kenya. Children aged <5 years admitted to KCH with clinically defined syndromic severe or very severe pneumonia were recruited. Respiratory samples were taken and tested for 15 virus targets, using real-time polymerase chain reaction. Unadjusted odds ratios were used to estimate the association between demographic and clinical characteristics and HCoV positivity. Results: From 2007 to 2019, we observed 11,445 pneumonia admissions, of which 314 (3.9%) tested positive for at least one of the HCoV types surveyed in the study. There were 129 (41.1%) OC43, 99 (31.5%) 229E, 74 (23.6%) NL63 positive cases and 12 (3.8%) cases of HCoV to HCoV coinfection.  Among HCoV positive cases, 47% (n=147) were coinfected with other respiratory virus pathogens. The majority of HCoV cases were among children aged <1 year (66%, n=208), though there was was no change in the proportion infected by age. HCoV-OC43 was predominant of the three HCoV types throughout the surveillance period. Evidence for seasonality was not identified. Conclusions: Overall, 4% of paediatric pneumonia admissions were associated with three endemic HCoVs, with a high proportion of cases co-occurring with another respiratory virus, no clear seasonal pattern, and with the age-distribution of cases following that of pneumonia admissions (i.e. highest in infants). These observations suggest, at most, a small severe disease contribution of endemic HCoVs in this tropical setting and offer insight into their potential future burden and epidemiological characteristics.

Molecular Epidemiology of Human Rhinovirus From 1-Year Surveillance Within a School Setting in Rural Coastal Kenya

Background

Human rhinovirus (HRV) is the most common cause of the common cold but may also lead to more severe respiratory illness in vulnerable populations. The epidemiology and genetic diversity of HRV within a school setting have not been previously described. The objective of this study was to characterize HRV molecular epidemiology in a primary school in a rural location of Kenya.

Methods

Between May 2017 and April 2018, over 3 school terms, we collected 1859 nasopharyngeal swabs (NPS) from pupils and teachers with symptoms of acute respiratory infection in a public primary school in Kilifi County, coastal Kenya. The samples were tested for HRV using real-time reverse transcription polymerase chain reaction. HRV-positive samples were sequenced in the VP4/VP2 coding region for species and genotype classification.

Results

A total of 307 NPS (16.4%) from 164 individuals were HRV positive, and 253 (82.4%) were successfully sequenced. The proportion of HRV in the lower primary classes was higher (19.8%) than upper primary classes (12.2%; P < .001). HRV-A was the most common species (134/253; 53.0%), followed by HRV-C (73/253; 28.9%) and HRV-B (46/253; 18.2%). Phylogenetic analysis identified 47 HRV genotypes. The most common genotypes were A2 and B70. Numerous (up to 22 in 1 school term) genotypes circulated simultaneously, there was no individual re-infection with the same genotype, and no genotype was detected in all 3 school terms.

Conclusions

HRV was frequently detected among school-going children with mild acute respiratory illness symptoms, particularly in the younger age groups (<5-year-olds). Multiple HRV introductions were observed that were characterized by considerable genotype diversity.

Genomic epidemiology and evolutionary dynamics of respiratory syncytial virus group B in Kilifi, Kenya, 2015-17

Respiratory syncytial virus (RSV) circulates worldwide, occurring seasonally in communities, and is a leading cause of acute respiratory illness in young children. There is paucity of genomic data from purposively sampled populations by which to investigate evolutionary dynamics and transmission patterns of RSV. Here we present an analysis of 295 RSV group B (RSVB) genomes from Kilifi, coastal Kenya, sampled from individuals seeking outpatient care in nine health facilities across a defined geographical area (∼890 km²), over two RSV epidemics between 2015 and 2017. RSVB diversity was characterized by multiple virus introductions into the area and co-circulation of distinct genetic clusters, which transmitted and diversified locally with varying frequency. Increase in relative genetic diversity paralleled seasonal virus incidence. Importantly, we identified a cluster of viruses that emerged in the 2016/17 epidemic, carrying distinct amino-acid signatures including a novel nonsynonymous change (K68Q) in antigenic site ∅ in the Fusion protein. RSVB diversity was additionally marked by signature nonsynonymous substitutions that were unique to particular genomic clusters, some under diversifying selection. Our findings provide insights into recent evolutionary and epidemiological behaviors of RSVB, and highlight possible emergence of a novel antigenic variant, which has implications on current prophylactic strategies in development.

Epidemiological and evolutionary dynamics of influenza B virus in coastal Kenya as revealed by genomic analysis of strains sampled over a single season

The genomic epidemiology of influenza B virus (IBV) remains understudied in Africa despite significance to design of effective local and global control strategies. We undertook surveillance throughout 2016 in coastal Kenya, recruiting individuals presenting with acute respiratory illness at nine outpatient health facilities (any age) or admitted to the Kilifi County Hospital (<5 years old). Whole genomes were sequenced for a selected 111 positives; 94 (84.7%) of B/Victoria lineage and 17 (15.3%) of B/Yamagata lineage. Inter-lineage reassortment was detected in ten viruses; nine with B/Yamagata backbone but B/Victoria NA and NP segments and one with a B/Victoria backbone but B/Yamagata PB2, PB1, PA, and MP segments. Five phylogenomic clusters were identified among the sequenced viruses; (i), pure B/Victoria clade 1A (n = 93, 83.8%), (ii), reassortant B/Victoria clade 1A (n = 1, 0.9%), (iii), pure B/Yamagata clade 2 (n = 2, 1.8%), (iv), pure B/Yamagata clade 3 (n = 6, 5.4%), and (v), reassortant B/Yamagata clade 3 (n = 9, 8.1%). Using divergence dates and clustering patterns in the presence of global background sequences, we counted up to twenty-nine independent IBV strain introductions into the study area (∼900 km2) in 2016. Local viruses, including the reassortant B/Yamagata strains, clustered closely with viruses from neighbouring Tanzania and Uganda. Our study demonstrated that genomic analysis provides a clearer picture of locally circulating IBV diversity. The high number of IBV introductions highlights the challenge in controlling local influenza epidemics by targeted approaches, for example, sub-population vaccination or patient quarantine. The finding of divergent IBV strains co-circulating within a single season emphasises why broad immunity vaccines are the most ideal for influenza control in Kenya.

Surveillance of endemic human coronaviruses (HCoV-NL63, OC43 and 229E) associated with pneumonia in Kilifi, Kenya

Introduction: Human coronaviruses (HCoVs) circulate endemically in human populations, often with seasonal variation. We describe the long-term patterns of paediatric disease associated with three of these viruses, HCoV-NL63, OC43 and 229E, in coastal Kenya. Methods: Continuous surveillance of pneumonia admissions was conducted at the Kilifi county hospital (KCH) located in the northern coastal region of Kenya. Children aged <5 years admitted to KCH with clinically defined syndromic severe or very severe pneumonia were recruited. Respiratory samples were taken and tested for 15 virus targets, using real-time polymerase chain reaction. Unadjusted odds ratios were used to estimate the association between demographic and clinical characteristics and HCoV positivity. Results: From 2007 to 2019, we observed 11,445 pneumonia admissions, of which 314 (3.9%) tested positive for at least one HCoV type. There were 129 (41.1%) OC43, 99 (31.5%) 229E, 74 (23.6%) NL63 positive cases and 12 (3.8%) cases of HCoV to HCoV coinfection.  Among HCoV positive cases, 47% (n=147) were coinfected with other respiratory virus pathogens. The majority of HCoV cases were among children aged <1 year (66%, n=208), though there was no age-dependence in the proportion testing positive. HCoV-OC43 was predominant of the three HCoV types throughout the surveillance period. Evidence for seasonality was not identified. Conclusions: Overall, 4% of paediatric pneumonia admissions were associated with three endemic HCoVs, with a high proportion of cases co-occurring with another respiratory virus, with no clear seasonal pattern, and with the age-distribution of cases following that of pneumonia admissions (i.e. highest in infants). These observations suggest, at most, a small severe disease contribution of endemic HCoVs in this tropical setting and offer insight into the potential future burden and epidemiological characteristics of SARS-CoV-2.

Transmission and evolutionary dynamics of human coronavirus OC43 strains in coastal Kenya investigated by partial spike sequence analysis, 2015-16

Human coronavirus OC43 (HCoV-OC43) is a major contributor to seasonal outbreaks of acute respiratory illness (ARI). The origins of locally circulating HCoV-OC43 strains and characteristics of their genetic diversity are unknown for most settings despite significance to effective HCoV control strategies. Between December 2015 and June 2016, we undertook ARI surveillance in coastal Kenya in nine outpatients and one inpatient health facility (HF). Ninety-two patient samples tested HCoV-OC43 positive and forty (43.5%) were successfully sequenced in spike (S) gene region (2,864 long, ∼70%). Phylogenetic analysis confirmed co-circulation of two distinct HCoV-OC43 clades that closely clustered with genotype G (n = 34, 85%) and genotype H (n = 6, 15%) reference strains. Local viruses within the same clade displayed low genetic diversity yielding identical sequences in multiple HF. Furthermore, the newly sequenced Kenyan viruses showed close phylogenetic relationship to other contemporaneous sampled strains (2015-16) including those originating from distant places (e.g. USA and China). Using a genetic similarity threshold of 99.1 per cent at nucleotide level, the HCoV-OC43 strains sampled globally between 1967 and 2019 fell into nine sequence clusters. Notably, some of these clusters appeared to have become extinct, or occurred only sporadically in a few geographical areas while others persisted globally for multiple years. In conclusion, we found that HCoV-OC43 strains spread rapidly both locally and across the globe with limited genetic evolution in the spike gene. Full-genome sequences that are spatio-temporally representative are required to advance understanding of the transmission pathways of this important human respiratory pathogen.

Surveillance of respiratory viruses among children attending a primary school in rural coastal Kenya

Background: Respiratory viruses are primary agents of respiratory tract diseases. Knowledge on the types and frequency of respiratory viruses affecting school-children is important in determining the role of schools in transmission in the community and identifying targets for interventions. Methods: We conducted a one-year (term-time) surveillance of respiratory viruses in a rural primary school in Kilifi County, coastal Kenya between May 2017 and April 2018. A sample of 60 students with symptoms of ARI were targeted for nasopharyngeal swab (NPS) collection weekly.  Swabs were screened for 15 respiratory virus targets using real time PCR diagnostics. Data from respiratory virus surveillance at the local primary healthcare facility was used for comparison. Results: Overall, 469 students aged 2-19 years were followed up for 220 days. A total of 1726 samples were collected from 325 symptomatic students; median age of 7 years (IQR 5-11). At least one virus target was detected in 384 (22%) of the samples with a frequency of 288 (16.7%) for rhinovirus, 47 (2.7%) parainfluenza virus, 35 (2.0%) coronavirus, 15 (0.9%) adenovirus, 11 (0.6%) respiratory syncytial virus (RSV) and 5 (0.3%) influenza virus.  The proportion of virus positive samples was higher among lower grades compared to upper grades (25.9% vs 17.5% respectively; χ 2 = 17.2, P -value <0.001). Individual virus target frequencies did not differ by age, sex, grade, school term or class size. Rhinovirus was predominant in both the school and outpatient setting. Conclusion: Multiple respiratory viruses circulated in this rural school population.  Rhinovirus was dominant in both the school and outpatient setting and RSV was of notably low frequency in the school. The role of school children in transmitting viruses to the household setting is still unclear and further studies linking molecular data to contact patterns between the school children and their households are required.

Genetic characterization of influenza A(H3N2) viruses circulating in coastal Kenya, 2009-2017

Background

Influenza viruses evolve rapidly and undergo immune driven selection, especially in the hemagglutinin (HA) protein. We report amino acid changes affecting antigenic epitopes and receptor‐binding sites of A(H3N2) viruses circulating in Kilifi, Kenya, from 2009 to 2017.

Methods

Next‐generation sequencing (NGS) was used to generate A(H3N2) virus genomic data from influenza‐positive specimens collected from hospital admissions and health facility outpatients presenting with acute respiratory illness to health facilities within the Kilifi Health and Demographic Surveillance System. Full‐length HA sequences were utilized to characterize A(H3N2) virus genetic and antigenic changes.

Results

From 186 (90 inpatient and 96 outpatient) influenza A virus‐positive specimens processed, 101 A(H3N2) virus whole genomes were obtained. Among viruses identified in inpatient specimens from 2009 to 2015, divergence of circulating A(H3N2) viruses from the vaccine strains A/Perth/16/2009, A/Texas/50/2012, and A/Switzerland/9715293/2013 formed 6 genetic clades (A/Victoria/208/2009‐like, 3B, 3C, 3C.2a, 4, and 7). Among viruses identified in outpatient specimens from 2015 to 2017, divergence of circulating A(H3N2) viruses from vaccine strain A/Hong Kong/4801/2014 formed clade 3C.2a, subclades 3C.2a2 and 3C.2a3, and subgroup 3C.2a1b. Several amino acid substitutions were associated with the continued genetic evolution of A(H3N2) strains in circulation.

Conclusions

Our results suggest continuing evolution of currently circulating A(H3N2) viruses in Kilifi, coastal Kenya and suggest the need for continuous genetic and antigenic viral surveillance of circulating seasonal influenza viruses with broad geographic representation to facilitate prompt and efficient selection of influenza strains for inclusion in future influenza vaccines.

Respiratory syncytial, parainfluenza and influenza virus infection in young children with acute lower respiratory infection in rural Gambia

Respiratory viral infections contribute significantly to morbidity and mortality worldwide, but representative data from sub-Saharan Africa are needed to inform vaccination strategies. We conducted population-based surveillance in rural Gambia using standardized criteria to identify and investigate children with acute lower respiratory infection (ALRI). Naso- and oropharyngeal swabs were collected. Each month from February through December 2015, specimens from 50 children aged 2–23 months were randomly selected to test for respiratory syncytial (RSV), parainfluenza (PIV) and influenza viruses. The expected number of viral-associated ALRI cases in the population was estimated using statistical simulation that accounted for the sampling design. RSV G and F proteins and influenza hemagglutinin genes were sequenced. 2385 children with ALRI were enrolled, 519 were randomly selected for viral testing. One or more viruses were detected in 303/519 children (58.4%). RSV-A was detected in 237 and RSV-B in seven. The expected incidence of ALRI associated with RSV, PIV or influenza was 140 cases (95% CI, 131–149) per 1000 person-years; RSV incidence was 112 cases (95% CI, 102–122) per 1000 person-years. Multiple strains of RSV and influenza circulated during the year. RSV circulated throughout most of the year and was associated with eight times the number of ALRI cases compared to PIV or IV. Gambian RSV viruses were closely related to viruses detected in other continents. An effective RSV vaccination strategy could have a major impact on the burden of ALRI in this setting.

Human metapneumovirus prevalence and patterns of subgroup persistence identified through surveillance of pediatric pneumonia hospital admissions in coastal Kenya, 2007-2016

BACKGROUND

Human metapneumovirus (HMPV) is an important respiratory pathogen that causes seasonal epidemics of acute respiratory illness and contributes significantly to childhood pneumonia. Current knowledge and understanding on its patterns of spread, prevalence and persistence in communities in low resource settings is limited.

METHODS

We present findings of a molecular-epidemiological analysis of nasal samples from children < 5 years of age admitted with syndromic pneumonia between 2007 and 2016 to Kilifi County Hospital, coastal Kenya. HMPV infection was detected using real-time RT-PCR and positives sequenced in the fusion (F) and attachment (G) genes followed by phylogenetic analysis. The association between disease severity and HMPV subgroup was assessed using Fisher's exact test.

RESULTS

Over 10 years, 274/6756 (4.1%) samples screened were HMPV positive. Annual prevalence fluctuated between years ranging 1.2 to 8.7% and lowest in the recent years (2014-2016). HMPV detections were most frequent between October of one year to April of the following year. Genotyping was successful for 205/274 (74.8%) positives revealing clades A2b (41.0%) and A2c (10.7%), and subgroups B1 (23.4%) and B2 (24.9%). The dominance patterns were: clade A2b between 2007 and 11, subgroup B1 between 2012 and 14, and clade A2c in more recent epidemics. Subgroup B2 viruses were present in all the years. Temporal phylogenetic clustering within the subgroups for both local and global sequence data was seen. Subgroups occurring in each epidemic season were comprised of multiple variants. Pneumonia severity did not vary by subgroup (p = 0.264). In both the F and G gene, the sequenced regions were found to be predominantly under purifying selection.

CONCLUSION

Subgroup patterns from this rural African setting temporally map with global strain distribution, suggesting a well-mixed global virus transmission pool of HMPV. Persistence in the local community is characterized by repeated introductions of HMPV variants from the global pool. The factors underlying the declining prevalence of HMPV in this population should be investigated.

Human rhinovirus spatial-temporal epidemiology in rural coastal Kenya, 2015-2016, observed through outpatient surveillance

Background: Human rhinovirus (HRV) is the predominant cause of upper respiratory tract infections, resulting in a significant public health burden. The virus circulates as many different types (168), each generating strong homologous, but weak heterotypic, immunity. The influence of these features on transmission patterns of HRV in the community is understudied. Methods: Nasopharyngeal swabs were collected from patients with symptoms of acute respiratory infection (ARI) at nine out-patient facilities across a Health and Demographic Surveillance System between December 2015 and November 2016. HRV was diagnosed by real-time RT-PCR, and the VP4/VP2 genomic region of the positive samples sequenced. Phylogenetic analysis was used to determine the HRV types. Classification models and G-test statistic were used to investigate HRV type spatial distribution. Demographic characteristics and clinical features of ARI were also compared. Results: Of 5,744 NPS samples collected, HRV was detected in 1057 (18.4%), of which 817 (77.3%) were successfully sequenced. HRV species A, B and C were identified in 360 (44.1%), 67 (8.2%) and 390 (47.7%) samples, respectively. In total, 87 types were determined: 39, 10 and 38 occurred within species A, B and C, respectively. HRV types presented heterogeneous temporal patterns of persistence. Spatially, identical types occurred over a wide distance at similar times, but there was statistically significant evidence for clustering of types between health facilities in close proximity or linked by major road networks. Conclusion: This study records a high prevalence of HRV in out-patient presentations exhibiting high type diversity. Patterns of occurrence suggest frequent and independent community invasion of different types. Temporal differences of persistence between types may reflect variation in type-specific population immunity. Spatial patterns suggest either rapid spread or multiple invasions of the same type, but evidence of similar types amongst close health facilities, or along road systems, indicate type partitioning structured by local spread.

Human rhinovirus spatial-temporal epidemiology in rural coastal Kenya, 2015-2016, observed through outpatient surveillance

Background: Human rhinovirus (HRV) is the predominant cause of upper respiratory tract infections, resulting in a significant public health burden. The virus circulates as many different types (~160), each generating strong homologous, but weak heterotypic, immunity. The influence of these features on transmission patterns of HRV in the community is understudied. Methods: Nasopharyngeal swabs were collected from patients with symptoms of acute respiratory infection (ARI) at nine out-patient facilities across a Health and Demographic Surveillance System between December 2015 and November 2016. HRV was diagnosed by real-time RT-PCR, and the VP4/VP2 genomic region of the positive samples sequenced. Phylogenetic analysis was used to determine the HRV types. Classification models and G-test statistic were used to investigate HRV type spatial distribution. Demographic characteristics and clinical features of ARI were also compared. Results: Of 5,744 NPS samples collected, HRV was detected in 1057 (18.4%), of which 817 (77.3%) were successfully sequenced. HRV species A, B and C were identified in 360 (44.1%), 67 (8.2%) and 390 (47.7%) samples, respectively. In total, 87 types were determined: 39, 10 and 38 occurred within species A, B and C, respectively. HRV types presented heterogeneous temporal patterns of persistence. Spatially, identical types occurred over a wide distance at similar times, but there was statistically significant evidence for clustering of types between health facilities in close proximity or linked by major road networks. Conclusion: This study records a high prevalence of HRV in out-patient presentations exhibiting high type diversity. Patterns of occurrence suggest frequent and independent community invasion of different types. Temporal differences of persistence between types may reflect variation in type-specific population immunity. Spatial patterns suggest either rapid spread or multiple invasions of the same type, but evidence of similar types amongst close health facilities, or along road systems, indicate type partitioning structured by local spread.

Whole genome analysis of local Kenyan and global sequences unravels the epidemiological and molecular evolutionary dynamics of RSV genotype ON1 strains

The respiratory syncytial virus (RSV) group A variant with the 72-nucleotide duplication in the G gene, genotype ON1, was first detected in Kilifi in 2012 and has almost completely replaced circulating genotype GA2 strains. This replacement suggests some fitness advantage of ON1 over the GA2 viruses in Kilifi, and might be accompanied by important genomic substitutions in ON1 viruses. Close observation of such a new virus genotype introduction over time provides an opportunity to better understand the transmission and evolutionary dynamics of the pathogen. We have generated and analysed 184 RSV-A whole-genome sequences (WGSs) from Kilifi (Kenya) collected between 2011 and 2016, the first ON1 genomes from Africa and the largest collection globally from a single location. Phylogenetic analysis indicates that RSV-A circulation in this coastal Kenya location is characterized by multiple introductions of viral lineages from diverse origins but with varied success in local transmission. We identified signature amino acid substitutions between ON1 and GA2 viruses’ surface proteins (G and F), polymerase (L), and matrix M2-1 proteins, some of which were positively selected, and thereby provide an enhanced picture of RSV-A diversity. Furthermore, five of the eleven RSV open reading frames (ORFs) (G, F, L, N, and P) formed distinct phylogenetic clusters for the two genotypes. This might suggest that coding regions outside of the most frequently studied G ORF also play a role in the adaptation of RSV to host populations, with the alternative possibility that some of the substitutions are neutral and provide no selective advantage. Our analysis provides insight into the epidemiological processes that define RSV spread, highlights the genetic substitutions that characterize emerging strains, and demonstrates the utility of large-scale WGS in molecular epidemiological studies.