INTRODUCTION
Human respiratory syncytial virus (RSV) is the second most frequent cause of respiratory infections, and the main cause of bronchiolitis and pneumonia in infants and young children [Reference Borchers1]. The incidence of RSV is known to be age-dependent, with an estimated 60% of children infected with RSV in their first year of life [Reference Glezen2]. Almost all children are infected by the time they turn 2 years of age [Reference Sorce3], whereas fewer infections are detected in older children and adults [Reference Hogan4]. Immunity to RSV following recovery from infection is partial and temporary, thus reinfection throughout early childhood is common [Reference Meng5]. Due to the lack of an effective vaccine, substantial health and economic costs are associated with RSV epidemics due to hospitalizations and treatment of RSV infection-related illnesses [Reference Slovic6].
RSV is a non-segmented, negative-sense single-strand RNA enveloped virus belonging to the family Paramyxoviridae [Reference Johnson7]. There are two major RSV subgroups (A and B) categorized on the basis of antigenic and genetic variability. Each subgroup is further categorized into genotypes based on the nucleotide sequence variation within the second hypervariable region of heavily glycosylated G glycoprotein (G protein). There are 12 genotypes for RSV-A (GA1–7, SAA1, NA1–2, and ON1–2) and 20 genotypes for RSV-B (GB1–4, SAB1–4, URU1–2, and BA1–10) [Reference Cui8, Reference Hirano9]. Depending on the genotypes, there are differences in the attachment G glycoprotein which interacts with host cell receptors [Reference Johnson7]. The RSV attachment G glycoprotein is responsible for virus binding to the host cell surface receptor and is a target of human neutralizing antibodies, together with the fusion (F) glycoprotein [Reference Sorce3]. The G protein shows the highest genetic and antigenic variability among all the RSV structural proteins. There is abundant evidence of accumulating amino acid changes in its hypervariable regions over time [Reference Johnson7, Reference Cane and Pringle10], making its study relevant in vaccine development strategies [Reference Otieno11].
The dynamics of RSV circulation are further demonstrated with the emergence of new genotypes that spread rapidly worldwide and replace previously circulating genotypes [Reference Slovic6]. Such changes were shown after the emergence of the BA genotype [Reference Trento12]. Among the RSV subgroup B, almost all strains detected after 2005 belong to the BA genotype [Reference Trento13]. Coincidently, a similar duplication event in RSV subgroup A led to the emergence of the ON1 genotype [Reference Eshaghi14], and the importance of epidemiological prevalence has also been reported for the ON1 genotype [Reference Kim15, Reference Pierangeli16]. It is assumed that these specific genetic characteristics provide a selective evolutionary advantage for these viruses [Reference Prifert17].
RSV coinfection with other respiratory viruses has been previously studied, in addition to its simple genetic characterization and annual distribution. The relevance to clinical severity in itself and coinfection with other respiratory viruses has also been investigated. For example, several studies revealed that the human rhinovirus (hRV) is the most common virus that coinfects with RSV. A greater number of viruses are not necessarily synonymous with greater disease severity [Reference Martinez-Roig18], and coinfection occurs frequently with RSV and hRV, as well as RSV and human bocavirus (hBoV). Overall, viral coinfection does not present greater severity, but can present mixed clinical features [Reference Calvo19]. However, some studies showed contradictory results. Children coinfected with RSV and other viruses presented more frequently with pneumonia than those with single RSV infection [Reference Asner20]. A previous study showed that during coinfections, one virus could block another simply by being the first to infect the available host cells [Reference Pinky and Dobrovolny21].
In this study, we investigated the molecular epidemiologic characteristics of RSV and its coinfection with other respiratory viruses in Gyeonggi Province in South Korea.
METHODS
Subjects and data collection
This study was conducted by the Korea National Institute of Health (KNIH) as part of the Korea Influenza and Respiratory Viruses Surveillance System (KINRESS). As part of the KINRESS, we collected a total of 4028 respiratory specimens from local hospitals, which are primary hospitals that first level of health system in Gyeonggi Province, South Korea from 2009 to 2014. Epidemiologic data and respiratory specimens from all patients with influenza-like illnesses were sent from local hospitals to the Gyeonggi Province Institute of Health and Environment (GIHE). Swab specimens were immediately placed at 4 °C in 2-ml cryovials containing 1·5 ml of cold viral transport medium (VTM; Difco, USA), and sent to the laboratory. Upon receipt, the specimens were processed immediately for virus detection, identification, and characterization. Aliquots of samples were also stored at −80 °C for additional analysis. Case reports and laboratory analysis data were entered into a KINRESS database (the Korea National Institute of Health, South Korea). The samples were collected for respiratory virus diagnosis, and written informed consent was obtained from the patients, their parents or legal guardians [Reference Kim15]. This study was approved by the Korea National Institute of Health Institutional Review Boards (Approval Nos. 2010-03EXP-1-R, 2011–06EXP-01-C, 2012-08EXP-06-3C, 2013-08EXP-03-5C and 2014-08EXP-08-6C-A).
Virus detection and subsequent molecular analysis
Viral RNA was extracted from 140 µl of each respiratory specimen using QIAamp Viral RNA Mini Kits (Qiagen GmbH, Hilden, Germany). Seven Multiplex-PCR assays were performed for screening 16 respiratory viral pathogens including RSV, influenza (IFV)-A and B, adenovirus (ADV), hBoV, parainfluenza virus (PIV) 1–3, human metapneumovirus (hMPV), human coronaviruses (hCoV) (229E, OC43, NL63), and hRV [Reference Kim15]. Confirmed RSV-positive samples were further screened for subgroup (A/B) [Reference Kim15] and genotyped by amplifying and sequencing the second hypervariable region of the G protein. The G gene-specific primer set for sequence analysis was as follows: forward primer G(151–173) F: CTGGCAATGATAATCTCAACTTC and reverse primer F(3–22) R: CAACTCCATTGTTATTTGCC [Reference Antoniassi da Silva22]. The cDNA was prepared using the viral RNA extraction method employed by the routine respiratory virus test. The reaction mixture contained 5 µl of RNA, which was mixed with a final concentration of 10 mM dNTPs, 20 µM random primer, 1× RT buffer, 200 U of Superscript III reverse transcriptase (Invitrogen, CA, USA), 40 U of RNase-OUT RNase inhibitor (Invitrogen), 25 mM MgCl2, 0·1 mM dithiothreitol (DTT), and RNase-free water in a final volume of 20 µl. The mixture was then incubated at 25 °C for 5 min, 50 °C for 60 min, and 72 °C for 5 min to terminate cDNA synthesis. Next, 5 µl of cDNA were added to a PCR mixture containing 1 µl of SP-Taq DNA polymerase (2·5 U/μl) (Cosmo Genetech, Seoul, South Korea), 36 µl of distilled water, 5 µl of 10× PCR buffer, 1 µl of 10 mM dNTPs, and 1 µl each of the forward and reverse primers (both 20 µM) for the G gene. Primary denaturation was conducted at 95 °C for 10 min, which was followed by 35 cycles of PCR where each cycle comprised denaturation for 40 s at 95 °C, annealing for 40 s at 54 °C, and elongation for 1 min at 72 °C, with a final extension cycle of 5 min at 72 °C. The PCR products were separated by electrophoresis using 1% agarose gel and visualized using 1× SYBR Safe DNA Gel Stain (Invitrogen) [Reference Kim15]. BigDye Terminator ver.3·1 (Applied Biosystem, Foster City, CA, USA) was utilized for the sequencing reaction, and nucleotide sequence analysis was performed using a 3730 DNA Analyzer (Applied Biosystem). Multiple nucleotide sequences were aligned and edited with ClustalW ver.1·8 software. Phylogenetic analysis was executed using the neighbor-joining method with a bootstrap value of 1000 replicates for testing statistical significance of the tree topology using MEGA ver.6·06 software [Reference Yoshihara23]. The representative 14 sequences were submitted to Genbank and assigned accession numbers of KY773693-KY773706. Network analysis was visualized using Cytoscape ver.2·6·1 software.
Statistical analysis
The categorical variables were compared using the two-tailed Chi-squared test and multivariate analysis, and multiple logistic regression analysis was applied to estimate the odds ratio (OR) and 95% confidence interval (CI). Statistical analyses were performed using the R.3·0·1 tool. P-values <0.05 were considered to indicate statistical significance.
RESULTS
Prevalence of respiratory viral cases
Among 4028 analyzed samples, 1635 samples were positive for respiratory viral infections of various types (Table 1). The annual percentage of cases that was positive for respiratory infections were 18·0% (75/416) in 2009, 26·9% (179/665) in 2010, 48·8% (335/686) in 2011, 51·7% (474/917) in 2012, 47·1% (345/732) in 2013, and 37·1% (227/612) in 2014. Most of these infections were due to IFV (overall: 22·8%, type A: 14·3%, type B: 8·5%), hRV (14·2%), and ADV (10·0%). RSV viral infection was detected in ~4·5% (n = 183/4028) of the patients enrolled.
Bold values are target virus value (RSV) in our study.
Epidemiological and clinical characteristics of patients with viral respiratory infections
The epidemiological and clinical characteristics of patients with viral respiratory infections are summarized in Tables 2 and 3; 50·4% of the patients with viral infections were female (PIV (50·5%), hCoV (54·5%), hMPV (57·5%), and IFV (52·9%)). Among the respiratory infections, hBoV and RSV infected the youngest patients (median age = 2 years), whereas IFV was more prevalent in older children (median age = 7 years) than any other viruses. Clinical symptoms such as fever (overall = 81·3%: pediatric patients, 72·7% vs. adult, 8·5%), cough (overall = 66·5%: pediatric patients, 58·5% vs. adult, 8·0%), nasal discharge (overall = 59·1%: pediatric patients, 52·2% vs. adult, 6·8%), and phlegm (overall = 47·8%: pediatric patients, 43·8% vs. adult, 4·0%) were frequently identified among the infected patients. However, several patients might have had conditions such as febrile seizure, and stomachache that were not listed as clinical symptoms. The most common pre-existing conditions among the patients were asthma (overall = 0·7%: pediatric patients, 0·6% vs. adult, 0·1%) and hypertension (overall = 0·9%: pediatric patients, 0% vs. adult, 0·9%). For example, patients infected with IFV had a high prevalence of asthma (n = 8, 30·8%) and hypertension (n = 8, 25·8%); additionally, these patients showed a high prevalence of antibiotic usage (140/210, 66·7%). Pediatric patients mean below 19 years old and adult means over 20 years old.
* Clinical conditions including febrile seizure and stomachache.
† Case was defined based WHO surveillance [43].
‡ Case with acute respiratory illness except influenza.
Annual incidence of RSV and RSV subgroups in patients with respiratory viruses
Among 4028 samples, 1635 samples had positive results for respiratory infections; of these 183 (4·5%) were RSV-positive, as confirmed by multiplex-PCR assays. Using real-time one-step RT–PCR to distinguish between RSV subgroups A and B, we categorized 141 RSV-positive samples as subgroup A (80·6%), 33 as subgroup B (18·8%), and 1 as RSV-A/B coinfection (0·6%). RSV-A was the major subgroup between January 2009 and November 2014, with January 2010 to December 2010 as the only exception. RSV-B (n = 8, 88·9%) was found in the majority of RSV cases from January 2013 to December 2013, while it was detected at only a slightly higher rate (n = 12, 57·1%) than RSV-A infections (n = 9, 42·9%). Most cases of RSV-A (n = 112, 79·4%) were detected from January 2011 to December 2012, whereas most cases of RSV-B (n = 12, 36·4%) were detected from January 2013 to December 2013.
Phylogenetic analysis of RSV and annual distribution of RSV subgroup A and B genotypes
We conducted viral genotyping based on the nucleic acid sequence of the second hypervariable region of the G protein. The target region was sequenced and aligned using MEGA ver.6·06. This phylogenetic analysis included 40 representative RSV-A isolates, 14 RSV-A reference sequences, 31 RSV-B isolates, and 26 RSV-B reference sequences. Of the RSV subgroup A (n = 141) and RSV subgroup B (n = 33) samples, 131 RSV-A (92·9%) and 31 RSV-B (93·9%) were sequenced successfully. RSV-A samples clustered in the genotypes of ON1 (n = 61, 43·3%), NA1 (n = 66, 46·8%), GA1 (n = 1, 0·7%), and GA5 (n = 3, 2·1%), whereas RSV-B samples clustered in the genotypes of BA9 (n = 29, 87·9%) and BA10 (n = 2, 6·1%) (Fig. 1). The ON1 and NA1 genotypes (total: n = 127, 96·9%) constituted the majority of RSV-A genotypes, while the BA9 genotype (n = 29, 87·9%) constituted the majority of RSV-B genotypes. In contrast, GA1 and GA5 genotypes (n = 4, 3·1%) in RSV-A, and BA10 genotype (n = 2, 6·1%) in RSV-B were observed less frequently. From 2009 to 2012, the NA1 genotype was consistently identified in the study subjects. The ON1 genotype was newly detected from 2011 to 2014. Between 2011 and 2012, both ON1 and NA1 genotypes were identified; however, the NA1 genotype was not identified from 2013 to 2014. These results suggest that the ON1 genotype of RSV newly emerged and was first isolated in December 2011. Then, it gradually dominated and replaced the NA1 genotype. For the RSV-B genotype, BA9 was observed continuously, except in 2009. There were no other notable genetic changes detected in RSV-B genotypes during the study period (Table 4, Fig. 2).
General information and clinical characteristics of patients with RSV subgroups and genotypes
Table 5 summarizes the demographic and clinical characteristics of patients associated with RSV subgroups and genotypes. The ON1 genotype (n = 31, 50·8%), RSV subgroup B (n = 19, 55·9%), BA9 genotype (n = 15, 51·7%), and BA10 genotype (n = 2, 100%) had more females than males. All RSV subgroups and genotypes were associated with the young patient group (median age of below 3 years) and the GA5 genotype with the youngest group (median age of 1·0 years). Most of the cases were 0–2 years of age (n = 116, 63·4%); however, two older patients (1 in the 50–65 year range and 1 >65 years) were infected with RSV subgroup B (BA9 genotype) (Fig. 3a ). We performed network analysis to address the correlation between RSV genotype and clinical symptoms (Fig. 3b ). Fever, cough, nasal discharge, and phlegm were the most common clinical symptoms; however, the ON1 genotype was associated with slightly more nasal discharge and phlegm than the NA1 genotype. Subjects infected with the BA9 genotype (RSV-B) had labored breathing, which is generally regarded as a serious symptom (Table 5).
Multiple logistic regression analysis of association of clinical characteristics in patients with RSV and RSV subgroups and genotypes
We performed multiple logistic regression analysis to understand the association between clinical features and RSV by dividing into pediatric patients (below 19 years) and adult (over 20 years) as a separate group. Demographical information (age and gender) and 16 clinical symptoms were included in the statistical analysis. OR values >1 indicate a positive association. Clinical symptoms such as cough (OR = 2·8, 95% CI 1·6–5·1), wheezing (OR = 2·8, 95% CI 1·7–4·4), and vomiting (OR = 2·2, 95% CI 1·2–3·9) were significantly higher in RSV-infected groups compared with non-RSV-infected groups on pediatric patients (P < 0·05). No statistically significant ORs were observed in the RSV subgroup A and BA9 genotype. In the ON1 genotype, only phlegm (OR = 11·8, 95% CI 3·8–46·7) was statistically significant (P < 0·05) and the NA1 genotype showed a significant association with gender (males, OR = 2·4, 95% CI 1·1–5·4) and chills (OR = 5·1, 95% CI 1·1–27·1). RSV subgroup B showed a significant association with nasal obstruction (OR = 4·6, 95% CI 1·2–20·0), but no other symptoms showed a significant association (Table 6). However, No statistically significant ORs were observed in the RSV-infected groups compared with non-RSV-infected groups on adult (data not shown) and there are only 3 data on adult in RSV subgroups and genotypes that is too little to analyze.
Statistically significant values are in bold.
* Statistically significant.
† ND [not determined]: missing data.
‡ ND [not determined]: one data point.
Coinfection with RSV and other respiratory viruses
We investigated simultaneous infections with RSV and other respiratory viruses using a multiplex-PCR assay. Among all tested respiratory viruses, hRV (n = 25, 47·2%) and ADV (n = 12, 22·6%) were detected most frequently with RSV. Most coinfections were detected from January 2011 to December 2012 (n = 45, 84·9%) proportionate with RSV-positive cases. Among all RSV-positive cases (n = 183), RSV single-infection was observed with the most frequency (n = 137, 74·5%), followed by mono-coinfections (n = 39, 21·2%) and dual-coinfections (n = 7, 3·8%) with RSV. Due to the number of cases in our study, only hRV coinfection with RSV (n = 21) could be evaluated with RSV single-infection (n = 138) in demographic and clinical comparisons between patients with RSV single- and coinfections. In RSV single-infections, fever (n = 117, 85·4%), cough (n = 115, 83·9%), nasal discharge (n = 100, 73·0%), and phlegm (n = 88, 64·2%) were the major clinical symptoms. In cases with hRV and RSV coinfection, fever (n = 20, 95·2%), cough (n = 18, 85·7%), nasal discharge (n = 17, 81·0%), and phlegm (n = 17, 81·0%) were the major clinical symptoms (Table 7). We performed a two-tailed χ 2 test to test the statistical difference between two independent groups [Reference Yoshihara23]. No significant association was detected between single (n = 137) and overall coinfection (n = 46) or between single (n = 137) and RSV and hRV coinfection (n = 21) (P > 0·05) (data not shown).
DISCUSSION
In this study, we investigated the distributions of all respiratory viruses, the RSV subgroups, and the coinfections of RSV with other respiratory viruses using real-time one-step RT–PCR. The RSV genotypes were determined by sequencing the second hypervariable region of the G protein, and the results were correlated with the demographic and clinical characteristics from consented clinical information from patients.
RSV evolved and changed continuously after an RSV presumed outbreak reported in 1941 [Reference Adams24]. For six consecutive years in Gyeonggi Province located in South Korea, RSV-A and B cocirculated; RSV-A predominated during most seasons, but RSV-B predominated in 2010 and 2013 (Table 4). This tendency was consistent with the results from previous studies in northern Italy [Reference Esposito25] and Japan [Reference Dapat26], but differed somewhat from the results from other countries in East Asia [Reference Cui8, Reference Cui27, Reference Yu28]. The difference appears to be due to various geographic regions and periods. Consistent with previous reports, higher proportions of infected patients were male [Reference Cui8, Reference Esposito25, Reference Yu28], and most RSV infections were detected in patients <2 years of age [Reference Gimferrer29]. In terms of seasonality, RSV epidemics commonly span the months from October to March, peaking between November and January [Reference Thorburn30], and our study showed an analogous pattern. However, different seasonal patterns were noted in tropical regions such as Vietnam, with peaks in the hot and dry season (July to September), suggesting a role of climate parameters such as temperature and relative humidity [Reference Yoshihara23].
After the first discovery of the ON1 genotype in Ontario, Canada in 2010 [Reference Eshaghi14], the ON1 genotype was reported in many other countries, including Japan [Reference Tsukagoshi31], China [Reference Cui32], Germany [Reference Prifert17], Italy [Reference Pierangeli16], India [Reference Choudhary33], and South Africa [Reference Valley-Omar34]. In our study, the ON1 genotype emerged in October 2011. In South Korea, the first ON1 genotype was discovered in August 2011 [Reference Lee35], which was the next season after Eshaghi et al. first discovered the novel ON1 genotype in Canada [Reference Eshaghi14, Reference Kim15]. Other previous studies in South Korea [Reference Kim15, Reference Lee35] used clinical samples collected from the Korean Influenza and Respiratory Viruses Surveillance System, which was incorporated nationwide with the Institute of Health and Environment. Therefore, the first ON1 genotype discovered in August 2011 may be from clinical specimens collected from another province in South Korea. From 2009 to 2010, most RSV subgroup A infections were clustered in the NA1 genotype; from 2011 to 2012, the NA1 genotype and ON1 genotypes coexisted but demonstrated opposite predominance: 2011 showed genotype distributions of 92% NA1 and 6% ON1, whereas 2012 showed distributions of 21% NA1 and 79% ON1. Furthermore, no cases of NA1 genotype (subgroup A) were detected from January 2013 to November 2014. Based on the pattern of disease epidemics, we surmise that the NA1 genotype (ancestor of ON1) was replaced with the ON1 genotype during 2011 to 2012. Such substitutions were observed in studies conducted in other countries [Reference Yoshihara23, Reference Esposito25].
The results of many studies [Reference Eshaghi14, Reference Fall36] demonstrate that the main characteristic of the ON1 genotype is the presence of a 72-nucleotide duplication insertion in the G protein, and this characteristic may have a positive evolutionary effect on RSV. For example, genetic variations may play a crucial biological role for enhancing the efficiency of viral attachment to the cell receptors or for faster viral replication during pathogenesis [Reference Yoshihara23], which may elicit strong resistance to herd immunity [Reference Cui27]. The BA genotype of RSV-B gained a 60-nucleotide duplication insertion in the G protein, a change that occurred since the BA genotype was first isolated in 1999 (Buenos Aires, Argentina) [Reference Trento12]. The GA1 genotype was newly isolated in September 2014, but reports showed that it was isolated previously in 2009. In addition, isolation of the GA1 genotype was reported in Senegal, Africa in 2015 [Reference Fall36] (Table 4).
Demographic and clinical information is important to complement the limitations of advanced molecular technology, as it can take several days to identify a respiratory virus genotype in an infected patient. Therefore, the relevance of clinical information and RSV-positive infections has been investigated continuously through various studies [Reference Pierangeli16, Reference Yoshihara23, Reference Esposito25, Reference Yu28, Reference Fall36]. We performed two analyses with respect to the clinical relevance of RSV-positive infections. We first performed a network analysis using simple frequencies, and then conducted a multiple logistic regression analysis. In RSV-infected patients, clinical symptoms such as fever, cough, nasal discharge, and phlegm were major signs. Serious symptoms such as labored breathing (n = 1) and constricted chest (n = 1) were reported in patients with the BA9 genotype in RSV-B, however, the low frequency prevented a robust analysis (Table 5, Fig. 3b ). On the other hand, multiple logistic regression analysis on pediatric patients showed that RSV-positive infections were associated with a 2·8-fold increased risk of cough and wheezing. Furthermore, the ON1 genotype was associated with an 11·8-fold increased risk of phlegm, and the NA1 genotype was associated with a 5·1-fold increased risk of chills and RSV subgroup B was associated with a 4·6-fold increased risk of nasal obstruction (Table 6). However, these results had limitations because the demographic and clinical information was not fully recorded in some patients and our study was performed through local hospital in Gyeonggi Province and these hospitals were all primary hospital. Therefore, most of clinical symptoms were mild and so few patients had serious symptoms. Further studies will be necessary to overcome these limitations and understand the exact associations of clinical symptoms and respiratory virus genotypes.
Many studies have reported contradictory results about coinfection and clinical severity between RSV subgroup and genotype [Reference Martinez-Roig18, Reference Calvo19, Reference Pinky and Dobrovolny21, Reference Yu28, Reference Tran37–Reference Skjerven40]. hRV was the most commonly detected virus (47·2%) with RSV among all respiratory viruses in Gyeoggi Province (Table 7). These results were consistent with other studies in China, Spain, and Norway [Reference Martinez-Roig18, Reference Calvo19, Reference Yu28, Reference Skjerven40]. No significant associations were detected between coinfection and the severity of clinical symptoms between RSV single-infection and overall coinfection, and between RSV single-infection and hRV and RSV coinfection (data not shown); these results were consistent with those reported in studies conducted in Vietnam [Reference Tran37], UK [Reference Cebey-López38], and Spain [Reference Martinez-Roig18, Reference Calvo19]. However, other studies have shown that coinfections with RSV and respiratory viruses such as hRV, hMPV, and PIV type 3 increase clinical severity [Reference Míguez, Iftimi and Montes39]. Recently, experimental studies were performed to determine the mechanism underlying coinfection and clinical severity. According to Pinky and Dobrovolny [Reference Pinky and Dobrovolny21], one virus can block another during coinfection by being the first to infect the available host cells. hRV, which is a fast-growing virus, reduces the replication of the remaining viruses. The slowest growing virus is suppressed in the presence of another virus, and the effect of viral coinfections on clinical outcomes is no more severe than a single-virus infection, and can even result in a less severe clinical impact. Skjerven et al. [Reference Skjerven40] reported that viral genomic load shows a positive correlation with disease severity. In other words, clinical severity depends on the viral genomic load rather than coinfection. However, other studies reported no correlation between viral load and clinical symptoms [Reference Lee41, Reference Franz42]. Further study is warranted to investigate this clinical aspect in future.
We investigated the molecular epidemiological characteristics of RSV in patients with viral respiratory infections. The association between clinical symptoms in patients and the molecular epidemiology of RSV infections will be useful for developing treatment therapies and vaccines against RSV.
ACKNOWLEDGEMENT
This study was supported by the Intramural Research Fund (grant no. 4851-300) of the Korea National Institute of Health.
DECLARATION OF INTEREST
None.