Emerg Infect DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention18507910260030808-027110.3201/eid1406.080271DispatchGlobal Distribution of Novel Rhinovirus GenotypeGlobal Distribution of Novel Rhinovirus GenotypeNovel Rhinovirus GenotypeBrieseThomas*RenwickNeil*VenterMarietjieJarmanRichard G.GhoshDhrubaa§KöndgenSophieShresthaSanjaya K.#HoeghA. Mette**CasasInmaculada††AdjogouaEdgard Valerie‡‡Akoua-KoffiChantal‡‡MyintKhin SawWilliamsDavid T.§§ChidlowGlenys¶¶van den BergRiaCalvoCristina##KochOrienkaPalaciosGustavo*KapoorVishal*VillariJoseph*DominguezSamuel R.***HolmesKathryn V.***HarnettGerry¶¶SmithDavid¶¶MackenzieJohn S.§§EllerbrokHeinzSchweigerBrunhildeSchønningKristian**ChadhaMandeep S.§LeendertzFabian H.MishraA.C.§GibbonsRobert V.HolmesEdward C.†††‡‡‡LipkinW. Ian*Columbia University, New York, New York, USAUniversity of Pretoria/NHLS Tswhane Academic Division, Pretoria, South AfricaArmed Forces Research Institute of Medical Sciences, Bangkok, ThailandNational Institute of Virology, Pune, IndiaRobert Koch-Institut, Berlin, GermanyWalter Reed AFRIMS Research Unit Nepal, Katmandu, NepalHvidovre University Hospital, Hvidovre, DenmarkCentro Nacional de Microbiologia, Instituto de Salud Carlos III, Majadahonda, Madrid, SpainInstitut Pasteur Côte d’Ivoire, Abidjan, Côte d’IvoireCurtin University of Technology, Perth, Western Australia, AustraliaPathWest Laboratory, Nedlands, Western Australia, AustraliaSevero-Ochoa Hospital, Leganés, Madrid, SpainUniversity of Colorado Denver School of Medicine, Aurora, Colorado, USAPennsylvania State University, University Park, Pennsylvania, USANational Institutes of Health, Bethesda, Maryland, USAAddress for correspondence: Thomas Briese, Center for Infection and Immunity, Mailman School of Public Health, Columbia University, 722 West 168th St, 18th Floor, New York, NY 10032, USA; email: thomas.briese@columbia.edu62008146944947

Global surveillance for a novel rhinovirus genotype indicated its association with community outbreaks and pediatric respiratory disease in Africa, Asia, Australia, Europe, and North America. Molecular dating indicates that these viruses have been circulating for at least 250 years.

Keywords: picornavirusrhinovirusHRV-Cmultiplex MassTag PCRlower respiratory tract infectionchildhood pneumoniadispatch

Acute respiratory illness (ARI) is the most frequent infectious disease of humans. Ordinary upper respiratory tract infections are usually self-limited; nevertheless, they result in major economic impact through loss of productivity and strain on healthcare systems. Lower respiratory tract infections (LRTIs) are among the leading causes of death in children <5 years of age worldwide, particularly in resource-poor regions (1). Streptococcus pneumoniae and Haemophilus influenzae are important bacterial causes of ARI, although their impact is expected to decline with increasing vaccine coverage. Collectively, however, viruses dominate as causative agents in ARI. Viruses frequently implicated in ARI include influenza virus, respiratory syncytial virus, metapneumovirus, parainfluenza virus, human enterovirus (HEV), and human rhinovirus (HRV).

HRVs are grouped taxonomically into Human rhinovirus A (HRV-A) and Human rhinovirus B (HRV-B), 2 species within the family Picornaviridae (International Committee on Taxonomy of Viruses database [ICTVdb]; http://phene.cpmc.columbia.edu). These nonenveloped, positive-sense, single-stranded RNA viruses have been classified serologically and on the basis of antiviral susceptibility profile, nucleotide sequence relatedness, and receptor usage (2). Phylogenetic analyses of viral protein VP4/VP2 and VP1 coding regions indicate the presence of 74 serotypes in genetic group A and 25 serotypes in genetic group B (2).

Isolated in the 1950s from persons with upper respiratory tract symptoms (2,3), HRVs have become known as the common cold virus because they are implicated in ≈50% of upper respiratory tract infections (4). Large community surveys, including the Virus Watch studies of the 1960–1970s (5), have shed light on some aspects of HRV biology and epidemiology. HRVs were also observed in LRTIs soon after their recognition (3), and data supporting a causative association have accumulated over the past decade (6,7). HRVs have also been implicated in exacerbations of asthma and chronic bronchitis and are increasingly reported in LRTIs of infants, elderly persons, and immunocompromised patients (4).

The Study

The advent of broad-range molecular assays, including multiplex PCR and microarray systems, promises new insights into the epidemiology and pathogenesis of respiratory disease (8,9), given that a laboratory diagnosis is not routinely achieved for a substantial portion of respiratory specimens from symptomatic patients. We recently described the application of a multiplex PCR method for microbial surveillance wherein primers are attached to tags of varying mass that serve as digital signatures for their genetic targets. Tags are cleaved from primers and recorded by mass spectroscopy, enabling a sensitive, inexpensive, and highly multiplexed microbial detection. We used the multiplex MassTag PCR system (10) to investigate respiratory samples that had tested negative during routine diagnostic assessment. This previous study yielded pathogen candidates in approximately one third of cases, and in 8 cases identified a novel genetic clade of picornaviruses divergent from the previously characterized clades, including HRV-A and-B (8). To assess whether this novel clade circulates outside New York state, where it was discovered in cases of influenzalike illness (ILI), we investigated respiratory specimens from Africa, Asia, Australia, and Europe. In most studies (Africa, Asia), the collecting laboratories performed MassTag PCR, and inert mass-tagged amplification products were sent for analysis by mass spectrometry (MS); in other instances (Europe, Australia), inactivated nasopharyngeal swabs or aspirates were sent to New York for MassTag PCR and MS analysis.

Samples in South Africa were collected through a program for comprehensive surveillance of causes of respiratory illness in hospitalized children in the Pretoria area. MassTag PCR was applied to 58 specimens collected during the 2006 season from symptomatic children in their first year of life with no diagnosis available from previous clinical laboratory evaluation. Analysis of amplification products by MS yielded positive signal for HEV/HRV in 14 (24%) samples. Independent amplification and sequence analysis of VP4/2 coding sequence (8) in both laboratories showed sequences that matched the novel genotype in 4 (29%) samples obtained from patients with LRTIs and respiratory distress (Table 1, Figure). Samples collected in Côte d’Ivoire, West Africa, were from symptomatic persons living in the vicinity of Taϊ National Park. This location was the most remote of our study; residents have limited contact with other human populations. In this location, 2 (10%) HRV-A were identified in the 52 samples available for analysis (Table 1, Figure).

Molecular diagnosis of ARI from 7 countries by using MassTag PCR* and VP4/2 sequencing†
CountrySeason(s)SamplesPicornavirus positive
TotalNovel cladeHRV-AHRV-BHEV% MaleAge range (mean/median)
South Africa200658144631710.4–30 mo (5.6/3)
Côte d’Ivoire2006522020010022–28 y (25/25)
Nepal2005–0680174751560.25–56 y (8.5/3)
India20075063300834–36 mo (17.8/18)
Australia20062211001004–6 mo (5/5)
Denmark20077075101571–8 mo (2.9/2)
Spain2003–200614‡146530861–96 mo (23.2/15.5)

*See (10).
†ARI, acute respiratory illness; HRV, human rhinovirus; HEV, human enterovirus.
‡With previous HRV diagnosis.

Phylogenetic analysis of VP4/2 coding region of viruses identified in association with acute respiratory illness (ARI) in South Africa, Côte d’Ivoire, Nepal, India, Western Australia, Denmark, and Spain (sequences deposited in GenBank under accession nos. EU697825–83). Phylogeny of VP4/2 nucleotide sequence (401 nt) was reconstructed by neighbor-joining analysis applying a Jukes-Cantor model; the scale bar indicates nucleotide substitutions per site. Included for reference are sequences belonging to the novel genotype identified in New York State (NY-003, –028, –041, –042, –060, and –074 [8]), similar viruses reported recently (QPM [11]; NAT01 and NAT045 [12]; and 024, 025, 026 [13]), and selected human rhinovirus A (HRV-A) serotypes (GenBank accession numbers for reference sequences are indicated in parentheses); HRV-B serotypes; human enterovirus C (HEV-C) viruses human coxsackievirus A1 and A24 (CV-A1, and CV-A24, respectively); human poliovirus 2 (PV-2); HEV-B viruses human echovirus 5 (E-5), human coxsackievirus B3 (CV-B3), and swine vesicular disease virus (CV-B5); HEV-D viruses human enterovirus 68 and 70 (EV-68, EV-70); porcine enterovirus B virus porcine enterovirus 9 (PEV-9); and HEV-A viruses human coxsackievirus A16 (CV-A16) and human enterovirus 71 (EV-71). SA, South Africa; IC, Côte d’Ivoire; HC, Nepal; IN, India; AUS, Australia; DK, Denmark; SO, Spain.

In Nepal, viruses of the novel genotype were identified in specimens collected during ILI surveillance or outbreaks of respiratory disease. Samples from ILI surveillance activities were collected in Kathmandu and Bharatpur. Outbreak samples were collected in the summer months from camps of >100,000 refugees from Bhutan located in Jhapa, southeast Nepal. Samples represented all age groups and were collected from December 2005 through July 2006. The novel genotype was identified by independent molecular typing in both laboratories in 4 (5%) samples (Table 1, Figure). In India, samples from 50 children with ARI, submitted for routine laboratory analysis during the 2007 season, were evaluated by MassTag PCR. Independent molecular typing in both laboratories indicated the novel genotype in 3 (6%) samples (Table 1, Figure).

Additional sample sets were obtained through main diagnostic laboratories in Western Australia, Denmark, and Spain, representing random respiratory specimens submitted for laboratory analysis. In 1 sample available from Western Australia, the novel genotype was identified in a preterm infant with undiagnosed, wheezy LRTI. The novel genotype was also found in 5 (7%) of 70 samples from Denmark and in 6 (43%) of 14 samples with previously diagnosed HRV infection from Spain (Table 1, Figure).

The 5% overall frequency of the novel genotype across our study samples, representing 34% of all detected picornavirus infections, and its observed global distribution, led us to analyze the accumulating sequence data for insights into their history. Rates of evolutionary change and the Time to the Most Recent Common Ancestor (TMRCA) of the novel clade were estimated by using the Bayesian Markov Chain Monte Carlo approach (BEAST package [14]; ), applying a relaxed molecular clock with an uncorrelated lognormal distribution of rates, a GTR + I + Γ4 model of nucleotide substitution (determined by MODELTEST [15]), and exponential population growth. Statistical uncertainty in each parameter estimate is expressed as 95% highest probability density (HPD) values. The estimated mean rate of evolutionary change was 6.6 × 10–4 substitutions/site/y (95% HPD = 0.3–14.6 × 10–4 substitutions/site/y; 38 dated samples collected over 32 mo (8,16) (S.R. Dominguez et al., unpub. data). Under this rate the mean TMRCA was estimated at 1,800 y, although with wide variance caused by the short sequence available (95% HPD = 279–5,201 y). Despite the inherent sampling error, this analysis suggests that this third clade of rhinovirus has been circulating for >250 years. The diversity observed within the novel clade and its genetic distance from other HRV/HEV were comparable to those seen for HRV-A, -B, or the HEV species (Table 2).

Percentage of intraspecies and interspecies conservation of VP4/2 nucleotide sequence*
VirusesHEV-AHEV-BHEV-CPV†HEV-DHRV-AHRV-BNew clade
HEV-A7261636363596160
HEV-B75646459596159
HEV-C757162616561
PV8160606261
HEV-D83596161
HRV-A806163
HRV-B8060
New clade75

*HEV, human entero
virus; PV, poliovirus; HRV, human rhinovirus.
†PV may be moved by the International Committee on Taxonomy of Viruses into HEV-C.

Conclusions

A clade of picornaviruses recently discovered in New York State is globally distributed and is found in association with community outbreaks of ARI and severe LRTIs of infants. These viruses contribute both to a substantial proportion of previously undiagnosed respiratory illness and to diagnosed, but nontyped cases of HRV infection. Similar viruses were recently characterized also in Queensland, Australia (11); California, USA (12); Hong Kong Special Administrative Region, People’s Republic of China (13); and Germany (16). Our findings indicate the need for further investigation into this third (HRV-C) group of rhinoviruses with emphasis on epidemiology, pathogenesis, and strategies to prevent and ameliorate disease caused by HRV infection.

Suggested citation for this article: Briese T, Renwick N, Venter M, Jarman RG, Ghosh D, Köndgen S, et al. Global distribution of novel rhinovirus genotype. Emerg Infect Dis [serial on the Internet] 2008 Jun [date cited]. Available from http://www.cdc.gov/eid/content/14/6/944.htm

Acknowledgments

We thank Ashlee N. Bennett and Jeffrey Hui for technical assistance.

This work was supported by National Institutes of Health awards AI062705, AI051292, AI059576, HL083850, and AI57158 (Northeast Biodefense Center–Lipkin), the South African National Health Laboratory Service Research Awards, award PI060532 by Fondo de Investigaciones Sanitarias, Instituto de Salud Carlos III, the Robert Koch-Institut, and the Max-Planck-Society. Support by the Ivorian Ministries of the Environment and Forests, of Research and of Health, and the Swiss Research Center, Abidjan, is gratefully acknowledged.

Dr Briese is associate director of the Center for Infection and Immunity and associate professor of epidemiology at Columbia University’s Mailman School of Public Health. His research focuses on the molecular epidemiology of viruses, virus-host interactions, and innovative methods for pathogen detection and diagnosis.

ReferencesBryce J, Boschi-Pinto C, Shibuya K, Black RE WHO estimates of the causes of death in children.Lancet 2005;365:114752 10.1016/S0140-6736(05)71877-815794969Turner RB, Couch RB Rhinoviruses. In: Knipe DM, Howley PM, editors. Fields virology. Philadelphia: Lippincott, Williams & Wilkins; 2007 p. 895–909.Ketler A, Hamparian VV, Hilleman MR Characterization and classification of ECHO 28-rhinovirus-coryzavirus agents.Proc Soc Exp Biol Med 1962;110:8213114455433Hayden FG Rhinovirus and the lower respiratory tract.Rev Med Virol 2004;14:1731 10.1002/rmv.40614716689Fox JP, Cooney MK, Hall CE The Seattle virus watch. V. Epidemiologic observations of rhinovirus infections, 1965–1969, in families with young children.Am J Epidemiol 1975;101:12243164769Papadopoulos NG, Bates PJ, Bardin PG, Papi A, Leir SH, Fraenkel DJ, Rhinoviruses infect the lower airways.J Infect Dis 2000;181:187584 10.1086/31551310837165Mosser AG, Vrtis R, Burchell L, Lee WM, Dick CR, Weisshaar E, Quantitative and qualitative analysis of rhinovirus infection in bronchial tissues.Am J Respir Crit Care Med 2004;171:64551 10.1164/rccm.200407-970OC15591468Lamson D, Renwick N, Kapoor V, Liu Z, Palacios G, Ju J, MassTag polymerase-chain-reaction detection of respiratory pathogens, including a new rhinovirus genotype, that caused influeza-like illness in New York State during 2004–2005.J Infect Dis 2006;194:1398402 10.1086/50855117054069Chiu CY, Rouskin S, Koshy A, Urisman A, Fischer K, Yagi S, Microarray detection of human parainfluenzavirus 4 infection associated with respiratory failure in an immunocompetent adult.Clin Infect Dis 2006;43:e716 10.1086/50789616983602Briese T, Palacios G, Kokoris M, Jabado O, Liu Z, Renwick N, Diagnostic system for rapid and sensitive differential detection of pathogens.Emerg Infect Dis 2005;11:310315752453Drummond AJ, Rambaut A BEAST: Bayesian evolutionary analysis by sampling trees.BMC Evol Biol 2007;7:214 10.1186/1471-2148-7-21417996036Posada D, Crandall KA MODELTEST: testing the model of DNA substitution.Bioinformatics 1998;14:8178 10.1093/bioinformatics/14.9.8179918953Renwick N, Schweiger B, Kapoor V, Liu Z, Villari J, Bullmann R, A recently identified rhinovirus genotype is associated with severe respiratory tract infection in children in Germany.J Infect Dis 2007;196:175460 10.1086/52431218190255McErlean P, Shackelton LA, Lambert SB, Nissen MD, Sloots TP, Mackay IM Characterisation of a newly identified human rhinovirus, HRV-QPM, discovered in infants with bronchiolitis.J Clin Virol 2007;39:6775 10.1016/j.jcv.2007.03.01217482871Kistler A, Avila PC, Rouskin S, Wang D, Ward T, Yagi S, Pan-viral screening of respiratory tract infections in adults with and without asthma reveals unexpected human coronavirus and human rhinovirus diversity.J Infect Dis 2007;196:81725 10.1086/52081617703411Lau SK, Yip CC, Tsoi HW, Lee RA, So LY, Lau YL, 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 2007;45:365564 10.1128/JCM.01254-0717804649