Molecular evolution of the fusion protein gene in human respiratory syncytial virus subgroup A

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

OBJECTIVES

This study investigated the prevalence, genetic diversity, and evolution of human respiratory syncytial virus (HRSV) in Barcelona from 2013 to 2023.

METHODS

Respiratory specimens from patients with RTI suspicion at Hospital Universitari Vall d'Hebron were collected from October 2013 to May 2023 for laboratory-confirmation of respiratory viruses. Next-generation sequencing was performed in randomly-selected samples with Illumina technology. Phylogenetic analyses of whole genome sequences were performed with BEAST v1.10.4. Signals of selection and evolutionary pressures were inferred by population dynamics and evolutionary analyses. Mutations in major surface proteins were genetic and structurally characterised, emphasizing those within antigenic epitopes.

RESULTS

Analyzing 139,625 samples, 5.3% were HRSV-positive (3008 HRSV-A, 3882 HRSV-B, 56 HRSV-A and -B, and 495 unsubtyped HRSV), with a higher prevalence observed in the paediatric population. Pandemic-related shifts in seasonal patterns returned to normal in 2022-2023. A total of 198 whole-genome sequences were obtained for HRSV-A (6.6% of the HRSV-A positive samples) belonging to GA2.3.5 lineage. For HRSV-B, 167 samples were sequenced (4.3% of the HRSV-B positive samples), belonging to GB5.0.2, GB5.0.4a and GB5.0.5a. HRSV-B exhibited a higher evolution rate. Post-SARS-CoV-2 pandemic, both subtypes showed increased evolutionary rates and decreased effective population size initially, followed by a sharp increase. Analyses indicated negative selective pressure on HRSV. Mutations in antigenic epitopes, including S276N and M274I in palivizumab-targeted site II, and I206M, Q209R, and S211N in nirsevimab-targeted site Ø, were identified.

DISCUSSION

Particularly in the context of the large-scale use in 2023-2024 season of nirsevimab, continuous epidemiological and genomic surveillance is crucial.

Characterization of a Panel of Monoclonal Antibodies Targeting the F-Protein of the Respiratory Syncytial Virus (RSV) for the Typing of Contemporary Circulating Strains

Respiratory syncytial virus (RSV) is the most common cause of upper and lower respiratory tract infections in infants and young children. Virus-specific monoclonal antibodies (mAbs) can be used for diagnosis, prophylaxis, and research of RSV pathogenesis. A panel of 16 anti-RSV mAbs was obtained from mice immunized by RSV strain Long. Half of them had virus-neutralizing activity. According to Western blot all of these mAbs effectively bound native oligomeric (homodimeric and homotrimeric) forms of the RSV fusion (F) protein. Only five of the mAbs interacted with the monomeric form, and only one of these possessed neutralizing activity. None of these mAbs, nor the commercial humanized neutralizing mAb palivizumab, reacted with the denaturated F protein. Thus, interaction of all these mAbs with F protein had clear conformational dependence. Competitive ELISA and neutralization assays allowed the identification of nine antigenic target sites for the interaction of mAb with the F protein. Five partially overlapping sites may represent a complex spatial structure of one antigenic determinant, including one neutralizing and four non-neutralizing epitopes. Four sites (three neutralizing and one non-neutralizing) were found to be distinct. As a result of virus cultivation RSV-A, strain Long, in the presence of a large amount of one of the neutralizing mAbs, an escape mutant with a substitution, N240S, in the F protein, was obtained. Thus, it was shown for the first time that position 240 is critical for the protective effect of an anti-RSV antibody. To assess the ability of these mAbs to interact with modern RSV strains circulating in St. Petersburg (Russia) between 2014 and 2022, 73 RSV-A and 22 RSV-B isolates were analyzed. Six mAbs were directed to conserved epitopes of the F protein as they interacted most efficiently with both RSV subtypes in a fixed cell-ELISA and could be used for diagnostic assays detecting RSV.

Prevalence, variation, and transmission patterns of human respiratory syncytial virus from pediatric patients in Hubei, China during 2020-2021

Human respiratory syncytial virus (RSV) is a severe threat to children and a main cause of acute lower respiratory tract infections. Nevertheless, the intra-host evolution and inter-regional diffusion of RSV are little known. In this study, we performed a systematic surveillance in hospitalized children in Hubei during 2020-2021, in which 106 RSV-positive samples were detected both clinically and by metagenomic next generation sequencing (mNGS). RSV-A and RSV-B groups co-circulated during surveillance with RSV-B being predominant. About 46 high-quality genomes were used for further analyses. A total of 163 intra-host nucleotide variation (iSNV) sites distributed in 34 samples were detected, and glycoprotein (G) gene was the most enriched gene for iSNVs, with non-synonymous substitutions more than synonymous substitutions. Evolutionary dynamic analysis showed that the evolutionary rates of G and NS2 genes were higher, and the population size of RSV groups changed over time. We also found evidences of inter-regional diffusion from Europe and Oceania to Hubei for RSV-A and RSV-B, respectively. This study highlighted the intra-host and inter-host evolution of RSV, and provided some evidences for understanding the evolution of RSV.

M2-2 gene as a new alternative molecular marker for phylogenetic, phylodynamic, and evolutionary studies of hRSV

The human Respiratory Syncytial Virus (hRSV) is the main causative agent of acute respiratory infections (ARI), such as pneumonia and bronchiolitis. One of the factors that lead to success in viral replication is the interaction of the M2-2 protein with the ribosomal complex. This interaction is responsible for the phase change of viral activity, acting as an inhibitor or inducer of viral replication, according to the concentration of mRNA. Based on the importance of M2-2 gene and protein have to viral physiology, we performed here evaluations of genetic diversity, phylogenetic reconstructions, phylodynamics, and selection test. Our results suggested an alternative way of classifying this virus in clades A and B, based on a new phylogenetic marker, the M2-2 gene. Therefore, our study is the first one to investigate the dynamics of the evolutionary diversification process of hRSV from the perspective of the M2-2 viral gene. In our study was also identified that the M2-2 gene is under the effect of purifying selection originated by population genetic bottlenecks. Therefore, the M2-2 gene demonstrated an interesting potential to be applied in evolutionary studies involving hRSV, recovering phylogenetic signals and traits of natural selection under the evolution of this virus.

Detailed Evolutionary Analyses of the F Gene in the Respiratory Syncytial Virus Subgroup A

We performed evolution, phylodynamics, and reinfection-related antigenicity analyses of respiratory syncytial virus subgroup A (RSV-A) fusion (F) gene in globally collected strains (1465 strains) using authentic bioinformatics methods. The time-scaled evolutionary tree using the Bayesian Markov chain Monte Carlo method estimated that a common ancestor of the RSV-A, RSV-B, and bovine-RSV diverged at around 450 years ago, and RSV-A and RSV-B diverged around 250 years ago. Finally, the RSV-A F gene formed eight genotypes (GA1-GA7 and NA1) over the last 80 years. Phylodynamics of RSV-A F gene, including all genotype strains, increased twice in the 1990s and 2010s, while patterns of each RSV-A genotype were different. Phylogenetic distance analysis suggested that the genetic distances of the strains were relatively short (less than 0.05). No positive selection sites were estimated, while many negative selection sites were found. Moreover, the F protein 3D structure mapping and conformational epitope analysis implied that the conformational epitopes did not correspond to the neutralizing antibody binding sites of the F protein. These results suggested that the RSV-A F gene is relatively conserved, and mismatches between conformational epitopes and neutralizing antibody binding sites of the F protein are responsible for the virus reinfection.

Genetic diversity and molecular evolution of human respiratory syncytial virus A and B

Human respiratory syncytial viruses (RSVs) are classified into two major groups (A and B) based on antigenic differences in the G glycoprotein. To investigate circulating characteristics and phylodynamic history of RSV, we analyzed the genetic variability and evolutionary pattern of RSVs from 1977 to 2019 in this study. The results revealed that there was no recombination event of intergroup. Single nucleotide polymorphisms (SNPs) were observed through the genome with the highest occurrence rate in the G gene. Five and six sites in G protein of RSV-A and RSV-B, respectively, were further identified with a strong positive selection. The mean evolutionary rates for RSV-A and -B were estimated to be 1.48 × 10^–3 and 1.92 × 10^–3 nucleotide substitutions/site/year, respectively. The Bayesian skyline plot showed a constant population size of RSV-A and a sharp expansion of population size of RSV-B since 2005, and an obvious decrease 5 years later, then became stable again. The total population size of RSVs showed a similar tendency to that of RSV-B. Time-scaled phylogeny suggested a temporal specificity of the RSV-genotypes. Monitoring nucleotide changes and analyzing evolution pattern for RSVs could give valuable insights for vaccine and therapy strategies against RSV infection.

Evolutionary analysis of human respiratory syncytial virus collected in Myanmar between 2015 and 2018

We studied genetic variation in the second hypervariable region (HVR) of the G gene of human respiratory syncytial virus (HRSV) from 1701 nasal swab samples collected from outpatients with acute respiratory infections at two general hospitals in the cities Yangon and Pyinmana in Myanmar from 2015 to 2018. HRSV genotypes were characterized using phylogenetic trees constructed using the maximum likelihood method. Time-scale phylogenetic tree analyses were performed using the Bayesian Markov chain Monte Carlo method. In total, 244 (14.3%) samples were HRSV-positive and were classified as HRSV-A (n = 84, 34.4%), HRSV-B (n = 158, 64.8%), and co-detection of HRSV-A/HRSV-B (n = 2, 0.8%). HRSV epidemics occurred seasonally between July (1.9%, 15/785) and August (10.5%, 108/1028), with peak infections in September (35.8%, 149/416) and October (58.2%, 89/153). HRSV infection rate was higher in children ≥1 year of age than in those <1 year of age (70.5% vs. 29.5%). The most common HRSV symptoms in children were cough (80%-90%) and rhinorrhea (70%-100%). The predominant genotypes were ON1for HRSV-A (78%) and BA9 for HRSV-B (64%). Time to the most recent common ancestor was 2014 (95% highest posterior density [HPD], 2012-2015) for HRSV-A ON1 and 2009 (95% HPD, 2004-2012) for HRSV-B BA9. The mean evolutionary rate (substitutions/site/year) for HRSV-B (2.12 × 10^-2, 95% HPD, 8.53 × 10^-3-3.63 × 10^-2) was slightly higher than that for HRSV-A (1.39 × 10^-2, 95% HPD, 6.03 × 10^-3-2.12 × 10^-2). The estimated effective population size (diversity) for HRSV-A increased from 2015 to 2016 and declined in mid-2018, whereas HRSV-B diversity was constant in 2015 and 2016 and increased in mid-2017. In conclusion, the dominant HRSV-A and HRSV-B genotypes in Myanmar were ON1 and BA9, respectively, between 2015 and 2018. HRSV-B evolved slightly faster than HRSV-A and exhibited unique phylogenetic characteristics.

Genetic characterization of G protein in respiratory syncytial virus ON-1 genotype in Tehran

Aim: We investigated the genetic characterization of the respiratory syncytial virus (RSV) ON-1 genotypes and their different lineages based on the G gene among children <2 years of age presenting with acute respiratory tract infections in Tehran, Iran. Materials & methods: A phylogenetic tree from the Iranian samples and ON-1 strains of various parts of the world were constructed. The amino acid composition of the RSV G protein of the ON-1 genotype was mapped. Results: Human RSV ON-1 genotypes from the Iranian samples clustered in three lineages. The most common amino acid substitutions were as follows: X218Q, I240S, L289P, Y304H and L310P. Conclusion: Continuing molecular epidemiological surveys in other regions of Iran will provide deeper insight into the nature of this replacement of the dominant RSV genotype from GA2 to ON-1 in Iran.

Evolutionary dynamics of influenza B strains detected from paediatric acute respiratory infections in central Vietnam

Influenza virus B belongs to the family Orthomyxoviridae with segmented negative-sense RNA genomes. Since 1970s, influenza B has diverged intoVictoria and Yamagata, which differs in antigenic and evolutionary characteristics. Yet, molecular-epidemiological information of influenza B from developing nations is limited. In central Vietnam, influenza A subtype-specific circulation pattern and clinical characteristics were previously described. However, molecular evolutionary characteristics of influenza B has not been discussed to date. We utilized the influenza B positives obtained from paediatric ARI surveillance during 2007-2013. Influenza B HA and NA genes were amplified, sequenced, and phylogenetic/molecular evolutionary analysis was performed using Maximum Likelihood and Bayesian MCMC. Phylodynamics analysis was performed with Bayesian Skyline Plot (BSP). Furthermore, we performed selection pressure analysis and estimated N-glycosylation sites. In the current study, overall positive rate for influenza B was 3.0%, and Victoria lineage immediately became predominant in post-A/H1N1pdm09 period. The noticeable shift in Victoria lineage WHO Group occurred. With respect to the evolutionary rate (substitutions/site/year), Victoria lineage HA gene was evolving faster than Yamagata lineage (2.43 × 10^-3 vs 2.00 × 10^-3). Furthermore, the evolutionary rate of Victoria Group 5 was greater than Group 1. BSP presented the rapid growth in Effective Population Size (EPS) of Victoria lineage occurred soon after the 1st A/H1N1pdm09 case was detected whereas the EPS of Yamagata lineage was stable for both genes. N-glycosylation pattern between lineages and among WHO Groups were slightly different, and HA gene had a total of 6 amino acid substitutions under positive section pressure (4 for Victoria and 2 for Yamagata). The current results highlight the importance of Victoria lineage in post-A/H1N1pdm09 period. Difference in evolutionary characteristics and phylodynamics may indicate lineage and WHO Group-specific evolutionary dynamics. It is necessary to further continue the molecular-epidemiological surveillance in local setting to gain a better understanding of local evolutionary characteristics of influenza B strains.

Molecular Evolution of the Fusion Protein (F) Gene in Human Respirovirus 3

To elucidate the evolution of human respirovirus 3 (HRV3), we performed detailed genetic analyses of the F gene (full-length) detected from hundreds of HRV3 strains obtained from various geographic regions. First, we performed time-scaled evolutionary analyses using the Bayesian Markov chain Monte Carlo method. Then, we performed analyses of phylodynamics, similarity, phylogenetic distance, selective pressure, and conformational B-cell epitope with the F-protein structural analyses. Time-scaled phylogenetic tree showed that the common ancestor of HRV3 and bovine respirovirus 3 diverged over 300 years ago and subdivided it into three major clusters and four subclusters during the most recent 100 years. The overall evolutionary rate was approximately 10^-3 substitutions/site/year. Indigenous similarity was seen in the present strains, and the mean phylogenetic distance were 0.033. Many negative selection sites were seen in the ectodomain. The conformational epitopes did not correspond to the neutralizing antibody binding sites. These results suggest that the HRV3 F gene is relatively conserved and restricted in this diversity to preserve the protein function, although these strains form many branches on the phylogenetic tree. Furthermore, HRV3 reinfection may be responsible for discordances between the conformational epitopes and the neutralizing antibody binding sites of the F protein. These findings contribute to a better understanding of HRV3 virology.

Molecular characterization of the glycoprotein and fusion protein in human respiratory syncytial virus subgroup A: Emergence of ON-1 genotype in Iran

HRSV is a principle cause of infant hospitalization, childhood wheezing and a common pathogen in the elderly. Limited information exists regarding HRSV genotypes in Iran. In order to better understand HRSV strain diversity, we performed an in-depth evaluation of the genetic variability of the HRSV F protein detected in children under two years of age that, presented with acute respiratory symptoms during 2015-2016 in Tehran. A total of 180 nasopharyngeal swabs were evaluated. The HRSV positive samples were genotyped for G and F gene sequences using RT-PCR and sequencing methods. Phylogenetic analysis was performed using the neighbor-joining and maximum likelihood methods. Genetic and antigenic characteristics of the F gene, nucleotide and amino acids in significant positions and immune system binding regions, as well as the p-distance, positive/negative selection site, linear epitopes and glycosylation sites were investigated in all selected sequences. Among the 83 HRSV positive samples, the Fifty-five cases were successfully sequenced. All of them were classified as subgroup A and belonged to the ON-1 genotype, which possessed 72-nt duplication in the G gene. This study is the first report on the emergence of ON-1 in Iran. ON-1 Iranian sequences clustered in three lineages according to virus fusion (F) gene variations. F gene sequence analysis showed that all genetic changes in the isolates from Iran were base substitutions and no deletion/insertions were identified. The low dN/dS ratio and lack of positively selected sites showed that the fusion genes found in the strains from Iran are not under host selective pressure. Continuing and long-term molecular epidemiological surveys for early detection of circulating and newly emerging genotypes are necessary to gain a better understanding of their epidemic potential.

Persistence and continuous evolution of the human respiratory syncytial virus in northern Taiwan for two decades

The study aimed to characterize the molecular epidemiology, phylogenetic relationship, and population dynamics of the G protein gene in clinical respiratory syncytial virus (RSV) strains isolated from northern Taiwan. We analyzed a total of 160 and 116 G protein gene sequences of RSV-A and RSV-B representative strains, respectively, from 804 clinical viral stocks collected between July 2000 and June 2016. Population dynamic patterns of the RSV G protein gene were analyzed using Bayesian inference through the Markov chain Monte Carlo framework. A phylogenetic analysis revealed that RSV-A from Taiwan could be categorized into GA2, GA5, and GA7 lineages. GA2 of RSV-A could be further divided into NA1, NA2, NA4, and ON1 clades. These RSV-A lineages has been replaced over time, whereas RSV-B strains from Taiwan continually evolved from a single lineage with significant time-dependent waves. Four putative positive selection sites were observed in both RSV-A and RSV-B. The Bayesian skyline plot revealed that the local population dynamics of RSV were associated with lineage displacement events. Both circulating subtypes and population dynamics represented a unique local pattern. Our results affirm the necessity of continuing molecular surveillance of RSV to attain a more comprehensive understanding of epidemics.

Detailed genetic analyses of the HN gene in human respirovirus 3 detected in children with acute respiratory illness in the Iwate Prefecture, Japan

We performed detailed genetic analyses of the partial hemagglutinin-neuraminidase (HN) gene in 34 human respirovirus 3 (HRV3) strains from children with acute respiratory illness during 2013-2015 in Iwate Prefecture, Japan. In addition, we performed analyses of the evolutionary timescale of the gene using the Bayesian Markov chain Monte Carlo (MCMC) method. Furthermore, we analyzed pairwise distances and performed selective pressure analyses followed by linear B-cell epitope mapping and N-glycosylation and phylodynamic analyses. A phylogenetic tree showed that the strains diversified at around 1939, and the rate of molecular evolution was 7.6 × 10^-4 substitutions/site/year. Although the pairwise distances were relatively short (0.03 ± 0.018 [mean ± standard deviation, SD]), two positive selection sites (Cys544Trp and Leu555Ser) and no amino acid substitutions were found in the active/catalytic sites. Six epitopes were estimated in this study, and three mouse monoclonal antibody binding sites (amino acid positions 278, 281, and 461) overlapped with two epitopes belonging to subcluster C3 strains. Bayesian skyline plot analyses indicated that subcluster C3 strains have been increasing from 2004, whereas subcluster C1 strains have declined from 2004. Based on these results, Iwate strains were divided into two subclusters and each subcluster evolved independently. Moreover, our results suggested that some predicted linear epitopes (epitopes 3 and 5) are candidates for an HRV3 vaccine motif. To better understand the details of the molecular evolution of HRV, further studies are needed.

GENETIC AND ANTIGENIC CHARACTERISTICS OF RESPIRATORY SYNCYTIAL VIRUS STRAINS ISOLATED IN ST. PETERSBURG IN 2013-2016

Antigenic and genetic characteristics of Russian RSV isolates are presented for the first time. Of the 69 strains isolated in St. Petersburg, 93% belonged to the RSV-A antigenic group. The antigenic variations in the F-protein RSV were analyzed using a panel from 6 monoclonal antibodies by the method of micro-cultural ELISA. Depending on the decrease in the effectiveness of interaction with monoclonal antibodies (relative to the reference strain Long), RSV-A isolates were divided into 4 antigenic subgroups. The results of 24 isolates sequencing showed that more than 60% of them had substitutions in significant F-protein sites compared to the ON67-1210A reference strain of the current RSV genotype ON1/GA2. The most variable were the signal peptide and antigenic site II. When comparing the results of ELISA and sequencing, it was not possible to identify any specific key substitutions in the amino acid sequence of the F-protein that affect the interaction of the virus with antibodies. The nucleotide sequence of the F-gene from 19 of the 24 characterized isolates was close to that of ON67-1210A reference virus and was significantly different from RSV-A Long and A2 viruses. A separate group consisted of 5 strains, in which the F-protein structure was approximated to RSV Long.

Molecular evolution of the fusion protein (F) gene in human respiratory syncytial virus subgroup B

In this study, we examined the molecular evolution of the fusion protein (F) gene in human respiratory syncytial virus subgroup B (HRSV-B). First, we performed time-scale evolution analyses using the Bayesian Markov chain Monte Carlo (MCMC) method. Next, we performed genetic distance, linear B-cell epitope prediction, N-glycosylation, positive/negative selection site, and Bayesian skyline plot analyses. We also constructed a structural model of the F protein and mapped the amino acid substitutions and the predicted B-cell epitopes. The MCMC-constructed phylogenetic tree indicated that the HRSV F gene diverged from the bovine respiratory syncytial virus gene approximately 580years ago and had a relatively low evolutionary rate (7.14×10^-4substitutions/site/year). Furthermore, a common ancestor of HRSV-A and -B diverged approximately 290years ago, while HRSV-B diverged into three clusters for approximately 60years. The genetic similarity of the present strains was very high. Although a maximum of 11 amino acid substitutions were observed in the structural model of the F protein, only one strain possessed an amino acid substitution located within the palivizumab epitope. Four epitopes were predicted, although these did not correspond to the neutralization sites of the F protein including the palivizumab epitope. In addition, five N-glycosylation sites of the present HRSV-B strains were inferred. No positive selection sites were identified; however, many sites were found to be under negative selection. The effective population size of the gene has remained almost constant. On the basis of these results, it can be concluded that the HRSV-B F gene is highly conserved, as is the F protein of HRSV-A. Moreover, our prediction of B-cell epitopes does not show that the palivizumab reaction site may be recognized as an epitope during naturally occurring infections.

Genetic analyses of the fusion protein genes in human parainfluenza virus types 1 and 3 among patients with acute respiratory infections in Eastern Japan from 2011 to 2015

PURPOSE

To genetically explore the fusion protein gene (F) in human parainfluenza virus type 1 (HPIV1) and type 3 (HPIV3) strains, we analysed them in patients with acute respiratory infections in Eastern Japan from 2011 to 2015.

METHODOLOGY

We constructed phylogenetic trees based on the HPIV and HPIV3 F gene using the maximum likelihood method and conducted P-distance and selective pressure analyses. We also predicted the linear epitopes of the protein in the prototype strains. Furthermore, we mapped the amino acid substitutions of the proteins.

RESULTS

Nineteen strains of HPIV1 and 53 strains of HPIV3 were detected among the clinical acute respiratory infection cases. The phylogenetic trees indicated that the HPIV1 and HPIV3 strains were classified into clusters II and III and cluster C, respectively. The P-distance values of the HPIV1 and HPIV3 F genes were <0.03. Two positive selection sites were inferred in the HPIV1 (aa 8 and aa 10), and one positive selection site was inferred in the HPIV3 (aa 108), but over 10 negative selection sites were inferred. Four epitopes were predicted for the HPIV1 prototype strains, while five epitopes were predicted for the HPIV3 prototype strain. A positive selection site (aa 108) or the HPIV3 F protein was involved in the predicted epitope. Additionally, we found that an amino acid substitution (R73K) in the LC76627 HPIV3 strain presumably may affect the resistance to neutralization by antibodies.

CONCLUSION

The F gene of HPIV1 and HPIV3 was relatively well conserved in the eastern part of Japan during the investigation period.

Evolution and Transmission of Respiratory Syncytial Group A (RSV-A) Viruses in Guangdong, China 2008-2015

Respiratory syncytial viruses (RSVs) including subgroups A (RSV-A) and B (RSV-B) are an important cause of acute respiratory tract infections worldwide. RSV-A include major epidemic strains. Fundamental questions concerning the evolution, persistence and transmission of RSV-A are critical for disease control and prevention, yet remain unanswered. In this study, we generated 64 complete G gene sequences of RSV-A strains collected between 2008 and 2015 in Guangdong, China. Phylogenetic analysis was undertaken by incorporating 572 publicly available RSV-A sequences. Current data indicate that genotypes GA1, GA4, and GA5 are endemic with limited epidemic activity. In contrast, the GA2 genotype which likely originated in 1980 has spread rapidly and caused epidemics worldwide. By analyzing GA2 genotype sequences across epidemic seasons within Guangdong, we find that RSV-A epidemics in Guangdong are caused by a combination of virus importation and local persistence, although the magnitude of the latter is likely overestimated due to infrequent sampling in other regions. Our results provide new insights into RSV-A evolution and transmission at global and local scales and highlights the rapid and wide spread of genotype GA2 compared to other genotypes. In order to control RSV transmission and outbreak, both local persistence and external introduction should be taken into account when designing optimal strategies.