key: cord-0745225-2y7eh9ao
authors: Naughtin, Monica; Sareth, Rith; Sentilhes, Anne-Charlotte; Vong, Sirenda; Joffret, Marie-Line; Cornillot, Emmanuel; Deubel, Vincent; Delpeyroux, Francis; Frutos, Roger; Buchy, Philippe
title: Genetic diversity of human rhinoviruses in Cambodia during a three-year period reveals novel genetic types
date: 2015-07-29
journal: Infect Genet Evol
DOI: 10.1016/j.meegid.2015.07.030
sha: 6b98d0fa40066d7036ebe9f283a20162a66e1a3b
doc_id: 745225
cord_uid: 2y7eh9ao

Acute respiratory viral infections are a major cause of morbidity during early childhood in developing countries. Human rhinoviruses are the most frequent cause of upper respiratory tract infections in humans, which can range in severity from asymptomatic to clinically severe disease. In this study we collected 4170 nasopharyngeal swabs from patients hospitalised with influenza-like illness in two Cambodian provincial hospitals between 2007 and 2010. Samples were screened for 18 respiratory viruses using 5 multiplex PCRs. A total of 11.2% of samples tested positive for human rhinoviruses (HRV). VP4/2 and VP1 regions were amplified and sequenced to study the distribution of rhinoviruses genotypes and species in Cambodia during this three-year period. Five novel genotypes, 2 species A, 2 species B and 1 species C were identified based on VP1 sequences. Co-infections with other viruses were demonstrated.

Human rhinoviruses are highly prevalent respiratory picornaviruses of the Enterovirus genus. They display a single stranded positive sense RNA genome of approximately 7200 nucleotides encoding a single polyprotein which is cleaved by viral proteases into 11 individual proteins. Human rhinoviruses are amongst the most frequent causes of upper respiratory tract infections, and are significantly associated with several lower respiratory tract infections including acute wheezing, pneumonia, and bronchiolitis, and cause about one half of asthma exacerbations (Busse et al., 2010) . One hundred classical serotypes have been defined by 20-fold or greater decrease in antigen cross-reactivity and form two genetically distinct species, A and B. There has recently been a renewed interest in rhinoviruses following the discovery of the highly prevalent and diverse species C group, using novel molecular-based approaches (McErlean et al., 2007; Briese et al., 2008; Kiang et al., 2008; Lee et al., 2007; Lau et al., 2007) .

Species B viruses are generally considered to be rare, with an average prevalence rate of around 7% of total HRV strains (Briese et al., 2008; Miller et al., 2009a) . Several studies suggest that species C viruses are a major cause of wheezing and exacerbation of asthma in children, and may have different clinical outcomes to other species [Lau et al., 2007; Miller et al., 2009b) , however, other studies find no significant differences (Piotrowska et al., 2009; McCulloch et al., 2014) .

The recent whole-genome sequencing of all species A and B, and several species C prototype strains provides a very useful framework for the genetic analysis and classification of clinical specimens, from partial or whole sequences (Palmenberg et al., 2009) . Recently, new criteria have been proposed for the classification of species C viruses, which has simplified the identification of novel genotypes, or ''genetic types.'' These criteria classify 33 genetic types, and 28 provisional types which require confirmation by additional VP1 gene sequencing . Within A and B species, there is generally a strong adherence to phylogenetic positioning (branching order) between the different structural protein genes, with a few exceptions (Savolainen et al., 2004; McIntyre et al., 2010) . Species C viruses cluster in two phylogenetic groups on the 5' untranslated region (5UTR), one with species A (designated in some studies as Ca) and the other as a 

Infection, Genetics and Evolution j o u r n a l h o m e p a g e : w w w . e l s e v i e r . c o m / l o c a t e / m e e g i d separate C cluster (designated Cc) (Palmenberg et al., 2010; Huang et al., 2009 ). There is a consensus view that partial genome sequencing of at least two regions (VP4/2, VP1 or 3CD genes) of the HRV genome is sufficient for genetic type classification based on nucleotide homology Arden and Mackay, 2010) . Several recent studies support classification guidelines based on pairwise nucleotide homology distance thresholds, based on studies across the genome, and these data are supported by serotypical distinctions in A and B species. The proposed inter-serotype pairwise nucleotide divergence threshold for A and B species is 13% in the VP1 gene and 10% in the VP4/2 gene McIntyre et al., 2010 McIntyre et al., , 2013 Savolainen et al., 2002a; Wisdom et al., 2009a) . These criteria also apply to the inter-genotype distance for species C viruses, which cannot be classified by serology due to their inability to grow in tissue culture McIntyre et al., 2010) . As a consequence, novel HRV strains were described in the past few years, particularly novel species C viruses (Briese et al., 2008; Renwick et al., 2007; Miller et al., 2009b; McIntyre et al., 2010; Wisdom et al., 2009b; Kistler et al., 2007) . Potentially novel species A viruses are rarely reported, however, recently the strain A101 strain has been observed in several studies worldwide .

Recombination, which occurs frequently in enteroviruses, was presumed to be uncommon in rhinoviruses, however, several studies have identified that recombination frequently occurs in the 5UTR region Huang et al., 2009; Tapparel et al., 2009) . Several studies have recently compared the location of ''breakpoints'' of recombination, to reveal typical recombination hotspots in the 5UTR Huang et al., 2009) . Recombination events in the coding regions seem to be rare Tapparel et al., 2009) .

In this work we have evaluated the incidence of acute lower respiratory tract illness (ALRI)-associated rhinoviruses infection during a 3-year hospital-based survey in Cambodia in two distinct geographical locations, and examined the genetic diversity of these viruses.

4170 nasopharyngeal respiratory samples were collected from patients of all ages hospitalised in Cambodia's Kampong Cham and Takeo provinces referral hospitals during the period between June 2007 and December 2009. The case definition of ALRI was as follows; for patients aged <5 years, ALRI was defined as an illness of <10 days duration that consisted of cough or difficult breathing, plus tachypnea. For patients aged between 5 and 14 years, case definition included the above symptoms as well as a fever >38°C. For patients 15 years or older ALRI was defined as a fever >38°C, plus chest pain or tachypnea or auscultatory crackles (Mardy et al., 2009) . A severe case for patients <5 years of age was defined as ALRI in addition to: respiratory rate P60/mn (age <2 months), respiratory rate P50/mn (age <11 months), respiratory rate P40/mn (aged between 1 and 5 years), saturation <93%, cardiac frequency >180 bpm (age <1 year), cardiac frequency >140 bpm (aged between 1 and 5 years), clinical respiratory distress, according to World Health Organisation (WHO) guidelines (WHO, 2005) .

Samples were screened for the presence of 18 common and novel viral respiratory pathogens, including HRSV, human metapneumovirus (HMPV), human bocavirus (HBoV), influenza A and B viruses, human coronaviruses OC43, 229E, HKU1, and NL63, severe acute respiratory syndrome-associated coronavirus (SARS-CoV), parainfluenza viruses 1 to 4, adenoviruses, human rhinoviruses, and enteroviruses, by multiplex reverse transcriptase PCR (RT-PCR)/PCR, as previously reported (Buecher et al., 2010; Arnott et al., 2011 Arnott et al., , 2013a Arnott et al., , 2013b . The RT-PCR step was amplifying an approximately 549 bp fragment of the human rhinoviruses VP4/VP2 and 5 0 -non-coding region (Savolainen et al., 2002b) and the subsequent semi-nested PCR was producing a 450 bp amplicon (Bellau-Pujol et al., 2005) .

Samples positive for rhinoviruses were amplified in the 5UTR and VP4/2 regions using the primers VP4/2F and VP4/2R (Supplementary Table 1 ) (Lau et al., 2007; Huang et al., 2009) . Certain samples could not be amplified using these primers and were amplified separately in the 5UTR and VP4/2 region (Supplementary Table 1 ). Amplification was performed with the Qiagen one-step RT-PCR kit, as described previously (Huang et al., 2009) . For VP1 amplification, RNA was reverse transcribed with random hexamer primers (Promega, United Kingdom) using the superscript III reverse transcriptase (Invitrogen). PCR was then conducted using Gotaq polymerase (Promega), and species-specific primers which have been previously described (Supplementary Table 1 ) (Tapparel et al., 2009) . For VP1 region amplification we used primers VP1 1.2F and VP1 1.18R (and nested PCR VP1 1.55F and VP1.18R) for species A, VP1 1.57F and VP1 2.19R (and nested PCR VP1 1.57F and VP1 1.18R for species B), VP1 P1.160F and VP1 2.68R (and nested VP1 P1.161F and VP1 2.69R) for species C. A touch-down program (annealing temp of 55-1°C per cycle for 10 cycles, then 30 cycles at 45°C) was used for both primary and nested PCRs. DNA products were separated by agarose gel electrophoresis and extracted using the QIAquick Gel Extraction Kit (Qiagen, Germany). Products were ligated into TOPO-TA vectors (Invitrogen, Carlsbad, CA). At least two plasmids were sequenced for each clinical sample. At least two recombinant clones were sequenced on an ABI3730 XL automatic DNA sequencer by Macrogen (Seoul, South Korea).

Alignments were conducted with Muscle, in SeaView 4.2.6 (Gouy et al., 2010) . Owing to the very high level of genome saturation detected, no Maximum-Likelihood analysis was conducted for building trees. Trees were therefore built using a neighbour-joining approach. The robustness of nodes was assessed with 500 bootstrap replicates.

Genetic similarity and difference matrices were constructed using Muscle in SeaView 4.2.6 (Gouy et al., 2010) . Sequence polymorphism, subsequent tests and statistical tests were performed using functions from the DnaSP 5.10.01 package (Librado and Rozas, 2009 ). The numbers of haplotype (alleles) were analysed. Sequence polymorphism was investigated by calculation of several measures. Nucleotide diversity, Pi (p), the average number of nucleotide differences per site between two sequences was calculated, using the Jukes and Cantor correction. Theta (Watterson's mutation parameter) was calculated for the whole sequence from S. Eta (g) is the total number of mutations, and S is the number of polymorphic sites. Ka (the number of non-synonymous substitutions per non-synonymous site), and Ks (the number of synonymous substitutions per synonymous site) for any pair of sequences were calculated. Tajima's D test was used for statistical testing the hypothesis that all mutations are selectively neutral.

Statistical tests D ⁄ and F ⁄ were used for testing the hypothesis that all mutations are selectively neutral.

The project (SISEA) was approved by the National Ethics Committee of Cambodia. All patients/parents of sick children who participated provided written informed consent.

The sequences described in this study were deposited in Genbank under the following accession numbers: VP1: KF034075, KF034076, KF034078 to KF034081, KF034083, LN623979 to LN623997, KM986111, KM986112 VP 2/4: KF034009 to KF034074 5 0 UTR sequences use for recombination were deposited under accession numbers KF033904 to KF034008.

The SISEA project (Surveillance of Infection in South East Asia) was established for surveillance of respiratory viruses in Cambodia, Vietnam, and Laos. Over a three-year period (2007-2009), we collected a total of 4170 respiratory samples from in-patients hospitalised with acute lower respiratory tract illness (ALRI) at two provincial hospitals in Cambodia (Takeo, southern Cambodia and Kampong Cham, central-north Cambodia).

Amongst the samples taken from Takeo hospital there were 977 children and 2279 adults. From the Kampong Cham hospital there was 584 children and 227 adults. RNA was extracted from nasopharyngeal (NPS) samples, and amplified by real time PCR to detect the presence of 18 viruses, as described in previous studies (Buecher et al., 2010) . From 4170 patient samples, 1105 samples tested positive for one or more virus. HRV was detected in 455 samples, representing the most frequently detected virus (10.9% of all patients, and 41.2% of patients positive for viral detections), consistent with evidence that rhinoviruses are usually highly prevalent (Miller et al., 2009a; Huang et al., 2009; Buecher et al., 2010; Follin et al., 2009; Loens et al., 2006; Vong et al., 2013; Guerrier et al., 2013) . Two hundred sixty-one of these patients were male, 194 female. Two hundred thirty-seven patients were under the age of 2 years (52%), with 273 (60%) under the age of 5, and 80 patients over the age of 50 (17.6%). As previously reported, ALRI in general and HRV infections in particular were observed throughout the year, without clear seasonality (Guerrier et al., 2013) . In addition to typical symptoms of ALRI recorded for all patients, supplementary comprehensive clinical data was available from 94 patients, which allowed us report on typical symptoms associated with disease (Table 1) , and also to form a case definition to identify severe rhinoviral infection. This definition included the following criteria; no bacterial infection (sputum-negative) or viral co-infection. For children >5 years, or adults, was the inclusion of at least two of the following; a systolic blood pressure of 690 mmHg, a cardiac frequency of P120 bpm, a respiratory rate of P30/mn, saturation of 693%, or a temperature of 635°C or P40°C. The criteria for children 65 years was a respiratory rate of either P60/mn (aged < 2 months), P50/mn (aged between 2 and 12 months), P40/mn (aged between 1 and 5 years), oxygen saturation of <93%, a cardiac frequency of >180 bpm (aged <1 year), or >140 bpm (aged between 1 and 5 years), and signs of severe respiratory distress. From this analysis 16 samples fulfilled the criteria for severe rhinoviral infection (approximately 17%).

In order to evaluate the genetic diversity of HRV circulating in the considered human population and the potential involvement of variants or novel types in the occurrence of the severe forms recorded, samples were analysed by PCR. PCR amplification was targeted on the VP1 gene and the VP4/2 gene region. A total of 28 VP1 sequences were obtained and genotyped ( Table 2) . The VP1 PCR failed on a number of samples, suggesting sequence divergence from the primer sequences used. Out of these 28 viruses 21 belonged to species A, 5 were members of species B and 2 corresponded to species C. Altogether, 5 viruses met the pairwise distance thresholds defined by McIntyre et al. (2013) to be recognised as novel types. 2 belonged to species A (i.e. D386 and K590 with 19.7% and 16.6% distance, respectively) and 2 others belonged to species B (i.e. K606 and P128 with 17% and 22.9% distance, respectively). The last sample corresponding to a new type, i.e. R005, belonged to species C and displayed a pairwise distance to the closest relative of 30.9%.

The VP4/2 sequence was also considered for genotyping the viruses. Sixty VP4/2 sequences were obtained which were distributed into 31, 11 and 18 sequences for species A, B and C, respectively (Table 3) . We also detected one enterovirus (EV68), which was previously recognised as HRV87 until it was reclassified as a variant of EV68 (Ishiko et al., 2002) . Viruses which were defined as clinically severe are highlighted in Table 3 . Sixteen viruses in species A displayed pairwise distances at nucleic acid level to the closest relative ranging from 10% to 19.9%. With respect to species B, 10 viruses out of 11 displayed pairwise distances ranging from 15.8% to 18.8%. The only virus from species C showing a VP4/2 pairwise distance above the 10% threshold to the closest relative was R005 (29.1%) which was identified as a novel type based on VP1 sequence distances. Interestingly, the VP4/2 closest relative, i.e. HRV-C15, was different than the VP1 closest relative, i.e. HRV-C23. However, when considering the pairwise distances at the protein level, a different feature was seen. Several viruses for which the nucleotide pairwise distance with the closest relative was clearly above the 10% threshold displayed a protein sequence 100% identical to that of the same closest virus. Interestingly, several of these putative new types based on nucleotide distance, ,displayed amino-acid distances lower than sample 

Owing to the variations observed between the nucleic acid and protein distances, nucleotide polymorphism and genome saturation were assessed to determine the solidity of a phylogenetic analysis. The ratio of observed transition (mutation from a purine nucleotide to another or a pyrimidine nucleotide to another) versus transversion (mutation for a purine nucleotide to pyrimidine one or from a pyrimidine nucleotide to purine one) was calculated as an indicator of genome saturation. Data are summarised in Supplementary Table 2. The ratio ranged from 0.714 to 1.086 for the coding sequences depending on the gene and virus species.

Owing to the presence of indels in the VP1 sequences, calculation was done on indel-free segments. For the 5'UTR, the transition/transversion ratio was 0.515. For reference, the transition/transversion ratio was calculated for the gene E of the dengue DENV-1 virus (Duong et al., 2013) with a value of 14.49. The rate of polymorphic sites (S) ranged from 44.3% to 71.3 5% depending on the sequence while H ranged from 42.8 to 101.8. The DENV-1 E gene used as a reference displayed a rate of polymorphic sites of 12.7% and a H value of 40.1. The ratio between parsimony informative sites (Pa) and singletons (Si) was strongly biased towards the former with a rate of Pa ranging from 88.5% to 96.1% for the coding regions while it was 84.3% for the 5'UTR. For reference, the values of Pa and Si were very close for the DENV-1 gene E, although higher for Si, with rates of 48.7% and 51.3%, respectively. A strong bias was also observed in the rate of synonymous replacement (Ns) which ranged from 86.2% to 96.4%, the rate for the DENV-1 E gene being in the same range with 84.9%. Accordingly, a very low value of the Ka/Ks was observed for all sequences. No bias was observed on the individual rate of A, T, G and C bases.

Phylogenetic analysis was conducted on VP1 and VP4/2 ( Fig. 1a and b ). Owing to the high level of genome saturation observed when analysing the transition/transversion ratio, the topology of trees based on protein alignments was also determined for the VP1 and VP4/2 proteins in order to assess congruent topologies between nucleic acids and proteins. The high level of genome saturation could not allow for running relevant evolution models for maximum likelihood analysis and trees were therefore built using the distance neighbour-joining method. For both VP1 and VP4/2 nucleic acids and protein trees were perfectly congruent (Supplementary Figs. 1 and 2) .

Several studies have examined recombination and found that the region between the VP4/2 and 5UTR could be a ''hotspot'' for recombination events (Huang et al., 2009; Tapparel et al., 2009; McIntyre et al., 2010) . A similar analysis was conducted and the putative recombination breakpoints for the viruses from species A and C observed in this study were found to be located in the same area as previously described (Huang et al., 2009 ) (Data not shown). 

This work is one of the largest respiratory virus surveillance studies, with more than 4000 clinical samples, providing evidence of the major rank occupied by HRV amongst respiratory viruses causing influenza-like illness. Rhinoviruses are frequently concurrently detected with other respiratory viruses due to their high circulation in the community. However, some studies have proposed that despite the frequency of co-detections, rhinoviruses are statistically the least likely to be involved in coinfection and may actually reduce the probability of another virus infecting the host (Duong et al., 2013) . HRV/HRV co-detections have Table 3 Species classification of clinical detections based on VP4/2 sequences.

Strains which were obtained from clinically severe hospitalised patients are highlighted in grey. Novel types identified based on VP1 distances are shown in bold italic. occasionally been reported Renwick et al., 2007; Huang et al., 2009) , although, interestingly no co-detections from the same species have been reported. In this work we report one HRV/HRV co-detection, and significantly both viruses detected from one patient were species A (HRV78 and HRV81). The frequency of HRV co-detections is difficult to assess given that often only a few clones are sequenced and therefore if one virus is present in a lesser proportion than the other, it might not be detected. In this work we analysed 455 viruses and we found out that 35 only out of these 455 were co-detected with other viruses (7.7% of total). This percentage is even lower than that reported in previous studies, i.e. 23% of total [40] which were concluding that rhinoviruses were statistically the least likely respiratory virus to be involved in co-infection (Greer et al., 2009) . This further highlights the independent role of rhinoviruses as major aetiological contributors to severe respiratory disease.

Few studies to date have been able to link clinical and virological data to identify strains with increased virulence. The viruses from clinically severe patients in this study comprised 6 members of species A, including one sample resembling the A101 novel genetic type , 1 virus from species B, 3 viruses from species C and 1 enterovirus (EV-68). Previous studies have suggested that species C viruses may be associated with increased disease severity (Lauinger et al., 2013) . However, no identifiable clustering of patterns of the clinically severe viruses could be seen in this work. The reported clinically severe viruses were spread in different genotypic groups, indicating that severity might not be an evolutionary pattern strictly related to viruses or sequences analysed in this work.

Since the characterisation of the prototype strains by serology in the 1960's there have been very few reports of novel viruses from A or B species. A recent study sequenced the entire genome of a virus designated A101, which has a nucleotide identity across the genome of only 76.6% with the nearest species A virus (HRV-53) . Another study conducted in Thailand found species A to be the most common species detected regardless of age groups (Fry et al., 2011) . Following the standards developed by McIntyre et al. (2013) five novel genotypes were identified in this work with respect to the VP1 gene. Two viruses from species A (i.e. D386 and K590), two viruses from species B (i.e. K606 and P128) and one virus from species C (i.e. R005). The divergence at the nucleic acid level was high enough above the threshold set by McIntyre et al. (2013) to consider these viruses actual new genotypes. This conclusion is strongly supported for three viruses, i.e. D386, P128 and R005, by the high divergence also observed at the protein level. This protein divergence was however very limited for the last two VP1 variants (i.e. K606 and K590). A similar situation is encountered when using VP4/2 as a target region for genotyping. Most of the putative variants identified on nucleotide distance display a very limited divergence with the closest relative at protein level (less than 3%) and in several cases there is a perfect identity. The best example found in this work is perhaps the novel VP1 type, R005, which displays the highest divergence to the closest relative of all novel types based on the VP1 gene and protein sequences and the highest divergence to the closest relative on the VP4/2 gene sequence but 100% VP4/2 protein sequence identity. In addition, the closest relative to R005 is not the same according to VP1 or VP4/2.

The polymorphism and genome saturation analyses conducted in this work suggest that HRV viruses might be exposed to two opposite trends: a very high rate of mutations generating genome saturation and a strong negative selective pressure resulting in highly conserved proteins. It is therefore difficult when in presence of conflicting nucleic-acid and amino-acid distance data to formally assign these viruses to a novel type. We suggest that the typology standards developed by McIntyre et al. (2013) should be completed with protein thresholds into order to assign unambiguously these conflicting detections to specific types.

The rhinoviruses studied in this work can be clearly discriminated into three genetic distinct taxa, i.e. the already recognised species A, B and C. This is especially true for the VP1 gene which displays indels. However, the relative position of the detected viruses within each taxon and their phylogeny cannot be established. Both the high rate of polymorphism and level of saturation observed can easily result in the occurrence of reverse mutations and blurred phylogeny. The clear genetic clustering and the high negative selective pressure observed on both the VP1 and VP4/2 genes suggest that the three species described might correspond to virus groups having undergone independent evolution and specialisation.

This work represents the insights into HRV genetic diversity in South East Asia and provides the first comprehensive view of the HRV genetic diversity and the association between disease severity and infecting genotype in Cambodia during a 3-year period. This work also provides the first demonstration of the very high level of genome saturation and nucleotidic polymorphism of HRV indicating that typing and classification standards should comprise protein distance thresholds. Finally, it shows that much more research is needed to understand the complex dynamics of this virus.

Newly identified human rhinoviruses: molecular methods heat up the cold viruses

A study of the genetic variability of human respiratory syncytial virus (HRSV) in Cambodia reveals the existence of a new HRSV group B genotype

Human bocavirus amongst an all-ages population hospitalised with acute lower respiratory infections in Cambodia

Genetic variability of human metapneumovirus amongst an all age population in Cambodia between

Development of three multiplex RT-PCR assays for the detection of 12 respiratory RNA viruses

Global distribution of novel rhinovirus genotype

Use of a multiplex PCR/RT-PCR approach to assess the viral causes of influenza-like illnesses in Cambodia during three consecutive dry seasons

Role of viral respiratory infections in asthma and asthma exacerbations

Genetic diversity and lineage dynamic of dengue virus serotype 1 (DENV-1) in Cambodia

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

Human rhinovirus infections in rural Thailand: epidemiological evidence for rhinovirus as both pathogen and bystander

SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building

Do rhinoviruses reduce the probability of viral co-detection during acute respiratory tract infections?

Acute viral lower respiratory tract infections in Cambodian children: clinical and epidemiologic characteristics

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

Human rhinovirus 87 identified as human enterovirus 68 by VP4-based molecular diagnosis

Assay for 5 0 noncoding region analysis of all human rhinovirus prototype strains

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

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

Patient characteristics and severity of human rhinovirus infections in children

A diverse group of previously unrecognized human rhinoviruses are common causes of respiratory illnesses in infants

DnaSP v5: a software for comprehensive analysis of DNA polymorphism data

Detection of rhinoviruses by tissue culture and two independent amplification techniques, nucleic acid sequence-based amplification and reverse transcription-PCR, in children with acute respiratory infections during a winter season

Influenza activity in Cambodia during

Severity of rhinovirus infection in hospitalized adults is unrelated to genotype

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

Proposals for the classification of human rhinovirus species A, B and C into genotypically assigned types

Analysis of genetic diversity and sites of recombination in human rhinovirus species C

Human rhinovirus C associated with wheezing in hospitalised children in the Middle East

A novel group of rhinoviruses is associated with asthma hospitalizations

Analysis of the complete genome sequences of human rhinovirus

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

Rhinoviruses are a major cause of wheezing and hospitalization in children less than 2 years of age

Full-genome sequence and analysis of a novel human rhinovirus strain within a divergent HRV-A clade

A recently identified rhinovirus genotype is associated with severe respiratory-tract infection in children in Germany

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

Phylogenetic analysis of rhinovirus isolates collected during successive epidemic seasons

Genetic clustering of all 102 human rhinovirus prototype strains: serotype 87 is close to human enterovirus 70

Proposals for the classification of human rhinovirus species C into genotypically-assigned types

New respiratory enterovirus and recombinant rhinoviruses among circulating picornaviruses

Acute lower respiratory infections in P5 year-old hospitalized patients in Cambodia, a low-income tropical country: clinical characteristics and pathogenic etiology

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

Genetics, recombination and clinical features of human rhinovirus species C (HRV-C) infections; interactions of HRV-C with other respiratory viruses

The SISEA project was funded by the Agence Française pour le Développement. This work was also supported by the Agence Nationale pour la Recherche (ANR 09 MIEN 019), the Fondation pour la Recherche Médicale (FRM DMI20091117313).

Supplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.meegid.2015.07. 030.