key: cord-291930-n7wq09rq
authors: Arden, K.E.; Faux, C.E.; O’Neill, N.T.; McErlean, P.; Nitsche, A.; Lambert, S.B.; Nissen, M.D.; Sloots, T.P.; Mackay, I.M.
title: Molecular characterization and distinguishing features of a novel human rhinovirus (HRV) C, HRVC-QCE, detected in children with fever, cough and wheeze during 2003
date: 2010-01-27
journal: J Clin Virol
DOI: 10.1016/j.jcv.2010.01.001
sha: 
doc_id: 291930
cord_uid: n7wq09rq

BACKGROUND: Human rhinoviruses (HRVs) are associated with more acute respiratory tract infections than any other viral group yet we know little about viral diversity, epidemiology or clinical outcome resulting from infection by strains, in particular the recently identified HRVs. OBJECTIVES: To determine whether HRVC-QCE was a distinct HRV-C strain, by determining its genome and prevalence, by cataloguing genomic features for strain discrimination and by observing clinical features in positive patients. STUDY DESIGN: Novel real-time RT-PCRs and retrospective chart reviews were used to investigate a well-defined population of 1247 specimen extracts to observe the prevalence and the clinical features of each HRV-QCE positive case from an in- and out-patient pediatric, hospital-based population during 2003. An objective illness severity score was determined for each HRVC-QCE positive patient. RESULTS: Differences in overall polyprotein and VP1 binding pocket residues and the predicted presence of a cis-acting replication element in 1B defined HRVC-QCE as a novel HRV-C strain. Twelve additional HRVC-QCE detections (1.0% prevalence) occurred among infants and toddlers (1–24 months) suffering mild to moderate illness, including fever and cough, who were often hospitalized. HRVC-QCE was frequently detected in the absence of another virus and was the only virus detected in three (23% of HRVC-QCE positives) children with asthma exacerbation and in two (15%) toddlers with febrile convulsion. CONCLUSIONS: HRVC-QCE is a newly identified, genetically distinct HRV strain detected in hospitalized children with a range of clinical features. HRV strains should be independently considered to ensure we do not overestimate the HRVs in asymptomatic illness.

Acute respiratory tract infections (ARTIs), a major cause of human morbidity and hospitalization, are most frequently viral in origin. The human rhinoviruses (HRVs) have for over a decade been associated, more than any other individual factor, with asthma 1-3 and chronic obstructive pulmonary disease exacerbations 4 and the majority of cold and flu-like illnesses (CFLIs). 5, 6 It has been shown that HRVs infect both upper and lower respiratory tract tissues. 7, 8 Classical serotypes (HRV-1 to 100) or "strains", the term we use to indicate a molecular equivalent, form two accepted species within the heavily populated genus Enterovirus, family Picornaviridae. Despite the considerable economic burden of HRV infections, 9 deliberate and comprehensive investigation of strain-specific clinical data is lacking.

In 2006 we proposed that a clade of divergent sequences belonged to a novel viral subgroup, HRV-A2, 10 which we later redefined 11 as a putative third species, HRV-C, in agreement with Lau et al. 12 To date, nine complete polyprotein sequences of HRV-C strains have been determined from Australia (HRVC-QPM 10 HRVC-NY074 14 ) , Hong Kong (HRVC-C024, HRVC-C025 and HRVC-C026 12 ) and Shanghai (HRVC-N4 and HRVC-N10 15 ).

Many HRV strains have distinct genetic, antigenic, immunogenic and epidemiologic profiles which are similar in nature to those of other respiratory virus species including human respiratory syncytial virus (HRSV). The prospect of variation in sensitivity to antivirals is also an important reason to investigate each novel strain comprehensively, as is the chance that infection by specific HRV species could be associated with more severe illness or particular pre-existing conditions. Disappointingly, the possibility that discrete HRV strains exist as distinct viral entities is not accounted for in asymptomatic control populations because strain identification is not routinely conducted. As a consequence elevated prominence is, by default, assigned to HRVs as if to a single virus.

We sought to genomically characterize HRV-QCE, an apparently distinct HRV-C strain whose novel subgenomic sequences were first detected in a child with exacerbated asthma and to define distinct molecular features among the HRV-Cs that are useful in discriminating strains for future studies and observe the clinical features amongst HRVC-QCE positive patients.

Clinical specimens (n = 1247; three additions were made to the previously described HRVC-QPM population 10 ) were predominantly nasopharyngeal aspirates (NPA; 93%) from patients (59% male) aged 1 day to 80 years (mean 9.2 years, median 1.3 years) presenting to Queensland hospitals with symptoms of ARTI during 2003. Extracts were stored and tested as previously described. 10 A clinical illness severity scoring system was applied as previously described. 11 This study was approved by the Royal Children's Hospital and University of Queensland Medical Research Ethics Committees (specific approvals #2006/097 and #2008/020 and #2007000004).

Primers for genome deduction were first designed using early HRVC-QCE sequences in 1A/1B (EU131891 and EU155152). All oligonucleotide sequences are available upon request. Additional primers supplemented the process. 10, 12, 16 To confirm our primary sequence belonged to a distinct virus and was not the result of mixed HRV template, a second round of RT-PCR and sequencing of overlapping fragments was performed on RNA from the original patient specimen using newly designed HRVC-QCE-specific primers. Discrepant sequence results were confirmed by additional rounds of RT-PCR and sequencing. Two specific diagnostic real-time RT-PCRs (RT-rtPCRs) were performed on the Rotor-Gene TM 3000 (Corbett Research) for diagnostic screening.

Polyprotein protease cleavage sites were sought using the Net-PicoRNA server. 17 Cis-acting replication element (cre) sequences 18 were sought using the RNAalifold server and structures were predicted using MFold. 19 

The polyprotein encoding sequence (GenBank #GQ323774) of HRVC-QCE, a new member of the species Human rhinovirus C, spanned 6423 nt; shorter than all HEVs and HRVs except HRV-NY074. 14, 20, 21 It had a G + C content of 43% and 2140 predicted amino acids terminating in an HRV-C-specific isoleucine. 11 The HRVC-QCE polyprotein shared 54% average amino acid identity with HRV-A strains, 50% with HRV-Bs and 79% with other HRV-Cs. Ten protease cleavage sites were predicted to reduce the HRVC-QCE polyprotein into typical picornavirus structural (VP1-4) and non-structural proteins (2A-C and 3A-D 17, 22, 23 ). HRVC-QCE shared a translation initiation site (MGAQVS) with most other HRVs and three motifs crucial for RNA polymerase binding (YGDD, TFLKR and SIRWT) shared by HRVs, and HEVs 24,25 but lacked two of three exposed VP1 motifs previously conserved among HRV polyproteins 26 were absent in HRV-QCE. Phylogeny inferred from a discontinuous alignment of 10 amino acids in VP3 and VP1 which comprise the predicted intercellular adhesion molecule (ICAM)-1 footprint 27 in major group HRV strains, placed HRVC-QCE amongst other HRV-Cs on a branch with HRV-1 (a minor group HRV which employs a different receptor) but distinct from the known major (employing the ICAM-1 molecule as receptor) or remaining minor group HRV strains (Fig. 1) . Key residues on the VP1 BC loop involved in receptor contact with the VLDLR 19 are missing in HRVC-QCE and there are differences in the HI loop compared to the classical strains. Antigenic site A 28 is absent in HRVC-QCE and the sequence at site B is unique to HRVC-QCE. As with other HRV-C strains identified to date, 18 the conserved loop sequence motif constituting a cre, R 1 NNNA 1 A 2 R 2 NNNNNNR 3 was identified (GCUCAAGCAAAUCA) in 1B of HRVC-QCE in the context of an appropriate predicted RNA stem-loop structure (not shown).

All classical HRV-A strains and 18/25 HRV-B strains are susceptible to the antiviral compound, pleconaril. Among the 25 key contact residues lining the antiviral binding pocket, there were 11 or 12 HRVC-QCE-specific differences compared to HRV-A or HRV-B consensus sequences, 29 respectively (Fig. 2) ; 11 differences exist between HRV-As and Bs and two differences between the consensus sequence of HRV-B sensitive versus resistant strains. Changes at two HRVC-QCE residues historically important for identifying naturally resistant strains were identified as Tyr 152 (in resistant strains) to Phe 152 and Val 191 to Thr 191 . His 245 , conserved in nearly all other classical HRV strains, varied in HRVC-QCE and among other HRV-Cs. Other conserved HRV residues including Ile/Leu 106 and Tyr/Phe/Ala 150, differed in HRVC-QCE (Fig. 2) .

To examine whether HRVC-QCE was a variant of a previously identified HRV strain, we first compared the full polyprotein against public sequence databases. The nearest matches were HRVC-C026 (93.5% identity, 2003 identical amino acids) and HRVC-QPM (92.9%, 1991 amino acids). The inferred amino acid sequences from complete 1D regions of eight HRVC-QCE-positive extracts (variants 001-004, 006, 009, 010, 012; intra-strain amino acid identity, 99.5%; GenBank accession number GQ323774 and bankit numbers, 1304337, 1304336, 1304343, 1304341, 1304344, 1304335, 1304320) were compared to those from known HRV and HEV serotypes and to previously described HRVC-QPM variants (n = 17; intra-strain identity of 99.8%; Fig. 3 ). The amino acid identity between variants of the HRVC-QPM and HRVC-QCE strains (91.2%) was lower than among variants of the two strains (≥99.5%). At least 93% polyprotein-length pairwise identity is also shared by 80% (n = 60) of antigenically distinct and officially recognised HRV-A strains (60 strains occur in at least one pairing in which sequence identity is ≥93%; HRV-8 and 95, HRV-1A and 1B and HRV-82 and sc were not considered discrete viruses). The HRVC-QCE VP1 (272 aa) is 9-24 amino acids shorter than that of other HRV serotypes. The HRVC-QCE residues aligning with the 10 residue ICAM-1 footprint differed by one from those of the most similar strain (HRVC-C026). In VP1, two of 6 amino acids from antigenic site B, four of 20 from the EF loop, one of 3 from the FG loop, two of 32 from the GH loop, three of 8 from the HI loop, one of 24 constituting the binding pocket and two of 48 nt in the 1B cre differed from the most similar HRV-C strain.

Since DQ875928 = NY60). The complete polyprotein of most of these strains has not been described and for the one that has (HRVC-NAT001), no additional variants were sought and only preliminary clinical studies were performed. HRVC-QCE RNA was detected in 1.0% (n = 13) of 1247 specimen extracts collected from January to December 2003 using a 1B targeted RT-rtPCR and confirmed by an assay targeting the 2C region. No control patients were sampled in 2003. HRVC-QCE prevalence was bimodal with the major peak in summer when 3.5% of specimens were positive (7.6% of 119 specimens from February were positive representing 76.9% of all HRVC-QCE positives and 16.1% of all virus detections (n = 56) that month). HRVC-QCE was the only HRV detected in January. The remaining 23.1% of HRVC-QCE positives were detected in winter (0.8% of specimens, peaking at 3.8% of all viral detections in June). No other viruses were detected in 61.5% (n = 8) of HRVC-QCE positives (Table 1 ) leaving five extracts in which another viral sequence was identified. Respiratory picornaviruses were the most frequently detected virus at 40.9% of detections (n = 351) followed by HBoV (11.8%, n = 101), HRSV (8.3%, n = 71) and HMPV (7.7%, n = 66). No viruses were found in 46.0% (n = 574) of extracts. Isolation of HRVC-QCE in culture was not attempted because of the age of the clinical specimens.

Chart review of the 13 HRVC-QCE positive cases revealed presentation with lower respiratory tract illnesses (6/13; 46.2%), most commonly involving wheeze and administration of bronchodilators (Table 1) . Antibiotics were not usually given to wheezers compared to non-wheezers. Eleven of the 13 cases (84.6%) were admitted to hospital and HRVC-QCE was the sole detection in 61.5% (8/13) of patients. The average severity score among single HRVC-QCE detections was higher (1.6) than that for co-detections (1.2) and was higher in summer (1.6) than in winter (1.0). Average scores were highest among lower respiratory tract presentations (1.8) and scores of 0 were obtained from an upper respiratory tract illness and a cystic fibrosis patient with an HRSV co-detection. There were no adult cases, despite adults comprising 15.8% of the screened population. Most (84.6%) patients were 2 years of age or younger.

We have described the detection, polyprotein coding region and distinguishing genetic features of a newly identified and evidently genetically distinct HRV-C strain. We have summarised the clinical features of a predominantly pediatric population from whom HRVC-QCE variants were detected. The relatively shortened genome with truncations and sequence diversity in the 1D capsid region may affect receptor usage, antigenicity and sensitivity to antivirals. Phe 152 in VP1 is a distinctive feature of pleconaril resistant clinical HRV-B strains, 30 and was present in the VP1 of HRVC-QCE as was Thr 191 , which is found in both resistant and sensitive strains. The predicted HRV-C-specific cre motif in 1B, the terminal isoleucine, the comparatively elevated G + C content and other features we have previously defined, 11 identified HRVC-QCE as a candidate for the proposed new species HRV-C.

HRVC-QCE appears to be a distinct global pathogen however no consensus molecular criteria exist which define strains or provide identity thresholds below which a clinical detection may be defined as a novel strain nor above which it may best be called a variant of a previously described strain. 11 The HRV literature notes infections which apparently have not elicited a cross-protective immune response as demonstrated by patients infected by distinct strains of the same or different species such over a defined period which share 62-93% amino acid identity in VP4. 12, 31, 32 These biological data together with strain-determining genetic criteria permit the proposal of a practical molecular threshold, for example ≤93% amino acid identity, that could be used as a surrogate to define the biological distinctiveness of HRVs that cannot be cultured. Strain-defining criteria are especially useful to determine whether it is more accurate to consider "the HRVs" as a single group or as a collection of related but mostly antigenically discrete viruses, as is the norm for HRSV or HMPV for example. Applying strain-defining criteria to control populations may very likely redefine our concept of the HRVs in asymptomatic illness, since the proportion of specimens positive for any distinct HRV strain in controls will be less than that which is positive for the HRV super-group when considered as a whole.

Approximately half of HRVC-QCE detections were from patients with lower respiratory tract symptoms and 85% of HRVC-QCE positives were detected from hospitalized children. HRVC-QCE was the sole virus detected in over half of all cases.

The paucity of sequence data, strain-defining criteria and clinical investigations of the novel HRVs or of individual classical strains parallel our limited knowledge of the extent and clinical impact of HRV diversity. The most recent significant healthcare impact attributed to HRVs was the 2009 H1N1v influenza pandemic when HRVs were the most common viruses detected among patients qualifying, but laboratory-negative for, pandemic virus testing. 33 Similar findings were described during the SARS-CoV outbreak. 34 A lack of viral sequence data weakens our differential diagnostic and overall detection capabilities which limits our ability to define epidemiological features as well as the occurrence of strain-, or species-specific clinical outcomes and potentially differentiating genetic features. Due to the retrospective nature of our study, no temporally comparable control population was available for testing. This is a limitation among most other HRV studies as is the current inability to culture the newly identified HRVs. When a control population is included however, strain typing is usually not conducted and so a higher proportion of asymptomatic illness is attributed to all HRVs as a total group. Additionally, because we and others cannot fulfil Koch's postulates for the currently unculturable HRV-Cs, the clinical observations cannot be specifically associated with the PCR-determined presence of HRVC-QCE RNA. Our findings highlight the need to better understand the genomics, taxonomy, epidemiology and clinical impact of the newly identified HRV strains to accurately quantify the scope of their clinical impact as distinct respiratory viruses.

None.

Personal exposure to nitrogen dioxide (NO2) and the severity of virus-induced asthma in children

Community study of role of viral infections in exacerbations of asthma in 9-11 year old children

Rhinovirus and respiratory syncytial virus in wheezing children requiring emergency care

Overview of virus-induced airway disease

Picornavirus infections: a primer for the practitioner

Viral etiology of common cold in children

Localization of human rhinovirus replication in the upper respiratory tract by in situ hybridization

Basal cells of differentiated bronchial epithelium are more susceptible to rhinovirus infection

The cost of community-managed viral respiratory illnesses in a cohort of healthy preschool-aged children

Characterisation of a newly identified human rhinovirus, HRV-QPM, discovered in infants with bronchiolitis

Distinguishing molecular features and clinical characteristics of a putative new rhinovirus species, human rhinovirus C (HRV C)

Clinical features and complete genome characterization of a distinct human rhinovirus genetic cluster, probably representing a previously undetected HRV species, HRV-C, associated with acute respiratory illness in children

Pan-viral screening of respiratory tract infections in adults with and without asthma reveals unexpected human coronavirus and human rhinovirus diversity

Sequencing and analyses of all known human rhinovirus genomes reveals structure and evolution

Evidence of recombination and genetic diversity in human rhinoviruses in children with acute respiratory infection

Frequency and dynamics of recombination within different species of human enteroviruses

Cleavage site analysis in picornaviral polyproteins: discovering cellular targets by neural networks

The cis-acting replication elements define human enterovirus and rhinovirus species

X-ray structure of a minor group human rhinovirus bound to a fragment of its cellular receptor protein

Recombination in circulating human enterovirus B: independent evolution of structural and non-structural genome regions

Enterovirus 68 is associated with respiratory illness and shares biological features with both the enteroviruses and the rhinoviruses

Human rhinovirus 2: complete nucleotide sequence and proteolytic processing signals in the capsid protein region

The complete nucleotide sequence of a common cold virus: human rhinovirus 14

Sequence analysis of human rhinoviruses in the RNA-dependent RNA polymerase coding region reveals within-species variation

Complete genomic sequencing shows that polioviruses and members of human enterovirus species C are closely related in the noncapsid coding region

Alignment of capsid protein VP1 sequences of all human rhinovirus prototype strains: conserved motifs and functional domains

Molecular characterization of human rhinovirus field strains isolated during surveillance of enteroviruses

Structure of human rhinovirus serotype 2 (HRV2)

VP1 sequencing of all human rhinovirus serotypes: insights into genus phylogeny and susceptibility to antiviral capsid-binding compounds

Insights into the genetic basis for natural phenotypic resistance of human rhinoviruses to pleconaril

Recurrent human rhinovirus infections in infants with refractory wheezing

Clinical and molecular epidemiology of human rhinovirus C in children and adults in Hong Kong reveals a possible distinct human rhinovirus C subgroup

A variety of respiratory viruses found in symptomatic travellers returning from countries with ongoing spread of the new influenza A(H1N1)v virus strain

SARS surveillance during emergency public health response

MEGA3: integrated software for molecular evolutionary genetic analysis and sequence alignment

We are grateful to Patrick Woo and Susanna Lau for primer and HRV strain sequences and to Daniel Gerlach for bioinformatics support. We thank Pathology Queensland Central for the provision of clinical specimens.Funding. This work was supported by an Australian NHMRC Project grant 455905.