Skip to main content

Etiology and clinical characteristics of SARS-CoV-2 and other human coronaviruses among children in Zhejiang Province, China 2017–2019



A novel coronavirus (SARS-CoV-2) emerging has put global public health institutes on high alert. Little is known about the epidemiology and clinical characteristics of human coronaviruses infections in relation to infections with other respiratory viruses.


From February 2017 to December 2019, 3660 respiratory samples submitted to Zhejiang Children Hospital with acute respiratory symptoms were tested for four human coronaviruses RNA by a novel two-tube multiplex reverse transcription polymerase chain reaction assays. Samples were also screened for the occurrence of SARS-CoV-2 by reverse transcription-PCR analysis.


Coronavirus RNAs were detected in 144 (3.93%) specimens: HCoV-HKU1 in 38 specimens, HCoV-NL63 in 62 specimens, HCoV-OC43 in 38 specimens and HCoV-229E in 8 specimens. Genomes for SARS-CoV-2 were absent in all specimens by RT-PCR analysis during the study period. The majority of HCoV infections occurred during fall months. No significant differences in gender, sample type, year were seen across species. 37.5 to 52.6% of coronaviruses detected were in specimens testing positive for other respiratory viruses. Phylogenic analysis identified that Zhejiang coronaviruses belong to multiple lineages of the coronaviruses circulating in other countries and areas.


Common HCoVs may have annual peaks of circulation in fall months in the Zhejiang province, China. Genetic relatedness to the coronaviruses in other regions suggests further surveillance on human coronaviruses in clinical samples are clearly needed to understand their patterns of activity and role in the emergence of novel coronaviruses.


Human coronaviruses (HCoVs) are common viruses that cause mild to severe respiratory infections. To date, there are seven CoVs have been confirmed to infect human, including two α-CoVs (HCoV-229E and HCoV-NL63) and five β -CoVs [HCoV-HKU1, HCoV-OC43, severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2, and Middle East respiratory syndrome coronavirus (MERS-CoV)] [1]. The ongoing pandemic of COVID-19 offers the most recent example of the tragedy when highly pathogenic coronaviruses emerge in the human population [2,3,4]. As of 23 March 2021, approximately 25 million cases have been reported to the WHO accounting for a 2.2% case fatality rate. Over the last two decades, the emergence of Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV) has been associated with severe atypical pneumonia and caused 774 deaths [5]. In 2012, a man in Saudi Arabia who infected by Middle East Respiratory Syndrome Coronavirus (MERS-CoV) died of acute pneumonia and renal failure [6]. The continuing spread of MERS-CoV resulted in unprecedented outbreaks and high case-fatality rates.

HCoV-NL63, HCoV-HKU1, HCoV-229E, and HCoV-OC43 circulate worldwide and cause a range of respiratory symptoms [7]. These four HCoVs result in 30% of common colds and acute respiratory infections such as pneumonia and bronchiolitis, especially in immunocompromised children [8, 9]. Lower respiratory tract infections are considered as the most common cause of death among children and are responsible for 15% mortality rate of children under 5 years of age globally [10, 11]. HCoVs pose serious threats to children health as they are associated with acute upper or lower respiratory tract infections [12, 13]. Although HCoVs are distributed globally, the predominant species may vary by region or year. Moreover, a notable feature of COVID-19 pandemic is the relative paucity of published report that describes the potential impact of SARS-CoV-2 on children with acute respiratory infection [14]. However, the prevalence and clinical profiles of HCoVs have been largely underdetermined and no studies described circulation of SARS-CoV-2 and four common HCoVs across multiple years in Zhejiang province, China.

Considering the discovery of novel SARS-CoV-2, we retrospectively investigated the role of SARS-CoV-2 and common human coronaviruses (HCoV-229E, HCoV-OC43, HCoV-NL63 and HCoV-HKU1) in childhood respiratory infections in Zhejiang province, China, and results were compared to data for other viral respiratory pathogens including influenza A and B, parainfluenza, respiratory syncytial virus, rhinovirus, and adenovirus.

Materials and methods

Ethical approval statement

The studies involving human participants were reviewed and approved by Zhejiang Provincial Center for Disease Control and Prevention of ethics committee. Children’s parents/guardians have been notified the purpose of the study and the right to keep information confidential. Children’s parents/guardians provided their written informed consent to participate in this study.

Sample collection

A total of 3660 retrospective respiratory samples were collected from children (age, ≤ 18 years) with acute respiratory tract infections who admitted to Children’s Hospital of Zhejiang University from 2017 to 2019. Written informed consent was obtained from the guardians of all participants before specimens and data collection. Specimens were collected in accordance with Zhejiang Provincial CDC guidance and contained nasopharyngeal swab (n = 26), sputum (n = 3590) and bronchoalveolar lavage (n = 41). 2138 samples were from male and the remaining were from female. Clinical specimens were initially kept in sterile EP tube with 5 mL viral transport medium (Becton, Dickinson and Company, NJ) and then transferred to Zhejiang Provincial CDC within the same day for routine monitoring. All clinical samples were stored at − 80 °C prior to blinded analysis. Clinical data including symptoms, history of illness, severity of infection, demographic data, drug usage and laboratory examination were recorded for each patient and reviewed retrospectively when patients confirmed with HCoV infections.

Nucleic acid extraction

Viral nucleic acid was extracted from 200 μL of clinical specimens using QIAGEN Cador Pathogen 96 QIAcube HT Kit (Qiagen, Hilden, German) and following manufacturer’s instructions. The extracts were eluted into 50 μL of DNase- and RNase-free water and stored at − 80 °C until used as a template for PCR assay.

Two-tube multiplex reverse transcription PCR (RT-PCR) assay

A two-tube multiplex RT-PCR assay (a two-tube assay) was performed for the identification of sixteen respiratory virus using virus-specific primers as described previously [15]. Briefly, one tube was added with nine pairs of chimeric primers for the detection of nine viruses (HRV, FluA, FluB, 229E, OC43, HKU1, PIV1, PIV3 and ADV), and the other tube was added with eight pairs of chimeric primers to detect another eight viruses (RSVA, RSVB, NL63, PIV2, PIV4, H3N2, H1N1 and HBoV). The primer sequences, target genes, amplicon sizes and primer working concentrations were summarized in Additional file 1: Table 1. Qiagen one-step RT-PCR Kit (Qiagen, Hilden, German) was used for amplification. A final volume of 25 μL PCR mixture containing 2 μL of template, 1 μL of primer-mix (mixture was prepared according to the primer working concentration in Table 1), 1 μL of enzyme-mix, 1 μL of deoxinucleotidetriphospate (dNTP), 5 μL of 5X buffer, 1 μL of Taq polymerase was subjected to multiplex RT-PCR screening. The amplification cycle was initial denaturation at 50 °C for 30 min, 95 °C for 15 min, followed by 10 cycles of 95 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s; 10 cycles of 95 °C for 30 s, 65 °C for 30 s, 72 °C for 30 s; 25 cycles of 95 °C for 30 s, 48 °C for 30 s, 72 °C for 30 s and one final extension at 72 °C for 3 min. The amplified fragment (15 μL) was separated on the QIAxcel automatic electrophoresis using QIAxcel DNA High-Resolution kit and confirmed as specific virus types following reference marker in the system.

Table 1 Distribution of 16 respiratory virus in children according to gender and age

The specificity of the two-tube assay was based on the automatic electrophoresis. All positive controls were observed and separated clearly from the other viral targets and no primer dimer was observed in either tube. The expected size of each virus type/subtype-specific amplicon was listed in Additional file 1: Table 1. The detection sensitivity in tube 1 and tube 2 was 2000 and 200 copies per reaction, respectively. The detection limit for HRV, PIV2, PIV3, RSVA, HBoV and ADV was 20 copies per reaction and the remaining 10 viruses had a detection limit of 200 copies per reaction.

Detection of SARS-CoV-2 by RT-qPCR

SARS-CoV-2 nucleic acid detection kit (Daan, China) was applied to detect the ORF1ab gene and N gene of SARS-CoV-2 according to protocols as described previously [16]. Primers and probes used for SARS-CoV-2 were listed in Additional file 1: Table 1. The reaction was performed in a 25 μL reaction mixture containing 17 μL of reaction mixture A, 3 μL of reaction mixture B and 5 μL of RNA. The RT-qPCR reactions were performed with initial one cycle at 50 °C for 15 min and 95 °C for 15 min, following by 45 cycles of 94 °C for 15 s and 55 °C for 45 s. All RT-qPCRs were performed on the 7500 Fast Real-Time PCR system (Applied Biosystems) and all samples were analyzed in duplicate.

Sequencing of spike (S) gene

The S genes of HKU1, OC43, NL63 and 229E positive samples were amplificated using PrimeScript™ One Step RT-PCR Kit Ver.2 (TaKaRa 055RA). The volume of PCR mixture was 25 μL, containing buffer 12.5 μL, enzyme-mix 1 μL, forward primer 0.5 μL, reverse primer 0.5 μL, RNA 5 μL, and Rnase free water 5.5 μL. Their forward and reverse primers were listed in Additional file 1: Table 2. The amplification cycle was reverse transcription at 50 °C for 30 min and initial denaturation at 94 °C for 2 min, followed by 40 cycles of 94 °C for 45 s, 50–56 °C (depending on species) for 30 s, 72 °C for 45 s, and one final extension at 72 °C for 10 min at the end of 40 cycles. Amplified PCR products were sequenced by Sanger sequencing method (TsingKe Biological Technology).

Table 2 Viral detection rate according to sample type

Phylogenetic analysis

DNA sequences were concatenated manually and aligned by Clustal W and downloaded to MEGA 7.0 to generate a phylogenetic tree constructed by neighbor-joining algorithm and jukes-cantor model with 1000 bootstrap. Genome sequences of other coronaviruses were retrieved from GeneBank and included in phylogenetic analysis for comparison.

Statistical analysis

The statistical software SPSS (Statistics 20, IBM, Armonk, NY) was used for data processing and statistical analysis. Fisher’s Chi-square test was performed to test for statistically significant associations between prevalence of human coronaviruses to season, age, patient gender, other respiratory viruses, and clinical sample type. P values of < 0.05 were considered as significant. Analyses of variance (ANOVAs) were also performed to test the null hypothesis that there was no difference between clinical manifestations and HCoVs infection. If a null hypothesis was rejected at the 0.05 level, a Tukey’s multiple-comparison test was used to identify differences in clinical characteristics.


Prevalence of four human coronaviruses in Zhejiang Province, China

To survey human coronaviruses (HCoVs) in Zhejiang province, nasopharyngeal swabs (n = 26), sputum (n = 3593) and bronchoalveolar-lavage fluid (n = 41) specimens were collected periodically between February 2017 and December 2019 from children with respiratory symptoms. Overall, all clinical specimens were negative for SARS-CoV-2 RNA by RT-PCR analysis. 38 (1.04%) were positive for HCoV-HKU1, 62 (1.69%) for HCoV-NL63, 8 (0.22%) for HCoV-229E and 38 (1.04%) for HCoV-OC43 (Table 1). Comparison of the percentage of positive samples for multiple respiratory viruses was used to infer prevalence (Fig. 1). Respiratory syncytial viruses (RSVA + RSVB, 14.51%), human rhinovirus (HRV, 13.74%), parainfluenza (PIVs 1–4, 10.96%) and influenza (FluA + FluB, 6.61%) were the four dominant viruses in Zhejiang area, which were detected more frequently than HCoVs (P < 0.001).

Fig. 1

The detection rates of respiratory viruses

Circulation trend

Overall, HCoVs dominant circulation season was in the transition of summer to fall, especially August to October (Fig. 2). The percent of HCoV positive tests varied seasonally throughout each year and the percent positive also varied annually by HCoV species (Fig. 3a). The detection rates of HCoVs were 2.18%, 3.94%, 4.51%, 5.22% of winter, spring, summer and autumn, respectively (Fig. 3b). HCoV-HKU1 demonstrated a distinct peak each of the three years, with a relatively large peak in Spring 2018. For HCoV- NL63, in contrast to the pattern described for HCoV-HKU1, the peak detection frequency occurred in the third season. 51.61% of positive HCoV- NL63 samples were detected between September to November. HCoV-OC43 also demonstrated a consistent detection frequency between seasons, with the testing peak occurring each year in spring. The exception was HCoV-229E, as 87.50% of detectable HCoV-229E were found in the spring and summer of 2019. To be noted, HCoV-229E detected only in spring and summer through the study period, and no HCoV-229E was detected in 2018.

Fig. 2

Monthly distribution of total HCoVs (orange), HCoV-NL63 (blue), HCoV-OC43 (yellow), HCoV-HKU1 (green) and HCoV-229E (purple) infections in hospitalized children with respiratory tract infections from February 2017 to December 2019

Fig. 3

Circulation trends of HCoV infections. a Yearly distribution of the four HCoV infections. b Seasonal distribution of the four HCoV infections

Clinical profiles associated with HCoV infections

Almost all HCoV-positive patients (143/144, 99.31%) had lower respiratory tract infections. The incidence of lower respiratory tract infections in patients with HCoVs (3.94%) infection was significantly lower than that of patients infected with RSV (14.51%), PIV (10.96%), or HRV (13.85%). In addition, only one patient infected with HCoV-NL63 was diagnosed with an upper respiratory tract infection from nasopharyngeal swab in this study. The positive rate of HCoV infections was 3.85% (1/26) for nasopharyngeal swab, 4.88% (2/41) for bronchoalveolar lavage and 3.93% (141/3590) for sputum (P > 0.05) (Table 2). The four investigated HCoV subtypes were detected in sputum samples, while only HCoV-NL63 infections were positive in nasopharyngeal swab and bronchoalveolar lavage specimens.

Fever, cough, expectoration and sore throat were the common clinical presentations with HCoV infections, among which cough (47.22%), sore throat (40.97%) and fever (23.61%) were the top three manifestations (Fig. 4a). Difficulty in breathing and respiratory failure were seen in patients who infected with HCoV-OC43. Slightly more patients with HCoV infections were males but no significant difference was in the sex distribution between the four species (P > 0.2) (Table 1). Age distribution of patients differed between HCoV species and infants (< 1 year old) with weak respiratory immunity were most susceptible to four types of HCoVs infections (Fig. 4b). To be noted, HCoV-229E had a higher detection proportion in children between 2 and 18 years old compared to HCoV-OC43, HCoV-NL63, and HCoV-HKU1.

Fig. 4

Clinical profiles of HCoV infections based on clinical manifestation (a) and age (b)

Coinfection of respiratory viruses

High rates of coinfection with other respiratory viruses were observed for all coronaviruses, with the coinfection rates of HCoV-HKU1, HCoV-NL63, HCoV-OC43 and HCoV-229E as 44.74% (17/38), 43.55% (27/62), 52.63% (20/38) and 37.50% (3/8), respectively (Table 3). Among the 146 positive detections, 54.11% reported a single HCoV species detection only, 1.37% reported two HCoV species, and 44.52% detected another respiratory virus. The most common HCoV co-infections were HCoV-HKU1 with HCoV-229E (2 specimens, 1.37%). Of all samples testing positive for coronavirus, 40% samples were coinfected with HRV. The most common co-detected non-HCoVs were human adenovirus (40.71%), human rhinovirus/bocavirus (37.57%), and human metapneumovirus (32.43%).

Table 3 Coinfection of 16 respiratory virus

Phylogenetic analysis of HCoV strains

Phylogenetic trees were further constructed based on the sequencing of HCoV spike genes. A representative subset of HCoV-OC43 (n = 13), HCoV-NL63 (n = 8) and HCoV-NL63 (n = 8) were included in the phylogenetic analysis and compared to the HCoVs from clinical sources worldwide. A high level of genetic diversity was observed among those HCoVs. The OC43 coronaviruses were clustered into genotype G (9, 69.23%) and genotype B (4, 30.77%) and related to the sequences detected in Malaysia and Beijing, respectively (Fig. 5a). No significant differences were observed between HCoV-OC43 strains identified in different years. Similarly, HCoV-NL63 strains in this study were clustered into genotype B (5, 62.5%) and genotype C (3, 37.5%) and shown close identity to the viruses isolated from China and USA (Fig. 5b). Moreover, the phylogenetic tree of the HCoV-HKU1 isolates shown the existence of two major clusters that contained both Hong Kong and France strains of HCoV-HKU1 (Fig. 5c). All HCoV-HKU1 subtypes isolated from 2017 were grouped into genotype B.

Fig. 5

Phylogenetic analysis based on nucleotide sequences of PCR products corresponding to the spike genes of HCoV-OC43, HCoV-NL63 and HCoV-HKU1. a Phylogenetic trees of HCoV-OC43 S gene (4.7 kb) constructed with neighbor-joining algorithm. Red circles indicated Zhejiang HCoV-OC43 strains. b Phylogenetic trees of HCoV-NL63 S gene (4.5 kb) constructed with neighbor-joining algorithm. Blue triangles indicated Zhejiang HCoV-NL63 strains. c Phylogenetic trees of HCoV-HKU1 S gene (4.6 kb) constructed with neighbor-joining algorithm. Green rhombus indicated Zhejiang HCoV- HKU1 strains


Emerging and reemerging pathogens represent a serious burden to global health. The recent emergence of SARS-CoV-2 clearly demonstrated the epidemic potential of coronaviruses. Here we demonstrated that HCoV-OC43, HCoV-NL63, HCoV-HKU1 and HCoV-229E were specifically associated with lower respiratory tract disease in pediatric population, accounting for 3.93% of all admissions for acute respiratory infections. Only one patient was associated with upper respiratory traction infection from nasopharyngeal swab samples. The newly discovered SARS-CoV-2 RNA was absent in all tested samples, suggesting no introduction of this highly pathogenic coronavirus into the Zhejiang population before 2020. RSV (14.51%) was the predominant respiratory virus in our pediatric population. Our results aligned with studies conducted in other countries whereby RSV was the predominant viral agent among young children [17,18,19,20]. All the HCoVs detected have been subtyped. HCoV-NL63 (1.69%) was the most prevalent coronaviruses in our study consistent with reports in Hong Kong [21], Japan [22], Brazil [23] and Kenya [24], but some studies demonstrated the occurrence of HCoV-OC43 was similar or even higher than that of HCoV-NL63 in United States [9], United Kingdom [25] and Guangdong province [26].

In general, common HCoVs are reported to have annual peaks of circulation in winter months and few detections in summer months in temperate regions [9, 22, 25, 27]. Our data shown a completely different trend as we found HCoVs had a summer-autumn peak of activity, suggesting seasonality of HCoVs in subtropical regions was not restricted to the winter. Individual HCoV subtypes shown variable circulation throughout each year. During the sampling period of this study, HCoV-OC43 was the major circulating strain in 2017, whereas HCoV-NL63 predominated in the following two years. This was expected, given previous observations of annual variation in the circulation of HCoVs in Japan [22], United States [9] and Thailand [28]. HCoV-229E was only detected in spring and summer with a comparatively low detection frequency. This reaffirmed with the report in United Kingdom indicating that HCoV-229E was lacking a discernible seasonality [25]. However, continuous surveillance of HCoV infections over a long period is required to delineate their circulation trend more precisely.

This study investigated children between age below 18 years old. The burden of HCoV infections in children was high, especially among infants ((< 1 year old) because of their immature immune system. The age distribution of patients varied between HCoV species. Similar to the epidemiological feature in United States [9], HCoV-OC43 detections were most prevalent in infants, while HCoV-229E infections were most common among elder children. In concordance with a previous study [17], we have also found that boys were more vulnerable (4.02%) to HCoV infections than girls (3.87%). Specimen types are important for the diagnosis of respiratory viral infections. Although the overall detection rate of HCoVs was similar between sample types, sputum was superior to bronchoalveolar lavage and nasopharyngeal swab as all HCoV subtypes were detected from sputum. This finding was consistent with previous reports that sputum was more sensitive than nasopharyngeal swabs for respiratory viral detection [29, 30]. Although more than half one-third of the HCoV positive samples were co-detected with other viruses, our data support that HCoVs lead to a substantial burden of respiratory tract infections in children.

Phylogenetic analysis based on the spike gene confirmed previously published results and illustrated the genetic diversity of HCoV strains. Recombination occurred in HCoV frequently. The primary circulating genotypes of HCoV-OC43 in our study were B and G, while none of the genotype A, C, D, E and F were detected. This finding was consistent with Guangzhou strains, where their dominant circulating genotypes were B and G [31]. The phylogenic tree based on spike gene of HCoV-NL63 revealed the genotypes circulating in Zhejiang area were B and C, which was different from Kenya as their circulating genotypes were A and B [32]. However, this database was too small to reflect the actual variation level of HCoVs. Further research is needed to apply whole genome sequencing as the primary sequencing method for the genetic characterization of HCoVs.

To best of our knowledge, this comprehensive survey is the first to describe the clinical patterns and genetic characterization of the four major HCoVs and SARS-CoV-2 in Zhejiang province during a multi-year period. However, some limitations still existed to this survey. We only evaluated respiratory symptoms among hospitalized children, while the data from asymptomatic children or patients with mild respiratory symptoms were not included. Additionally, we assessed HCoV percent positivity based on the submission of specimen collection. Aggregate data might contain multiple specimens from the same patient, potentially affecting the demographic characteristics of our total study populations. Therefore, our data probably had some bias to represent the overall health burden of HCoVs infections of the general pediatric population in Zhejiang province, China.


Taken together, this study revealed that human coronaviruses were a significant cause of acute respiratory infections among hospitalized children, especially infants. Individual HCoVs may show variable circulation from year to year. As the global COVID-19 outbreak continues to grow, more investigations and rigorous surveillance are warranted to ascertain the circulation patterns and molecular epidemiology of human coronaviruses in China.

Availability of date and materials

The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.



Human coronaviruses


Ribonucleic acid


Severe acute respiratory syndrome coronavirus


Severe acute respiratory syndrome coronavirus 2


Middle east respiratory syndrome coronavirus


Human coronavirus NL63


Human coronavirus OC43


Human coronavirus HKU1


Human coronavirus 229E


Polymerase chain reaction


Respiratory syncytial viruses


Respiratory syncytial virus A


Respiratory syncytial virus B


Parainfluenza virus 1


Parainfluenza virus 2


Parainfluenza virus 3


Parainfluenza virus 4


Human rhinovirus


Influenza virus A


Influenza virus B


FluAInfluenza virus H3N2


Influenza virus H1N1


Human bocavirus




  1. 1.

    Viruses CSGotICoTo. The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2. Nat Microbiol. 04 2020;5(4):536–544.

  2. 2.

    Zhu N, Zhang D, Wang W, et al. A novel Coronavirus from patients with pneumonia in China, 2019. N Engl J Med. 2020;382(8):727–33.

    CAS  Article  Google Scholar 

  3. 3.

    Hui DS, Azhar E, Madani TA, et al. The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health—the latest 2019 novel coronavirus outbreak in Wuhan. China Int J Infect Dis. 2020;91:264–6.

    CAS  Article  Google Scholar 

  4. 4.

    Han Q, Lin Q, Jin S, You L. Coronavirus 2019-nCoV: a brief perspective from the front line. J Infect. 2020;80(4):373–7.

    CAS  Article  Google Scholar 

  5. 5.

    Drosten C, Günther S, Preiser W, et al. Identification of a novel coronavirus in patients with severe acute respiratory syndrome. N Engl J Med. 2003;348(20):1967–76.

    CAS  Article  Google Scholar 

  6. 6.

    Zaki AM, van Boheemen S, Bestebroer TM, Osterhaus AD, Fouchier RA. Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia. N Engl J Med. 2012;367(19):1814–20.

    CAS  Article  Google Scholar 

  7. 7.

    Su S, Wong G, Shi W, et al. Epidemiology, genetic recombination, and pathogenesis of Coronaviruses. Trends Microbiol. 2016;24(6):490–502.

    CAS  Article  Google Scholar 

  8. 8.

    van der Hoek L, Pyrc K, Berkhout B. Human coronavirus NL63, a new respiratory virus. FEMS Microbiol Rev. 2006;30(5):760–73.

    Article  Google Scholar 

  9. 9.

    Killerby ME, Biggs HM, Haynes A, et al. Human coronavirus circulation in the United States 2014–2017. J Clin Virol. 2018;04(101):52–6.

    Article  Google Scholar 

  10. 10.

    Klig JE, Shah NB. Office pediatrics: current issues in lower respiratory infections in children. Curr Opin Pediatr. 2005;17(1):111–8.

    Article  Google Scholar 

  11. 11.

    Griffin MR, Walker FJ, Iwane MK, et al. Epidemiology of respiratory infections in young children: insights from the new vaccine surveillance network. Pediatr Infect Dis J. 2004;23(11 Suppl):S188-192.

    Article  Google Scholar 

  12. 12.

    Jevšnik M, Uršič T, Zigon N, Lusa L, Krivec U, Petrovec M. Coronavirus infections in hospitalized pediatric patients with acute respiratory tract disease. BMC Infect Dis. 2012;12:365.

    Article  Google Scholar 

  13. 13.

    Dominguez SR, Robinson CC, Holmes KV. Detection of four human coronaviruses in respiratory infections in children: a one-year study in Colorado. J Med Virol. 2009;81(9):1597–604.

    Article  Google Scholar 

  14. 14.

    Williams PCM, Howard-Jones AR, Hsu P, et al. SARS-CoV-2 in children: spectrum of disease, transmission and immunopathological underpinnings. Pathology. 2020;52(7):801–8.

    CAS  Article  Google Scholar 

  15. 15.

    Li J, Qi S, Zhang C, et al. A two-tube multiplex reverse transcription PCR assay for simultaneous detection of sixteen human respiratory virus types/subtypes. Biomed Res Int. 2013;2013:327620.

    PubMed  PubMed Central  Google Scholar 

  16. 16.

    Zhu N, Wang W, Liu Z, et al. Morphogenesis and cytopathic effect of SARS-CoV-2 infection in human airway epithelial cells. Nat Commun. 2020;11(1):3910.

    CAS  Article  Google Scholar 

  17. 17.

    Khor CS, Sam IC, Hooi PS, Quek KF, Chan YF. Epidemiology and seasonality of respiratory viral infections in hospitalized children in Kuala Lumpur, Malaysia: a retrospective study of 27 years. BMC Pediatr. 2012;12:32.

    Article  Google Scholar 

  18. 18.

    Shafik CF, Mohareb EW, Yassin AS, et al. Viral etiologies of lower respiratory tract infections among Egyptian children under five years of age. BMC Infect Dis. 2012;12:350.

    Article  Google Scholar 

  19. 19.

    Kwofie TB, Anane YA, Nkrumah B, Annan A, Nguah SB, Owusu M. Respiratory viruses in children hospitalized for acute lower respiratory tract infection in Ghana. Virol J. 2012;9:78.

    Article  Google Scholar 

  20. 20.

    Calderaro A, De Conto F, Buttrini M, et al. Human respiratory viruses, including SARS-CoV-2, circulating in the winter season 2019–2020 in Parma, Northern Italy. Int J Infect Dis. 2021;102:79–84.

    CAS  Article  Google Scholar 

  21. 21.

    Chiu SS, Chan KH, Chu KW, et al. Human coronavirus NL63 infection and other coronavirus infections in children hospitalized with acute respiratory disease in Hong Kong, China. Clin Infect Dis. 2005;40(12):1721–9.

    Article  Google Scholar 

  22. 22.

    Matoba Y, Abiko C, Ikeda T, et al. Detection of the human coronavirus 229E, HKU1, NL63, and OC43 between 2010 and 2013 in Yamagata, Japan. Jpn J Infect Dis. 2015;68(2):138–41.

    Article  Google Scholar 

  23. 23.

    Cabeça TK, Granato C, Bellei N. Epidemiological and clinical features of human coronavirus infections among different subsets of patients. Influenza Other Respir Viruses. 2013;7(6):1040–7.

    Article  Google Scholar 

  24. 24.

    Sipulwa LA, Ongus JR, Coldren RL, Bulimo WD. Molecular characterization of human coronaviruses and their circulation dynamics in Kenya, 2009–2012. Virol J. 2016;13:18.

    Article  Google Scholar 

  25. 25.

    Gaunt ER, Hardie A, Claas EC, Simmonds P, Templeton KE. Epidemiology and clinical presentations of the four human coronaviruses 229E, HKU1, NL63, and OC43 detected over 3 years using a novel multiplex real-time PCR method. J Clin Microbiol. 2010;48(8):2940–7.

    CAS  Article  Google Scholar 

  26. 26.

    Liu P, Shi L, Zhang W, et al. Prevalence and genetic diversity analysis of human coronaviruses among cross-border children. Virol J. 2017;14(1):230.

    Article  Google Scholar 

  27. 27.

    Vabret A, Mourez T, Gouarin S, Petitjean J, Freymuth F. An outbreak of coronavirus OC43 respiratory infection in Normandy, France. Clin Infect Dis. 2003;36(8):985–9.

    Article  Google Scholar 

  28. 28.

    Dare RK, Fry AM, Chittaganpitch M, Sawanpanyalert P, Olsen SJ, Erdman DD. Human coronavirus infections in rural Thailand: a comprehensive study using real-time reverse-transcription polymerase chain reaction assays. J Infect Dis. 2007;196(9):1321–8.

    CAS  Article  Google Scholar 

  29. 29.

    Falsey AR, Formica MA, Walsh EE. Yield of sputum for viral detection by reverse transcriptase PCR in adults hospitalized with respiratory illness. J Clin Microbiol. 2012;50(1):21–4.

    CAS  Article  Google Scholar 

  30. 30.

    Jeong JH, Kim KH, Jeong SH, Park JW, Lee SM, Seo YH. Comparison of sputum and nasopharyngeal swabs for detection of respiratory viruses. J Med Virol. 2014;86(12):2122–7.

    CAS  Article  Google Scholar 

  31. 31.

    Zhang SF, Tuo JL, Huang XB, et al. Epidemiology characteristics of human coronaviruses in patients with respiratory infection symptoms and phylogenetic analysis of HCoV-OC43 during 2010–2015 in Guangzhou. PLoS ONE. 2018;13(1):e0191789.

    Article  Google Scholar 

  32. 32.

    Kiyuka PK, Agoti CN, Munywoki PK, et al. Human Coronavirus NL63 molecular epidemiology and evolutionary patterns in rural coastal Kenya. J Infect Dis. 2018;217(11):1728–39.

    CAS  Article  Google Scholar 

Download references


We greatly appreciate the support from the pathology department of Zhejiang Children Hospital for providing clinical specimens.


This work was funded by Major National Science and Technology Projects in the 13th Five-year Plan (2017ZX10103008-002), Health leading Talents Program of Zhejiang Province (CX-9), China Postdoctoral Science Foundation (2020T130104ZX) and Zhejiang Postdoctoral Science Foundation (ZJ2020027).

Author information




LF and LS drafted the manuscript. LS, SY and DZ involved sample collection, clinical data collection and virus detection. YZ, LF and LS participated in spike gene amplification. YZ and LF interpreted the surveillance and virological data. YZ designed the whole project. HM and YC provided important guidance. All authors reviewed and revised the first and final drafts of this manuscript. All authors read and agreed the final manuscript.

Corresponding author

Correspondence to Lei Fang.

Ethics declarations

Ethics approval and consent to participate

This study was ethically approved by Zhejiang Provincial Center for Disease Control and Prevention of ethics committee. Written informed consent was obtained from the guardians of all participants before the sample and data collection.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

Supplementary table 1. Sequences of oligonucleotides used for molecular analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhang, Y., Su, L., Chen, Y. et al. Etiology and clinical characteristics of SARS-CoV-2 and other human coronaviruses among children in Zhejiang Province, China 2017–2019. Virol J 18, 89 (2021).

Download citation


  • Human coronavirus
  • Hospitalized children
  • Phylogenic analysis
  • Acute respiratory infection
  • SARS-CoV-2