Transmission patterns and evolution of respiratory syncytial virus in a community outbreak identified by genomic analysis

Exploring RSV transmission patterns in different age groups in the United States

Respiratory syncytial virus (RSV) infections are a significant public health concern for pediatric populations and older adults, with seasonal winter outbreaks in the United States (US). Little is known about the timing of RSV epidemics across age groups and the relative contribution of within-group and between-group transmission of RSV in each age group. The lack of understanding of age-specific RSV transmission patterns limits our ability to inform vaccination policies. In this study, we examine the timing and transmission patterns of RSV epidemics across different age groups in 12 US states from 2018 to 2024. We found that infants under 1 year and young children aged 1-4 years experienced the earliest epidemic timing, while the elderly group had the latest. Using a semi-mechanistic age-structured spatiotemporal model, we further showed that between-group transmission greatly contributes (>50%) to the burden of RSV hospitalizations for children under 1, school-age children aged 5-17, and adults aged 18-64. By contrast, incidence in the elderly group (above 65 years) was primarily driven by transmission within the age group. Our findings indicate that distinct age groups play unique roles in propagating RSV epidemics in the US, with age-specific transmission patterns that can guide more effective RSV vaccination policies.

The genomic evolutionary dynamics and global circulation patterns of respiratory syncytial virus

Respiratory syncytial virus (RSV) is a leading cause of acute lower respiratory tract infection in young children and the second leading cause of infant death worldwide. While global circulation has been extensively studied for respiratory viruses such as seasonal influenza, and more recently also in great detail for SARS-CoV-2, a lack of global multi-annual sampling of complete RSV genomes limits our understanding of RSV molecular epidemiology. Here, we capitalise on the genomic surveillance by the INFORM-RSV study and apply phylodynamic approaches to uncover how selection and neutral epidemiological processes shape RSV diversity. Using complete viral genome sequences, we show similar patterns of site-specific diversifying selection among RSVA and RSVB and recover the imprint of non-neutral epidemic processes on their genealogies. Using a phylogeographic approach, we provide evidence for air travel governing the global patterns of RSVA and RSVB spread, which results in a considerable degree of phylogenetic mixing across countries. Our findings highlight the potential of systematic global RSV genomic surveillance for transforming our understanding of global RSV spread.

Potential Effects on Elderly People From Nirsevimab Use in Infants

Nirsevimab therapy has the potential to revolutionize infant respiratory syncytial virus (RSV) prophylaxis. But other populations suffering RSV, such the elderly or those over 60, may also be protected by using this novel antibody in the infant group. It is true that some studies link the use of nirsevimab to a reduction in the virus's ability to spread by lowering the viral load in infants as a result of the drug's long half-life. However, this protective effect may not be very significant because RSV transmission in the elderly typically comes from other elderly people or from school-aged children. Furthermore, RSV may be transmitted at any time of the year and not just during the period of nirsevimab protection due to its existence in human reservoirs. The reasons made here show that, even though nirsevimab treatment in infants may protect the elderly, this benefit would be limited and testimonial. Therefore, immunizing the elderly with currently licensed and developing vaccines should be a priority.

Genomic characterisation of respiratory syncytial virus: a novel system for whole genome sequencing and full-length G and F gene sequences

To advance our understanding of respiratory syncytial virus (RSV) impact through genomic surveillance, we describe two PCR-based sequencing systems, (i) RSVAB-WGS for generic whole-genome sequencing and (ii) RSVAB-GF, which targets major viral antigens, G and F, and is used as a complement for challenging cases with low viral load. These methods monitor RSV genetic diversity to inform molecular epidemiology, vaccine effectiveness and treatment strategies, contributing also to the standardisation of surveillance in a new era of vaccines.

High-throughput sequencing approaches applied to SARS-CoV-2

High-throughput sequencing is crucial for surveillance and control of viral outbreaks. During the ongoing coronavirus disease 2019 (COVID-19) pandemic, advances in the high-throughput sequencing technology resources have enhanced diagnosis, surveillance, and vaccine discovery. From the onset of the pandemic in December 2019, several genome-sequencing approaches have been developed and supported across the major sequencing platforms such as Illumina, Oxford Nanopore, PacBio, MGI DNBSEQTM and Ion Torrent. Here, we share insights from the sequencing approaches developed for sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) between December 2019 and October 2022.

Nationwide and long-term molecular epidemiologic studies of mumps viruses that circulated in Japan between 1986 and 2017

In Japan, major mumps outbreaks still occur every 4–5 years because of low mumps vaccine coverage (30–40%) owing to the voluntary immunization program. Herein, to prepare for a regular immunization program, we aimed to reveal the nationwide and long-term molecular epidemiological trends of the mumps virus (MuV) in Japan. Additionally, we performed whole-genome sequencing (WGS) using next-generation sequencing to assess results from conventional genotyping using MuV sequences of the small-hydrophobic (SH) gene. We analyzed 1,064 SH gene sequences from mumps clinical samples and MuV isolates collected from 25 prefectures from 1986 to 2017. The results showed that six genotypes, namely B (110), F (1), G (900), H (3), J (41), and L (9) were identified, and the dominant genotypes changed every decade in Japan since the 1980s. Genotype G has been exclusively circulating since the early 2000s. Seven clades were identified for genotype G using SH sequence-based classification. To verify the results, we performed WGS on 77 representative isolates of genotype G using NGS and phylogenetically analyzed them. Five clades were identified with high bootstrap values and designated as Japanese clade (JPC)-1, -2, -3, -4, -5. JPC-1 and -3 accounted for over 80% of the total genotype G isolates (68.3 and 13.8%, respectively). Of these, JPC-2 and -5, were newly identified clades in Japan through this study. This is the first report describing the nationwide and long-term molecular epidemiology of MuV in Japan. The results provide information about Japanese domestic genotypes, which is essential for evaluating the mumps elimination progress in Japan after the forthcoming introduction of the mumps vaccine into Japan’s regular immunization program. Furthermore, the study shows that WGS analysis using NGS is more accurate than results obtained from conventional SH sequence-based classification and is a powerful tool for accurate molecular epidemiology studies.

Methods Combining Genomic and Epidemiological Data in the Reconstruction of Transmission Trees: A Systematic Review

In order to better understand transmission dynamics and appropriately target control and preventive measures, studies have aimed to identify who-infected-whom in actual outbreaks. Numerous reconstruction methods exist, each with their own assumptions, types of data, and inference strategy. Thus, selecting a method can be difficult. Following PRISMA guidelines, we systematically reviewed the literature for methods combing epidemiological and genomic data in transmission tree reconstruction. We identified 22 methods from the 41 selected articles. We defined three families according to how genomic data was handled: a non-phylogenetic family, a sequential phylogenetic family, and a simultaneous phylogenetic family. We discussed methods according to the data needed as well as the underlying sequence mutation, within-host evolution, transmission, and case observation. In the non-phylogenetic family consisting of eight methods, pairwise genetic distances were estimated. In the phylogenetic families, transmission trees were inferred from phylogenetic trees either simultaneously (nine methods) or sequentially (five methods). While a majority of methods (17/22) modeled the transmission process, few (8/22) took into account imperfect case detection. Within-host evolution was generally (7/8) modeled as a coalescent process. These practical and theoretical considerations were highlighted in order to help select the appropriate method for an outbreak.

Mohammed K, W. Lambisia A, R. de Laurent Z, W. Mburu M, M. Ong’era E, M. Morobe J, Otieno E, Abdou Azali H, Said Abdallah K, Diarra A, Ahmed Yahaya A, Borus P, Gumede Moeletsi N, Fred Athanasius D, Tsofa B, Bejon P, James Nokes D, Isabella Ochola-Oyier L. Detection of SARS-CoV-2 variant 501Y.V2 in Comoros Islands in January 2021

Background. Genomic data is key in understanding the spread and evolution of SARS-CoV-2 pandemic and informing the design and evaluation of interventions. However, SARS-CoV-2 genomic data remains scarce across Africa, with no reports yet from the Indian Ocean islands. Methods. We genome sequenced six SARS-CoV-2 positive samples from the first major infection wave in the Union of Comoros in January 2021 and undertook detailed phylogenetic analysis. Results. All the recovered six genomes classified within the 501Y.V2 variant of concern (also known as lineage B.1.351) and appeared to be from 2 sub-clusters with the most recent common ancestor dated 30 th Oct-2020 (95% Credibility Interval: 06 th Sep-2020 to 10 th Dec-2020). Comparison of the Comoros genomes with those of 501Y.V2 variant of concern from other countries deposited into the GISAID database revealed their close association with viruses identified in France and Mayotte (part of the Comoros archipelago and a France, Overseas Department). Conclusions. The recovered genomes, albeit few, confirmed local transmission following probably multiple introductions of the SARS-CoV-2 501Y.V2 variant of concern during the Comoros's first major COVID-19 wave. These findings demonstrate the importance of genomic surveillance and have implications for ongoing control strategies on the islands.

Functional Features of the Respiratory Syncytial Virus G Protein

Respiratory syncytial virus (RSV) is a major cause of serious lower respiratory tract infections in children <5 years of age worldwide and repeated infections throughout life leading to serious disease in the elderly and persons with compromised immune, cardiac, and pulmonary systems. The disease burden has made it a high priority for vaccine and antiviral drug development but without success except for immune prophylaxis for certain young infants. Two RSV proteins are associated with protection, F and G, and F is most often pursued for vaccine and antiviral drug development. Several features of the G protein suggest it could also be an important to vaccine or antiviral drug target design. We review features of G that effect biology of infection, the host immune response, and disease associated with infection. Though it is not clear how to fit these together into an integrated picture, it is clear that G mediates cell surface binding and facilitates cellular infection, modulates host responses that affect both immunity and disease, and its CX3C aa motif contributes to many of these effects. These features of G and the ability to block the effects with antibody, suggest G has substantial potential in vaccine and antiviral drug design.

Community factors associated with local epidemic timing of respiratory syncytial virus: A spatiotemporal modeling study

Respiratory syncytial virus (RSV) causes a large burden of morbidity in young children and the elderly. Spatial variability in the timing of RSV epidemics provides an opportunity to probe the factors driving its transmission, including factors that influence epidemic seeding and growth rates. Using hospitalization data from Connecticut, New Jersey, and New York, we estimated epidemic timing at the ZIP code level using harmonic regression and then used a Bayesian meta-regression model to evaluate correlates of epidemic timing. Earlier epidemics were associated with larger household size and greater population density. Nearby localities had similar epidemic timing. Our results suggest that RSV epidemics grow faster in areas with more local contact opportunities, and that epidemic spread follows a spatial diffusion process based on geographic proximity. Our findings can inform the timing of delivery of RSV extended half-life prophylaxis and maternal vaccines and guide future studies on the transmission dynamics of RSV.

Mathematical modelling of respiratory syncytial virus (RSV) in low- and middle-income countries: A systematic review

Background:

Due to high burden of respiratory syncytial virus (RSV) in low- and middle-income countries (LMIC), international funding organizations have prioritized the development of RSV vaccines. Mathematical models of RSV will play an important role in assessing the relative value of these interventions. Our objectives were to provide an overview of the existing RSV modelling literature in LMIC and summarize available results on population-level effectiveness and cost-effectiveness.

Methods:

We searched MEDLINE from 2000 to 2020 for English language publications that employed a mathematical model of RSV calibrated to LMIC. Qualitative data were extracted on study and model characteristics. Quantitative data were collected on key model input assumptions and base case effectiveness and cost-effectiveness estimates for various immunization strategies.

Findings:

Of the 283 articles reviewed, 15 met inclusion criteria. Ten studies used modelling techniques to explore RSV transmission and/or natural history, while eight studies evaluated RSV vaccines and/or monoclonal antibodies, three of which included cost-effectiveness analyses. Six studies employed deterministic compartmental models, five studies employed individual transmission models, and four studies used different types of cohort models. Nearly every model was calibrated to at least one middle-income country, while four were calibrated to low-income countries.

Interpretation:

The mathematical modelling literature in LMIC has demonstrated the potential effectiveness of RSV vaccines and monoclonal antibodies. This review has demonstrated the importance of accounting for seasonality, social contact rates, immunity from prior infection and maternal antibody transfer. Future models should consider incorporating individual-level risk factors, subtype-specific effects, long-term sequelae of RSV infections, and out-of-hospital mortality.

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.

Towards a unified classification for human respiratory syncytial virus genotypes

Since the first human respiratory syncytial virus (HRSV) genotype classification in 1998, inconsistent conclusions have been drawn regarding the criteria that define HRSV genotypes and their nomenclature, challenging data comparisons between research groups. In this study, we aim to unify the field of HRSV genotype classification by reviewing the different methods that have been used in the past to define HRSV genotypes and by proposing a new classification procedure, based on well-established phylogenetic methods. All available complete HRSV genomes (>12,000 bp) were downloaded from GenBank and divided into the two subgroups: HRSV-A and HRSV-B. From whole-genome alignments, the regions that correspond to the open reading frame of the glycoprotein G and the second hypervariable region (HVR2) of the ectodomain were extracted. In the resulting partial alignments, the phylogenetic signal within each fragment was assessed. Maximum likelihood phylogenetic trees were reconstructed using the complete genome alignments. Patristic distances were calculated between all pairs of tips in the phylogenetic tree and summarized as a density plot in order to determine a cutoff value at the lowest point following the major distance peak. Our data show that neither the HVR2 fragment nor the G gene contains sufficient phylogenetic signal to perform reliable phylogenetic reconstruction. Therefore, whole-genome alignments were used to determine HRSV genotypes. We define a genotype using the following criteria: a bootstrap support of ≥ 70 per cent for the respective clade and a maximum patristic distance between all members of the clade of ≤0.018 substitutions per site for HRSV-A or ≤0.026 substitutions per site for HRSV-B. By applying this definition, we distinguish twenty-three genotypes within subtype HRSV-A and six genotypes within subtype HRSV-B. Applying the genotype criteria on subsampled data sets confirmed the robustness of the method.

Primary and Repeated Respiratory Viral Infections Among Infants in Rural Nepal

BACKGROUND

Respiratory viruses cause significant morbidity and death in infants; 99% of such deaths occur in resource-limited settings. Risk factors for initial and repeated respiratory viral infections in young infants in resource-limited settings have not been well described.

METHODS

From 2011 to 2014, a birth cohort of infants in rural Nepal was enrolled and followed with weekly household-based active surveillance for respiratory symptoms until 6 months of age. Respiratory illness was defined as having any of the following: fever, cough, wheeze, difficulty breathing, and/or a draining ear. We tested nasal swabs of infants with respiratory illness for multiple respiratory viruses by using a reverse transcription polymerase chain reaction assay. The risk of primary and repeated infections with the same virus was evaluated using Poisson regression.

RESULTS

Of 3528 infants, 1726 (49%) had a primary infection, and 419 (12%) had a repeated infection. The incidences of respiratory viral infection in infants were 1816 per 1000 person-years for primary infections and 1204 per 1000 person-years for repeated infection with the same virus. Exposure to other children and male sex were each associated with an increased risk for primary infection (risk ratios, 1.13 [95% confidence interval (CI), 1.06-1.20] and 1.14 [95% CI, 1.02-1.27], respectively), whereas higher maternal education was associated with a decreased risk for both primary and repeated infections (risk ratio, 0.96 [95% CI, 0.95-0.98]). The incidence of subsequent infection did not change when previous infection with the same or another respiratory virus occurred. Illness duration and severity were not significantly different in the infants between the first and second episodes for any respiratory virus tested.

CONCLUSIONS

In infants in rural Nepal, repeated respiratory virus infections were frequent, and we found no decrease in illness severity with repeated infections and no evidence of replacement with another virus. Vaccine strategies and public health interventions that provide durable protection in the first 6 months of life could decrease the burden of repeated infections by multiple respiratory viruses, particularly in low-resource countries.

High-throughput sequencing (HTS) for the analysis of viral populations

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

Non-gradient and genotype-dependent patterns of RSV gene expression

Respiratory syncytial virus (RSV) is a nonsegmented negative-strand RNA virus (NSV) and a leading cause of severe lower respiratory tract illness in infants and the elderly. Transcription of the ten RSV genes proceeds sequentially from the 3’ promoter and requires conserved gene start (GS) and gene end (GE) signals. Previous studies using the prototypical GA1 genotype Long and A2 strains have indicated a gradient of gene transcription extending across the genome, with the highest level of mRNA coming from the most promoter-proximal gene, the first nonstructural (NS1) gene, and mRNA levels from subsequent genes dropping until reaching a minimum at the most promoter-distal gene, the polymerase (L) gene. However, recent reports show non-gradient levels of mRNA, with higher than expected levels from the attachment (G) gene. It is unknown to what extent different transcript stabilities might shape measured mRNA levels. It is also unclear whether patterns of RSV gene expression vary, or show strain- or genotype-dependence. To address this, mRNA abundances from five RSV genes were measured by quantitative real-time PCR (qPCR) in three cell lines and in cotton rats infected with RSV isolates belonging to four genotypes (GA1, ON, GB1, BA). Relative mRNA levels reached steady-state between four and 24 hours post-infection. Steady-state patterns were non-gradient and genotype-specific, where mRNA levels from the G gene exceeded those from the more promoter-proximal nucleocapsid (N) gene across isolates. Transcript stabilities could not account for the non-gradient patterns observed, indicating that relative mRNA levels more strongly reflect transcription than decay. Our results indicate that gene expression from a small but diverse set of RSV genotypes is non-gradient and genotype-dependent. We propose novel models of RSV transcription that can account for non-gradient transcription.

Genomic analysis of respiratory syncytial virus infections in households and utility in inferring who infects the infant

Infants (under 1-year-old) are at most risk of life threatening respiratory syncytial virus (RSV) disease. RSV epidemiological data alone has been insufficient in defining who acquires infection from whom (WAIFW) within households. We investigated RSV genomic variation within and between infected individuals and assessed its potential utility in tracking transmission in households. Over an entire single RSV season in coastal Kenya, nasal swabs were collected from members of 20 households every 3-4 days regardless of symptom status and screened for RSV nucleic acid. Next generation sequencing was used to generate >90% RSV full-length genomes for 51.1% of positive samples (191/374). Single nucleotide polymorphisms (SNPs) observed during household infection outbreaks ranged from 0-21 (median: 3) while SNPs observed during single-host infection episodes ranged from 0-17 (median: 1). Using the viral genomic data alone there was insufficient resolution to fully reconstruct within-household transmission chains. For households with clear index cases, the most likely source of infant infection was via a toddler (aged 1 to <3 years-old) or school-aged (aged 6 to <12 years-old) co-occupant. However, for best resolution of WAIFW within households, we suggest an integrated analysis of RSV genomic and epidemiological data.

White L, Poovorawan K, Soonthornworasiri N, Sukontamarn P, Chanthavilay P, F. Medley G, Pan-ngum W. Modeling household dynamics on Respiratory Syncytial Virus (RSV)

Respiratory Syncytial Virus (RSV) is the most common cause of respiratory tract infection in infants and children and shows increasing trend among elderly people worldwide. In many developing country settings, population and household structures have gone through some significant changes in the past decades, namely fewer births, more elderly population, and smaller household size but more RSV high-risk individuals. These dynamics have been captured in a mathematical model with RSV transmission dynamics to predict the disease burden on the detailed population for future targeted interventions. The population and disease dynamics model was constructed and tested against the hospitalization data for Acute Lower Respiratory Tract Infection due to RSV in rural Thai settings between 2005 and 2011. The proportion of extended families is predicted to increase by about 10% from 2005 to 2020, especially for those with elderly population, while the classic nuclear family type (with adults and children) will decline by about 10%. For RSV, infections from extended family type (approximately 60% of all household types) have majorly contributed to the force of infection (FOI). While the model predicted the increase of FOI from the extended family by 15% from 2005 to 2020, the FOI contributed by other household types would be either stable or decrease in the same time period. RSV incidence rate is predominantly high among babies (92.2%) and has been predicted to decrease slightly over time (from 940 to 864 cases per 100,000 population by 2020), while the incidence rates among children and elderly people may remain steadily low over the same period. However, the estimated incidence rates among elderly people were twice than those in children. The model predicts that approximately 60% of FOI for RSV will come from members of the extended family type. The incidence rate of RSV among children and elderly in extended families was about 20 times lower than that in infants and the trend is steady. Targeted intervention strategies, such as health education in some specific groups and targeted vaccination, may be considered, with the focus on extended family type. Target interventions on babies can lessen the transmission to children and elderly especially when transmission within households of extended family type is high.

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.