Epidemiological Characteristics and Whole-Genome Analysis of Respiratory Syncytial Virus in Jining City from February 2023 to December 2024

Epidemiological Characteristics and Whole-Genome Analysis of Respiratory Syncytial Virus in Jining City from February 2023 to December 2024 | Authorea try { document.documentElement.classList.add('js'); } catch (e) { } var _gaq = _gaq || []; _gaq.push(['_setAccount', 'G-8VDV14Y67G']); _gaq.push(['_trackPageview']); (function() { var ga = document.createElement('script'); ga.type = 'text/javascript'; ga.async = true; ga.src = ('https:' == document.location.protocol ? 'https://ssl' : 'http://www') + '.google-analytics.com/ga.js'; var s = document.getElementsByTagName('script')[0]; s.parentNode.insertBefore(ga, s); })(); Skip to main content Preprints Collections Wiley Open Research IET Open Research Ecological Society of Japan All Collections About About Authorea FAQs Contact Us Quick Search anywhere Search for preprint articles, keywords, etc. Search Search ADVANCED SEARCH SCROLL This is a preprint and has not been peer reviewed. Data may be preliminary. 7 August 2025 V1 Latest version Share on Epidemiological Characteristics and Whole-Genome Analysis of Respiratory Syncytial Virus in Jining City from February 2023 to December 2024 Authors : Huixin Dou , Yongjian Jia , Yajuan Jiang , Xuezhen Shi 0000-0002-6709-0202 , Xiaoyu Wang , Ying Yue , Wei Liu , Bo-Yan Jiao [email protected] , and Feifei Authors Info & Affiliations https://doi.org/10.22541/au.175456810.02664656/v1 189 views 133 downloads Contents Abstract Supplementary Material Information & Authors Metrics & Citations View Options References Figures Tables Media Share Abstract Background: Respiratory syncytial virus (RSV) is a major viral pathogen causing respiratory tract infections in children, often leading to bronchiolitis, pneumonia, and even death. This study aimed to investigate the epidemiological patterns and whole-genome characteristics of RSV circulating in Jining City between February 2023 and December 2024. Methods: From February 2023 to December 2024, a total of 5,042 throat swab samples were collected from influenza-like illness (ILI) cases at two sentinel hospitals in Jining. RSV was detected using reverse transcription quantitative real-time PCR (RT-qPCR). RSV-positive samples were subjected to whole-genome sequencing. Phylogenetic trees were constructed based on the whole genome and G gene sequences, and antigenic variation in viral proteins was analyzed. Results: RSV positivity was 1.98% (100/5042), with higher rates in children under 5 years of age. RSV activity peaked in April–May 2023 and December 2023–January 2024. A total of 29 RSV genomes were sequenced, including 18 RSV-A and 11 RSV-B. RSV-A was dominant during the first peak, and RSV-B during the second. Most RSV-A strains belonged to clade AD.3 and genotype ON1; RSV-B strains clustered into clades B.D.E.1, B.D.4.1.1, and B.D.E.2, all within genotype BA9. Mutations were identified in antigenic epitopes of the G and F proteins, including amino acid substitutions and changes in glycosylation and phosphorylation sites. Conclusions: RSV-A and RSV-B co-circulated in Jining in an alternating pattern, with evidence of ongoing antigenic evolution. Continued RSV surveillance and accelerated vaccine development are essential. 1. Introduction Respiratory syncytial virus (RSV), a membrane-enveloped RNA virus belonging to the genus Orthopneumovirus within the family Paramyxoviridae , exists in both spherical (10–350 nm) and filamentous (60–200 nm) forms(1-3). RSV infection can cause a wide range of clinical symptoms, including fever, rhinorrhea, cough, dyspnea, wheezing, bronchiolitis, pneumonia, and even death(4, 5). Globally, RSV is responsible for an estimated 33 million cases, 3.6 million hospitalizations, and approximately 100,000 deaths annually, representing a significant disease burden(6). RSV encodes several structural and non-structural proteins, including NS1, NS2, N, P, M, SH, G, F, M2-1, M2-2, and L. Among these, the G (attachment) and F (fusion) glycoproteins play critical roles in viral entry and host cell interaction and contain major antigenic epitopes(1, 2). Based on antigenic and genetic variability in the G gene, RSV is classified into two antigenic subtypes: RSV-A and RSV-B(7, 8). Further genotyping based on the G gene sequences has led to the identification of multiple genotypes: GA1–GA7, SAA1, NA1–NA4, and ON1–ON2 for RSV-A; and BA1–BA14, GB1–GB5, SAB1–SAB4, URU1–URU2, NZB1–NZB2, BA-CCA, BA-CCB, BA-C, CBB, and CB1 for RSV-B(2, 9, 10). With the development of next-generation sequencing (NGS), whole-genome sequencing has been increasingly applied to study RSV evolution and clade classification at a higher resolution (11-15). In China, RSV is one of the leading causes of respiratory illness, with significant regional and seasonal variation in prevalence. Typically, RSV circulates during the winter and spring seasons, although outbreaks have also been observed in early summer in certain years. RSV-A and RSV-B often alternate in dominance between epidemic seasons, with RSV-A generally being more prevalent(16, 17). Jining is a major city in northern China with a population exceeding 8.2 million(18, 19). As a regional hub, it plays an important role in infectious disease surveillance and has been the focus of several viral epidemiology studies(20-22). In this study, we analyzed the epidemiological characteristics of RSV in Jining from February 2023 to December 2024 and performed evolutionary analyses on 29 RSV whole-genome sequences. The findings provide important insights into RSV circulation and molecular evolution in this region, contributing valuable data for RSV surveillance and early warning systems. 2. Materials and Methods 2.1 Study Participants From February 2023 to December 2024, throat swab samples were collected within 3 days of symptom onset from patients with influenza-like illness (defined as fever ≥38 °C accompanied by cough or sore throat) at Jining No.1 People’s Hospital and Rencheng District Maternal and Child Health Hospital. Swabs were placed in viral transport medium produced by Youkang Biotechnology Co., Ltd. (Beijing, China; catalog no. MT0901). A total of 5,042 samples were collected, including 4,450 from outpatients and 592 from hospitalized patients. 2.2 Nucleic Acid Detection A 200 μL aliquot of each throat swab sample was used for nucleic acid extraction, performed using the SSNP-9600A automated extraction system and purification kit (catalog no. SDKF60101) from Jiangsu Shuoshi Biotechnology Co., Ltd. RSV RNA was detected using a fluorescent RT-qPCR kit (catalog no. SKY-8237) from Shenzhen Shengkeyuan Biotechnology Co., Ltd. with the following cycling conditions: 50 °C for 15 min; 95 °C for 3 min; followed by 45 cycles of 95 °C for 5 s and 55 °C for 30 s. RSV subtype (A/B) was determined using a dual detection RT-qPCR kit (catalog no. SKY-8240) from the same manufacturer under identical cycling conditions. 2.3 Whole-Genome Sequencing Whole-genome sequencing was performed on 29 RSV-positive samples with Ct values ≤30. Among the 29 cases, 16 were male and 13 female; age distribution was as follows: 17 aged 1–<5 years, 6 aged 5–<15 years, 3 aged 25–<60 years, and 3 aged ≥60 years; 21 samples were collected in 2023 and 8 in 2024. Reverse transcription was performed using the RSV Whole-Genome Capture Kit (catalog no. BK-RSVAB024) from Hangzhou Boyi Technology Co., Ltd. with a 15 μL reaction mix (3 μL RT Super Mix, 12 μL RNA template). Conditions: 25 °C for 10 min, 50 °C for 10 min, 85 °C for 5 min, 4 °C hold. Amplification was conducted with 2.5 μL RT product, 12.5 μL HotStart HiFi Mix, 4 μL primer pool (A1/A2/B1/B2), and 6 μL nuclease-free water. Cycling conditions: 98 °C for 30 s; 5 cycles of 98 °C for 15 s and 63 °C for 5 min; followed by 25–30 cycles of 98 °C for 15 s and 65 °C for 4 min; final extension at 72 °C for 2 min. Purification was performed using AMPure XP beads (Beckman Coulter, catalog no. A63880), and quantification was done using a Qubit 3 fluorometer (Invitrogen). cDNA fragmentation was conducted using the Nextera XT DNA Library Preparation Kit (Illumina, catalog no. 15032350), and adapter ligation was carried out using Nextera XT Index Kit v2 Set A (catalog no. 15052163) with the following cycling conditions: 72 °C for 3 min, 95 °C for 30 s; then 12 cycles of 95 °C for 10 s, 55 °C for 30 s, and 72 °C for 30 s; final extension at 72 °C for 5 min. Libraries were sequenced using the Illumina NextSeq 2000 platform with a 300-cycle NextSeqTM 1000/2000 P1 reagent cartridge (catalog no. 20049920). 2.4 Sequence Alignment and Analysis Sequence assembly and analysis were performed using CLC Genomics Workbench 24 (QIAGEN). RSV whole-genome sequences were aligned and compared using BLAST against reference sequences in the NCBI database. Representative early ON1 and BA9 strains from China—RSV-A isolate BJ/40180 (GenBank: MF614946.1) and RSV-B strain RSVB/BCH-Y/2016 (GenBank: KY924878.1)—as well as recent RSV genomes were downloaded from GenBank. Multiple sequence alignments were performed using MEGA 7.0.14, and maximum-likelihood phylogenetic trees were constructed with 1,000 bootstrap replicates based on the whole genome and G gene sequences. Similarity analysis of RSV genes and encoded proteins was performed using MEGA 7.0.14 and MegAlign (DNASTAR 7.0.1). G protein amino acid variations were analyzed using MEGA 7.0.14. 2.5 Statistical Analysis All statistical analyses were conducted using SPSS version 25.0. Group comparisons were performed using the chi-square (χ²) test. A p-value < 0.05 was considered statistically significant. 3. Results 3.1 Epidemiological Characteristics of RSV Among the 5,042 influenza-like illness (ILI) samples collected from February 2023 to December 2024, 100 tested positive for RSV, yielding a positivity rate of 1.98% (100/5042). The infection rate was highest among children under 1 year and those aged 1 to <5 years, while lower rates were observed in the 15–<25 and 25–<60 age groups. The differences in RSV positivity across age groups were statistically significant (χ² = 108.53, P < 0.001). The RSV positivity rate was higher in males (2.41%, 60/2485) compared to females (1.56%, 40/2557), and the gender difference was statistically significant (χ² = 4.685, P = 0.03). However, no significant difference was observed between outpatients and inpatients (χ² = 0.01, P = 0.93) (Supplementary Table 1, Figure 1). Seasonal variation in RSV activity was evident. Two distinct peaks were observed: one in April–May 2023 and another in December 2023–January 2024, during which RSV positivity rates were significantly higher compared to other months (χ² = 74.68, P < 0.001). RSV-A was the predominant subtype during the April–May 2023 peak, whereas RSV-B dominated in December 2023–January 2024 (Supplementary Table 2, Figure 1). 3.2 Whole-Genome Sequencing of RSV Whole-genome sequencing was successfully performed on 29 RSV-positive samples, yielding complete genomes ranging from 14,590 to 15,221 base pairs in length. Among them, 18 sequences were classified as RSV-A and 11 as RSV-B. Of the 21 sequences obtained in 2023, 16 were RSV-A and 5 were RSV-B. In contrast, the 8 sequences from 2024 included 2 RSV-A and 6 RSV-B strains. All RSV-A strains belonged to the ON1 genotype, and all RSV-B strains were classified as BA9 genotype. 3.3 Sequence Similarity Analysis Among the 18 ON1 genotype RSV-A sequences from Jining, whole-genome nucleotide identity ranged from 98.33% to 99.99%. The nucleotide similarity for the G and F genes ranged from 94.61% to 100% and 98.61% to 100%, respectively. At the amino acid level, the G and F proteins shared 90.13%–100% and 98.95%–100% identity, respectively. Compared with an earlier Chinese ON1 strain, RSV-A isolate BJ/40180 (GenBank: MF614946.1)(23), whole-genome nucleotide identity ranged from 98.84% to 99.07%, G and F gene nucleotide identity ranged from 96.41% to 97.04% and 99.07% to 99.48%, and G and F protein amino acid identity ranged from 93.31% to 96.18% and 99.48% to 99.83%, respectively (Supplementary Table 3 and Figures 2A, 2B, 3C, 3D). For the 11 BA9 genotype RSV-B sequences, whole-genome nucleotide similarity ranged from 98.46% to 99.98%. The G and F gene nucleotide identity ranged from 95.70% to 100% and 98.43% to 100%, while the G and F protein amino acid identity ranged from 91.77% to 100% and 98.95% to 100%, respectively. Compared with the early Chinese BA9 reference strain RSVB/BCH-Y/2016(24), whole-genome nucleotide identity ranged from 97.92% to 98.18%. G and F gene nucleotide identity ranged from 96.02% to 97.69% and 97.97% to 98.49%, while G and F protein amino acid identity ranged from 93.35% to 95.89% and 98.26% to 99.13%, respectively (Supplementary Table 4 and Figures 2E, 2F, 3G, 3H). 3.4 Phylogenetic Analysis Phylogenetic analysis of the RSV-A G gene was performed using 18 sequences obtained from Jining, along with representative sequences of globally circulating RSV-A genotypes from recent years. All 18 Jining strains clustered within the ON1 lineage, confirming their classification as ON1 genotype. Notably, these sequences were distributed across six distinct sublineages within the ON1 clade, suggesting substantial intra-genotypic diversification. This finding indicates that multiple ON1 variants co-circulated in Jining during the study period, contributing to local RSV-A transmission (Figure 3 and Supplementary Table 5). To further characterize RSV-A evolutionary dynamics, whole-genome clade assignments were determined using the Nextstrain platform (https://clades.nextstrain.org/). Of the 18 RSV-A genomes, 17 were classified as clade A.D.3 and one as A.D.5.2. These results were corroborated by a maximum-likelihood phylogenetic tree constructed from whole-genome sequences, which yielded identical clade groupings (Figure 5 and Supplementary Table 5.). Similarly, phylogenetic analysis of the G gene from 11 RSV-B strains revealed that all isolates clustered within the BA9 lineage, indicating a single prevailing genotype in Jining. However, the sequences fell into three distinct subclusters within the BA9 clade, reflecting ongoing genetic diversification among local RSV-B strains (Figure 4 and Supplementary Table 5). Clade classification of the RSV-B whole genomes using the Nextstrain platform showed that four sequences belonged to clade B.D.E.1, two to B.D.E.1.2, three to B.D.4.1.1, and two to B.D.E.2. These assignments were consistent with results from whole-genome phylogenetic reconstruction, which also placed six sequences in B.D.E.1, three in B.D.4.1.1, and two in B.D.E.2 (Figure 6 and Supplementary Table 5). 3.5 Amino Acid Variation Analysis Comparative analysis with the early Chinese ON1 reference strain BJ/40180 revealed that all 18 RSV-A ON1 genomes from Jining harbored amino acid substitutions in multiple proteins, including N, P, M, G, F, M2-1, M2-2, and L. The G glycoprotein, which mediates host cell attachment and monoclonal antibody recognition, contains a conserved CX3C motif responsible for binding to the host receptor CX3CR1 and facilitating viral adhesion(1, 2). This motif remained intact in all 18 Jining ON1 strains. The G and F proteins are the primary targets for virus-neutralizing antibodies(25). Critical antigenic epitopes of the G protein are located at residues 153–197 and the C-terminal 84–85 residues(26-28). Among the ON1 sequences, a single strain carried an N178G mutation within the 153–197 region. However, extensive variability was observed in the C-terminal region, where 28 unique amino acid substitutions were identified (Table 1). The F protein comprises six major antigenic sites: site Ø (aa 62–69, 196–209), site I (aa 380–400), site II (aa 255–276), site III (aa 46–54, 305–310), site IV (aa 427–437), and site V (aa 148–194)(25). Three ON1 strains exhibited an S276N substitution located in site II. Within the p27 domain (aa 110–136), which is critical for membrane fusion(29), observed mutations included R113K, Y117H, N121S (in 9 strains), and T122A (Table 1). The L, P, N, and M2-1 proteins constitute the ribonucleoprotein (RNP) complex, which is essential for viral RNA transcription and replication(30, 31). All known functional residues of M2-1 required for RNA binding (S58, S61, K92, K150, R151, K159) and interaction with P (R126, T130, Y134, S137, Q144, R151, P153, T160, N163) remained fully conserved(32, 33). One ON1 strain carried an E239D substitution in the P protein, located at a critical interface for N protein binding(34), which may affect RNP complex formation. Residues in the N protein mediating binding to P (R259, L272, Y288, F300, F320) and to viral RNA (K170, D175, R184, R185) were conserved in all sequences(35). The L protein contains three enzymatic domains—RNA-dependent RNA polymerase (RdRp), methyltransferase (Mtase), and polyribonucleotidyltransferase (PRNTase). Catalytic residues for RdRp (D700, G810, D811, N812)(36-38), Mtase (K1831, D1936, K1973, E2004)(39), and PRNTase (G1264, T1267, H1338, R1339)(40, 41) were fully conserved (Table 1). The M, P, and F proteins cooperate in the assembly of virus-like particles (VLPs)(42, 43). Critical regions on the P protein required for interaction with M (aa 98–109, 115–125, and 131–151) remained unchanged(44). One ON1 strain exhibited an E49D substitution within the 39–57 region of P, which is necessary for VLP formation and known to require dephosphorylation(45); however, this mutation is unlikely to affect phosphorylation status. Mutations in the M protein included N27D and M73L within the N-terminal region (aa 1–200), which mediates self-interaction with the C-terminal domain(46). Phosphorylation sites S95 and T205, essential for M-mediated nuclear transport(47), were conserved across all ON1 strains (Table 1). For the 11 RSV-B BA9 strains, comparison with the early Chinese reference strain RSVB/BCH-Y/2016 revealed amino acid changes in NS1, NS2, N, P, SH, G, F, M2-1, M2-2, and L proteins. The CX3C motif in the G protein remained conserved. Within the key G protein antigenic region (aa 153–197), N176S (2 strains) and S189N (1 strain) substitutions were detected. Additionally, 25 mutations and a Q311 premature stop codon were observed in the C-terminal region. In the F protein, all 11 strains harbored I206M and Q209R mutations within site Ø, while six strains carried an S389P substitution in site I. All strains exhibited a K191R substitution, and six strains carried S190N mutations in site V. Within the p27 domain, H114Y was consistently observed, while L125P was present in two strains (Table 2). Mutations affecting the M2-1–P interaction interface included an A160T substitution in M2-1 (present in all strains) and a T92A mutation in the P protein (one strain), which may influence RNP complex assembly. All critical residues in the N protein for interaction with P and viral RNA were conserved. The RdRp, Mtase, and PRNTase catalytic residues in L and P–M interaction domains remained unchanged (Table 2). NS1 and NS2 proteins suppress host type I interferon (IFN) responses. The α3 helix region (aa 119–139) of NS1, essential for IFN suppression(48, 49), and NS2 residues 24–26 and 32–33, known to inhibit IFN-β expression(50), were conserved in all BA9 strains (Table 2). 3.6 Analysis of Glycosylation and Phosphorylation Sites N-linked glycosylation site prediction using NetNGlyc 1.0 revealed notable alterations in RSV-encoded proteins among strains circulating in Jining. In the RSV-A ON1 subset, 12 out of 18 strains acquired a novel glycosylation site at asparagine 299 (N299) of the G glycoprotein. This site resides within the immunologically critical 84–85 C-terminal region, and the acquisition of N299 glycosylation may potentially modulate epitope conformation and reduce recognition by host neutralizing antibodies. Among the 11 RSV-B BA9 strains, glycosylation changes were more diverse: two strains acquired an additional N69 glycosylation site in NS2; two gained N114, and five gained N256 glycosylation sites in the G protein. In contrast, loss of glycosylation sites was also observed, including N294 (in two strains), N308 (in all strains), and N167 in the L protein (five strains). Notably, the glycosylation sites N256, N294, and N308 in the G protein are located within or adjacent to the C-terminal antigenic domain and are implicated in modulating viral immune evasion (Table 3). Phosphorylation site prediction using NetPhos 3.1 identified several substitutions that may alter protein phosphorylation patterns, particularly within known antigenic regions. In RSV-A ON1 strains, three exhibited a gain-of-function S284 phosphorylation site in the G protein, while phosphorylation potential at T241 (3 strains), S260, S283, S291, Y297, and S299 (12 strains) was reduced. All these residues map to the immunodominant C-terminal epitope region of the G protein, suggesting a potential role in immune modulation. In the F glycoprotein, nine strains exhibited an additional S121 phosphorylation site within the p27 domain, which is critical for membrane fusion. Moreover, 11 ON1 strains demonstrated loss of T92 phosphorylation in the P protein, a residue located within the P protein’s 90–110 domain, which mediates interaction with M2-1(Table 3). In RSV-B BA9 strains, two sequences showed a novel phosphorylation site at S176 within the G protein’s 153–197 antigenic epitope. Additionally, phosphorylation gain was noted at T252 (7 strains), S289 (1 strain), and Y294 (1 strain), all localized within the C-terminal epitope region. Conversely, a substantial reduction in phosphorylation potential was observed at multiple sites within this region: S189, T253, T254, T264, S265, T274, S275, Y285, T288, S295, T300, and T310—several of which were conserved losses across most or all strains. These extensive modifications within the G protein’s antigenic domain underscore the dynamic evolution of post-translational modification (PTM) landscapes, potentially contributing to altered antigenicity and immune escape. In the F glycoprotein of BA9 strains, a reduction in phosphorylation at S389—a site located within antigenic site I—was observed in six strains, which may impact antibody recognition or membrane fusion efficiency. All 11 strains displayed phosphorylation at Y114 within the p27 domain, further supporting the immunomodulatory relevance of this region. Furthermore, gain of phosphorylation at T160 in the M protein was detected in all BA9 strains. This residue is critical for M2-1 interaction with the P protein, suggesting that this PTM may influence RNP complex formation. Finally, one BA9 strain exhibited loss of T92 phosphorylation in the P protein, again located in the M2-1 interaction interface, and potentially affecting complex stability and viral transcription dynamics (Table 3). 4. Discussion Viral infections remain a leading cause of morbidity and mortality worldwide, underscoring their persistent burden on global health systems and their impact on vulnerable populations(51-58). RSV exhibits strong seasonality in China, with peak activity typically occurring in winter and spring, and occasional outbreaks in early summer in certain years(59, 60). The ON1 and BA9 genotypes currently represent the predominant circulating RSV lineages in China(60, 61). In this study, two distinct RSV epidemic peaks were identified in Jining between February 2023 and December 2024: one in April–May 2023 and another in December 2023–January 2024. These findings are consistent with prior reports indicating RSV circulation during both winter–spring and early summer periods in China, including a nationwide RSV surge in April–May 2023(16). The first peak was primarily caused by RSV-A, while the second was dominated by RSV-B, reflecting the typical genotype alternation pattern previously observed across China(59). The higher positivity rate in children under five aligns with the known age distribution of RSV infection in China(59). Interestingly, while some studies have found no significant difference in RSV positivity between males and females(17), our data revealed a significantly higher positivity rate in males in Jining during the study period, suggesting potential sex-related or strain-specific differences in local RSV transmission. A total of 29 complete RSV genomes were obtained, including 18 RSV-A ON1 and 11 RSV-B BA9 strains, consistent with national trends. The persistence of ON1 and BA9 as dominant genotypes underscores the importance of continued genomic surveillance and the urgent need for genotype-specific vaccine development. Notably, the RSV-A ON1 genomes in Jining were classified into two clades (A.D.3 and A.D.5.2), while the BA9 genomes were distributed among four clades (B.D.E.1, B.D.E.1.2, B.D.4.1.1, and B.D.E.2). These findings highlight the genetic complexity of RSV circulation in Jining, with multiple clade variants co-circulating within the same genotype. Among RSV genes, the G gene is known to be the most variable(28, 62). Consistent with this, our analysis revealed the lowest nucleotide and amino acid similarity for the G gene among both ON1 and BA9 sequences, indicating continuous antigenic drift. Compared to early reference strains (BJ/40180 for ON1 and RSVB/BCH-Y/2016 for BA9), 29 amino acid substitutions were observed in the G protein antigenic regions of ON1 strains, and 27 substitutions plus one premature stop codon (Q311*) were detected in BA9 strains. In contrast, only a single antigenic mutation was found in the F protein of ON1 strains, while five were identified in BA9 strains. These G and F protein variations may alter epitope structure, impair host immune recognition, and contribute to reinfection. The p27 domain of the F protein plays a critical role in membrane fusion. Four amino acid substitutions were detected in the p27 domain of ON1 strains, and two in BA9 strains, which may impact F protein function. In addition, the E239D mutation in the P protein (detected in one ON1 strain) maps to a critical interface for N protein binding. Similarly, the A160T mutation in M2-1 and T92A mutation in P—both detected in BA9 strains—are located within the P–M2-1 interaction domain and may affect ribonucleoprotein (RNP) complex stability. The E49D mutation in the P protein (ON1 strain) was located in the critical VLP-forming region (aa 39–57) and may influence particle assembly. Post-translational modifications (PTMs), including glycosylation and phosphorylation, are essential for RSV protein function and immune evasion. N-linked glycosylation at residues such as N27, N70, N116, N126, and N500 on the F protein can alter immunogenicity and fusion activity(63-65). The G protein contains multiple glycosylation sites that contribute to immune escape(66-69). In our study, N299 glycosylation was newly acquired in 12 ON1 G proteins, and N256, N294, and N308 glycosylation changes were found in BA9 strains—all located within G protein antigenic regions—suggesting a role in immune modulation and escape. Phosphorylation also plays a pivotal role in regulating RSV replication and host interaction. The P protein contains more than 40 phosphorylation sites, including S54, T105, T108, S116, S117, S119, S143, S156, T160, T188, S203, T210, S232, and S237(70-73). S54 phosphorylation regulates P–M binding and membrane translocation of RNP complexes(70, 74), while T108 is required for P–M2-1 complex formation(75). Phosphorylation at S143, T160, and T210 modulates transcriptional activity; S116, S117, S119, S232, and S237 are required for M2-2-mediated inhibition of viral RNA synthesis(76); and sites such as T105, T188, S203, and T210 are critical for viral replication(77). The phosphorylation status of M protein residues S95 and T205 influences oligomerization and viral production(47, 78). M2-1 phosphorylation at S58 and S61 regulates RNA binding and transcriptional activity(32, 79). Y88 phosphorylation in the N protein affects its oligomerization and polymerase activity(80). Phosphorylation can also alter antigenicity(81, 82), and our data show extensive phosphorylation changes within the G and F protein epitopes of ON1 and BA9 strains, potentially contributing to immune evasion. Moreover, phosphorylation changes at P and M2-1 interaction sites may impact RNP assembly and viral replication. In conclusion, this study provides a comprehensive analysis of RSV epidemiology and genomic variation in Jining, China, from February 2023 to December 2024. The findings offer valuable insights into genotype distribution, clade diversity, and evolving antigenic features of RSV circulating in northern China. These results emphasize the need for continuous genomic surveillance, timely public health response, and the acceleration of genotype-specific RSV vaccine development to mitigate the impact of future outbreaks. Availability of Data and Materials Data is provided within the supplementary information files. Consent for publication Not applicable. Funding This work was supported by Science and Technology Development Funds for Shandong Medical and Health (202112060725, 202412060875, 202412060648), Shandong Province Medical Staff Science and Technology Innovation Program (SDYWZGKCJHLH202467) Acknowledgements Not applicable. Author Contributions WL、BJ conceived and designed the experiments. HD、YJ、YJ performed the experiments. XW and YY responsible for epidemiological investigation and sample collection，HD、YJ、YJ、FH、XS、YY、WL、BJ analyzed the data. HD、YJ、YJ、FH、XW、BJ interpreted the data. WL、BJ wrote the manuscript. All authors contributed to the article and approved the submitted version. Competing Interests The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. \received DD MMMM YYYY \acceptedDD MMMM YYYY Reference 1. Duan Y, Liu Z, Zang N, Cong B, Shi Y, Xu L, Jiang M, Wang P, Zou J, Zhang H, Feng Z, Feng L, Ren L, Liu E, Li Y, Zhang Y, Xie Z. 2024. Landscape of respiratory syncytial virus. Chin Med J (Engl) 137:2953-2978.2. Sanz-Munoz I, Sanchez-de Prada L, Castrodeza-Sanz J, Eiros JM. 2024. Microbiological and epidemiological features of respiratory syncytial virus. Rev Esp Quimioter 37:209-220.3. Jeffree CE, Rixon HW, Brown G, Aitken J, Sugrue RJ. 2003. Distribution of the attachment (G) glycoprotein and GM1 within the envelope of mature respiratory syncytial virus filaments revealed using field emission scanning electron microscopy. Virology 306:254-67.4. Chaiut W, Sapbamrer R, Dacha S, Sudjaritruk T, Parwati I, Sumarpo A, Malasao R. 2023. Characteristics of Respiratory Syncytial Virus Infection in Hospitalized Children Before and During the COVID-19 Pandemic in Thailand. J Prev Med Public Health 56:212-220.5. Li Y, Zhai Y, Lin Y, Lu C, He Z, Wu S, Yang C, Zhou R, Liu W, Chen D. 2023. Epidemiology of respiratory syncytial virus in hospitalized children with community-acquired pneumonia in Guangzhou: a 10-year study. J Thorac Dis 15:967-976.6. Li Y, Wang X, Blau DM, Caballero MT, Feikin DR, Gill CJ, Madhi SA, Omer SB, Simoes EAF, Campbell H, Pariente AB, Bardach D, Bassat Q, Casalegno JS, Chakhunashvili G, Crawford N, Danilenko D, Do LAH, Echavarria M, Gentile A, Gordon A, Heikkinen T, Huang QS, Jullien S, Krishnan A, Lopez EL, Markic J, Mira-Iglesias A, Moore HC, Moyes J, Mwananyanda L, Nokes DJ, Noordeen F, Obodai E, Palani N, Romero C, Salimi V, Satav A, Seo E, Shchomak Z, Singleton R, Stolyarov K, Stoszek SK, von Gottberg A, Wurzel D, Yoshida LM, Yung CF, Zar HJ, Respiratory Virus Global Epidemiology N, Nair H, et al. 2022. Global, regional, and national disease burden estimates of acute lower respiratory infections due to respiratory syncytial virus in children younger than 5 years in 2019: a systematic analysis. Lancet 399:2047-2064.7. Anderson LJ, Hierholzer JC, Tsou C, Hendry RM, Fernie BF, Stone Y, McIntosh K. 1985. Antigenic characterization of respiratory syncytial virus strains with monoclonal antibodies. J Infect Dis 151:626-33.8. Sullender WM. 2000. Respiratory syncytial virus genetic and antigenic diversity. Clin Microbiol Rev 13:1-15, table of contents.9. Peret TC, Hall CB, Schnabel KC, Golub JA, Anderson LJ. 1998. Circulation patterns of genetically distinct group A and B strains of human respiratory syncytial virus in a community. J Gen Virol 79 ( Pt 9):2221-9.10. Bashir U, Alam MM, Sadia H, Zaidi SS, Kazi BM. 2015. Correction: Molecular Characterization of Circulating Respiratory Syncytial Virus (RSV) Genotypes in Gilgit Baltistan Province of Pakistan during 2011-2012 Winter Season. PLoS One 10:e0145599.11. Pangesti KNA, Abd El Ghany M, Walsh MG, Kesson AM, Hill-Cawthorne GA. 2018. Molecular epidemiology of respiratory syncytial virus. Rev Med Virol 28.12. Dong X, Edwards S, Deng YM, Dapat C, Hirankitti A, Wordsworth R, Whitney P, Baird R, Freeman K, Daley AJ, Barr IG. 2025. An Improved Rapid and Sensitive Long Amplicon Method for Nanopore-Based RSV Whole-Genome Sequencing. Influenza Other Respir Viruses 19:e70106.13. Wei X, Wang L, Li M, Qi J, Kang L, Hu G, Gong C, Wang C, Wang Y, Huang F, Gao GF. 2024. Novel imported clades accelerated the RSV surge in Beijing, China, 2023-2024. J Infect 89:106321.14. Sondlane H, Ogunbayo A, Donato C, Mogotsi M, Esona M, Hallbauer U, Bester P, Goedhals D, Nyaga M. 2024. Whole genome molecular analysis of respiratory syncytial virus pre and during the COVID-19 pandemic in Free State province, South Africa. Virus Res 347:199421.15. Tramuto F, Maida CM, Randazzo G, Guzzetta V, Santino A, Li Muli R, Costantino C, Graziano G, Amodio E, Mazzucco W, Vitale F. 2024. Whole-Genome Sequencing and Genetic Diversity of Human Respiratory Syncytial Virus in Patients with Influenza-like Illness in Sicily (Italy) from 2017 to 2023. Viruses 16.16. Li BS, Shi Y, Geng MJ, Guo YQ, Lin F, Zhang YP, Li ZJ, Wang LP. 2024. [Epidemiological characteristics of human respiratory syncytial virus infection in febrile respiratory syndrome cases in China from 2010 to 2020]. Zhonghua Yu Fang Yi Xue Za Zhi 58:1311-1317.17. Cui AL, Xia BC, Zhu Z, Xie ZB, Sun LW, Xu J, Xu J, Li Z, Zhao LQ, Long XR, Yu DS, Zhu B, Zhang F, Mu M, Xie H, Cai L, Zhu Y, Tian XL, Wang B, Gao ZG, Liu XQ, Ren BZ, Han GY, Hu KX, Zhang Y. 2024. [Epidemiological characteristics of human respiratory syncytial virus among acute respiratory infection cases in 16 provinces of China from 2009 to 2023]. Zhonghua Yu Fang Yi Xue Za Zhi 58:945-951.18. Yin Q, Liu W, Jiang Y, Feng Q, Wang X, Dou H, Liu Z, He F, Fan Y, Jiao B, Jiao B. 2024. Comprehensive genomic analysis of the SARS-CoV-2 Omicron variant BA.2.76 in Jining City, China, 2022. BMC Genomics 25:378.19. Jiang Y, Dou H, Wang X, Song T, Jia Y, Yue Y, Li L, He F, Kong L, Wu Z, Huang X, Liang Y, Jiao B, Jiao B. 2024. Analysis of seasonal H3N2 influenza virus epidemic characteristics and whole genome features in Jining City from 2018 to 2023. J Med Virol 96:e29846.20. Wu C, Fu Z, Xie C, Zhao J, He F, Jiao B, Jiao B. 2024. Epidemiological Characteristics and Genotypic Features of Rotavirus and Norovirus in Jining City, 2021-2022. Viruses 16.21. Fu Z, Wu C, Yu Y, Zhao J, Dong Y, Liu W, Dou H, Shi X, Cai C, Jiao B, Liu T, Jiao B. 2025. Genomic characterization of norovirus and Sapovirus from a diarrhea outbreak in a school linked to heavy rainfall. Front Microbiol 16:1570161.22. Dou H, Chen C, Song T, Sun X, He F, Jia Y, Wang X, Jiang Y, Yue Y, Huang S, Yan S, Jiao B, Jiao B. 2025. Epidemiological characteristics and genomic analysis of respiratory adenovirus in Jining City from February 2023 to July 2024. BMC Genomics 26:369.23. Cui G, Zhu R, Deng J, Qian Y. 2017. Two Complete Sequences for Genomes of Human Respiratory Syncytial Virus Genotype ON1 Field Strains from Beijing, China. Genome Announc 5.24. Xu L, Gao H, Zeng J, Liu J, Lu C, Guan X, Qian S, Xie Z. 2018. A fatal case associated with respiratory syncytial virus infection in a young child. BMC Infect Dis 18:217.25. Huang J, Diaz D, Mousa JJ. 2019. Antibody Epitopes of Pneumovirus Fusion Proteins. Front Immunol 10:2778.26. Jones HG, Ritschel T, Pascual G, Brakenhoff JPJ, Keogh E, Furmanova-Hollenstein P, Lanckacker E, Wadia JS, Gilman MSA, Williamson RA, Roymans D, van ’t Wout AB, Langedijk JP, McLellan JS. 2018. Structural basis for recognition of the central conserved region of RSV G by neutralizing human antibodies. PLoS Pathog 14:e1006935.27. Cane PA, Thomas HM, Simpson AF, Evans JE, Hart CA, Pringle CR. 1996. Analysis of the human serological immune response to a variable region of the attachment (G) protein of respiratory syncytial virus during primary infection. J Med Virol 48:253-61.28. Anderson LJ, Jadhao SJ, Paden CR, Tong S. 2021. Functional Features of the Respiratory Syncytial Virus G Protein. Viruses 13.29. Neal HE, Barrett CT, Edmonds K, Moncman CL, Dutch RE. 2024. Examination of respiratory syncytial virus fusion protein proteolytic processing and roles of the P27 domain. J Virol 98:e0163924.30. Cao D, Gao Y, Liang B. 2021. Structural Insights into the Respiratory Syncytial Virus RNA Synthesis Complexes. Viruses 13.31. Collins PL, Hill MG, Cristina J, Grosfeld H. 1996. Transcription elongation factor of respiratory syncytial virus, a nonsegmented negative-strand RNA virus. Proc Natl Acad Sci U S A 93:81-5.32. Tanner SJ, Ariza A, Richard CA, Kyle HF, Dods RL, Blondot ML, Wu W, Trincao J, Trinh CH, Hiscox JA, Carroll MW, Silman NJ, Eleouet JF, Edwards TA, Barr JN. 2014. Crystal structure of the essential transcription antiterminator M2-1 protein of human respiratory syncytial virus and implications of its phosphorylation. Proc Natl Acad Sci U S A 111:1580-5.33. Gao Y, Cao D, Pawnikar S, John KP, Ahn HM, Hill S, Ha JM, Parikh P, Ogilvie C, Swain A, Yang A, Bell A, Salazar A, Miao Y, Liang B. 2020. Structure of the Human Respiratory Syncytial Virus M2-1 Protein in Complex with a Short Positive-Sense Gene-End RNA. Structure 28:979-990 e4.34. Tran TL, Castagne N, Bhella D, Varela PF, Bernard J, Chilmonczyk S, Berkenkamp S, Benhamo V, Grznarova K, Grosclaude J, Nespoulos C, Rey FA, Eleouet JF. 2007. The nine C-terminal amino acids of the respiratory syncytial virus protein P are necessary and sufficient for binding to ribonucleoprotein complexes in which six ribonucleotides are contacted per N protein protomer. J Gen Virol 88:196-206.35. Esneau C, Raynal B, Roblin P, Brule S, Richard CA, Fix J, Eleouet JF, Galloux M. 2019. Biochemical characterization of the respiratory syncytial virus N(0)-P complex in solution. J Biol Chem 294:3647-3660.36. Shirai T, Sada M, Saraya T, Kurai D, Sunagawa S, Ishii H, Kimura H. 2022. Detailed Analyses of Molecular Interactions between Favipiravir and RNA Viruses In Silico. Viruses 14.37. Noton SL, Aljabr W, Hiscox JA, Matthews DA, Fearns R. 2014. Factors affecting de novo RNA synthesis and back-priming by the respiratory syncytial virus polymerase. Virology 462-463:318-27.38. Sutto-Ortiz P, Eleouet JF, Ferron F, Decroly E. 2023. Biochemistry of the Respiratory Syncytial Virus L Protein Embedding RNA Polymerase and Capping Activities. Viruses 15.39. Sutto-Ortiz P, Tcherniuk S, Ysebaert N, Abeywickrema P, Noel M, Decombe A, Debart F, Vasseur JJ, Canard B, Roymans D, Rigaux P, Eleouet JF, Decroly E. 2021. The methyltransferase domain of the Respiratory Syncytial Virus L protein catalyzes cap N7 and 2’-O-methylation. PLoS Pathog 17:e1009562.40. Gilman MSA, Liu C, Fung A, Behera I, Jordan P, Rigaux P, Ysebaert N, Tcherniuk S, Sourimant J, Eleouet JF, Sutto-Ortiz P, Decroly E, Roymans D, Jin Z, McLellan JS. 2019. Structure of the Respiratory Syncytial Virus Polymerase Complex. Cell 179:193-204 e14.41. Liang B, Li Z, Jenni S, Rahmeh AA, Morin BM, Grant T, Grigorieff N, Harrison SC, Whelan SPJ. 2015. Structure of the L Protein of Vesicular Stomatitis Virus from Electron Cryomicroscopy. Cell 162:314-327.42. Meshram CD, Baviskar PS, Ognibene CM, Oomens AGP. 2016. The Respiratory Syncytial Virus Phosphoprotein, Matrix Protein, and Fusion Protein Carboxy-Terminal Domain Drive Efficient Filamentous Virus-Like Particle Formation. J Virol 90:10612-10628.43. Mandviwala AS, Munje AK, Huckriede ALW, Arankalle VA, Patil HP. 2025. Immunogenicity evaluation of respiratory syncytial virus prefusogenic-F based virus-like-particles consisting of G and M proteins in mice. Vaccine 56:127203.44. Bajorek M, Galloux M, Richard CA, Szekely O, Rosenzweig R, Sizun C, Eleouet JF. 2021. Tetramerization of Phosphoprotein is Essential for Respiratory Syncytial Virus Budding while its N Terminal Region Mediates Direct Interactions with the Matrix Protein. J Virol 95.45. Meshram CD, Oomens AGP. 2019. Identification of a human respiratory syncytial virus phosphoprotein domain required for virus-like-particle formation. Virology 532:48-54.46. Trevisan M, Di Antonio V, Radeghieri A, Palu G, Ghildyal R, Alvisi G. 2018. Molecular Requirements for Self-Interaction of the Respiratory Syncytial Virus Matrix Protein in Living Mammalian Cells. Viruses 10.47. Ghildyal R, Teng MN, Tran KC, Mills J, Casarotto MG, Bardin PG, Jans DA. 2022. Nuclear Transport of Respiratory Syncytial Virus Matrix Protein Is Regulated by Dual Phosphorylation Sites. Int J Mol Sci 23.48. Merritt TN, Pei J, Leung DW. 2023. Pathogenicity and virulence of human respiratory syncytial virus: Multifunctional nonstructural proteins NS1 and NS2. Virulence doi:10.1080/21505594.2023.2283897:2283897.49. Chatterjee S, Luthra P, Esaulova E, Agapov E, Yen BC, Borek DM, Edwards MR, Mittal A, Jordan DS, Ramanan P, Moore ML, Pappu RV, Holtzman MJ, Artyomov MN, Basler CF, Amarasinghe GK, Leung DW. 2017. Structural basis for human respiratory syncytial virus NS1-mediated modulation of host responses. Nat Microbiol 2:17101.50. Pei J, Wagner ND, Zou AJ, Chatterjee S, Borek D, Cole AR, Kim PJ, Basler CF, Otwinowski Z, Gross ML, Amarasinghe GK, Leung DW. 2021. Structural basis for IFN antagonism by human respiratory syncytial virus nonstructural protein 2. Proc Natl Acad Sci U S A 118.51. Chen X, Ye H, Li S, Jiao B, Wu J, Zeng P, Chen L. 2017. Severe fever with thrombocytopenia syndrome virus inhibits exogenous Type I IFN signaling pathway through its NSs invitro. PLoS One 12:e0172744.52. Duan X, Li S, Holmes JA, Tu Z, Li Y, Cai D, Liu X, Li W, Yang C, Jiao B, Schaefer EA, Fusco DN, Salloum S, Chen L, Lin W, Chung RT. 2018. MicroRNA 130a Regulates both Hepatitis C Virus and Hepatitis B Virus Replication through a Central Metabolic Pathway. J Virol 92.53. Chen Y, Jiao B, Yao M, Shi X, Zheng Z, Li S, Chen L. 2017. ISG12a inhibits HCV replication and potentiates the anti-HCV activity of IFN-alpha through activation of the Jak/STAT signaling pathway independent of autophagy and apoptosis. Virus Res 227:231-239.54. Shi X, Jiao B, Chen Y, Li S, Chen L. 2017. MxA is a positive regulator of type I IFN signaling in HCV infection. J Med Virol 89:2173-2180.55. Yuan Y, Jiao B, Qu L, Yang D, Liu R. 2023. The development of COVID-19 treatment. Front Immunol 14:1125246.56. Li Y, Li S, Duan X, Chen Y, Jiao B, Ye H, Yao M, Chen L. 2016. Interferon-Stimulated Gene 15 Conjugation Stimulates Hepatitis B Virus Production Independent of Type I Interferon Signaling Pathway In Vitro. Mediators Inflamm 2016:7417648.57. Jiao B, Shi X, Chen Y, Ye H, Yao M, Hong W, Li S, Duan X, Li Y, Wang Y, Chen L. 2017. Insulin receptor substrate-4 interacts with ubiquitin-specific protease 18 to activate the Jak/STAT signaling pathway. Oncotarget 8:105923-105935.58. Li H, Liu T, Shi X, Du H, Cai C, Yang D, Qu L, Dou H, Jiao B, Jiao B. 2025. Mechanisms and therapeutic potential of pharmacological agents targeting inflammasomes. Biomed Pharmacother 189:118164.59. Luo HJ, Huang XB, Zhong HL, Ye CX, Tan X, Zhou K, Yuan L, Zhang SF, Zhu X, Lin CJ, Wang WJ, Xu L, Cao KY. 2020. Epidemiological characteristics and phylogenic analysis of human respiratory syncytial virus in patients with respiratory infections during 2011-2016 in southern China. Int J Infect Dis 90:5-17.60. Jiang M, Wang F, Xu Y, Han Z, Guo Q, Sun Y, Zhu R, Chen D, Zhou Y, Qu D, Cao L, Qu W, Lu F, Zhao L. 2025. The Changed Endemic Pattern of Respiratory Syncytial Virus Infections in Pediatric Patients With the Relaxation of NPIs in Beijing, China, 2022-2023. J Med Virol 97:e70348.61. Li F, Zhu Y, Li Q, Guan X, Zhang H, Li C, Zhang M, Li L, Fu Y, Duan Y, Huang L, Xie Z, Chen X. 2024. Genomic characteristics of human respiratory syncytial virus from children in China during 2017-2020. Arch Virol 169:219.62. Johnson PR, Jr., Olmsted RA, Prince GA, Murphy BR, Alling DW, Walsh EE, Collins PL. 1987. Antigenic relatedness between glycoproteins of human respiratory syncytial virus subgroups A and B: evaluation of the contributions of F and G glycoproteins to immunity. J Virol 61:3163-6.63. Tian D, Battles MB, Moin SM, Chen M, Modjarrad K, Kumar A, Kanekiyo M, Graepel KW, Taher NM, Hotard AL, Moore ML, Zhao M, Zheng ZZ, Xia NS, McLellan JS, Graham BS. 2017. Structural basis of respiratory syncytial virus subtype-dependent neutralization by an antibody targeting the fusion glycoprotein. Nat Commun 8:1877.64. Jacobs L, Leemans A, Stobbelaar K, Fransen A, Cos P, Delputte P. 2024. Evaluating the Impact of N-Glycan Sequon Removal in the p27 Peptide on RSV F Protein Immunogenicity and Functionality. Viruses 16.65. Leemans A, Boeren M, Van der Gucht W, Martinet W, Caljon G, Maes L, Cos P, Delputte P. 2019. Characterization of the role of N-glycosylation sites in the respiratory syncytial virus fusion protein in virus replication, syncytium formation and antigenicity. Virus Res 266:58-68.66. Aziz F, Farooqui N, Abbas T, Javaid M, Rafaqat W, Zhamalbekova A, Ali SA, Ali S, Abid SH. 2024. Phylogenetic and phylodynamic analysis of respiratory syncytial virus strains circulating in children less than five years of age in Karachi-Pakistan. Infect Genet Evol 126:105694.67. Gorman JJ, McKimm-Breschkin JL, Norton RS, Barnham KJ. 2001. Antiviral activity and structural characteristics of the nonglycosylated central subdomain of human respiratory syncytial virus attachment (G) glycoprotein. J Biol Chem 276:38988-94.68. Rawling J, Melero JA. 2007. The use of monoclonal antibodies and lectins to identify changes in viral glycoproteins that are influenced by glycosylation: the case of human respiratory syncytial virus attachment (G) glycoprotein. Methods Mol Biol 379:109-25.69. Garcia-Beato R, Martinez I, Franci C, Real FX, Garcia-Barreno B, Melero JA. 1996. Host cell effect upon glycosylation and antigenicity of human respiratory syncytial virus G glycoprotein. Virology 221:301-9.70. Asenjo A, Rodriguez L, Villanueva N. 2005. Determination of phosphorylated residues from human respiratory syncytial virus P protein that are dynamically dephosphorylated by cellular phosphatases: a possible role for serine 54. J Gen Virol 86:1109-1120.71. Mazumder B, Adhikary G, Barik S. 1994. Bacterial expression of human respiratory syncytial viral phosphoprotein P and identification of Ser237 as the site of phosphorylation by cellular casein kinase II. Virology 205:93-103.72. Asenjo A, Mendieta J, Gomez-Puertas P, Villanueva N. 2008. Residues in human respiratory syncytial virus P protein that are essential for its activity on RNA viral synthesis. Virus Res 132:160-73.73. Navarro J, Lopez-Otin C, Villanueva N. 1991. Location of phosphorylated residues in human respiratory syncytial virus phosphoprotein. J Gen Virol 72 ( Pt 6):1455-9.74. Asenjo A, Gonzalez-Armas JC, Villanueva N. 2008. Phosphorylation of human respiratory syncytial virus P protein at serine 54 regulates viral uncoating. Virology 380:26-33.75. Asenjo A, Calvo E, Villanueva N. 2006. Phosphorylation of human respiratory syncytial virus P protein at threonine 108 controls its interaction with the M2-1 protein in the viral RNA polymerase complex. J Gen Virol 87:3637-3642.76. Asenjo A, Villanueva N. 2016. Phosphorylation of the human respiratory syncytial virus P protein mediates M2-2 regulation of viral RNA synthesis, a process that involves two P proteins. Virus Res 211:117-25.77. Beavis AC, Tran KC, Barrozo ER, Phan SI, Teng MN, He B. 2021. Respiratory Syncytial Virus Phosphoprotein Residue S156 Plays a Role in Regulating Genome Transcription and Replication. J Virol 95:e0120621.78. Bajorek M, Caly L, Tran KC, Maertens GN, Tripp RA, Bacharach E, Teng MN, Ghildyal R, Jans DA. 2014. The Thr205 phosphorylation site within respiratory syncytial virus matrix (M) protein modulates M oligomerization and virus production. J Virol 88:6380-93.79. Cartee TL, Wertz GW. 2001. Respiratory syncytial virus M2-1 protein requires phosphorylation for efficient function and binds viral RNA during infection. J Virol 75:12188-97.80. Basse V, Wang Y, Rodrigues-Machado C, Henry C, Richard CA, Leyrat C, Galloux M. 2025. Regulation of respiratory syncytial virus nucleoprotein oligomerization by phosphorylation. J Biol Chem 301:108256.81. Min S, Mohallem R, Aryal UK, Kinzer-Ursem TL, Rochet JC. 2023. Effects of Neighboring Phosphorylation Events on the Affinities of pT181-Tau Antibodies. Anal Chem 95:18241-18248.82. Scherf M, Koy C, Rower C, Neamtu A, Glocker MO. 2023. Characterization of Phosphorylation-Dependent Antibody Binding to Cancer-Mutated Linkers of C(2)H(2) Zinc Finger Proteins by Intact Transition Epitope Mapping-Thermodynamic Weak-Force Order Analysis. J Am Soc Mass Spectrom 34:171-181. Figure Legend Figure 1. Epidemiological characteristics of RSV in Jining, China (February 2023 – December 2024). (A) RSV positivity rates across different age groups. The x-axis represents age groups, and the y-axis indicates positivity rates (%). (B) RSV positivity rates by sex. The x-axis represents sex (male vs. female), and the y-axis indicates positivity rates (%). (C) RSV positivity rates among inpatients and outpatients. The x-axis represents patient setting (inpatient vs. outpatient), and the y-axis indicates positivity rates (%). (D) Monthly distribution of RSV positivity rates. The x-axis represents calendar months, and the y-axis indicates positivity rates (%). Figure 2. Similarity analysis of RSV genomes and encoded proteins. (A) Nucleotide similarity among the 18 RSV-A ON1 genomes from Jining. (B) Amino acid similarity of encoded proteins among the 18 RSV-A ON1 sequences from Jining. (C) Nucleotide similarity between the 18 RSV-A ON1 genomes from Jining and the early Chinese ON1 reference strain BJ/40180. (D) Amino acid similarity of encoded proteins between the 18 RSV-A ON1 sequences from Jining and BJ/40180. (E) Nucleotide similarity among the 11 RSV-B BA9 genomes from Jining. (F) Amino acid similarity of encoded proteins among the 11 RSV-B BA9 sequences from Jining. (G) Nucleotide similarity between the 11 RSV-B BA9 genomes from Jining and the early Chinese BA9 reference strain RSVB/BCH-Y/2016. (H) Amino acid similarity of encoded proteins between the 11 RSV-B BA9 sequences from Jining and RSVB/BCH-Y/2016. Figure 3. Phylogenetic analysis of RSV-A G gene sequences. The phylogenetic tree was constructed using MEGA version 7.0.14 with the maximum-likelihood method and 1,000 bootstrap replicates. Brown branches represent RSV-A G gene sequences from Jining collected in 2023; green branches represent sequences collected in 2024. Figure 4. Phylogenetic analysis of complete RSV-A genomes. The phylogenetic tree was constructed using MEGA version 7.0.14 with the maximum-likelihood method and 1,000 bootstrap replicates. Brown branches represent RSV-A complete genome sequences from Jining collected in 2023; green branches represent sequences collected in 2024. Figure 5. Phylogenetic analysis of RSV-B G gene sequences. The phylogenetic tree was constructed using MEGA version 7.0.14 with the maximum-likelihood method and 1,000 bootstrap replicates. Brown branches represent RSV-B G gene sequences from Jining collected in 2023; green branches represent sequences collected in 2024. Figure 6. Phylogenetic analysis of complete RSV-B genomes. The phylogenetic tree was constructed using MEGA version 7.0.14 with the maximum-likelihood method and 1,000 bootstrap replicates. Brown branches represent RSV-B complete genome sequences from Jining collected in 2023; green branches represent sequences collected in 2024. Supplementary Material File (table 1.docx) Download 17.83 KB File (table 2.docx) Download 17.19 KB File (table 3.docx) Download 19.24 KB Information & Authors Information Version history V1 Version 1 07 August 2025 Copyright This work is licensed under a Non Exclusive No Reuse License. Keywords antigenic variation epidemiology genetics respiratory syncytial virus virus classification Authors Affiliations Huixin Dou Centers for Disease Control and Prevention Global Health China View all articles by this author Yongjian Jia Centers for Disease Control and Prevention Global Health China View all articles by this author Yajuan Jiang Centers for Disease Control and Prevention Global Health China View all articles by this author Xuezhen Shi 0000-0002-6709-0202 Baoding No 1 Central Hospital Department of Pathology View all articles by this author Xiaoyu Wang Centers for Disease Control and Prevention Global Health China View all articles by this author Ying Yue Centers for Disease Control and Prevention Global Health China View all articles by this author Wei Liu Centers for Disease Control and Prevention Global Health China View all articles by this author Bo-Yan Jiao [email protected] Centers for Disease Control and Prevention Global Health China View all articles by this author Feifei Northern Arizona University View all articles by this author Metrics & Citations Metrics Article Usage 189 views 133 downloads .FvxKWukQNSOunydq8rnd { width: 100px; } Citations Download citation Huixin Dou, Yongjian Jia, Yajuan Jiang, et al. Epidemiological Characteristics and Whole-Genome Analysis of Respiratory Syncytial Virus in Jining City from February 2023 to December 2024. Authorea . 07 August 2025. DOI: https://doi.org/10.22541/au.175456810.02664656/v1 If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download. For more information or tips please see 'Downloading to a citation manager' in the Help menu . Format Please select one from the list RIS (ProCite, Reference Manager) EndNote BibTex Medlars RefWorks Direct import Tips for downloading citations document.getElementById('citMgrHelpLink').addEventListener('click', function() { popupHelp(this.href); return false; }); $(".js__slcInclude").on("change", function(e){ if ($(this).val() == 'refworks') $('#direct').prop("checked", false); $('#direct').prop("disabled", ($(this).val() == 'refworks')); }); View Options View options PDF View PDF Figures Tables Media Share Share Share article link Copy Link Copied! Copying failed. Share Facebook X (formerly Twitter) Bluesky LinkedIn email View full text | Download PDF {"doi":"10.22541/au.175456810.02664656/v1","type":"Article"} Now Reading: Share Figures Tables Close figure viewer Back to article Figure title goes here Change zoom level Go to figure location within the article Download figure Toggle share panel Toggle share panel Share Toggle information panel Toggle information panel Go to previous graphic Go to next graphic Go to previous table Go to next table All figures All tables View all material View all material xrefBack.goTo xrefBack.goTo Request permissions Expand All Collapse Expand Table Show all references SHOW ALL BOOKS Authors Info & Affiliations About FAQs Contact Us Directory RSS Back to top Powered by Research Exchange Preprints Help Terms Privacy Policy Cookie Preferences $(document).ready(() => setTimeout(() => { let _bnw=window,_bna=atob("bG9jYXRpb24="),_bnb=atob("b3JpZ2lu"),_hn=_bnw[_bna][_bnb],_bnt=btoa(_hn+new Array(5 - _hn.length % 4).join(" ")); $.get("/resource/lodash?t="+_bnt); },4000)); (function(){function c(){var b=a.contentDocument||a.contentWindow.document;if(b){var d=b.createElement('script');d.innerHTML="window.__CF$cv$params={r:'9ffe98462cb43fe2',t:'MTc3OTQ4MTk5NQ=='};var a=document.createElement('script');a.src='/cdn-cgi/challenge-platform/scripts/jsd/main.js';document.getElementsByTagName('head')[0].appendChild(a);";b.getElementsByTagName('head')[0].appendChild(d)}}if(document.body){var a=document.createElement('iframe');a.height=1;a.width=1;a.style.position='absolute';a.style.top=0;a.style.left=0;a.style.border='none';a.style.visibility='hidden';document.body.appendChild(a);if('loading'!==document.readyState)c();else if(window.addEventListener)document.addEventListener('DOMContentLoaded',c);else{var e=document.onreadystatechange||function(){};document.onreadystatechange=function(b){e(b);'loading'!==document.readyState&&(document.onreadystatechange=e,c())}}}})();