key: cord-0875451-m2b0huwk
authors: Jia, Ran; Lu, Lijuan; Li, Shu; Liu, Pengcheng; Xu, Menghua; Cao, Lingfeng; Su, Liyun; Xu, Jin
title: Human rhinoviruses prevailed among children in the setting of wearing face masks in Shanghai, 2020
date: 2022-03-14
journal: BMC Infect Dis
DOI: 10.1186/s12879-022-07225-5
sha: d05d42019a427b7e2c5207eb7bfb894417e79127
doc_id: 875451
cord_uid: m2b0huwk

BACKGROUND: Human rhinovirus (HRV) is the predominant etiological agent of the common cold in children and adults. A recent study showed that the inhibitory effect of face masks on viral shedding of HRV was less prominent than that on other respiratory viruses. Considering that most Chinese people have worn face masks in public area since the outbreak of coronavirus disease 2019, we aimed to find out whether HRV prevailed among children in 2020 and demonstrate the details of the epidemiological features of HRV under such a special circumstance. METHODS: We summarized the incidences of various respiratory virus infections in patients who visited the Children’s Hospital of Fudan University during 2018–2020, and genotyped HRV positive nasopharyngeal specimens collected from 316 inpatients and 72 outpatients that visited the hospital in 2020. RESULTS: There was a major prevalence of HRV among children in the latter half of 2020, with a clear seasonality that HRV-As prevailed in summer while HRV-Cs in autumn. HRV-As were more prone to cause severe lower respiratory tract infections (LRTI), while HRV-Cs were closely associated with childhood wheezing. The predominant genotypes were A11, A28, A47, A82, A101, C40 and C43. Notably, A21, A82 and A101 took up larger proportions in severe cases than in non-severe cases. CONCLUSIONS: Our findings described a major prevalence of HRVs among children in 2020, which highlight the unique transmitting pattern of HRV and help to narrow the targets for antiviral strategies.

Human rhinoviruses (HRVs) are the leading cause of upper respiratory tract infections (URTIs) since its first isolation in the 1950s [1] . HRVs also cause pneumonia hospitalization in vulnerable people such as children, the elderly and those with underlying diseases. HRVassociated diseases pose great socio-economic burdens to the country annually [2] . However, given that HRVinfected people are usually manifest self-limited and mild symptoms or even asymptomatic, HRVs have long been afforded little attention and no antivirals or vaccines have been approved for HRVs up to now [3] .

HRVs belong to the Picornaviridae family, and are single-stranded, positive-sense RNA viruses, indicating that it contains the sense strand of RNA as their genome which can be readily translated into proteins [4] . The genome is approximately 7,200 base pair (bp), including a single open reading frame (ORF) (~ 6500 bp), a 5' untranslated region (UTR) (~ 650 bp) and a 3′ UTR (~ 50 bp) [5] . About 100 serotypes which were culturable in vitro were classified into HRV-As and HRV-Bs based on the similarity of partial genetic sequences in the 1990s. Afterwards, at the beginning of the 2000s, researchers identified at least 50 more new HRV strains which couldn't be cultured and were classified into a unique species now named as HRV-C [2] . So far, more than 160 HRV genotypes have been identified [6] .

After the outbreak of coronavirus disease 2019 (COVID- 19) , most people have developed the habit of wearing face masks in public area in order to inhibit the transmission of respiratory pathogens. However, Leung et al. quantified the amount of respiratory viruses in exhaled breath of participants with acute respiratory illnesses and found that wearing medical face masks significantly reduced the RNA level of influenza viruses and coronaviruses (OC43 and NL63) in respiratory droplets or aerosols, but not in HRVs [7] , suggesting that the inhibiting effect of face masks may be less effective in HRV transmission. Hence, we conducted this research to further figure out whether HRVs could still spread among children in spite of the popularization of face masks and meanwhile demonstrate the details of the epidemiological features of HRVs. The findings in this study will expand the knowledge of HRV epidemiology and arouse people's attention to HRV's unique transmission pattern under such a special background.

A total of 316 nasopharyngeal aspirates from inpatients with lower respiratory tract infection (LRTI) hospitalized in the Children's Hospital of Fudan University in Shanghai from June 2020 to November 2020 were collected in this study. All the inpatients were diagnosed with LRTI supported by symptoms and radiographic changes and were defined as HRV positive after routine screening for common respiratory viruses including respiratory syncytial virus (RSV), adenovirus (AdV), influenza A and B viruses (IAV and IBV), parainfluenza virus type 1 (PIV-1), PIV-2, PIV-3, human rhinoviruses (HRV) and human metapneumovirus (MPV). For HRV screening, RNA from respiratory samples were extracted using a magnetic beads-based nucleic acid extraction system NP968-C (Tianlong Technology, China) according to the manufacturer's instruction. Then a one-step real time quantitative polymerase chain reaction (RT-qPCR) kit (Land medical, China) with primers targeting the 5′UTR (263 bp) of HRVs was used to detect HRV RNA. The remaining viruses and mycoplasma were detected using an immunofluorescence assay kit (Diagnostic Hybrids, USA). Briefly, nasopharyngeal aspirates were centrifuged and the cell pellet was fixed in acetone. A mixture of fluorescein-labeled monoclonal antibodies directed against the target viruses were added onto the cells, followed by an incubation of 30 min at 37 °C. A Mounting Fluid containing glycerol was added onto the stained cells and then a coverslip was placed on the prepared cells. The cells were examined using a fluorescence microscope (Nikon, Japan). Isolation and culture of bacterial and fungal pathogens were carried out according to the routine microbiology examination and diagnosis. Bacterial and fungal strains were identified using VIETEK automated bacterial analyzer (France) or MALDI-TOF/MS mass spectrometry (Bruck, France).

A total of 703 nasopharyngeal swabs from outpatients with URTI who visited the hospital during June 2020 to November 2020 were collected randomly and screened for HRV by RT-qPCR. The randomization was done as follows: first, one staff member covered all the information of the patients on the swabs with a blank tag paper. Then another staff member was asked to choose the swabs randomly to avoid biases in patients' gender, age, and illness.

LRTIs are illnesses that affect the respiratory system below the throat. The severity-based classification of the patients was performed by experienced clinicians according to the World Health Organization (WHO)'s latest definition of severe LRTI cases [8, 9] . Briefly, a child of any age with danger signs (e.g. cyanosis, seizures, lethargic/unconscious, unable to drink/breastfeed, respiratory failure) were defined as severe LRTI cases [8] [9] [10] . All experiments in the study were carried out in accordance with relevant guidelines and regulations. The study was reviewed and approved by the Ethics Committee of the Children's Hospital of Fudan University on Feb 2020 (Approval Number: 202027).

For genotyping, the extracted RNA were reverse transcribed and amplified using a nested RT-PCR strategy. HRV molecular subtyping was performed using primers targeting the VP4/VP2 regions (540 bp) of HRVs as reviewed in a previously published paper [11] . To increase both the sensitivity and efficiency of genotyping, we used a modified nested PCR method [12] . Briefly, the reverse transcription and the first amplification step were performed using a one-step RT-PCR kit (Rui'an Biotechnology, China) with outer primers: VP-OS (5′-CCG GCC CCT GAA TGY GGC TAA-3′) and VP-OAS (5′-ACA TRT TYTSNCCAAANAYDCCCAT-3′). The second amplification step was performed using a Premix Taq kit (Takara, Japan) with inner primers: VP-IS (5′-ACC RAC TAC TTT GGG TGT CCGTG-3′) and VP-IAS (5′-TCWG-GHARY TTC CAMCACCANCC-3′) [11, 13] . The amplification products were sequenced by Sangon Biotech Co., Ltd., China, followed by subjection to phylogenetic analysis using MEGA software.

Proportions for categorical variables were compared using the χ 2 test or Fisher's exact test. Independent group t-test was used for the comparison of means for continuous variables that were normally distributed. The Mann-Whitney U test was used for continuous variables not normally distributed. All statistical analyses were performed using GraphPad Prism software. Two-sided p-values of less than 0.05 were considered statistically significant.

A total of 4481 patients were tested positive for at least one of the following viruses as RSV, AdV, IAV, IBV, PIV-1, PIV-2, PIV-3, MPV and HRV in the Children's Hospital of Fudan University from June 2018, when the hospital started the HRV test for patients, to December 2020. Most respiratory viruses were barely detected after the outbreak of COVID-19 ( Fig. 1A) , but HRVs showed a remarkable increase in the middle of 2020 in June 2020 among the 3 years. The detection rate of HRV (HRV positive patients/total patients tested for the virus) also reached the climax (30.7%, 94/306) in June 2020. There was a mild increase in HRV infection in September 2020. PIV-3s and RSVs also increased gradually but were much less prominent than HRVs. The proportion of HRVs in the total virus positive cases was 52.5% in 2020, which was much higher than that in 2019 (26.7%) (Fig. 2B) . Collectively, these data indicated that there was a major prevalence of HRV in the year of 2020.

A total of 316 nasopharyngeal aspirates from HRV positive inpatients during the epidemics of HRVs in 2020 were collected and 82.3% (260/316) were successfully genotyped. We also randomly collected 703 nasopharyngeal swabs from outpatients with URTI during the same period of the inpatients, among which 10.2% (72/703) were determined as HRV positive and 90.3% (65/72) were successfully genotyped. The genetic variability of HRV genotypes in our data was very wide, as shown in the phylogenetic tree (Fig. 2) . A total of 29 HRV-A genotypes, 8 HRV-B genotypes and 22 HRV-C genotypes were detected in the patients. HRV-A was the most frequently detected species both in the inpatients (51.6%, 163/316) and the outpatients (45.8%, 33/72), followed by HRV-C (27.2%, 86/316 for inpatients; 33.3%, 24/72 for outpatients) and HRV-B (3.5%, 11/316 for inpatients; 11.1%, 8/72 for outpatients). The monthly distribution of patients revealed that HRV-As mainly prevailed in the summer (June to August) with A11, A47, A82 and A101 being the most frequent, while HRV-Cs quickly caught up in the autumn and peaked in September represented by C40 and C43. Notably, unlike most HRV-As, A28 mainly prevailed in the autumn (September to November) rather than summer (Fig. 3 ). HRVs infected more males than females both in the inpatients and the outpatients, with the male/female ratio being 1.34:1 and 1.4:1 respectively. 55.7% of total HRV-positive patients were infants under 1 year of age, and no discernable differences was found in the proportions of HRV species among different age groups (Fig. 4A) . Also, the Phylogenetic trees of VP4/VP2 gene sequences of HRVs. Phylogenetic trees were generated from manually trimmed 540 bp fragments using the neighbor-joining method and branch supported with 1000 bootstrap iterations using MEGA software. Bootstrap values were shown on tree nodes. Study sequences were identified by accession number. Reference sequences from GenBank were identified by accession number (HRV type) with black dots. Purple branches, HRV-As; Green branches, HRV-Bs; Red branches, HRV-Cs predominant genotypes appeared to be similar among children of different ages (Fig. 4B) .

To find out the association between HRV and clinical characteristics, we collected the clinical information of patients including symptoms, co-infections and underlying diseases. Among the three species, HRV-Bs seemed to infect more females than males, although the number of HRV-Bs were small (Table 1) . Notably, HRV-Cs appeared to be the most frequently detected species in the 46 patients with wheezing (60.9%, 28/46). 11.7% (23/196) of the HRV-A positive patients were defined as severe LRTI, which was much higher than HRV-B (0%, 0/19) and HRV-C (4.5%, 5/110), suggesting that HRV-A is more prone to cause severe illness.

Afterwards we classified the patients into three groups according to their disease severity, including the outpatients with URTI, the inpatients with non-severe LRTI or severe LRTI (Table 2) . No significant differences were found in the age distribution among the three groups. Females were more likely to develop HRV-associated severe LRTI, with 64.5% (20/31) severe cases being girls. In addition to longer hospitalizations, severe LRTI cases were more likely to have a cough than the other two groups. Notably, the percentages of underlying diseases increased in sequence from URTI (33.3%), non-severe LRTI (62.4%) to severe LRTI (87.1%).

The percentage of severe cases in HRV positive cases showed an increase in October and November 2020 (Fig. 5A) , despite that the detection rate of HRVs displayed a downward trend since September 2020 (Fig. 1A) , indicating that the detection rate and severity of HRV infections did not correspond completely. The genotypes detected in severe group included A11, A21, A28, A47, A82, A101, C40, C43 and C45, all of which were also detected in non-severe LRTI patients (Fig. 5B) . Notably, there were three genotypes which made up significantly larger proportions in severe LRTI cases than in non-severe LRTI cases, namely A21 (9.7% vs. 0.4%, p = 0.0004), A82 (16.1% vs. 5.3%, p = 0.0183) and A101 (22.6% vs. 7.7%, p = 0.0065). To be noted, A21 was the only genotype that was more frequently detected in severe LRTI cases (3 cases) than in non-severe LRTI cases (1 case), despite that the total Fig. 4 The distribution of HRV genotypes among patients of different ages. A The proportions of HRV species in different age groups from infants to adolescents. B The numbers of patients detected positive for the indicating genotypes were shown by age number of non-severe LRTI cases was almost ten times of the severe LRTI cases.

To further investigate the viral load-associated factors, we collected the Ct values of HRV positive samples. In our data, the patients co-infected with other respiratory viruses showed comparable viral loads with those with HRV mono-infection (Fig. 6A) . Also, the viral loads didn't seem to be correlated with the disease severity (Fig. 6B) . But HRV-Bs showed lower viral loads (higher Ct values) compared with the other two species (Fig. 6C) , which was in line with its lowest prevalence and least possibility to cause severe illness.

HRV infections were mainly transmitted via aerosols generated by coughing, sneezing and nose [14] , which is supposed to be effectively decreased by face masks. But the major HRV prevalence among children in 2020 indicates a weakened inhibitory effect of face masks [7] . Still, the unique transmitting pattern of HRV which enabled itself to escape from face masks deserves further investigation.

The majority of people wear disposable medical masks in public during the COVID-19 pandemic as WHO recommended, given that medical masks could help block large-particle droplets, splashes, sprays, or splatter that may contain viruses or bacteria [15] . But as we previously mentioned, the filtering effect of medical masks was insufficient to block HRV shedding [7] , and the increased HRV infection during the COVID-19 pandemic has been reported in various countries [16, 17] . N95 masks, which are class II medical devices, are designed to achieve a very close facial fit and very efficient filtration of airborne particles. Unlike medical masks, N95 masks could confer much better protection and have been proved to effectively block viruses like the influenza virus and HRV [18, 19] . Hence, it is advisable for HRV positive patients to wear N95 masks in order to reduce transmission.

Generally speaking, non-enveloped viruses (eg, HRV and AdV) are more heat-resistant and could survive Table 1 The clinical features of patients among the three HRV species HRV human rhinovirus, LRTI lower respiratory infection, RSV respiratory syncytial virus, AdV adenovirus *P < 0.05 was considered of significant difference among the three groups longer in a dry and acidic environment than enveloped respiratory viruses (eg, RSV, IAV, PIV, and CoV), which largely increases their chances of spread [20] . Also, it is very easy for children to touch contaminated surfaces/ objects (fomites). Children usually couldn't wash their hands timely and couldn't avoid close personal contact, which facilitate fomite-mediated viral transmission including HRV, enteroviruses, AdV, and rotavirus [21] . Hence, the spread of HRV among children might be attributed to both the reduced effect of face masks and children's uncontrolled behavior [22] . To be noted, RSV was reported to be the most common reason for LRTIassociated hospitalization in children less than 1 year of age, while HRV was reported to be the most common reason for LRTI-associated hospitalization in older children [20, 23] , which might be due to the limited independent activity of children under 1 year of age. Zhao et al. 's paper based on the respiratory samples of children in Shanghai during 2013-2015 shared some similar findings with ours, such as the age/gender preferences of HRV and the seasonality of HRV-C [24] . But HRVs were most frequently detected during winter in Zhao's paper but summer in ours. In Zhao's paper, the predominant genotypes included A78, A12, A89, A61, B70, C2, C6, C24 and C16, none of which were the main genotypes in our study. Moreover, we summarized the genotypes in papers focusing on various countries and concluded that the prevailing genotypes changed greatly with time and place [25] [26] [27] . But what these papers have in common was that HRV-As and HRV-Cs were the most frequently detected species and usually prevailed alternatively and seasonally. Considering of the substantial genetic diversity of HRVs, long-term and large-population-based studies are needed for a comprehensive understanding of HRV prevalence.

In consistence with our findings, several studies also reported that HRV-Cs are more commonly associated with early childhood asthma than the other two species [6, [28] [29] [30] , which might be attributed to the different cellular receptors. HRV-As and HRV-Bs use intercellular adhesion molecule 1 (ICAM-1) and the low-density-lipoprotein receptor (LDLR) for viral binding [31] [32] [33] , while Bochkov et al. found that HRV-Cs possibly use cadherinrelated family member 3 (CDHR3) for viral binding [34] . Notably, CDHR3 is a susceptibility locus for wheezing illness and early childhood asthma [34] . Hence, anti-childhood wheezing and subsequent asthma control strategies should pay more attention to HRV-Cs. It was reported that the patients co-infected with other respiratory viruses showed higher viral loads than those with HRV mono-infection [35] , but it is not the case in our data. Also, the viral loads wasn't correlated with the disease severity both in our study and other studies [36, 37] , while the rates of underlying diseases increased progressively with disease severity, suggesting that host factors bear important responsibility for the disease severity. Notably, the viral load in nasopharyngeal swabs may not reflect the viral load in the lower respiratory tract, and the relationship between viral load in lower respiratory tract and severe LRTI deserves further exploration. Lee et al. reported that the detection rate and severity of HRV infections did not correspond, and HRV-As and HRV-Cs were more likely to develop severe LRTI than HRV-Bs [38] , which is also the case in our data. A21 was more frequently detected in severe LRTIs than non-severe LRTIs and URTIs in our study, which is in line with the findings of a paper focusing on adults, although they didn't find specific site mutations in the sequences of A21 obtained from severe cases [39] . Whether there are particular A21 mutations that facilitate viral replication and host adaptation, especially in the lower respiratory tract tropism, deserves to be further demonstrated.

There is a growing understanding on the pathogenesis of viral and bacterial coinfections. For instance, viral infection in the respiratory tract could induce airway damage, promote bacterial adherence, decrease mucociliary clearance and impair the immune system, all of which facilitate bacterial co-infection [40, 41] . Conversely, primary bacterial infection may predispose to viral infections by facilitating viral propagation and infection within the respiratory system [40] . In terms of the host's factors, studies focusing on bacterial co-infections in COVID-19 patients found that advanced age and other comorbidities, such as chronic kidney disease, diabetes, and chronic heart disease, are associated with bacterial coinfections [42, 43] . But considering the small number of patients with bacterial co-infection (n = 22) in our study, we didn't analyze the risk factors of bacterial coinfection in HRV-positive patients, which is a limitation of the study. There are also other limitations in this paper. For example, our data only collected the samples from children in 2020, which makes us fail to compare the epidemiological features of HRV genotypes before and after the outbreak of COVID-19. Moreover, genetic analysis is needed to figure out whether there are meaningful site mutations in the prevailing HRV genotypes, such as A21, A82 and A101. More efforts are needed for better understanding of the individual and viral factors that contribute to more severe illnesses, so as to reduce the overall burden of respiratory illness.

Collectively, our findings described the details of the HRV prevalence among children in 2020, which is worthy of our reflection on the distinct transmitting pattern of HRVs. Moreover, our data suggested that the antiviral strategies to reduce HRV-related morbidity in high-risk children should focus on HRV-As and HRV-Cs. In a word, our findings add to the knowledge of the epidemiological features of HRV among children and underline the necessity to control HRV infection despite of the use of face masks. 

Rhinovirus biology, antigenic diversity, and advancements in the design of a human rhinovirus vaccine

The economic burden of non-influenza-related viral respiratory tract infection in the United States

The ABCs of rhinoviruses, wheezing, and asthma

Detection methods targeting the positive-and negative-sense RNA transcripts from plus-stranded RNA viruses

Rhinoviruses replicate effectively at lower airway temperatures

Wheezing rhinovirus illnesses in early life predict asthma development in high-risk children

Respiratory virus shedding in exhaled breath and efficacy of face masks

Integrated management of childhood illness: distance learning course. Geneva: World Health Organization

Early-life exposure to indoor air pollution or tobacco smoke and lower respiratory tract illness and wheezing in African infants: a longitudinal birth cohort study

Hospitalizations for severe lower respiratory tract infections

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

Detection of canine distemper virus (CDV) through one step RT-PCR combined with nested PCR

Human rhinovirus C infections mirror those of human rhinovirus A in children with community-acquired pneumonia

Susceptibility of an airborne common cold virus to relative humidity

Effectiveness of face masks in preventing airborne transmission of SARS-CoV-2. mSphere

Implications of the unexpected persistence of human rhinovirus/enterovirus during the COVID-19 pandemic in Canada. Influenza Other Respir Viruses

Increased risk of rhinovirus infection in children during the coronavirus disease-19 pandemic. Influenza Other Respir Viruses

Assessment of a respiratory face mask for capturing air pollutants and pathogens including human influenza and rhinoviruses

Potential utilities of mask-wearing and instant hand hygiene for fighting SARS-CoV-2

Rhinovirus and preschool wheeze

Transmission of SARS-CoV-2: a review of viral, host, and environmental factors

To mask or not to mask children to overcome COVID-19

Role of viral infections in the development and exacerbation of asthma in children

Genotypic diversity and epidemiology of human rhinovirus among children with severe acute respiratory tract infection in Shanghai

Molecular subtyping of human rhinovirus in children from three sub-Saharan African countries

Molecular epidemiology of human rhinovirus from 1-year surveillance within a school setting in rural coastal Kenya

Molecular epidemiology and clinical impact of rhinovirus infections in adults during three epidemic seasons in 11 European countries

Rhinovirus type in severe bronchiolitis and the development of asthma

Human rhinovirus species C infection in young children with acute wheeze is associated with increased acute respiratory hospital admissions

Rhinovirus C is associated with severe wheezing and febrile respiratory illness in young children

The major human rhinovirus receptor is ICAM-1

Members of the low density lipoprotein receptor family mediate cell entry of a minor-group common cold virus

Cadherin-related family member 3, a childhood asthma susceptibility gene product, mediates rhinovirus C binding and replication

A genome-wide association study identifies CDHR3 as a susceptibility locus for early childhood asthma with severe exacerbations

Viral loads and disease severity in children with rhinovirusassociated illnesses

Influenza and rhinovirus viral load and disease severity in upper respiratory tract infections

Initial high viral load is associated with prolonged shedding of human rhinovirus in allogeneic hematopoietic cell transplant recipients

Human rhinovirus species and season of infection determine illness severity

• fast, convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year submit your research ? Choose BMC

Viral and bacterial interactions in the upper respiratory tract

Viral-bacterial co-infections in the respiratory tract

Incidence of co-infections and superinfections in hospitalized patients with COVID-19: a retrospective cohort study

Empiric antibacterial therapy and community-onset bacterial coinfection in patients hospitalized with coronavirus disease 2019 (COVID-19): a multi-hospital cohort study

Publisher's Note

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

Not applicable. 

The authors have no competing interests.