Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Rhinovirus Genome Variation during Chronic Upper and Lower Respiratory Tract Infections

  • Caroline Tapparel ,

    caroline.tapparel@hcuge.ch

    Affiliations Laboratory of Virology, Division of Infectious Diseases and Division of Laboratory Medicine, University of Geneva Hospitals, Geneva, Switzerland, University of Geneva Medical School, Geneva, Switzerland

  • Samuel Cordey ,

    Contributed equally to this work with: Samuel Cordey, Thomas Junier

    Affiliations Laboratory of Virology, Division of Infectious Diseases and Division of Laboratory Medicine, University of Geneva Hospitals, Geneva, Switzerland, University of Geneva Medical School, Geneva, Switzerland

  • Thomas Junier ,

    Contributed equally to this work with: Samuel Cordey, Thomas Junier

    Affiliations Department of Genetic Medicine and Development, University of Geneva Medical School, Geneva, Switzerland, Swiss Institute of Bioinformatics, Geneva, Switzerland

  • Laurent Farinelli,

    Affiliation Fasteris SA, Plan-les-Ouates, Switzerland

  • Sandra Van Belle,

    Affiliations Laboratory of Virology, Division of Infectious Diseases and Division of Laboratory Medicine, University of Geneva Hospitals, Geneva, Switzerland, University of Geneva Medical School, Geneva, Switzerland

  • Paola M. Soccal,

    Affiliation Division of Pulmonary Medicine and Clinic of Thoracic Surgery, University of Geneva Hospitals, Geneva, Switzerland

  • John-David Aubert,

    Affiliation Division of Pulmonary Medicine, University of Lausanne Hospitals, Lausanne, Switzerland

  • Evgeny Zdobnov,

    Affiliations Department of Genetic Medicine and Development, University of Geneva Medical School, Geneva, Switzerland, Swiss Institute of Bioinformatics, Geneva, Switzerland

  • Laurent Kaiser

    Affiliations Laboratory of Virology, Division of Infectious Diseases and Division of Laboratory Medicine, University of Geneva Hospitals, Geneva, Switzerland, University of Geneva Medical School, Geneva, Switzerland

Abstract

Routine screening of lung transplant recipients and hospital patients for respiratory virus infections allowed to identify human rhinovirus (HRV) in the upper and lower respiratory tracts, including immunocompromised hosts chronically infected with the same strain over weeks or months. Phylogenetic analysis of 144 HRV-positive samples showed no apparent correlation between a given viral genotype or species and their ability to invade the lower respiratory tract or lead to protracted infection. By contrast, protracted infections were found almost exclusively in immunocompromised patients, thus suggesting that host factors rather than the virus genotype modulate disease outcome, in particular the immune response. Complete genome sequencing of five chronic cases to study rhinovirus genome adaptation showed that the calculated mutation frequency was in the range observed during acute human infections. Analysis of mutation hot spot regions between specimens collected at different times or in different body sites revealed that non-synonymous changes were mostly concentrated in the viral capsid genes VP1, VP2 and VP3, independent of the HRV type. In an immunosuppressed lung transplant recipient infected with the same HRV strain for more than two years, both classical and ultra-deep sequencing of samples collected at different time points in the upper and lower respiratory tracts showed that these virus populations were phylogenetically indistinguishable over the course of infection, except for the last month. Specific signatures were found in the last two lower respiratory tract populations, including changes in the 5′UTR polypyrimidine tract and the VP2 immunogenic site 2. These results highlight for the first time the ability of a given rhinovirus to evolve in the course of a natural infection in immunocompromised patients and complement data obtained from previous experimental inoculation studies in immunocompetent volunteers.

Introduction

Human rhinoviruses (HRV) usually cause self-limited upper respiratory tract (URT) illness. However, they are increasingly reported to be associated with complications, such as asthma [1], [2], [3], chronic obstructive pulmonary disease (COPD) exacerbations [4], pneumonia, and bronchiolitis in young children [5]. HRVs have the ability to infect the lower respiratory tract (LRT) [6], [7] and can cause chronic infections in immunocompromised hosts [8], [9]. Unlike enteroviruses, most rhinoviruses replicate optimally at lower airway temperatures, which is thought to explain their URT tropism [10]. Rhinovirus strains infecting the LRT may require specific adaptative genomic changes that have not yet been identified.

For most viral infections, disease severity is the result of a close interplay between viral and host factors. The host immune status is known to play a critical role for viral clearance and resolution of infection. Abnormal innate immune responses to HRV have been observed in individuals presenting HRV-induced COPD or asthma exacerbation [11], [12]. By contrast, reports are rare on the implication of viral determinants in HRV disease severity. Palmemberg et al proposed that the rhinovirus 5′UTR polypyrimidine tract may affect virulence [13]. At the serotype or species level, a clear association between the genotype and induced pathology remains to be demonstrated. HRVs were previously classified into two species, HRV-A and HRV-B. In 2006, a third species, HRV-C, was identified [14], [15], [16], [17], [18], [19] and appears to be highly prevalent and circulating worldwide [20], [21], [22]. Several studies have suggested that HRV-C types are prone to induce more severe illness in children [17], [21], [23], [24], [25].

Due to their error-prone RNA polymerase, RNA viruses are subject to constant evolution. The error rate of picornavirus RNA polymerases has been estimated to range between 10−3 and 10−4 errors/nucleotide/cycle of replication [26], [27]. This variability is a driving force for virus evolution. The number of genera, species, and serotypes illustrates this diversity (http://www.picornaviridae.com) with at least 74 HRV-A and 25 HRV-B serotypes identified, and 61 HRV-C types proposed [22].

While the evolution of rhinoviruses has been well characterized at the serotype or species level [13], [28], [29], [30], [31], [32], little is known about the diversity generated during the course of a natural human infection. Several studies have investigated the genomic variability of enteroviruses and vaccine-derived polioviruses during chronic infections in immunocompromised patients, but there is no such report for rhinoviruses [33], [34], [35], [36].

Our group analyzed recently HRV-39 genome evolution in experimental immunocompetent human infections [37]. We estimated an in vivo mutation frequency of 3.4×10−4 mutations/nucleotides over the entire open reading frame (ORF) during a five-day acute infection period, and identified regions of mutation hot spots in the viral capsid (VP1, VP2 and VP3), and 2C and 3C genes.

In the present study, we performed first a phylogenetic analysis to compare the relative distribution of HRV species or serotypes according to the respiratory site (URT versus LRT) and in protracted infection in hospital patients and lung transplant recipients. In the latter group, we evaluated the frequency of mutations and characterized mutation hot spot regions during the course of five natural human protracted infections, each caused by a different rhinovirus serotype. In one case, we followed intra-host rhinovirus genome variation at different time points in the URT and LRT over a period of 27 months using both classical and ultra-deep sequencing methods.

Materials and Methods

Ethics statement

Written informed consent was obtained from all individuals prior to study participation. The studies were approved by the Institutional Review Board of the University of Geneva Hospitals.

RNA extraction, RT and real-time PCR assay

Clinical specimens were extracted by the HCV Amplicor Specimen Preparation kit (Roche, Rotkreuz, Switzerland), TRIzol (Invitrogen, Carlsbad, CA, USA), or Easymag (bioMérieux, Geneva, Switzerland), according to the manufacturers' instructions. Reverse transcription (RT) and real-time PCR with the “Panenterhino/Ge/08” primers and probe assay were performed as previously described [38]. Picornavirus-positive samples were detected in patients enrolled in a cohort of lung transplant recipients (September 2008 to November 2010) [6], [39] and in hospital patients screened by the routine laboratory of the University of Geneva Hospitals (February to November 2009) (Table S1).

Sequencing

Complete genome sequences were obtained from the LRT and URT samples of five chronically-infected individuals collected at different time points during the course of infection. Fragments were amplified by PCR using degenerate primers designed to anneal highly-conserved sequences of the corresponding HRV reference strain. Overlapping fragments were assembled and, when necessary, specific non-degenerate primers were designed to fill the gaps between the original PCR products. All primers used in this study are listed in Table S2. PCR (primers 46 and P1.204) and semi-nested PCR (primers 47 and P1.204, Table S2) were performed for VP4/VP2 sequences.

PCR products were purified with the microcon columns (Millipore, Zug, Switzerland) before sequencing. Each PCR product was sequenced at least twice. Nucleotide changes were confirmed by a second PCR and sequencing. Chromatograms produced with the ABI Prism 3130XL DNA Sequencer (Applied Biosystems, PE Europe BV, Basel, Switzerland) were directly imported for proofreading with the Geneious Pro 5.0.3 software (Biomatters Ltd., Auckland, New Zealand). All sequences are available at GenBank under accession numbers HM347236-727; JF285163-176; JF285179; JF285186; JF285192,193,195,197,198,200,204,206-8,210,213-214,224-225, JF285228-307.

Phylogenetic analysis

Alignments were constructed using Muscle [40] with a maximum of 64 iterations. Multiple FastA was converted into PHYLIP format (for tree-building) with the EMBOSS program ‘seqret’ [41]. Trees were built with PhyML [42] using the GTR model BIONJ for the initial tree and optimized tree topology and branch lengths. Trees with fewer than 50 species used 16 rate categories and larger trees used 8. Transition/transversion ratio, proportion of invariant sites, and the shape parameter (alpha) of the gamma distribution were estimated. Tree processing (including rooting, computation of support values, and displaying) was done with the Newick Utilities [43].

Statistical tests

All statistical data analyses were carried out using the SAS/STAT software version 9.1.3 (SAS Institute Inc., Cary, NC, USA). Associations between HRV species, infection site, and study cohorts were evaluated using chi-square tests when possible, or Fisher's exact tests when some expected values were smaller than 5.

Full-length genome amplification for ultra-deep sequencing

URT and LRT full-length sequences were obtained by amplification of eight overlapping nested PCR products with forward and reverse primers (Table S2). RT-PCR assays were performed in duplicate to discard any mutations introduced by polymerase errors. We used Platinium Taq DNA Polymerase High Fidelity (Invitrogen) for PCRs, according to the manufacturer's instructions. All amplicons were purified with the microcon columns (Millipore) before sequencing. The eight overlapping PCRs were pooled at equimolar concentrations and used for ultra-deep sequencing for each individual sample.

Ultra-deep sequencing analysis

Samples.

Libraries were prepared according to the manufacturer's protocol (Illumina, Inc., San Diego, CA, USA) using indexed adapters designed by Fasteris (Fasteris SA, Plan-les-Ouates, Switzerland) (see Methods S1).

Genome analyzer run.

Libraries were pooled and sequenced on an Illumina Genome Analyzer GAIIx (Illumina) single-read channel for 76 cycles using a version 4 sequencing kit. We performed base-calling using Illumina pipeline RTA SCS.2.6 and CASAVA 1.6, which produced over 24 million pass filter reads attributed unambiguously or 1.85G bases.

Bioinformatics data analyses.

Mapping was performed using the MAQ software version 0.7.1 (http://maq.sourceforge.net/maq-man.shtml). An average of 95% reads were mapped on the “F78 URT 1.4 m” consensus sequence for each sample analyzed with a mean coverage of 21.2K reads per nucleotide. Mapping results were used to extract a list of single nucleotide polymorphisms (SNPs) throughout the genome (see Supplementary method: bioinformatics data analyses).

Mutation analysis along the whole HRV genome.

Density of mutations was represented with a Gaussian kernel density function using a smoothing band width of 0.1 kernel standard deviation. Curves are shown at the same scale and normalized so that a value of 1 is the highest density found over all genomes. Graphs, including kernel estimates, were produced with the R statistical package.

HRV ORF amino acid conservation plot.

The conservation plot was computed with the EMBOSS package Plotcon program using default parameters, except for the sliding window size, which was set to 30 aa. The same program was used in text mode to compute the mean and standard deviations.

Results

Prevalence of each rhinovirus species in protracted and LRT infections

We reported previously three cases of patients chronically infected with a unique rhinovirus strain [8]. Since then, we have identified two additional similar cases. All five patients were immunocompromised lung transplant recipients infected both in the URT and LRT (Table 1). Three were infected with a HRV-A species (HRV-A64, -A24 and -A9), and two with a HRV-B species (HRV-B3 and -B27).

thumbnail
Table 1. Baseline and main characteristics of lung transplant recipients with chronic rhinoviral infection.

https://doi.org/10.1371/journal.pone.0021163.t001

We genotyped strains found in nasopharyngeal specimens (URT), bronchoalveolar lavage fluids (LRT), and strains identified in cases with infections lasting more than three weeks (protracted infection) to assess whether members of the HRV-A, -B and -C species were equally represented in the LRT and/or in the case of protracted infection. Samples were collected from a prospective cohort of lung transplant recipients screened routinely from September 2008 to November 2010, as well as hospital patients screened at the routine laboratory of the University of Geneva Hospitals between February and November 2009 (Table S1). VP4/VP2 sequencing was performed on all positive cases whenever possible. Sequences were obtained for 50 HRV-positive specimens collected from the lung transplant recipient cohort, and 94 specimens collected from hospital patients (56 from children, 36 from adults, and 1 with age not available). Phylogenetic analysis was performed on all sequences obtained, as well as corresponding sequences of HRV-A and -B reference types and proposed HRV-C reference types [22]. Of the 144 sequenced samples, 127 represented different infectious episodes with 71 (55.5%) due to HRV-A types, 10 (7.8%) to HRV-B types, and 47 (36.7%) to HRV-C types (Figure 1 and Figures S1A, B, and C).

thumbnail
Figure 1. Repartition of upper respiratory tract, lower respiratory tract, and protracted infection episodes according to HRV species.

A. Percentage of each species associated with upper respiratory tract infection (URT), lower respiratory tract infection (LRT), and protracted infection (PI) among lung transplant recipients and hospital patients. B. Detailed repartition of URT, LRT and PI episodes in lung transplant recipients, and pediatric and adult hospital patients according to HRV species. Details of HRV-A, -B and -C types are shown in Figures S1A to C.

https://doi.org/10.1371/journal.pone.0021163.g001

No association was observed between a given rhinovirus species and the type of infection (URT, LRT, protracted infection) in the whole study population (Fisher's exact test, 4DF; P = 0.41). HRV-A and -C were present at similar frequencies among LRT infections (22.45% and 19.15%, respectively) in lung transplant recipients and hospital patients (Figure 1). We found also a good representation of different rhinovirus serotypes in all types of specimens (URT, LRT, protracted infection) (Figure S1). By contrast, there was a clear association between the populations studied (lung transplant recipients versus adult and pediatric hospital patients) and the type of infection (URT, LRT, protracted infection) (Fisher's exact test, 2DF; P<0.05). Six protracted infections were observed; five in lung transplant recipients and one in a 3-month-old infant. Finally, we observed an association between the infecting HRV species and patient age (Fisher's exact test, 2DF; P<0.05). The proportion of HRV-C positive samples (LRT, URT, protracted infection) was significantly higher in pediatric hospital patients than in adult hospital patients (51% versus 17%, respectively) (chi-Square 1DF; P<0.05) (Figure 1B).

Mutation frequency and mutation hot spots in chronic cases

Complete genome sequences were obtained at different time points in the URT and/or LRT for the five chronically-infected lung transplant recipients (Table 1). The number of substitutions per nucleotide and per day is reported in Table 1 for each patient, except for A46 for whom samples collected at different time points were not available. Calculated mutation frequencies ranged between 7.27×10−6 and 3.88×10−5. No overlapping mutations in the URT or LRT were found between the five patients. However, superimposition of all changes along a reference genome (HRV-A2) (Figure 2A) revealed that whereas synonymous changes were distributed along the genome without any apparent specific pattern, non-synonymous mutations were clustering mostly in the capsid genes VP2, VP3 and VP1, the first half of 2A genes, and in 3A genes. Cold spot regions were located in the viral capsid gene VP4, as well as in 2B, 2C, 3B, 3C, and most of 3D. Only three changes were found in the 5′UTR and one in the 3′UTR. The variation hot spots observed in the capsid genes were similar to the amino acid variability observed among HRV reference types (Figure 2B).

thumbnail
Figure 2. Synonymous and non-synonymous mutation hot spot regions.

A. The synonymous (upper panel), 5′UTR and non-synonymous (lower panel) changes observed between the first and last samples sequenced from patients P37, V38, R58 and F78, as well as between the URT and LRT samples from patient A46, were plotted together along the HRV-2 genome (GenBank accession number X02316) to map mutation hot and cold spot regions. B. The conservation plot was calculated based on an alignment of the amino acid sequence of 99 rhinovirus serotypes and clinical strains as previously published [13]. The average sequence similarity measure for the ORF was 0.114 (red line). sd: standard deviation.

https://doi.org/10.1371/journal.pone.0021163.g002

Intra-host rhinovirus genome variation according to time and respiratory site

Genetic variation was extensively studied for patient F78 as several samples were available throughout his entire 27-month infection period, including samples collected the same day in the URT and LRT. A phylogenetic tree (Figure 3A) performed with 11 complete genome sequences showed that the URT and LRT viral populations did not segregate for at least 26 months. Furthermore, a LRT subpopulation present in a sample collected after 15.4 months was the founder of all URT and LRT viral populations later collected. The non-synonymous changes that appeared and subsequently fixed were located in VP2 (3), VP3 (1), VP1 (6), 2C (1), 3A (1), and 3D (1). The mutation rate then slowed down shortly before the patient's death (4.3×107 and 1.3×108 viral copy/ml for the last URT and LRT, respectively), despite a high viral load at the end of our follow-up period. This may account for the overall slower mutation rate calculated for this patient when compared with those infected for a shorter period (Table 1). The last three samples collected simultaneously in both the URT and LRT were further analyzed by high throughput sequencing to compare variation at the minority population level (≥5% of the total population). Most minority variants present in the LRT after 26.6 months persisted in this site, but were absent in the corresponding URT viral populations (Figure 3B). A similar observation was made with the majority variants and demonstrates that the genomes of the URT and LRT viral populations varied separately, despite high replication levels at both sites. Of note, several samples from this patient (including sample F78 LRT 26.6 m) were successfully grown in cell culture. Finally, the sequences of the last URT (F78 URT 27.4 m) and LRT (F78 LRT 27.4 m) samples were compared to the first complete genome sequence (F78 URT 1.4 m) and changes specific to the URT (upper panel) and LRT (lower panel) were mapped along the HRV genome (Figure 3C). Specific synonymous changes were found in both the URT and LRT sequences. However, non-synonymous changes were observed only in the LRT where we identified one change in the 5′UTR polypyrimidine tract [44], one non-synonymous change in the VP2 immunogenic site 2, and one non-synonymous change in protein 2A [32], [45].

thumbnail
Figure 3. Intra-host genetic variation in patient F78 at the majority and minority population level.

A. Phylogenetic tree of the 11 complete genome sequences corresponding to samples collected over a 27.4-month infection period in the URT and LRT of patient F78, as well as HRV-A67 reference type (FJ445149) and 3 HRV-9 complete genome sequences (FJ445114,15,17). HRV-A 67 and FJ445117 were used as outgroups to root the tree and are not shown. B. Minority mutation repartition along the HRV genome as recorded for the last 3 URT and LRT samples collected from patient F78. Curves indicate minority mutation densities (estimated by a Gaussian kernel function) along the genome. C. Changes specific to the URT and LRT. The sequences of the last URT and LRT samples collected from patient F78 were compared to the first complete genome sequence (F78 URT 1.4 m) and changes specific to the URT (upper panel) and LRT (lower panel) are mapped along the HRV genome.

https://doi.org/10.1371/journal.pone.0021163.g003

Discussion

Based on our phylogenetic analysis, the respective frequency distribution of strains infecting the URT and LRT did not reveal any apparent correlation between a given HRV serotype or species and their ability to infect the LRT. In both lung transplant recipients and hospital patients, HRV-A and -C strains presented a similar propensity to infect the LRT. No HRV-B strains were found to be associated with LRT, but the relatively small number of infected samples precludes us from drawing firm conclusions. In this respect, we have previously described chronic LRT HRV-B infections [8]. In our study, protracted infections were all found in hosts with a very low level of immunity (mostly lung transplant recipients and a 3-month-old infant). These results suggest that host factors are determinants of disease outcome, rather than virus type.

Five lung transplant recipients were chronically infected with HRV during periods of time ranging from three to 27 months, thus allowing us to study the intra-host genetic variation during natural chronic infections in immunocompromised hosts. Complete Sanger-based and ultra-deep genome sequencing were performed at different time points and sites during infection. Mutation mapping along the HRV genome pointed out that synonymous changes were roughly spread along the entire ORF, whereas non-synonymous changes clustered mostly in the capsid VP2, VP3, and VP1 genes. These capsid genes are also the most variable during acute infections in immunocompetent hosts [37], as well as the most diverse among the 99 HRV reference types [13]. The VP2, VP3 and VP1 proteins form the outer shell of the virion and are under immune pressure, which explains their fast evolution. The localization of mutation hot spots in these genes in immunocompromised hosts may appear unexpected. However, although these patients are highly immunosuppressed at the time of transplantation, immunosuppressive therapy is then adapted to prevent graft rejection and to preserve a residual cellular and humoral immunity that may account for this observation. This residual immunity is also necessary for the control and recovery of viral and bacterial infections observed in these patients.

We calculated the mutation frequency occurring during these five chronic infections, each caused by a distinct HRV serotype and for various periods of time. Mutation frequency was different in each host and likely reflects different viral replication levels, different host environments, and different durations of infection. The calculated mutation frequency ranged between 7.27×10−06 and 3.88×10−05 change/nt/day. Recent data on HRV-39 genome evolution over a 5-day acute experimental infection in human immunocompetent volunteers revealed a mutation frequency of 3.4×10−4 change/nt (equivalent to 6.83×10−5 change/nt/day) [37]. Therefore, the HRV mutation frequency observed in some immunocompromised hosts infected over months was of the same order of magnitude as that observed in acute 5-day infections.

Extensive analysis of rhinovirus genome modification was performed in the case of one patient (F78) infected for more than two years with a HRV-A9 strain. Phylogenetic analysis of 11 complete genome sequences (including ultra-deep sequences performed at three different time points) revealed that the URT and LRT populations were phylogenetically indistinguishable for a prolonged period of time, either because of co-variation or, more likely, constant viral population mixing. Once potentially adapted to the host through the accumulation of non-synonymous changes mainly in the capsid gene VP1, the viral populations appeared to evolve more slowly. Finally, we observed signatures of putative adaptation to lower airway conditions after several months of infection. Indeed, the last two LRT populations studied presented specific changes in the 5′UTR polypyrimidine tract and two non-synonymous changes, one in the VP2 immunogenic site 2 and the second in protein 2A. These changes were absent in the URT population and, upon their appearance, the two populations remained separated as confirmed by ultra-deep sequencing analysis. Although the direct impact of these changes on virus growth ability at lower airway conditions remains to be demonstrated experimentally, the very high viral load observed at this stage suggests that both populations were adapted to their sites.

Taken together, our data suggest that immunocompromised patients cannot clear viral infections as immunocompetent individuals, and represent a potential reservoir for the emergence of new variants and inter-host transmission due to longer chronic viral infection. In addition, these patients may be co-infected by two viruses, thus opening the door to recombination, another putative driving force of rhinovirus evolution [30]. With the emergence of new therapies and progress in transplantation, the population of immunocompromised patients is constantly increasing. Our results suggest that this could accelerate the ability of viruses to adapt to the host, evolve, and propagate and may favor a more rapid emergence of new viral variants.

Supporting Information

Figure S1.

Repartition of protracted and lower respiratory tract infections among the HRV-A (panel A), HRV-B (panel B), and HRV-C (panel C) reference serotypes. VP4-VP2 cladograms of rhinoviruses isolated from lung transplant recipients (blue) and routinely screened hospital patients (red) (Table S1), as well as the 74 HRV-A (panel A), 25 HRV-B (panel B), and 61 proposed HRV-C reference types (black) (panel C). SV1 (GenBank accession number AY064708) was used as an outgroup. LRT and PI are highlighted by arrows and black lines respectively.

https://doi.org/10.1371/journal.pone.0021163.s001

(PDF)

Table S1.

Samples collected from the lung transplant recipient cohort and samples from hospital patients screened by the routine laboratory of the University Hospitals of Geneva. *Samples are designated by the first letter of their family name followed by the last two numbers of their birth year (eg. B48), the month and year of sample collection (eg. 0209), and the site of sampling (U, upper respiratory tract; L, lower respiratory tract).

https://doi.org/10.1371/journal.pone.0021163.s002

(XLS)

Table S2.

Primers used to amplify and sequence the genomes of the clinical strains described in Table 1 .

https://doi.org/10.1371/journal.pone.0021163.s003

(XLS)

Acknowledgments

We would like to thank Lara Turin, Chantal Gaille, and Ghislaine Wagner for technical assistance, Anouk Deruaz for help with statistical tests, and Rosemary Sudan for editorial assistance.

Author Contributions

Conceived and designed the experiments: CT LK. Performed the experiments: SC SVB. Analyzed the data: TJ EZ LF. Contributed reagents/materials/analysis tools: JA PS TJ EZ LF. Wrote the paper: CT. Revised the manuscript: SC LK.

References

  1. 1. Gern JE, Busse WW (1999) Association of rhinovirus infections with asthma. Clin Microbiol Rev 12: 9–18.
  2. 2. Johnston SL, Pattemore PK, Sanderson G, Smith S, Campbell MJ, et al. (1996) The relationship between upper respiratory infections and hospital admissions for asthma: a time-trend analysis. Am J Respir Crit Care Med 154: 654–660.
  3. 3. Johnston SL, Pattemore PK, Sanderson G, Smith S, Lampe F, et al. (1995) Community study of role of viral infections in exacerbations of asthma in 9–11 year old children. BMJ 310: 1225–1229.
  4. 4. Seemungal TA, Harper-Owen R, Bhowmik A, Jeffries DJ, Wedzicha JA (2000) Detection of rhinovirus in induced sputum at exacerbation of chronic obstructive pulmonary disease. Eur Respir J 16: 677–683.
  5. 5. Renwick N, Schweiger B, Kapoor V, Liu Z, Villari J, et al. (2007) A recently identified rhinovirus genotype is associated with severe respiratory-tract infection in children in Germany. J Infect Dis 196: 1754–1760.
  6. 6. Garbino J, Soccal PM, Aubert JD, Rochat T, Meylan P, et al. (2009) Respiratory viruses in bronchoalveolar lavage: a hospital-based cohort study in adults. Thorax 64: 399–404.
  7. 7. Papadopoulos NG, Bates PJ, Bardin PG, Papi A, Leir SH, et al. (2000) Rhinoviruses infect the lower airways. J Infect Dis 181: 1875–1884.
  8. 8. Kaiser L, Aubert JD, Pache JC, Deffernez C, Rochat T, et al. (2006) Chronic rhinoviral infection in lung transplant recipients. Am J Respir Crit Care Med 174: 1392–1399.
  9. 9. Milano F, Campbell AP, Guthrie KA, Kuypers J, Englund JA, et al. (2010) Human rhinovirus and coronavirus detection among allogeneic hematopoietic stem cell transplantation recipients. Blood 115: 2088–2094.
  10. 10. Papadopoulos NG, Sanderson G, Hunter J, Johnston SL (1999) Rhinoviruses replicate effectively at lower airway temperatures. J Med Virol 58: 100–104.
  11. 11. Mallia P, Message SD, Gielen V, Contoli M, Gray K, et al. (2011) Experimental rhinovirus infection as a human model of chronic obstructive pulmonary disease exacerbation. Am J Respir Crit Care Med 183: 734–742.
  12. 12. Message SD, Laza-Stanca V, Mallia P, Parker HL, Zhu J, et al. (2008) Rhinovirus-induced lower respiratory illness is increased in asthma and related to virus load and Th1/2 cytokine and IL-10 production. Proc Natl Acad Sci USA 105: 13562–13567.
  13. 13. Palmenberg AC, Spiro D, Kuzmickas R, Wang S, Djikeng A, et al. (2009) Sequencing and analyses of all known human rhinovirus genomes reveal structure and evolution. Science 324: 55–59.
  14. 14. Arden KE, McErlean P, Nissen MD, Sloots TP, Mackay IM (2006) Frequent detection of human rhinoviruses, paramyxoviruses, coronaviruses, and bocavirus during acute respiratory tract infections. J Med Virol 78: 1232–1240.
  15. 15. Kistler A, Avila PC, Rouskin S, Wang D, Ward T, et al. (2007) Pan-viral screening of respiratory tract infections in adults with and without asthma reveals unexpected human coronavirus and human rhinovirus diversity. J Infect Dis 196: 817–825.
  16. 16. Lamson D, Renwick N, Kapoor V, Liu Z, Palacios G, et al. (2006) MassTag polymerase-chain-reaction detection of respiratory pathogens, including a new rhinovirus genotype, that caused influenza-like illness in New York State during 2004–2005. J Infect Dis 194: 1398–1402.
  17. 17. Lau SK, Yip CC, Tsoi HW, Lee RA, So LY, et al. (2007) Clinical features and complete genome characterization of a distinct human rhinovirus (HRV) genetic cluster, probably representing a previously undetected HRV species, HRV-C, associated with acute respiratory illness in children. J Clin Microbiol 45: 3655–3664.
  18. 18. Lee WM, Kiesner C, Pappas T, Lee I, Grindle K, et al. (2007) A diverse group of previously unrecognized human rhinoviruses are common causes of respiratory illnesses in infants. PLoS ONE 2: e966.
  19. 19. McErlean P, Shackelton LA, Lambert SB, Nissen MD, Sloots TP, et al. (2007) Characterisation of a newly identified human rhinovirus, HRV-QPM, discovered in infants with bronchiolitis. J Clin Virol 39: 67–75.
  20. 20. Huang T, Wang W, Bessaud M, Ren P, Sheng J, et al. (2009) Evidence of recombination and genetic diversity in human rhinoviruses in children with acute respiratory infection. PLoS One 4: e6355.
  21. 21. Wisdom A, Leitch EC, Gaunt E, Harvala H, Simmonds P (2009) Screening respiratory samples for detection of human rhinoviruses (HRVs) and enteroviruses: comprehensive VP4-VP2 typing reveals high incidence and genetic diversity of HRV species C. J Clin Microbiol 47: 3958–3967.
  22. 22. Simmonds P, McIntyre C, Savolainen-Kopra C, Tapparel C, Mackay IM, et al. (2010) Proposals for the classification of human rhinovirus species C into genotypically assigned types. J Gen Virol 91: 2409–2419.
  23. 23. Khetsuriani N, Lu X, Teague WG, Kazerouni N, Anderson LJ, et al. (2008) Novel human rhinoviruses and exacerbation of asthma in children. Emerg Infect Dis 14: 1793–1796.
  24. 24. Miller EK, Edwards KM, Weinberg GA, Iwane MK, Griffin MR, et al. (2009) A novel group of rhinoviruses is associated with asthma hospitalizations. J Allergy Clin Immunol 123: 98–104 e101.
  25. 25. Tapparel C, L'Huillier AG, Rougemont AL, Beghetti M, Barazzone-Argiroffo C, et al. (2009) Pneumonia and pericarditis in a child with HRV-C infection: a case report. J Clin Virol 45: 157–160.
  26. 26. Drake JW (1999) The distribution of rates of spontaneous mutation over viruses, prokaryotes, and eukaryotes. Ann N Y Acad Sci 870: 100–107.
  27. 27. Harvala H, Simmonds P (2009) Human parechoviruses: biology, epidemiology and clinical significance. J Clin Virol 45: 1–9.
  28. 28. McIntyre CL, McWilliam Leitch EC, Savolainen-Kopra C, Hovi T, Simmonds P (2010) Analysis of genetic diversity and sites of recombination in human rhinovirus species C. J Virol 84: 10297–10310.
  29. 29. Tapparel C, Junier T, Gerlach D, Cordey S, Van Belle S, et al. (2007) New complete genome sequences of human rhinoviruses shed light on their phylogeny and genomic features. BMC Genomics 8: 224.
  30. 30. Tapparel C, Junier T, Gerlach D, Van-Belle S, Turin L, et al. (2009) New respiratory enterovirus and recombinant rhinoviruses among circulating picornaviruses. Emerg Infect Dis 15: 719–726.
  31. 31. Kistler AL, Webster DR, Rouskin S, Magrini V, Credle JJ, et al. (2007) Genome-wide diversity and selective pressure in the human rhinovirus. Virol J 4: 40.
  32. 32. Lewis-Rogers N, Bendall ML, Crandall KA (2009) Phylogenetic relationships and molecular adaptation dynamics of human rhinoviruses. Mol Biol Evol 26: 969–981.
  33. 33. Odoom JK, Yunus Z, Dunn G, Minor PD, Martin J (2008) Changes in population dynamics during long-term evolution of sabin type 1 poliovirus in an immunodeficient patient. J Virol 82: 9179–9190.
  34. 34. Hovi T, Lindholm N, Savolainen C, Stenvik M, Burns C (2004) Evolution of wild-type 1 poliovirus in two healthy siblings excreting the virus over a period of 6 months. J Gen Virol 85: 369–377.
  35. 35. Bailly JL, Chambon M, Henquell C, Icart J, Peigue-Lafeuille H (2000) Genomic variations in echovirus 30 persistent isolates recovered from a chronically infected immunodeficient child and comparison with the reference strain. J Clin Microbiol 38: 552–557.
  36. 36. Kew OM, Sutter RW, Nottay BK, McDonough MJ, Prevots DR, et al. (1998) Prolonged replication of a type 1 vaccine-derived poliovirus in an immunodeficient patient. J Clin Microbiol 36: 2893–2899.
  37. 37. Cordey S, Junier T, Gerlach D, Gobbini F, Farinelli L, et al. (2010) Rhinovirus genome evolution during experimental human infection. PLoS One 5: e10588.
  38. 38. Tapparel C, Cordey S, Van Belle S, Turin L, Lee WM, et al. (2009) New molecular detection tools adapted to emerging rhinoviruses and enteroviruses. J Clin Microbiol 47: 1742–1749.
  39. 39. Soccal PM, Aubert JD, Bridevaux PO, Garbino J, Thomas Y, et al. (2010) Upper and lower respiratory tract viral infections and acute graft rejection in lung transplant recipients. Clin Infect Dis 51: 163–170.
  40. 40. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792–1797.
  41. 41. Rice P, Longden I, Bleasby A (2000) EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet 16: 276–277.
  42. 42. Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52: 696–704.
  43. 43. Junier T, Zdobnov EM (2010) The Newick utilities: high-throughput phylogenetic tree processing in the UNIX shell. Bioinformatics 26: 1669–1670.
  44. 44. Anderson EC, Hunt SL, Jackson RJ (2007) Internal initiation of translation from the human rhinovirus-2 internal ribosome entry site requires the binding of Unr to two distinct sites on the 5′ untranslated region. J Gen Virol 88: 3043–3052.
  45. 45. Sherry B, Mosser AG, Colonno RJ, Rueckert RR (1986) Use of monoclonal antibodies to identify four neutralization immunogens on a common cold picornavirus, human rhinovirus 14. J Virol 57: 246–257.