Genetic diversity and epidemiology of human rhinovirus among children with severe acute respiratory tract infection in Guangzhou, China
Virology Journal volume 18, Article number: 174 (2021)
Human rhinovirus (HRV) is one of the major viruses of acute respiratory tract disease among infants and young children. This work aimed to understand the epidemiological and phylogenetic features of HRV in Guangzhou, China. In addition, the clinical characteristics of hospitalized children infected with different subtype of HRV was investigated.
Hospitalized children aged < 14 years old with acute respiratory tract infections were enrolled from August 2018 to December 2019. HRV was screened for by a real-time reverse-transcription PCR targeting the viral 5′UTR.
HRV was detected in 6.41% of the 655 specimens. HRV infection was frequently observed in children under 2 years old (57.13%). HRV-A and HRV-C were detected in 18 (45%) and 22 (55%) specimens. All 40 HRV strains detected were classified into 29 genotypes. The molecular evolutionary rate of HRV-C was estimated to be 3.34 × 10–3 substitutions/site/year and was faster than HRV-A (7.79 × 10–4 substitutions/site/year). Children who experienced rhinorrhoea were more common in the HRV-C infection patients than HRV-A. The viral load was higher in HRV-C detection group than HRV-A detection group (p = 0.0148). The median peak symptom score was higher in patients with HRV-C infection as compared to HRV-A (p = 0.0543), even though the difference did not significance.
This study revealed the molecular epidemiological characteristics of HRV in patients with respiratory infections in southern China. Children infected with HRV-C caused more severe disease characteristics than HRV-A, which might be connected with higher viral load in patients infected with HRV-C. These findings will provide valuable information for the pathogenic mechanism and treatment of HRV infection.
Human rhinoviruses (HRV) are the most commonly encountered respiratory viruses and the most frequent causes of acute respiratory infections in young children and infants. HRV is generally associated with common cold and mild upper respiratory infections, but may also lead to more severe lower respiratory tract illnesses, such as pneumonia, bronchiolitis and asthma [1, 2].
HRV belongs to the genus Enterovirus and family Picornaviridae, and currently has more than 100 serotypes, classified into three species: HRV-A, B and C. Species HRV-A and HRV-B were discovered in the 1950s , but HRV-C was identified in 2006 using novel molecular-based techniques [4, 5]. There are more than 100 distinct serotypes of HRVs identified according to their capsid proteins. All HRV species have been identified in throughout the year in many regions, and HRV-A and HRV-C appear to be the predominant species detected in patients with acute respiratory infection . Recent studies suggest that the illness severity differs among HRV species [7, 8], and HRV-C has been shown more frequently associated with severe asthma attacks and lower respiratory tract infections compared with other HRV species [9,10,11]. Also, recently studies have shown that HRV-C is possibly more virulent and cause more severe illness [8, 12,13,14].
In addition, the association of HRV species with clinical severity is not well understood. In this study, we investigated the epidemiological, evolution, and clinical characteristics of HRV infections in children with acute respiratory tract infections. Furthermore, we studied the impact of HRV species and nasopharyngeal viral load on the disease severity of acute respiratory tract infections.
Study population and sample collection
A total of 655 nasopharyngeal swab specimens were collected from hospitalized infants or children with acute respiratory tract infection from Guangdong Panyu Maternal and Child Health Hospital between August 2018 and December 2019. Specimens were only taken from individuals primary diagnosis of viral infection and with \(\leqslant\) 3 days of fever (temperature \(\geqslant\) 37.5 °C), cough, sputum, throat sore, dyspnea and/or other acute respiratory tract infection symptoms. Nasopharyngeal swabs were collected within 24 h after admission by medical professionals, and all the samples were placed in viral transportation medium and stored at − 80 °C. Demographic information and clinical characteristics were recorded for each patient. The degree of disease severity of each HRV infection patient was estimated and scored according to a severity scoring system as described previously .
Detection of respiratory viruses
Total viral nucleic acids were extracted from the viral transportation medium using the QIAamp MinElute Virus Spin Kit (Qiagen, Valencia, CA) according to the manufacturer’s instructions. Reverse transcription was performed using RevertAid First Strand cDNA Synthesis kit (Invitrogen Life Technologies, Carlsbad, CA) with random primers. The cDNA was used for virus detection immediately or stored at -20˚C until further use. HRV infection was detected by using qRT-PCR with HRV specific primers and a probe, HRV forward primer (5′-TGGACAGGGTGTGAAGAGC-3′), reverse primer (5′-CAAAGTAGTCGGTCCCATCC-3′), and probe (FAM-TCCTCCGGCCCCTGAATG-BHQ1). Other viral infections were simultaneously screened, including human parainfluenza virus, influenza virus, respiratory syncytial virus, human coronaviruses (229E, NL63, HKU1, and OC43), metapneumovirus, adenovirus, and bocavirus.
Samples that tested positive for HRV were used to further determine genotypes, the VP4/VP2 and 5′UTR regions of HRV were amplified using nest-PCR, primers were synthesized from the published primer sequences . PCR was initiated at 95 °C for 10 s, followed by 35 cycles of 95 °C for 5 s, 55 °C for 30 s and 72 °C for 1 min, with a final extension at 72 °C for 10 min. Specimens from which amplification of the VP4/VP2 regions failed were defined as untyped. All sequencing was performed by Sangon Biotech Co., Ltd. (Shanghai, China) using ABI-PRISM 3730XL DNA sequencer (Applied Biosystems).
Phylogenetic analyses by neighbour-joining (NJ) and Bayesian Markov Chain Monte Carlo (MCMC) methods
The sequences obtained in this study were aligned with representative sequences retrieved from GenBank using Clustal W. The phylogenetic tree was constructed using the NJ method in MEGA 7.0 software , and the reliability of the tree was estimated with 1000 bootstrap replications.
Molecular evolutionary analysis was performed with Bayesian Markov Chain Monte Carlo (MCMC) method using BEAST ver.2.5.1 . The most suitable nucleotide substitution model (GTR + G) was selected using jModelTest 2.1.10 , and the datasets were analyzed with an uncorrelated lognormal relaxed clock model. Convergence was assessed using Tracer version 1.7.1, and it was accepted when the MCMC chain was run through enough steps to make the effective sample size (ESS) above 200 after a 10% burn-in. The maximum clade credibility (MCC) tree was constructed by Tree Annotate 2.5.1 after removing the first 10% of trees as burn-in, and the phylogenetic tree was visualized by FigTree v1.4.4. The uncertainties of the estimates were indicated by 95% highest posterior density (HPD) intervals.
Continuous variables were analyzed with the one-way analysis of variance (ANOVA) and Student’s t-test; chi-square was performed for ordinal or categorical data. Mann–Whitney U-test was used to compare severity scores between groups. A p value below 0.05 was considered statistically significant. Statistical analysis was performed with SPSS software (version 17.0; SPSS, Inc., Chicago, IL, USA).
Demographic characteristics and seasonal distribution of HRV infection
A total of 655 nasopharyngeal swab specimens were collected from children with acute respiratory illness. Of these, real-time PCR revealed that 42 (6.41%) of 655 hospitalized children were HRV positive, and there is no HRV-positive patients were co-infected with other respiratory viruses. As shown in Table 1, approximately 90.5% (38/42) of these children were under 5 years of age, and the majority of them were under 2 years old (57.13%). The age and gender distributions of the HRV-positive patients were not statistically significant. The prevalence of HRV infections during August 2018 to December 2019 in Guangzhou is shown in Fig. 1. Among the collected specimens, HRV-A was identified in 42.9% of samples (18/42), HRV-C in 52.4% (22/42) and HRV-B was not identified. In addition, there are two samples that were not successfully typed. HRV was prevalent mainly in the winter, especially in October to November.
Phylogenetic analysis of HRV strains using the neighbor joining method
The VP4/VP2 region (approximately 420 bases) was sequenced for 40 real-time PCR positive samples. In order to confirm the classification of HRV types of the 40 sequences, phylogenetic analysis of the VP4/VP2 region sequences was performed. As shown in Fig. 2, the HRV-A strains for each genotype were as follows: A15 (1); A16 (1); A21 (1); A24 (1); A29 (1); A36 (1); A40 (1); A41 (1); A46 (1); A47 (1); A57 (1); A61 (3); A78 (1); A80 (1); A101 (1). The present HRV-C strains in each genotype were as follows: HRV-C2 (7); HRV-C3 (1); HRV-C5 (1); HRV-C6 (1); HRV-C13 (1); HRV-C14 (1); HRV-C17 (1); HRV-C20 (1); HRV-C21 (1); HRV-C31 (1); HRV-C35 (1); HRV-C41 (2); HRV-C46 (3); HRV-C48 (1).
The phylogenetic analysis of the 5’UTR sequences was also performed. As shown in Fig. 3, HRV-B and majority of HRV-A sequences formed a distinct clade, while the tree formed a clade that includes intermixed HRV-A and HRV-C strains. All HRV-positive samples were grouped into 18 genotypes, and one sample (3903-GD-CHN-2018) was not typed. Therefore, our observations indicate that the main species of HRV circulating in Guangzhou from 2018 to 2019 were HRV-A and HRV-C, whereas the prevalent genotype was HRV-C2. Moreover, multiple HRV serotypes could be detected within a time period.
Estimation of the time to the most recent common ancestor (tMRCA) and molecular evolutionary rate for HRV strains using the Bayesian MCMC method
The phylogenetic trees with Bayesian Markov Chain Monte Carlo (MCMC) method as shown in Fig. 4 were constructed to estimate the time-scaled evolution of HRV VP4/VP2 region for HRV global strains. The time-scaled Maximum Clade Credibility (MCC) tree revealed that the times to the most recent common ancestor (tMRCAs) of HRV-A was around 1310 years ago (95% highest probability density (HPD): 590–3396) and 1370 years ago (95% HPD: 590–3962) for HRV-C.
The molecular evolutionary rate of HRV-A strains was estimated to be 7.79 × 10–4 substitutions/site/year (95% HPD: 2.87 × 10–5-1.89 × 10–4), while that of HRV-C was estimated to be 3.34 × 10–3 substitutions/site/year (95% HPD: 1.63 × 10–3 to 5.26 × 10–3). The rate of HRV-C strains was significantly faster than that of HRV-A strains.
Clinical and laboratory characteristics in children with HRV infection
To characterize the clinical manifestations caused by each HRV species, the clinical symptomsand laboratory tests of HRV-positive patients are listed in Table 2. All HRV-positive patients presented with severe acute respiratory infection (SARI). Cough (87.5%), rhinorrhoea (62.5%), fever (62.5%) and expectoration (55%) were the most common symptoms at presentation. In addition, HRV-C positive patients more often were expectoration (p = 0.026) as compared to HRV-A positive patients. In other clinical features, no statistically significant difference was recorded.
Viral load and severity of infection
To determine the relationship between the viral load of different HRV genotypes and the severity of the disease. We compared median symptom score and Ct values for each HRV species. As shown in Fig. 5, the Ct value of HRV-A ranged from 28.15 to 37.05 and HRV-C ranged from 23.1 to 32.47. The viral load of HRV-A was significantly higher than that of HRV-C (p = 0.0148). According to the severity scoring system, the median peak symptom score for the 18 patients with HRV-A infection was 6 (range 4–10.25), and that 22 patients with HRV-C infection was 8 (range 5.5–10). And the median peak symptom score was higher in patients with HRV-C infection as compared to HRV-A (p = 0.0543), although the difference did not reach significance.
HRV is the predominant pathogen identified in hospitalized infants and young children with acute respiratory tract infection . HRV although mostly associated with mild disease, can also be a cause of severe illness . In this study, we analyzed the epidemiology and genotypic diversity of HRV within hospitalized children. During the monitoring, 42 (6.41%) were positive for HRV, which is consistent with previous research conducted in Guangdong (5.47%) . However, the prevalence of HRV could be different in across regions and years [22,23,24].
The highest incidence of HRV infection occurred among the 0–1-year-old infants, and 71.4% (30/42) of total HRV occurred in children younger than 3 years old. This is comparable with the results found in China and other countries [13, 24,25,26]. These observations support the view that infants were more susceptible to HRV.
The phylogenetic analysis shows that only HRV-A and HRV-C are circulating in children in the population we analyzed. The proportion of the HRV species revealed in this study (HRV-A, 45% and HRV-C, 55%) is consistent with prior studies [2, 23, 27, 28]. The evolutionary rate of viral genes differs among respiratory viruses as previous reported [29, 30]. In the present study, the evolutionary rate of the VP4/VP2 region of HRV-A (3.34 × 10–3 substitutions/site/year) was faster than that reported for global HRV-C strains (6.6 × 10–4 substitutions/site/year) , but similar to that reported for Japan HRV-C strain(3.07 × 10–3 substitutions/site/year) . Therefore, the evolution rate of HRV strains varies from region to region, but the precise mechanisms are not known yet.
Clinical symptoms of HRV infected patients were similar to those previously reported cough, rhinorrhoea, fever and expectoration . In this study there were no significant differences observed in clinical symptoms between two species, except more frequent expectoration in patients with HRV-C species. Previous studies have reported that HRV-C causes more severe clinical manifestations compared with the other HRV types [8, 12, 34]. And some observations suggested that in HRV-C infected patients with pneumonia, higher mean viral load was reported . However, another research shown there was no significant difference in median peak viral load between HRV species in hospitalized patients . In the present analysis, we observed that HRV-C viral load in hospitalized children was significantly higher than HRV-A, indicate that association between HRV species and viral load was found. Interestingly, HRV-C infected patient exhibits higher symptom score in comparison to patient infected with HRV-A. Although statistically significant differences was not found, the disease was more severe in patients infected with HRV-C comparison with HRV-A. Together, our findings suggest that the pathophysiology of infection with HRV-C may differ from HRV-A. Previous studies have demonstrated that HRV-A and HRV-B utilize the intercellular adhesion molecule 1 (ICAM-1) and low-density lipoprotein receptor (LDLR) family members as cellular receptor [37, 38]. A unique feature of HRV-C is its use of the cadherin-related family member 3 (CDHR3) for host cell entry. Moreover, a single-nucleotide polymorphism (C529) in CDHR3 is associated with enhanced viral binding and promoting viral replication with a consequent increase in viral load . Further study will be needed to determine the mechanism of HRV-A and HRV-C resulting in different disease severities.
In conclusion, we analyzed the molecular characterization of HRV and indicated that HRV-A and HRV-C were predominant in HRV infections. The genetic variability analysis of HRV highlights the genetic complexity and rapid evolution of HRV that warrant continuous molecular surveillance. Clinical data analysis reveals HRV-C and high viral load were shown to be the important determinants of the severity of acute respiratory illnesses. Our study highlight the fact that it is necessary to study the mechanism of specific HRV genotype in causing severe disease.
Availability of data and materials
The datasets used and/or analyzed in the current study are available from the corresponding author upon reasonable request.
severe acute respiratory infections
real-time reverse transcription polymerase chain reaction
molecular evolutionary genetics analysis
statistical product and service solutions
intercellular adhesion molecule 1
low-density lipoprotein receptor
cadherin-related family member 3
Markov Chain Monte Carlo
effective sample size
maximum clade credibility
highest posterior density
one-way analysis of variance
times to the most recent common ancestor
Nair H, et al. Global and regional burden of hospital admissions for severe acute lower respiratory infections in young children in 2010: a systematic analysis. The Lancet. 2013;381:1380–90.
Rahamat-Langendoen JC, et al. The significance of rhinovirus detection in hospitalized children: clinical, epidemiological and virological features. Clin Microbiol Infect Off Publ Eur Soc Clin Microbiol Infect Dis. 2013;19:E435-442.
Price WH. The isolation of a new virus associated with respiratory clinical disease in humans. Proc Natl Acad Sci U S A. 1956;42(12):892–6.
Lamson D, Renwick N, Kapoor V, Liu Z, Palacios G, Ju J, Dean A, St George K, Briese T, Lipkin WI. MassTag polymerase-chain-reaction detection of respiratory pathogens, including a new rhinovirus genotype, that caused influenza-like illness in New York State during 2004–2005. J Infect Dis. 2006;194(10):1398–402.
Lau SK, et al. Clinical features and complete genome characterization of a distinct human rhinovirus (HRV) genetic cluster, probably representing a previously undetected HRV species, HRV-C, associated with acute respiratory illness in children. J Clin Microbiol. 2007;45:3655–64.
Arakawa M, et al. Molecular epidemiological study of human rhinovirus species A, B and C from patients with acute respiratory illnesses in Japan. J Med Microbiol. 2012;61:410–9.
Calvo C, et al. Role of rhinovirus C in apparently life-threatening events in infants, Spain. Emerg Infect Dis. 2009;15:1506–8.
Khetsuriani N, et al. Novel human rhinoviruses and exacerbation of asthma in children. Emerg Infect Dis. 2008;14:1793–6.
Linsuwanon P, et al. High prevalence of human rhinovirus C infection in Thai children with acute lower respiratory tract disease. J Infect. 2009;59:115–21.
Bizzintino J, et al. Association between human rhinovirus C and severity of acute asthma in children. Eur Respir J. 2011;37:1037–42.
Gern JE. The ABCs of rhinoviruses, wheezing, and asthma. J Virol. 2010;84:7418–26.
Cox DW, et al. Human rhinovirus species C infection in young children with acute wheeze is associated with increased acute respiratory hospital admissions. Am J Respir Crit Care Med. 2013;188:1358–64.
Iwane MK, et al. Human rhinovirus species associated with hospitalizations for acute respiratory illness in young US children. J Infect Dis. 2011;204:1702–10.
Lee WM, et al. Human rhinovirus species and season of infection determine illness severity. Am J Respir Crit Care Med. 2012;186:886–91.
Gern JE, et al. Relationships among specific viral pathogens, virus-induced interleukin-8, and respiratory symptoms in infancy. Pediatr Allergy Immunol. 2002;13:386–93.
Wisdom A, et al. Screening respiratory samples for detection of human rhinoviruses (HRVs) and enteroviruses: comprehensive VP4–VP2 typing reveals high incidence and genetic diversity of HRV species C. J Clin Microbiol. 2009;47:3958–67.
Kumar S, et al. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Bouckaert R, et al. BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput Biol. 2014;10:6.
Darriba D, et al. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772–772.
Drysdale SB, et al. Rhinovirus: not just the common cold. J Infect. 2017;74:S41–6.
Luo HJ, et al. 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 Off Publ Int Soc Infect Dis. 2020;90:5–17.
Zheng SY, et al. Epidemiological analysis and follow-up of human rhinovirus infection in children with asthma exacerbation. J Med Virol. 2018;90:219–28.
Zhao Y, et al. Genotypic diversity and epidemiology of human rhinovirus among children with severe acute respiratory tract infection in Shanghai, 2013–2015. Front Microbiol. 2018;9:1836.
Yan Y, et al. Clinical and epidemiological profiles including meteorological factors of low respiratory tract infection due to human rhinovirus in hospitalized children. Ital J Pediatr. 2017;43:23.
Zeng SZ, et al. Prevalence of human rhinovirus in children admitted to hospital with acute lower respiratory tract infections in Changsha, China. J Med Virol. 2014;86:1983–9.
Linder JE, et al. Human rhinovirus C: age, season, and lower respiratory illness over the past 3 decades. J Allergy Clin Immunol. 2013;131(69–77):e61-66.
Lauinger IL, et al. Patient characteristics and severity of human rhinovirus infections in children. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2013;58:216–20.
Jacobs SE, et al. Clinical and molecular epidemiology of human rhinovirus infections in patients with hematologic malignancy. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2015;71:51–8.
Mizuta K, et al. Detailed genetic analysis of hemagglutinin-neuraminidase glycoprotein gene in human parainfluenza virus type 1 isolates from patients with acute respiratory infection between 2002 and 2009 in Yamagata Prefecture, Japan. Virol J. 2011;8:533.
Peck KM, Lauring AS. Complexities of viral mutation rates. J Virol. 2018;92:e01031-17.
Briese T, et al. Global distribution of novel rhinovirus genotype. Emerg Infect Dis. 2008;14:944–7.
Kiyota N, et al. Genetic analysis of human rhinovirus species A to C detected in patients with acute respiratory infection in Kumamoto Prefecture, Japan 2011–2012. Infect Genet Evol J Mol Epidemiol Evol Genet Infect Dis. 2014;21:90–102.
Peltola V, et al. Rhinovirus infections in children: a retrospective and prospective hospital-based study. J Med Virol. 2009;81:1831–8.
Pierangeli A, et al. Molecular epidemiology and genetic diversity of human rhinovirus affecting hospitalized children in Rome. Med Microbiol Immunol. 2013;202:303–11.
Choi SH, et al. Clinical and molecular characterization of rhinoviruses A, B, and C in adult patients with pneumonia. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2015;63:70–5.
Piralla A, et al. Clinical severity and molecular typing of human rhinovirus C strains during a fall outbreak affecting hospitalized patients. J Clin Virol Off Publ Pan Am Soc Clin Virol. 2009;45:311–7.
Greve JM, et al. The major human rhinovirus receptor is ICAM-1. Cell. 1989;56:839–47.
Hofer F, et al. Members of the low density lipoprotein receptor family mediate cell entry of a minor-group common cold virus. Proc Natl Acad Sci U S A. 1994;91:1839–42.
Bochkov YA, et al. Cadherin-related family member 3, a childhood asthma susceptibility gene product, mediates rhinovirus C binding and replication. Proc Natl Acad Sci U S A. 2015;112:5485–90.
We appreciate the medical doctors, nurses and laboratory technicians at First Affiliated Hospital of Jinan University and Guangdong Panyu Maternal and Child Health Hospital for their support on epidemiological data and clinical sample collection.
This research has been supported by grants from the National Major Projects of Major Infectious Disease Control and Prevention, the Ministry of Science and Technology of the People’s Republic of China (No. 2017ZX10103011-007), the National Natural Science Foundation of China (Nos. 32070149, 31770178).
Ethics approval and consent to participate
The project was approved by The First Affiliated Hospital of Jinan University, and the committee’s reference number is KY-2019-056. We explained the details of our study to each subject and written informed consent was obtained from all participants prior to their inclusion in the study. Nasal samples and medical data were collected and analyzed anonymously.
Consent for publication
The authors have declared that no competing interests exist.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Li, W., Yu, B., Zhou, J. et al. Genetic diversity and epidemiology of human rhinovirus among children with severe acute respiratory tract infection in Guangzhou, China. Virol J 18, 174 (2021). https://doi.org/10.1186/s12985-021-01645-6