Rhinovirus C and Respiratory Exacerbations in Children with Cystic Fibrosis

To investigate a possible role for human rhinovirus C in respiratory exacerbations of children with cystic fibrosis, we conducted microbiologic testing on respiratory specimens from 103 such patients in São Paulo, Brazil, during 2006–2007. A significant association was found between the presence of human rhinovirus C and respiratory exacerbations.

To investigate a possible role for human rhinovirus C in respiratory exacerbations of children with cystic fi brosis, we conducted microbiologic testing on respiratory specimens from 103 such patients in São Paulo, Brazil, during 2006-2007. A signifi cant association was found between the presence of human rhinovirus C and respiratory exacerbations.
C ystic fi brosis (CF) is an autosomal inherited disease characterized by recurrent and chronic respiratory infections that ultimately lead to the need for a lung transplant early in life or to death (1). The role of bacterial infections in CF is well established, and most treatments focus on eradication or suppression of bacterial infections (mainly those caused by Pseudomonas aeruginosa) (1).
Respiratory viruses such as respiratory syncytial virus (RSV) and infl uenza also seem to cause early damage or increase the risk for respiratory exacerbations (2,3) in these patients. However, the role of newly described respiratory viruses is not well known. Infection of these patients with human rhinovirus (HRV), a member of the family Picornaviridae, has been described. Although some studies have suggested a substantial pathogenic role for these viruses (3), controversy still exists (4).
Recently a new clade of human rhinovirus, named rhinovirus C (5), was identifi ed through molecular methods. This new clade has been found throughout the world (6), and some studies have attributed severe respiratory infections in children to this agent (7,8). We investigated whether this agent played a role in the respiratory exacerbations of children and adolescents with CF who attended the Instituto da Crinaça, Hospital das Clínicas da Faculdade de Medicina da Universidade de São Paulo.

The Study
A total of 103 CF patients (49 girls, 54 boys; median age 8.9 years; age range 3.8 months-17.8 years) were enrolled in the study from September 6, 2006 through September 4, 2007. Nasopharyngeal aspirates or nasal mucus specimens for viral investigation, as well as sputum or oropharyngeal samples for microbiology culture, were collected during scheduled visits or unscheduled visits on 408 occasions, with a median ± SD of 4 ± 1.74 visits per patient (range 1-9 visits).
Clinical and lung function data were obtained at all visits. Exacerbation of respiratory disease was defi ned as the presence of >2 of the following signs or symptoms: fever, increase in the amount of secretion or cough intensity, change in sputum's color, worsening of dyspnea, loss of appetite, a decrease of forced expiratory volume in 1 s >10%, and weight loss.
Total nucleic acids were extracted from nasopharyngeal samples by using a QIAmp Viral RNA Mini Kit (QIA-GEN, Hamburg, Germany), according to manufacturer's instructions. Reverse transcription was conducted with High Capacity cDNA Archive Kit (Applied Biosystems, Foster City, CA, USA) by using 20 μL of the previously extracted RNA. Respiratory viruses were identifi ed by individual reverse transcription-PCRs or PCRs selective for RSV; infl uenza viruses A and B; human parainfl uenza viruses 1, 2, and 3; human coronavirus; human metapneumovirus; adenovirus; human bocavirus; picornavirus; and the β-actin gene (9). For picornavirus, we used the primer pair OL26-OL27, which include a portion of the 5′ noncoding region (NCR) common to all picornaviruses, in the same conditions previously described (10). After sequence amplifi cation, products were examined by capillary electrophoresis in an automated DNA sequencer (MegaBace, General Electric Healthcare-Amersham Biosciences, Little Chalfont, UK), and results were visualized through the MegaBace FragmentProfi ler software, which discriminates fragment sizes and fl uorescent intensities.
Samples in which picornavirus cDNA had been identifi ed were submitted to a TaqMan-based real-time PCR protocol (11) to identify HRV and enterovirus. Only samples in which rhinovirus had been identifi ed by real-time PCR were submitted to 5′ ΝCR sequencing with an ABI 377 automated sequencer (Applied Biosystems), and results were submitted to the GenBank database, accession nos. GU933027-GU933118. Sequencing of the 5′ NCR region has been shown to accurately discriminate among subtypes of rhinovirus, including genotypes C and A2 (12). All sequencing chromatograms obtained were edited manually to obtain contiguous fragments (contigs), by using Sequence Navigator software (Applied Biosystems). All sequences were screened at the National Center for Biotechnology Information website by using the Basic Local Alignment Search Tool (BLAST) (http://blast.ncbi.nlm.nih.gov/Blast. cgi). HRV genotype was confi rmed by phylogenetic analysis as described below. Sequences generated and a set of reference strains representative of HRV genotypes, available in the GenBank database, were aligned by using ClustalW (www.ebi.ac.uk/Tools/clustalw2). Minor manual adjustments were made to improve the alignment with BioEdit software (www.mbio.ncsu.edu/BioEdit/BioEdit.html). The sequence of echovirus 11 was used as the outgroup. Phylogenetic analysis was performed with PAUP* version 4b10 (http://paup.csit.fsu.edu). Neighbor-joining and maximumlikelihood trees were constructed on the basis of appropriate nucleotide substitution models determined by Modeltest v3.7 (University of Vigo, Vigo, Spain). Bootstrapping was assessed by using 1,000 replicates. Trees were visualized by using the TreeView program (www.taxonomy.zoology. gla.ac.uk/rod/treeview.html).
To account for correlations among samples from the same patient, we used binomial generalized linear models to identify the virologic variables associated with the main endpoints (respiratory exacerbation and hospital admission). Results were presented as odds ratios (ORs) and 95% confi dence intervals.
At least 1 respiratory virus was identifi ed in 203 (49.8%) of 408 samples; rhinovirus was the main identifi ed agent (139 samples, 34.1%). The results of virus identification, in relation to clinical status, are shown in Table 1. Co-infections were found in only 26 samples (6.4%); rhinovirus was the most frequent agent.
Sequencing was performed in 93 of 139 samples that were positive for rhinovirus and showed a predominance of genotype A (36 samples; 38.7%) (Figure 1). Rhinovirus subtypes A2 and C were identifi ed in 14 samples each. Isolates in 3 samples were identifi ed as HRV87, coxsackievirus, and echovirus. Therefore, 5 isolates were enteroviruses, which indicates that the TaqMan-based real-time PCR platform misidentifi ed them. Sequencing was not successful in 1 sample. A maximum-likelihood phylogenetic tree of our samples and reference HRV sequences was constructed ( Figure 2).
Patients were examined during acute exacerbation of respiratory disease on 142 occasions, and patients required hospital admission on 31 occasions. The identifi cation of respiratory viruses was not associated with pulmonary exacerbations. Because rhinovirus was the main agent identifi ed among stable patients, we verifi ed the effects of respiratory viruses, excluding rhinovirus from the analysis. A signifi cant association with respiratory exacerbation was found (OR 1.195, p = 0.010) ( Table 2).
In contrast, when looking at rhinovirus subtypes, we noticed that identifi cation of rhinovirus subtypes A2 or C was also signifi cantly associated with respiratory exacerbations (OR 1.213) ( Table 2). Identifi cation of infl uenza A was the only variable associated with an increase in the risk for hospital admission (OR 1.988, 95% confi dence interval 1.238-3.194, p = 0.004).

Conclusions
These new HRV genotypes were initially described in samples from patients with infl uenza-like illnesses (13). Recently, evidence has been increasing for the involvement of HRV-C in severe respiratory conditions such as bronchiolitis in infants (7) and in exacerbation of asthma (14). However, these studies that attempted to clarify the patho-genicity of these new species did not include control groups or nonsymptomatic persons.
In the study reported here, we obtained samples from patients during routine visits and exacerbations, which enabled us to identify a distinct role of different HRV subtypes. Our fi ndings, however, cannot be extrapolated to infants with CF, which were underrepresented in our study. These infants may be at greater risk of HRV-C respiratory infections as they are for with RSV infections (2). Previous studies of respiratory virus infections in CF patients provided confl icting results on the potential effect of rhinovirus. Smyth et al. (15), who evaluated respiratory exacerbations in CF patients, identifi ed rhinovirus as the leading agent, with exacerbations following a more benign course but resulting in greater use of intravenous antimicrobial drugs (15). Olesen et al. (4) studied 75 CF patients for 1 year and used PCR to investigate the presence of 7 different viruses in sputum or laryngeal aspirates. The group reported that viral infections did not reduce lung function or increase respiratory symptoms. Rhinovirus was by far the leading agent identifi ed throughout their study (4). A more recent study by Wat et al. (3), in which they used real-time nucleic acid sequence-based amplifi cation in conjunction with molecular markers to investigate the presence of 9 respiratory viruses in children with CF, described an association of viral infections with respiratory exacerbations, particularly those caused by infl uenza A, infl uenza B, and rhinovirus (3).
We report infections by the novel rhinovirus subtypes A2 and C in CF patients. A signifi cant association was found between the presence of these agents and respiratory exacerbations. Our fi ndings indicate the need for further investigation of HRV-C in CF patients.