key: cord-1014976-gxjm149n authors: da Costa Souza, Larissa; Bello, Edson José Monteiro; dos Santos, Eliane Maria; Nagata, Tatsuya title: Molecular and clinical characteristics related to rhinovirus infection in Brasília, Brazil date: 2021-01-06 journal: Braz J Microbiol DOI: 10.1007/s42770-020-00411-0 sha: 79ae950da2f803b8d35b7b99fda3210973368aaf doc_id: 1014976 cord_uid: gxjm149n INTRODUCTION: Human rhinovirus (HRV) is one of the most common human viral pathogens related to infections of the upper and lower respiratory tract, which can result in bronchiolitis and pneumonia. However, the relevance of HRV in human health was under-estimated for long time due to the absence of molecular targets for influenza and influenza-like syndrome surveillance in Brasília, Brazil. OBJECTIVES: The main objective of this study was analyze the clinical characteristics and outcomes of HRV infections in comparison with patients without HRV and other common respiratory viruses. MATERIALS AND METHODS: For this purpose, new specific primer sets were designed based on the high throughput sequencing analysis in previous study. These primers were used for HRV detection by RT-qPCR and Sanger sequencing of amplified cDNA of 5′ genomic region. The phylogenetic tree with representative HRV isolates was constructed using the Mega X software. Statistical analysis considering the patient profiles were performed using IBM SPSS program with non-parametric tests. RESULTS: The most prevalent virus in negative samples was rhinovirus (n = 40), including three rhinovirus species (rhinovirus A, B, and C). The odds ratio associated with HRV infection was 2.160 for patients younger than 2 years and 4.367 for people living in rural areas. The multiple analysis showed lower chance of patients with HRV presenting respiratory distress. CONCLUSION: In this study, it was reported the predominance of rhinoviruses in cases of respiratory illness for negative patients for the influenza and influenza-like syndrome surveillance, being rhinorrhea, the most significant symptom associated with the disease. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s42770-020-00411-0. Human rhinovirus (HRV) is one of the most common human viral pathogens. They are non-enveloped single-stranded R N A v i r u s , m e m b e r o f E n t e r o v i r u s g e n u s , o f Picornaviridae family [1, 2] . About half of common cold are caused by rhinoviruses, which have more than 160 serotypes described, in three species (rhinovirus A, B, and C) [3] . In addition, cold leads to economic problems related clinical morbidity, overloading the health system, also causing absences at work and in schools [2, 4] . Rhinoviruses cause not only upper respiratory tract infections but also lower respiratory tract infections, mainly in children and patients with chronic pulmonary disease or immunocompromised [5] . The most common symptoms of upper infection are sore throat, rhinorrhea, headache, cough, malaise, and lower fever. However, in some cases, there may be aggravation for rhinosinusitis and otitis media [6, 7] . Infections of the lower respiratory tract are often associated with bronchiolitis and pneumonia especially in children [8] [9] [10] . In Brazil, the implementation of the sentinel surveillance system began in 2000 [11] , so referral hospitals have been conducting active surveillance to detect respiratory viruses. Such surveillance includes notification and laboratory investigation of cases with the diagnostic criteria of influenza-like syndrome (ILS) and severe acute respiratory syndrome (SARS). This viral respiratory infection monitoring has resulted in important information about the circulation of community-acquired respiratory viruses. However, some viruses may not be identified by the absence of RT-qPCR target performed at central laboratories. The present study reports the molecular and clinical characteristics, and outcomes of HRV infections in comparison with patients without HRV and other common respiratory viruses. Detection of HRV by RT-qPCR The nasopharyngeal and tracheal secretion samples received in the laboratory in June, July, and August 2016 were tested for the respiratory virus diagnostic panel performed in Lacen-DF (Central Laboratory at Federal District, Brasilia, Brazil) by the standard RT-qPCR [12] (for influenza A virus, influenza B virus, human orthopneumovirus, human metapneumovirus, human mastadenovirus C, human respirovirus 1, human orthorubulavirus 2, and human respirovirus 3). The target genes of the standard RT-qPCR are shown in the Supplementary Table 1 . During transport, the samples were stored on ice and sent to the laboratory within 48 h after collection. Those that had negative results were separated and stored in freezer at − 70°C for subsequent DNA/RNA extraction. The metagenomic analyses by high throughput sequencing (HTS) of those samples as one pooled virome were previously studied [13] . Based on this study, HRV universal primers targeting conserved regions in 5' UTR (untranslated region) were newly designed for HRV detection (Table 1 ) using the Geneious R8.1 program (Biomatters, Auckland, New Zealand). The RNA was extracted from 145 samples using the Magna Pure LC Total Nucleic Acid Kit (Hoffmann-La Roche, Basel, Switzerland), and then all the 145 samples were projected to the RT-qPCR using GoTaq® Probe 1-Step RT-qPCR System (Promega, Madison, USA). Positivity for rhinovirus was defined with cycle threshold (CT) less than or equal to 40 (total of 45 cycles in the reaction) [12] . Human RNase P gene was used for endogenous control for the same sample set in separate tube. RT-PCR and sanger sequencing The HRV_Com_For and HRV_Com_Rev primers were newly designed (Table 1) targeting conserved regions based on virome analyses by HTS. This primer pair was used to amplify an approximately 540-bp fragment for amplicon sequencing, including part of 5' UTR and the VP4/VP2 protein gene of HRV A, B, and C. cDNA of the selected samples were synthesized using MMLV transcriptase (Thermo Fisher Scientific, Waltham, USA) with random hexamer primer and, then, amplified with LongAmp Taq DNA Polymerase (New England BioLabs, Ipswich, USA). The PCR product was purified and sequenced by the Sanger method at Macrogen Inc. (Seoul, South Korea). Sequence analysis The nucleotide (nt) sequences were determined using the Geneious R8.1 program (Biomatters, Auckland, New Zealand). Complete rhinovirus genome sequences were retrieved from GenBank as references and aligned with the sequences obtained in this study. This alignment was again trimmed to set the sequence size determined in this study (approximately 500 bases, the size has small variation due to indel mutation). The phylogenetic tree based on the 5' UTR and the VP4/VP2 protein gene was constructed using the Mega X software [14] . The tree was composed by eight nucleotide sequences from this study and other forty rhinovirus isolates, that were complete genome sequences available on databases. Phylogenetic tree was inferred using maximum likelihood method with Hasegawa-Kishino-Yano model as indicated by jModelTest analysis [15] . In order to obtain more precise nt identity in pairwise comparison, three more related virus sequences for each sequence in this study were retrieved by BlastN search.The pairwise comparison was then performed. Statistical analysis Medical records were obtained from compulsory notification sheet for ILS and SARS sent with the sample to the laboratory. Data were analyzed using the software IBM SPSS (Chicago, USA). Non-parametric tests were used as appropriate, for qualitative and quantitative variables. All p values were two tailed and a value of < 0.05 was considered significant. Rhinoviruses were identified as the main virus found in negative samples for the regular respiratory virus survey in Lacen- (Table 1) , which is most conserved among rhinoviruses. To identify HRV at species level, primers which target a region of approximately 540 bp of the 5' UTR-VP4/VP2 protein gene region were sequenced. Five samples were concluded of HRV-A, one of HRV-B and two of HRV-C ( Fig. 1) , out of eight sequenced samples of the patients in the Federal District during the study period. The other samples were not of sufficient quality for sequencing. All sequences were deposited in the databank: rhinovirus A BSB_8 (Accession number LC549198), rhinovirus A BSB_16 (LC549199), rhinovirus A BSB_26 (LC549200), rhinovirus A BSB_33 (LC549201), rhinovirus A BSB_38 (LC549202), rhinovirus B BSB_30 (LC549203), rhinovirus C BSB_2 (LC549204), and rhinovirus C BSB_15 (LC549205). It is interesting to note that although a small number of samples were sequenced, the genetic variability of HRV-A was very wide, positioning in five different clusters in the phylogenetic tree ( Fig. 1 ). Among them, only nt identities of 72.94-83.86% were observed ( Table 2) . On the other hand, two isolates of HRV-C were very similar, possessing 98.99% of nt identity. When identifying the location of the most related virus (in the databases), it can be observed that there is almost no relation in the geographical distribution ( Table 3 ). The most related rhinovirus isolates came from several continents such as America, Asia, Europe, and Oceania. This observation suggests that the dissemination of the rhinoviruses can happen worldwide in long distance and is different from the arbovirus geographical dissemination. Patients with (n = 40) and without (n = 105) rhinovirus were compared in relation to socio-demographic and clinical parameters in order to evaluate possible risk factors associated with the involvement of these viruses. It can be observed in Table 4 that patients with rhinovirus were significantly different from patients negative for the main respiratory viruses in relation to age (P = 0.042) and to the area of residence (P = 0.019). Patients younger than 2 years presented 2.160 times more likely to be infected by rhinovirus than patients with 2 years old or more. In this study, patients with rhinovirus were significantly younger than patients without rhinovirus (P = 0.042) ( Table 5 ). The median age for patients with the virus was 2.50 (interquartile range = 27.38) versus 26.00 (interquartile range = 50.13) years for patients without this virus. Thus, greater susceptibility of children to rhinovirus infection is reported. Figure 2 shows the difference between the groups with and without the virus in relation to the days of symptoms and age. It is observed that the days of symptoms did not differ significantly, but the age distribution of the patients was significantly lower in patients who presented rhinovirus respiratory symptoms (P = 0.042). Patients with rhinovirus were located significantly more in rural areas when compared to negative patients. The odds ratio for the rural area was 4.367 (1/0.229), which means that people resident in rural areas were 4.367 times more affected by rhinovirus when compared to urban dwellers. Regarding the main symptoms of patients with acute respiratory infection, it can be observed in Table 6 that patients with rhinovirus presented significantly more rhinorrhea in relation to the patients without this virus detected (P = 0.001). Patients with rhinovirus were 3.981 times more likely to presenting rhinorrhea in relation to negative patients. There was no statistically significant difference in relation to the other symptoms. A multiple logistic regression with 20 explanatory variables (age, days of symptoms, month, diagnostic criteria, sex, breed, zone, risk factor, vaccinated, hospitalization, medication, use of ventilatory support, fever, cough, sore throat, dyspnea, respiratory distress, myalgia, O 2 saturation < 95%, and rhinorrhea) was performed to analyze the variables associated with rhinovirus infection (response variable). In the model containing only the constant, there was 71% agreement with the classification table. When entering the explanatory variables, the concordance increased to 79.8%. The regression model was statistically significant (P = 0.005), with adequate adjustments (P = 0.862-Hosmer's test and Lemeshow) and R 2 = 0.391. To improve the model, a technique of variable selection was used by the backwise stepwise method. The best fit was obtained after 15 steps of exclusion of variables, remaining six explanatory variables. The new model was also statistically significant (P < 0.001), with optimal adjustment (P = 0.923-Hosmer and Lemeshow test) and R 2 = 0.319. Three significant explanatory variables were obtained: respiratory distress, rhinorrhea, and medication (Table 7) . Patients with rhinovirus presented 2.985 (1/0.335) times less chance of experiencing respiratory distress, 6.601 and 3.892 times more chance of having rhinorrhea and using medication (oseltamivir), respectively. Thus, it is observed that the main symptoms related to rhinovirus infection were the presence of rhinorrhea and absence of respiratory discomfort. The use of antiviral medication, oseltamivir, was also significantly more observed in patients with rhinovirus in the multiple analysis. To assess how respiratory symptoms were associated in patients with rhinovirus infection, a cluster analysis was performed with the 40 positive patients, using dendrogram with Ward method and measure of the quadratic Euclidean distance. It is observed that two groups of symptoms were mainly formed: fever, cough, and rhinorrhea in the first group and myalgia, O 2 saturation < 95%, respiratory distress, dyspnea, and sore throat in the second group (Fig. 3) . The symptoms of the first group may be associated with cases of common cold (upper airways infection), and those of the second group with cases of complications of infection in the lower respiratory tract. To investigate whether age was not a covariate of rhinorrhea, the Mantel-Haenszel test was performed, which evaluated the association between rhinorrhea symptom and age-corrected rhinovirus. It was observed that the association between rhinorrhea and rhinovirus occurred in the group of patients older than 2 years (P = 0.001), but not in the group younger than 2 years (P = 1,000) (Supplementary Table 2) ; that is, children with and without rhinovirus presented Since the beginning of the monitoring of respiratory virus circulation in the country, rhinovirus has been found with high frequency, either alone or co-detected with other respiratory viruses, mainly in the south and southeast regions [16, 17] . To date, however, the circulation of these viruses in the central region of the country has not been well studied. Therefore, the identification of these viruses is crucial to determine the etiologic agent of respiratory illness. In this case, the 5' UTR region is the most suitable to identify the virus by RT-qPCR, and this region has already been used by other laboratories with similar oligonucleotides [18, 19] . Rhinoviruses are reported to be the main cause of the common cold and spread from person to person via direct and indirect contact, mainly by infected respiratory secretions. The transmission by hands is important in this context [20] . The easy contagion of rhinovirus can explain its rapid spread, over long distances, as shown in Table 3 . HRV infections were associated with diseases such as asthma, chronic pulmonary disease, bronchiolitis, and pneumonia in children, the elderly, and immunocompromised. Thus, proper clinical management, timely diagnosis, and containment of outbreaks are becoming important [5, 16] . Rhinovirus infection was significantly associated with residency in rural areas. In this study, 7% came from rural areas, and among those infected by the virus, the percentage was 15%. A study in Brazil has already reported that children from urban areas of school age are more susceptible to respiratory inflammatory problems compared to those in rural regions [21] . Thus, the negative cases for rhinovirus and other common respiratory viruses may be associated with allergic inflammatory processes, which is more frequent in urban areas. With this, rural patients would present respiratory problems related mainly to viral infections, as observed in relation to HRV. In the multiple analysis, the relationship between rural residence and rhinovirus was not observed. Patients younger than 2 years were more subject to HRV infection and had a significantly lower proportion of cases with fever, and more comorbidities such as asthma than patients without HRV [16, 22, 23] . In this study, it was observed that children younger than 2 years were more susceptible to rhinovirus infection, but there was no association with comorbidities. The main symptom associated with HRV was rhinorrhea rather than low frequency of fever. The HRV infections predominated during autumn and winter with significant negative correlation between the number of HRV cases and the average temperature [16] . So, the result obtained in this study, with rhinovirus dominance in negative samples for the other common respiratory viruses, was expected since the samples were collected in late fall and winter (June, July, and August), when the virus circulation is high. The data here presented showed that patients with rhinovirus were significantly younger than patients without rhinovirus, consistent with other works that showed the detection rate of HRV infection in children was significantly higher than that in adults, mainly the HRV-C. Clinical manifestations associated with HRV-C appear to be more severe in children [24] and often associated with respiratory complications, like acute wheezing illness and bronchiolitis [25, 26] . In adults, HRV-A is more common and is associated with cases of viral pneumonia [27] . Fever and nasal congestion were more frequently observed in pediatric group than the adults which were most diagnosed with severe respiratory disease entity [28] . Figure 3 presented two groups of symptoms. This division may be associated with the main affected area, upper or lower respiratory tract. Rhinorrhea was the clinical manifestation significantly associated with HRV infection in this study. Other studies also show rhinorrhea as the main symptom associated with rhinovirus infection [29, 30] . The use of oseltamivir was significantly associated with HRV infection in the multiple analysis (Table 7) . In univariate analysis, this association was not verified (Table 4 ). However, a high percentage (37.5%) of patients with rhinoviruses who took the medication incorrectly were observed. Although rhinovirus and influenza symptoms can be similar and therefore confounded [31] , the use of oseltamivir in patients with HRV has already been reported, but in a much smaller percentage (2.1-2.5%) [22, 31] . This study showed the predominance of rhinoviruses in cases of respiratory illness for negative patients for the influenza and influenza-like syndrome surveillance. Rhinorrhea was the most significant symptom associated with the disease. All three virus species were identified in the Federal District, Brazil. A lower chance of presenting respiratory distress was also observed in the multiple analysis, although the symptoms of patients with HRV were associated with both upper airways infection and complications of infection in the lower respiratory tract (Fig. 3) . In the patients with HRV, 55.0 and 37.5% presented exacerbations such as dyspnea and respiratory distress, and 57.5% had to be hospitalized (Tables 4 and 6 ). As HRV is also identified in asymptomatic patients [32, 33] , more detailed studies need to be made to better understand the prevalence, risk factors, and outcomes related to rhinovirus infection in the population. The online version contains supplementary material available at https://doi.org/10.1007/s42770-020-00411-0. Authors' contributions All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by all the authors. The first draft of the manuscript was written by Larissa da Costa Souza and Tatsuya Nagata and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript. Funding This work was supported by grant from CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico) with the project number of 309670/2018-7. Availability of data and material The data for this study were made available by the Health Department of the Federal District through a project partnership and approved by an ethics committee. All genetic sequences of rhinovirus were deposited in the databank. Conflict of interest The authors declare that they have no conflict of interest. Ethics approval This research project was reviewed and approved by ethics committee of the Faculty of Health Sciences (University of Brasília), approval letter No. 3.052.443. Consent to participate This study obtained waiver of the informed consent form. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Sequencing and analyses of all known human rhinovirus genomes reveal structure and evolution The ABCs of rhinoviruses, wheezing, and asthma Analysis of the complete genome sequences of human rhinovirus The economic burden of non-influenza-related viral respiratory tract infection in the United States Human rhinoviruses Virological and serological analysis of rhinovirus infections during the first two years of life in a cohort of children Physiologic abnormalities in the paranasal sinuses during experimental rhinovirus colds Prospective genotyping of human rhinoviruses in children and adults during the winter of 2009-2010 Respiratory viruses in adults with community-acquired pneumonia Human rhinovirus C-associated severe pneumonia in a neonate O desafio da influenza: epidemiologia e organização da vigilância no Brasil. Boletim Eletrônico Epidemiológico da Secretaria de Vigilância em Saúde CDC protocol of realtime RTPCR for influenza A (H1N1) Human virome in nasopharynx and tracheal secretion samples MEGA X: molecular evolutionary genetics analysis across computing platforms jModelTest: phylogenetic model averaging Impact and seasonality of human rhinovirus infection in hospitalized patients for two consecutive years Viral etiology among the elderly presenting acute respiratory infection during the influenza season RespiFinder: a new multiparameter test to differentially identify fifteen respiratory viruses multi-centre study in the Netherlands examining laboratory ability to detect enterovirus 68, an emerging respiratory pathogen Hand-to-hand transmission of rhinovirus colds Inflammatory mediators in nasal lavage among school-age children from urban and rural areas in São Paulo Clinical characteristics and outcomes of human rhinovirus positivity in hospitalized children Association between human rhinovirus C and severity of acute asthma in children Human rhinovirus species and season of infection determine illness severity New respiratory enterovirus and recombinant rhinoviruses among circulating picornaviruses Genetics, recombination and clinical features of human rhinovirus species C (HRV-C) infections; interactions of HRV-C with other respiratory viruses Human rhinovirus C: age, season, and lower respiratory illness over the past 3 decades Molecular epidemiology and clinical features of rhinovirus infections among hospitalized patients in a medical center in Taiwan Prospective assessment of rhinovirus symptoms and species recurrence in children with and without an acute wheezing exacerbation Kalgoorlie Otitis Media Research Project Team (2012) Prevalence of and risk factors for human rhinovirus infection in healthy aboriginal and nonaboriginal Western Australian children Role of picornaviruses in flu-like illnesses of adults enrolled in an oseltamivir treatment study who had no evidence of influenza virus infection Code availability Not applicable.