Genome-wide diversity and selective pressure in the human rhinovirus
© Kistler et al; licensee BioMed Central Ltd. 2007
Received: 13 January 2007
Accepted: 03 May 2007
Published: 03 May 2007
The human rhinoviruses (HRV) are one of the most common and diverse respiratory pathogens of humans. Over 100 distinct HRV serotypes are known, yet only 6 genomes are available. Due to the paucity of HRV genome sequence, little is known about the genetic diversity within HRV or the forces driving this diversity. Previous comparative genome sequence analyses indicate that recombination drives diversification in multiple genera of the picornavirus family, yet it remains unclear if this holds for HRV.
To resolve this and gain insight into the forces driving diversification in HRV, we generated a representative set of 34 fully sequenced HRVs. Analysis of these genomes shows consistent phylogenies across the genome, conserved non-coding elements, and only limited recombination. However, spikes of genetic diversity at both the nucleotide and amino acid level are detectable within every locus of the genome. Despite this, the HRV genome as a whole is under purifying selective pressure, with islands of diversifying pressure in the VP1, VP2, and VP3 structural genes and two non-structural genes, the 3C protease and 3D polymerase. Mapping diversifying residues in these factors onto available 3-dimensional structures revealed the diversifying capsid residues partition to the external surface of the viral particle in statistically significant proximity to antigenic sites. Diversifying pressure in the pleconaril binding site is confined to a single residue known to confer drug resistance (VP1 191). In contrast, diversifying pressure in the non-structural genes is less clear, mapping both nearby and beyond characterized functional domains of these factors.
This work provides a foundation for understanding HRV genetic diversity and insight into the underlying biology driving evolution in HRV. It expands our knowledge of the genome sequence space that HRV reference serotypes occupy and how the pattern of genetic diversity across HRV genomes differs from other picornaviruses. It also reveals evidence of diversifying selective pressure in both structural genes known to interact with the host immune system and in domains of unassigned function in the non-structural 3C and 3D genes, raising the possibility that diversification of undiscovered functions in these essential factors may influence HRV fitness and evolution.
Human rhinoviruses (HRV) are the major cause of the common cold, accounting for as much as 80% of upper respiratory infections in the fall cold season (reviewed in ). In the United States, the common cold is estimated to account for approximately 1 billion upper respiratory infections per year, 22 million days of missed school, and $40 billion in direct and indirect costs due to lost work and productivity . Thus, despite typically presenting as a mild, self-limited upper respiratory infection, HRVs exact a significant health and economic burden on society in general. Moreover, recent evidence suggests that HRV infections may not always be mild or restricted to the upper respiratory tract. Results from in vitro and in vivo experimental studies have demonstrated that HRVs can both penetrate and damage bronchial epithelial cells in the lower respiratory tract [3–8]. HRV infections can cause acute bronchitis in healthy children and adults (especially the elderly), precipitate exacerbations in patients with asthma, chronic obstructive pulmonary disease, and cystic fibrosis, and can lead to fatal pneumonia in immunocompromised patients (reviewed in [9–12]).
Despite the ubiquity of HRV infections among healthy populations and their potentially severe clinical consequences in vulnerable populations, no preventive or curative therapies are currently available. Development of such therapies against HRV has in large part been hampered by the great diversity within the HRV genus, and the fact that multiple serotypes co-circulate during each cold season. This diversity has been traditionally characterized via a set of distinct types of phenotypic assays. Antisera neutralization studies performed in the 1960s to 1970s identified 102 distinct HRV serotypes . Subsequent drug susceptibility analysis divided these 102 HRV prototype strains into two major groupings, subgroup A (HRVA), with 77 serotypes, and subgroup B (HRVB), with 25 serotypes . A single serotype, HRV87, falls into neither of these two groups and is actually more similar to human enteroviruses (HEVs) than human rhinoviruses [15, 16]. Identification of two cellular receptors for HRV further divided these serotypes into 2 additional groups [17, 18]: the major cellular receptor (intracellular adhesion molecule 1, ICAM1) group, composed of 90 HRV serotypes [19, 20], and the minor cellular receptor (low density lipoprotein receptor, LDLR) group, made up of 11 HRV serotypes .
More recent molecular genetic analyses of a number of subgenomic regions of HRV have largely corroborated these phenotypic classifications of the HRVs [17, 22–29]. However, due to the paucity of available HRV genome sequences, it is unclear how well the diversity detected in these assays reflects the genome-wide diversity present among the characterized HRV serotypes. The genomes of only six HRV serotypes are publicly available (HRV2 , HRV16 , HRV1b , HRV14 [33, 34], HRV89 , and HRV39 ). These genome sequences represent only a small fraction of the HRV genomic sequence space, and provide limited insight into the genome-wide diversity within this genus, or how this diversity is generated and continues to propagate from year to year.
Here, we expand this set of 6 fully sequenced HRV genomes to a more representative set of 34 genomes through whole genome shotgun sequencing of 27 diverse HRV reference serotypes and a single clinical isolate of HRV associated with an outbreak of severe lower respiratory illness in an elder care facility in Santa Cruz, CA . We have used this larger and more diverse set of HRV genomes to analyze the genome-wide diversity in HRVs and to determine the selective pressure operating at each codon of the HRV genome. Mapping these selective pressure data onto available three dimensional HRV protein structures relative to known functional domains has provided insight into the underlying biology driving evolution of these HRV prototypes and serves as a springboard for future analyses of novel and currently circulating HRVs and the drugs developed to inhibit them.
Generation of a representative set of HRV genome sequences for analysis
In order to obtain an accurate picture of the genetic diversity and selective pressure across the HRV genome, our first task was to expand the set of 6 fully sequenced HRV serotypes to a larger set of HRV genomes that more fully captured the genetic diversity of the known set of 102 serotypes. Since the capsid region has been found to be the most variable portion of other fully sequenced picornavirus genomes [38, 39], we utilized previously generated capsid gene phylogenies of the 102 HRV serotypes [25, 26, 28] to identify an additional set of HRV serotypes that would prove most informative for our analysis. We identified 28 additional serotypes from across the HRV gene capsid phylogenies (Additional File 1, Figure S1) that yielded selective pressure results for the VP1 gene that were well-correlated with the results obtained from the full set of 102 HRV serotype VP1 gene sequences (Materials and Methods, Additional file 1, Figure S2). We thus focused our whole genome shotgun sequence analysis efforts on recovery of genome sequence from these 28 HRV serotypes. Combined with the 6 previously sequenced HRV genomes and the rhino/entero HRV87 genome, this provided a larger, more representative set of 35 HRV genomes for further analysis.
Consistent phylogenetic pattern observed at every locus of the HRV genome
With this expanded set of HRV genomes in hand, we next examined the agreement between the HRV genomic and subgenomic phylogenies. Prior comparative sequence analysis of two other picornaviruses, the human enteroviruses (HEVs) and the Foot-and-Mouth Disease viruses (FMDVs) have uncovered significant incongruences between the genomic and subgenomic phylogenies of these viruses that suggest that recombination plays a significant role in generating diversity in the picornavirus family [38, 40–42]. Comparison of the phylogenies of more extensively sequenced structural and non-structural subgenomic regions of the HRV genome have suggested that similar phylogenetic incongruences may be present in the HRV genome [25, 26, 28, 29]. However, more recent analysis of the prior set of 5 fully sequenced HRVA genomes and a review of the subgenomic data has cast doubt on these conclusions .
Pairwise sequence analysis shows consistent diversity across the genome
Average % pairwise identity detected among HRV nucleotide sequences
HRV ave (min, max)*
HRVA ave (min, max)**
HRVB ave (min, max)#
HRVAvsHRVB ave (min, max)$
Average % pairwise identity detected among HRV amino acid sequences
HRV ave (min, max)*
HRVA ave (min, max)**
HRVB ave (min, max)#
HRVAvsHRVB ave (min, max)$
Recombination scan predicts only small, scattered events in the HRV genome
To directly compare the type and frequency of recombination events in HRV relative to other members of the picornavirus family, we performed a genome-wide scan for recombination events among the fully sequenced HRV genomes (Materials and Methods). This analysis identified ten putative recombination events (Additional file 2, Table S1). However, in contrast to the large-scale single crossover events that have been previously detected between the structural and non-structural genes of HEV and FMDV genomes [38–44], all of the events detected in the HRV genomes were small in size (average length: 281 bp, range: 84–474 bp) and predicted to result from double crossover events localized mainly in the 5'NCR of the genome and a few distinct loci scattered throughout the coding region of the genome (Additional file 2, Table S1). Thus, the extent and scope of recombination predicted to have occurred in these representative HRV genomes is indeed quite different from that seen for HEVs and FMDVs.
Selective pressure across the human rhinovirus genome
We next investigated how HRV diversity might have arisen by analyzing the types of evolutionary forces acting on the HRV genome. We utilized the genome-based HRV phylogeny and the available genome sequences to compute the ratio of non-synonymous to synonymous changes (dN/dS) for each codon in the HRVA and HRVB genomes (Materials and Methods). Such calculations allowed us to create selective pressure profiles for the HRVA and HRVB genomes as a whole, providing an overview of the evolutionary landscape of the HRV genome (Figure 2D).
Overall, we detected similar selective pressure profiles for the HRVA and HRVB genomes (Figure 2D). Intriguingly, this selective pressure analysis reveals that a large proportion of the genome is under purifying selective pressure (82.65% for HRVA and 86.74% for HRVB), exhibiting codon-specific dN/dS ratios at the lower limits of detection (<0.06), despite the high level of genetic diversity we detected across the HRV genomes by scanning pairwise analysis. However, this purifying selective pressure is not distributed uniformly across the genome. It predominates in the central region of the genome that includes a set of non-structural genes (2A, 2B, 2C, 3A, and 3B) that interact with both viral factors and essential host cell factors during the viral replication cycle, and is also detectable across the majority of the 1A gene, which encodes the VP4 capsid protein that assembles on the interior side of the viral particle. Interrupting these regions of purifying selective pressure are two major clusters of residues with elevated dN/dS values: one in a subset of the structural genes (1B, 1C, and 1D) which lie on the outer surface of the viral capsid, and another in a pair of the non-structural genes (3C and 3D) which encode a protease and polymerase essential for viral replication.
Structure-function mapping of diversifying residues in structural genes
To gain insight into the functional significance of these clusters of diversifying selective pressure detected within the HRV genome, we next examined how the location of the clusters of diversifying residues correlated with previously characterized functional and structural domains within the HRV genome. We first focused on the diversifying structural genes and examined the location of diversifying capsid residues relative to three previously characterized functional domains of the HRV virion: the neutralizing immunogen (NIm) sites, the cellular receptor contacts, and the binding pocket of pleconaril, a potent capsid inhibitor of HRVs and HEVs .
In contrast, analysis of the selective pressure in the capsid residues within the pleconaril binding site revealed an overall paucity of diversifying selective pressure (Additional file 3, Table S2). However, one of the residues lining the pleconaril binding site in the VP1 gene (residue #191) has diversifying selective pressure detectable above background. Intriguingly, this residue corresponds to one of two residues in the binding pocket shared among naturally occurring pleconaril resistant HRVB serotypes. When mutated in a susceptible HRVB serotype, residue #191 has been shown to confer a 30-fold reduction in pleconaril susceptibility .
Structure-function mapping of diversifying residues in non-structural genes
Given the essential nature of the functions performed by the products of the non-structural genes, it was quite surprising to detect a cluster of diversifying selective pressure within the 3C and 3D genes of the HRV genome. The wealth of structural and functional observations concerning these two factors allowed for analysis of the correlation in location of diversifying residues relative to the structural and functional domains previously characterized in each of these two non-structural genes.
A number of diversifying residues were also observed in regions of the 3D protein for which functional data is lacking. This is the case for a large set of diversifying residues found to localize to the outer surface of the fingers subdomain of the polymerase (Figures 8A and 8B). The role that this large domain plays in polymerase activity is not completely resolved. Recent work has demonstrated at least one residue in this domain (the highly conserved G64) can influence polymerase fidelity [57–60]. However, because this residue lies distant from the diversifying residues we detect on the surface of the fingers subdomain, their possible functional significance is unclear. Taken together, these data indicate, that like the 3C protease, proximity to characterized functional domains of the 3D polymerase does fully explain the diversifying pressure detected in this essential viral factor.
Conservation of non-coding RNAs and essential structural elements
Like all members of the Picornaviridae family, HRVs possess a number of essential cis-acting RNA elements that are required for, or enhance viral replication . An essential cloverleaf structure and internal ribosomal entry site (IRES) have been identified in the 5' non-coding region of the genome, while a small hairpin RNA element that enhances replication has been found in the 3' non-coding region. An additional essential RNA structure, a small stem-loop cis-acting replication element (CRE) resides within the coding sequences of the Picornaviridae genomes.
In our analysis of 34 HRV genome sequences, evidence for conservation of each of these elements was detected at both the primary sequence and secondary structure level (Additional file 4, Data S3 and S4). While these structures have been inferred previously from phylogenetic comparisons of available HRV genomes , our analysis provides a robust HRV consensus structure for each element in the 5' and 3' non-coding region (Additional file 4, Data S3 and S4).
Here, we have addressed a gap in our understanding of the evolutionary forces driving diversification of HRV and deepened our understanding of HRV biology in a number of ways. First, we have augmented the set of 6 fully sequenced HRV serotypes to a more representative subset of 34 genomes from across the HRV phylogeny. Second, we have performed a comprehensive analysis of the genetic diversity and evolutionary pressures operating upon the HRV genus. We have found a uniform pattern of genetic variability across the genome that is unlikely to be driven by large-scale recombination events as has been observed among other genera of the picornavirus family. We have also obtained a molecular portrait of the HRV genomic evolutionary landscape, which has revealed clusters of diversifying residues in both structural and non-structural genes cast against a background of purifying selective pressure. Finally, we have provided insight into the possible functional relevance of the detected diversifying pressure in both the structural and non-structural genes of HRV through comparison of the overlap in these residues with structural and functional domains previously characterized in HRV.
Correlation in genetic and phenotypic subgroupings of HRV
Our results indicate that the 2 major genetic subgroups of HRV correlate directly with phenotypic groupings based on in vitro studies of HRV susceptibility to a set of early generation "pocket factor" binding drugs that interact with the capsid gene products of the virus . This puzzling correlation between pocket factor susceptibility and the genetic relationships of non-structural genes in HRV was first noted almost 20 years ago in the original drug susceptibility study when only a limited set of non-structural gene sequences were available . More recent subgenomic sequence analyses have largely corroborated these findings [25, 26, 28]. Here, we extend these results to every locus of the HRV genome.
In general, this observation has been somewhat difficult to understand since these drugs could not have shaped HRV evolution, given that they have not been commonly used to treat viral infections in general, or HRV infections in particular. Our results provide a possible explanation. Because there is a consistent level of sequence diversity across the HRV genome, each locus in the genome possesses a genetic relationship identical to that of the structural genes targeted by the drug. Thus, the correlation between genotype and drug susceptibility phenotype is easily detectable at each loci in the genome, regardless of its potential to interact directly with the drug.
Recombination and diversification in the HRV genome
Our analysis has also revealed a lack of significant recombination within the HRV genome that is surprising in light of the fact that multiple serotypes that utilize the same cellular receptor are known to co-circulate during each HRV season . Moreover, this is also quite distinct from what has been observed for other genera in the Picornaviridae family, where recombination has been proposed to play a significant role in genetic diversification (reviewed in ). Taken together, our results favor the possibility that genetic drift is likely to be the major driving force for diversification in the HRV genus. These conclusions extend and agree with the recent work of Simmonds . It would appear that the known HRV isolates act as independently segregating genomes, with little potential for inter-genome recombination, in contrast to the non-segregating, highly recombinant genomes such as HEV, FMDV, the teschoviruses, and bovine enteroviruses.
Furthermore, it has been hypothesized that there is a biological compatibility barrier for recombination among HRV serotypes, since experimental evidence has demonstrated recombinants from similarly diverged picornaviruses tend to be inviable (reviewed in ). It is also possible that there may be additional barriers related to the characteristics of HRV infection (intracellular partitioning, persistence time in the cell, viral titer, blocks to co-infection, etc) that might preclude the opportunity for recombination to occur. With a diverse array of HRV genome sequences in hand, such hypotheses can now be directly tested.
Purifying selective pressure dominates in the HRV genome
Despite a notoriously error-prone polymerase and a significant amount of genetic diversity across the HRV genome, our selective pressure analysis indicates that overall, the HRV genome is under strong pressure to preserve the amino acid sequences encoded within genome. This sort of profile is not unique to HRV, since a similar bias towards purifying selection has been detected in selective pressure analysis of the capsid region of FMDV field isolates . A preponderance of purifying selective pressure is particularly obvious for the central region of the genome encoding the non-structural P2 gene products (P2A protease, P2B 'viroporin', and P2C ATPase and membrane association factor) and the 3A and 3B gene products. Each of these viral gene products is known to proteolyze or to interact with essential cellular factors, which are highly conserved. Thus, it may be that the lifecycle of HRV and its requirement to interact with and inactivate a variety of host factors results in significant sequence constraints within this portion of the genome.
Although these results may appear to contradict recent studies demonstrating that at least one Picornaviridae family member, poliovirus, evolves through quasispeciation , they actually do not rule out a similar process occurring in HRV. Rather, our results reflect the overall selective pressure acting on the HRV genome derived from the consensus sequences generated from our shotgun assemblies, and we have not focused on the potential minority polymorphisms that may exist within the population of each of the HRV prototypes. Inspection of each of our shotgun assemblies does reveal high quality sequence polymorphisms in a minority of the shotgun reads throughout the assembled genomes (data not shown). However, a greater depth of sequencing for each isolate would be required to unambiguously address the extent of HRV quasispeciation.
Implications of diversifying selective pressure in the structural genes
Although we detected overlap with each of the functional domains found on the viral particle, the diversifying capsid residues overlap significantly only with previously identified antigenic sites from both the HRVA and HRVB genomes. This result is intriguing in light of the variability in genetic diversity and serotype diversity known to exist in some of the Picornaviridae family members, such as the FMDVs and HEVs. The FMDVs are similar to HRVs, in that over 100 distinct serotypes have also been identified . These observations suggest that the icosahedral viral particle of these picornaviruses is relatively flexible, and is able to accommodate a wide array of non-synonymous changes. However, this immunogenic diversity is not generally shared among the capsids of all Picornaviridae family members. In particular, poliovirus has only 3 characterized serotypes. Moreover, recent analysis of vaccine-derived poliovirus isolates indicates that many of the most frequent non-synonymous changes which develop in the capsid genes do not alter the immunogenicity of the virus, despite being present in antigenic determinants . It is unclear if these results are unique to poliovirus or extend to other picornaviruses.
This is particularly relevant for our analysis, since we were unable to explain all of the diversifying selective pressure by direct overlap with antigenic sites on the surface of the viral pentamer. While many of our diversifying residues map within close proximity to these NIms, it is unclear if diversification of sites proximal to NIms actually alters their antigenicity. Such questions are difficult to resolve at this time, since the known antigenic determinants of HRV have been identified through sequence analysis of HRVs able to escape neutralization of a limited set of monoclonal antibodies raised against only 2 of the 102 HRV serotypes [46–52]. Thus, a more complete understanding of the statistically significant proximity detected here between diversifying capsid residues and the NIms awaits more comprehensive characterization of additional distinct antigenic sites on the HRV capsid.
Although not statistically significant, a surprising amount of overlap was also detected between the diversifying capsid residues and the characterized HRV cellular receptor contacts. Whether diversification of in these residues actually alters the functionality of these domains in the capsid, or merely reflects as-yet undiscovered functions, or regions of the HRV capsid that are under immune surveillance is unclear from these observations. However, it has been established that important functional domains in viruses are not excluded from immune surveillance, and that mutations within antigenic targets that overlap functional domains can abolish antibody interaction with little or no impact on interactions required in the functional domain (reviewed in ). Whether such observations also apply to this set of diversifying residues requires a more comprehensive understanding of both the antigenic determinants of the HRV capsid as well as the binding affinities to the HRV cellular receptors across different HRV serotypes.
Implications of diversifying selective pressure in the non-structural genes
Perhaps one of the most surprising results from this analysis was the detection of clusters of diversifying residues within two non-structural genes that perform essential functions during viral replication. Why did we detect any diversifying residues in these genes? We attempted to investigate this question through similar mapping of the location of the diversifying residues onto available crystal structures of the 3C protease and 3D polymerase. As was observed for the diversifying capsid residues, the diversifying residues in both the 3C protease and 3D polymerase map to surface-exposed residues; however, here we observed less of a bias towards a particular location or functional domain on the surface of each of these factors. We did detect a large proportion of the diversifying residues in the 3C protease and 3D polymerase positioned in the vicinity of characterized domains that are likely to influence RNA/VPg primer binding (for 3C protease) or hypothesized oligomerization domain interactions, protein binding and/or the coordination of subdomain movements that have been hypothesized to influence catalytic activity (for 3D polymerase).
However, the remaining fraction of the diversifying residues within these non-structural genes map to regions in each of these factors for which functions have not yet been assigned. We have not detected a correlation between the 3C protease and 3D polymerase diversifying residues with MHC class I presenting peptides detectable in 3C and 3D. Likewise, we were also unable to detect any correlation between variation in electrostatic potential on the surface of the 3C protease and 3D polymerase, or significant covariation with any other diversifying residues in the genome. Thus, the role these diversifying residues may play in specific functions of the 3C protease and 3D polymerase, or in overall viral fitness, requires further exploration.
Such studies are particularly relevant given recent discoveries highlighting our incomplete knowledge of the functional domains within these two factors. Recently, a previously uncharacterized region of the poliovirus 3D polymerase lying outside the catalytic domain was shown to influence polymerase activity and thus fidelity [58, 59, 68]. Similarly, mutational analysis of the poliovirus 3C protein has recently uncovered a number of residues required for viral replication and VPg binding that happen localize outside the defined protease and RNA binding/VPg primer binding domains but in proximity with these unassigned diversifying residues, (C.E. Cameron, personal communication). Additional progress in structural analysis of the poliovirus 3CD precursor also indicates potential intersubunit (3C–3D) and intrasubunit (3D–3D) interactions in domains of the 3C and 3D subunits within close proximity to a number of the diversifying residues we have identified within regions of currently unassigned function . A complete understanding of the possible functional role that these diversifying residues may play in either of these individual factors or the active 3CD precursor awaits additional functional studies. The convergence of our results with these independent studies suggesting novel functional domains and interactions within the non-structural genes points to the utility of selective pressure analysis to uncover potentially important functional domains within a genome that may influence viability and overall fitness.
Conservation of essential non-coding RNA elements in the HRV genome
Analysis of RNA elements present in both the non-coding (5' cloverleaf and IRES, and 3' stem-loop element) and coding regions (CRE) of the HRV genome indicates conservation of both sequence and secondary structures in these regulatory elements in both HRVA and HRVB genomes. Although the consensus secondary structures among these elements appear similar to those generated based on a much smaller set of HRV genome sequences , subtle sequence variations can be detected between the HRVA and HRVB subgroup members, as well as within each of the subgroup members (Additional file 4, Data S3 and S4). Such differences are of particular interest as these elements have been shown to be essential for viral replication, translation, overall viability, and in the case of poliovirus, for pathogenicity and tissue tropism [72–75]. Comprehensive analyses of the functional implications and associated clinical implications of diversity in sequence and secondary structure of these regions of the HRV genome have not been performed. Correlations in variation of the known functions of these RNAs with the sequence variation and structural diversity found within this subset of HRVs will shed light on the role they play in viral growth and replication, and may further clarify the role non-coding regions in HRV pathogenesis.
Potential role for selective pressure analysis in drug development
To date, two drugs targeting conserved regions of the HRV genome have advanced to Phase III clinical trials. Pleconaril, a potent capsid inhibitor of HRVs and HEVs, binds to a surface-accessible hydrophobic pocket in the VP1 protein on the external face of the viral particle . Ruprintrivir targets the proteolytic active site of the 3C protein and exhibits broad inhibition of HRV growth in vitro .
Unfortunately, neither of these drugs has demonstrated sufficient symptom relief, or in the case of pleconaril, exhibited untoward interactions with other drugs. Thus, FDA approval was not granted for either of these potential therapies. Moreover, pleconaril treatment has been shown to give rise to drug resistant viruses at a low frequency . This has not been observed with rupritrivir. Such observations can be explained in the context of our selective pressure analysis. Inspection of our data for the residues targeted by these two drugs reveals only a single residue to possess diversifying selective pressure above background (Additional File 3, Tables S2 and Additional file 4, Table S3). This residue lies within the pleconaril binding site and corresponds to VP1 residue 191. Prior work identified this residue to be one of two residues that varied from the consensus valine in pleconaril susceptible HRV serotypes to leucine in resistant HRV serotypes . In fact, a V191L mutation engineered in a susceptible HRVB serotype was found to be sufficient to confer a 30-fold reduction in susceptibility to pleconaril .
Having identified the only residue known to yield pleconaril resistance, these results illustrate the potential utility of selective pressure analysis with respect to drug development. In early stages of drug development, selective pressure analysis combined with assays for drug efficacy and viral pathogenicity could prove valuable in de novo choice of drug targets. The diversifying potential of residues within or flanking drug binding sites could be evaluated in silico, and mutations in such residues could be engineered and assayed for drug binding, normal substrate binding, and viral growth. Ultimately, incorporating such analysis in the drug development pipeline may allow the avoidance of targets with high potential for drug resistance or increased virulence.
This analysis has closed a gap in our understanding of the genetic diversity and evolutionary pressures across the HRV genome. It has provided a deeper understanding of the similarities and differences between the genetic diversity present in HRV compared to other genera of the picornavirus family. These results have also raised several testable questions related to several domains of unknown function and HRV evolution itself. Ultimately, such knowledge may serve to elucidate the determinants of pathogenicity within the HRV genome and aid in the development of therapeutics to reduce or eliminate the clinical symptoms associated with this ubiquitous respiratory pathogen.
Isolation of RNA from low passage HRV prototype stocks
Low passage tissue culture supernatants from tissue culture cells infected with the HRV serotypes (indicated in Additional File 1, Figure S1 and Additional File 5, Table S4) were obtained from the California Department of Health Services (CaDHS). Supernatants were centrifuged briefly to pellet cellular debris, then passed through 0.2 μm filters, brought to 10 mM CaCl2, and incubated with 600 units of micrococcal nuclease (Fermentas) for 3 hours at 37°C. RNA was then isolated from the culture supernatants via Trizol:chloroform extraction, followed by isopropanol precipitation.
Amplification and shotgun sequencing of HRV prototype stock RNA
RNA isolated from HRV prototype culture supernatants was reverse transcribed, randomly amplified as previously described , and cloned into the pCR2.1 TOPO TA vector (Invitrogen) to generate plasmid libraries for each HRV serotype. The resulting libraries were transformed into bacteria. Plasmid DNA prepared from each library of transformants was sequenced using the Big Dye terminator v. 3.1 (Applied Biosystems) containing either -21 universal or -28 reverse primer and analyzed on an ABI 3730xl sequencer (Applied Biosystems).
Shotgun sequence analysis and assembly of HRV genomes
Approximately 7 Mb of DNA derived from 14,208 reads, with an average length of 500 bp, were shotgun sequenced to generate the initial HRV genome assemblies. Contaminating human and bacterial reads (60% of all reads) were identified and removed through BLAST analysis . A total of 8,278 viral reads were processed and assembled with the CONSED software suite . Overall, each genome assembly contained an average of 304 input viral reads, with an average read depth of 22, and average quality score of 86.4 (Additional file 5, Table S5). Specific PCR was performed to obtain sequences at the extreme 5'end and 3'end of each genome sequenced and to close any internal gaps. For the ends, a single high quality (minimum phred score of 20) sequencing read with at least 100 nucleotides of overlap with the shotgun assembly reads was required to consider each genome finished. For the internal gaps, a minimum of 2 high quality forward and reverse reads with overlap of at least 100 nucleotides with shotgun contigs were required to consider internal gaps closed. A shotgun sequence assembly derived from the previously sequenced HRV001b  was used to validate the quality of sequences obtained by these methods. The resulting shotgun assembly of HRV001b was 99.6% identical (6198 identities of 6223 nucleotides assembled) to the fully sequenced HRV001b present in NCBI (genbank identifier 221708).
Sequence alignment and phylogenetic analysis
Inferred amino acid sequence of the coding regions of the 34 complete HRV genomes were aligned using the CLUSTALW program . This alignment was then back-translated into nucleotide sequence and combined with alignments of the 5' and 3' non-coding regions, generated using CLUSTALW, to form the whole-genome nucleotide alignment used for analysis. Neighbor-joining phylogenetic trees were generated from the alignment using CLUSTALW with Kimura's correction for multiple base substitutions. Maximum likelihood trees were generated using baseml from the PAML  package and DNAML from the Phylip  package. Trees generated using neighbor-joining and maximum likelihood methods contained similar topologies, and differed only in computed branch lengths. The HKY85 model of nucleotide substitution was used, and the values of the transition/transversion rate and the alpha parameter in baseml were estimated through maximum likelihood calculation. Alignment positions with gaps were ignored in all cases.
Scanning average pairwise sequence identity plots were generated using a moving window of 100 nucleotides or 50 amino acids across the whole-genome nucleotide alignment and the corresponding amino acid translation in the coding region of the genome.
The genomic nucleotide alignment of the 34 complete HRV genomes was analyzed using RDP version 2 . Six automated recombination analysis algorithms were run: RDP, GENECONV , BOOTSCAN , MaxChi , Chimaera , and Sister Scanning . These algorithms were selected from the set of published recombination detection methods based on their ability to identify recombinant sequences, the associated breakpoints, and parental sequences. In computational and empirical comparative tests, no single method performed best under all conditions, and consistent results from more than one method was the best indicator of recombination [87, 89]. Resulting predictions of recombination events with p-values less than 0.05 were analyzed manually using all six methods. Events supported by evidence from more than one method were further characterized by manual analysis of bootstrapped phylogenetic trees of the relevant genomic locus to determine the genotypes involved in the recombination event.
Selective pressure analysis
Codon-based models of evolution of coding sequence allowing for variable selection pressure among sites in a maximum-likelihood framework were used to evaluate the selective pressure operating on each gene individually. Codon-substitution models [90, 91] were compared using likelihood ratio tests (LRT) to test for significant diversifying selection within each gene.
These codon-substitution models, allowing for variable ω (dN/dS) parameters among sites, were fit to the nucleotide alignment of the coding sequence of the genome. Model M1a, or the neutral model, incorporates a class of sites under purifying selection with ω0 < 1, and a second class of sites with ω1 = 1. Model M2a adds a third class of sites ω2 > 1, to allow for diversifying selection. Similarily, Model 7 incorporates a discrete beta distribution (10 classes) to model values of ω between 0 and 1, while Model 8 adds an additional parameter ω > 1. Likelihood ratio tests were performed between nested models (M1a versus M2a, or M7 versus M8) to calculate the significance of diversifying selection within a gene (Additional file 5, Table S6). An empirical Bayesian approach was then used to calculate the posterior probability that a site belongs to each of the ω site classes. This probability value was then used to compute an estimate of dN/dS for each site in the sequence. Maximum likelihood calculations on the substitution models were implemented using the codeml program from version 3.14 of the PAML package .
To ascertain how well the resulting dN/dS values computed from the subset of 34 reference genomes reflected the selective pressure present in the full set of 102 known HRV serotypes, we compared the dN/dS values computed for each residue in the VP1 gene of this set of HRVA and HRVB serotypes to the same dN/dS values obtained independently from the available VP1 sequences of all 102 HRV serotypes [25, 26]. Although the absolute value of the dN/dS ratios differed between the two sets, their relative rankings were well correlated (0.91 and 0.80, for HRVA and HRVB genomes, respectively; Additional file 1, Figure S2), with few potential false positives and false negatives detected. Thus, it appears that the relative rank, rather than absolute magnitude of the dN/dS values we have computed from this subset of HRV genomes accurately approximates the selective pressures detectable among the full set of 102 HRV reference serotypes.
Tests of heterogeneous synonymous substitution rates among sites were performed using the REL analysis implemented in the HYPHY  phylogenetic package. This method of analysis is very similar to that described above, but differs in codon models available, and in the modeling of site classes (REL site classes are modeled as N discrete classes, similar to model M3 in codonml). Analysis using the GY  model of codon evolution with six discrete classes of non-synonymous and synonymous mutation rates was used to determine the effects of variable dS across sites on the data. Although varying dS resulted in a lowered magnitude of a number of capsid residues in the smaller dataset of HRVB genomes, it did not significantly impact the per-residue dN/dS values for the HRVA genomes or confer any significant changes in the overall identity or localization of the 5% highest scoring dN/dS residues of the capsid genes (Additional file 5, Figure S9). Thus, for the sake of simplicity, dN/dS values discussed in the results section were those derived from the calculations described above assuming a homogeneous synonymous substitution rate.
Mapping dN/dS values onto 3-dimensional crystal structures
Viral pentamer structures were generated from the NCBI Protein Database (pdb) files of HRV2 (pdb id 1FPN), HRV14 (pdb id 4RHV), and HRV16 (pdb id 1AYM) using the Oligomer Generator utility from the VIPERdb website . Analysis of the 3C protease and 3D polymerase was performed using the HRV2 3C protease (pdb id 1CQQ), and HRV14 3D polymerase (pdb id 1XR5), respectively. The molecular structure visualization program, Chimera , was used to generate images of the viral proteins.
Calculations of the significance of the overlap in structure space between sets of dN/dS data were calculated using an average minimum distance between residues metric. Observed average minimum distance between two sets (A and B) of residues was calculated by taking the average of the minimum three-dimensional Cartesian distance from each residue of set A to the nearest residue from set B. In effect this is a measurement of how closely correlated the positions of set A are to any subset of the positions in set B. To calculate the significance of this observed distance, 100,000 iterations of this calculation were computed, randomizing the locations of the residues in set A for each calculation. The distribution of the resulting average minimum distance values was used to calculate a p-value for the significance of the observed value.
The GenBank accession numbers for the sequenced HRV genomes range from DQ473485–DQ473512.
This work was supported by a grant from the Sandler Program for Asthma Research, the Packard Foundation, the Doris Duke Charitable Foundation, the Howard Hughes Medical Institute, and a grant from the National Institutes of Health grant R21 AI057506. We are grateful to Lisa Cook, Donald Williams, Jim Eldred, and Matthew Hickenbotham for providing technical support with shotgun sequencing, and to D. Ganem, C. Chiu, K. Fischer, P. Tang, A. Urisman, and R. Andino for advice and comments.
- Heikkinen T, Jarvinen A: The common cold. Lancet 2003,361(9351):51-59. 10.1016/S0140-6736(03)12162-9PubMedView ArticleGoogle Scholar
- Fendrick AM, Monto AS, Nightengale B, Sarnes M: The economic burden of non-influenza-related viral respiratory tract infection in the United States. Arch Intern Med 2003,163(4):487-494. 10.1001/archinte.163.4.487PubMedView ArticleGoogle Scholar
- Gern JE, Galagan DM, Jarjour NN, Dick EC, Busse WW: Detection of rhinovirus RNA in lower airway cells during experimentally induced infection. Am J Respir Crit Care Med 1997,155(3):1159-1161.PubMedView ArticleGoogle Scholar
- Papadopoulos NG, Bates PJ, Bardin PG, Papi A, Leir SH, Fraenkel DJ, Meyer J, Lackie PM, Sanderson G, Holgate ST, et al.: Rhinoviruses infect the lower airways. J Infect Dis 2000,181(6):1875-1884. 10.1086/315513PubMedView ArticleGoogle Scholar
- Papadopoulos NG, Johnston SL: Rhinoviruses as pathogens of the lower respiratory tract. Can Respir J 2000,7(5):409-414.PubMedGoogle Scholar
- Papadopoulos NG, Sanderson G, Hunter J, Johnston SL: Rhinoviruses replicate effectively at lower airway temperatures. J Med Virol 1999,58(1):100-104. Publisher Full Text 10.1002/(SICI)1096-9071(199905)58:1<100::AID-JMV16>3.0.CO;2-DPubMedView ArticleGoogle Scholar
- Schroth MK, Grimm E, Frindt P, Galagan DM, Konno SI, Love R, Gern JE: Rhinovirus replication causes RANTES production in primary bronchial epithelial cells. Am J Respir Cell Mol Biol 1999,20(6):1220-1228.PubMedView ArticleGoogle Scholar
- Subauste MC, Jacoby DB, Richards SM, Proud D: Infection of a human respiratory epithelial cell line with rhinovirus. Induction of cytokine release and modulation of susceptibility to infection by cytokine exposure. J Clin Invest 1995,96(1):549-557.PubMedPubMed CentralView ArticleGoogle Scholar
- Hayden FG: Rhinovirus and the lower respiratory tract. Rev Med Virol 2004,14(1):17-31. 10.1002/rmv.406PubMedView ArticleGoogle Scholar
- Ghosh S, Champlin R, Couch R, Englund J, Raad I, Malik S, Luna M, Whimbey E: Rhinovirus infections in myelosuppressed adult blood and marrow transplant recipients. Clin Infect Dis 1999,29(3):528-532.PubMedView ArticleGoogle Scholar
- Ison MG, Hayden FG, Kaiser L, Corey L, Boeckh M: Rhinovirus infections in hematopoietic stem cell transplant recipients with pneumonia. Clin Infect Dis 2003,36(9):1139-1143. 10.1086/374340PubMedView ArticleGoogle Scholar
- Garbino J, Gerbase MW, Wunderli W, Deffernez C, Thomas Y, Rochat T, Ninet B, Schrenzel J, Yerly S, Perrin L, et al.: Lower respiratory viral illnesses: improved diagnosis by molecular methods and clinical impact. Am J Respir Crit Care Med 2004,170(11):1197-1203. 10.1164/rccm.200406-781OCPubMedView ArticleGoogle Scholar
- Hamparian VV, Colonno RJ, Cooney MK, Dick EC, Gwaltney JM Jr, Hughes JH, Jordan WS Jr, Kapikian AZ, Mogabgab WJ, Monto A, et al.: A collaborative report: rhinoviruses – extension of the numbering system from 89 to 100. Virology 1987,159(1):191-192. 10.1016/0042-6822(87)90367-9PubMedView ArticleGoogle Scholar
- Andries K, Dewindt B, Snoeks J, Wouters L, Moereels H, Lewi PJ, Janssen PA: Two groups of rhinoviruses revealed by a panel of antiviral compounds present sequence divergence and differential pathogenicity. J Virol 1990,64(3):1117-1123.PubMedPubMed CentralGoogle Scholar
- Blomqvist S, Savolainen C, Raman L, Roivainen M, Hovi T: Human rhinovirus 87 and enterovirus 68 represent a unique serotype with rhinovirus and enterovirus features. J Clin Microbiol 2002,40(11):4218-4223. 10.1128/JCM.40.11.4218-4223.2002PubMedPubMed CentralView ArticleGoogle Scholar
- Oberste MS, Maher K, Schnurr D, Flemister MR, Lovchik JC, Peters H, Sessions W, Kirk C, Chatterjee N, Fuller S, et al.: Enterovirus 68 is associated with respiratory illness and shares biological features with both the enteroviruses and the rhinoviruses. J Gen Virol 2004,85(Pt 9):2577-2584. 10.1099/vir.0.79925-0PubMedView ArticleGoogle Scholar
- Abraham G, Colonno RJ: Many rhinovirus serotypes share the same cellular receptor. J Virol 1984,51(2):340-345.PubMedPubMed CentralGoogle Scholar
- Uncapher CR, DeWitt CM, Colonno RJ: The major and minor group receptor families contain all but one human rhinovirus serotype. Virology 1991,180(2):814-817. 10.1016/0042-6822(91)90098-VPubMedView ArticleGoogle Scholar
- Greve JM, Davis G, Meyer AM, Forte CP, Yost SC, Marlor CW, Kamarck ME, McClelland A: The major human rhinovirus receptor is ICAM-1. Cell 1989,56(5):839-847. 10.1016/0092-8674(89)90688-0PubMedView ArticleGoogle Scholar
- Staunton DE, Merluzzi VJ, Rothlein R, Barton R, Marlin SD, Springer TA: A cell adhesion molecule, ICAM-1, is the major surface receptor for rhinoviruses. Cell 1989,56(5):849-853. 10.1016/0092-8674(89)90689-2PubMedView ArticleGoogle Scholar
- Hofer F, Gruenberger M, Kowalski H, Machat H, Huettinger M, Kuechler E, Blass D: Members of the low density lipoprotein receptor family mediate cell entry of a minor-group common cold virus. Proc Natl Acad Sci USA 1994,91(5):1839-1842. 10.1073/pnas.91.5.1839PubMedPubMed CentralView ArticleGoogle Scholar
- Binford SL, Maldonado F, Brothers MA, Weady PT, Zalman LS, Meador JW 3rd, Matthews DA, Patick AK: Conservation of amino acids in human rhinovirus 3C protease correlates with broad-spectrum antiviral activity of rupintrivir, a novel human rhinovirus 3C protease inhibitor. Antimicrob Agents Chemother 2005,49(2):619-626. 10.1128/AAC.49.2.619-626.2005PubMedPubMed CentralView ArticleGoogle Scholar
- Deffernez C, Wunderli W, Thomas Y, Yerly S, Perrin L, Kaiser L: Amplicon sequencing and improved detection of human rhinovirus in respiratory samples. J Clin Microbiol 2004,42(7):3212-3218. 10.1128/JCM.42.7.3212-3218.2004PubMedPubMed CentralView ArticleGoogle Scholar
- Horsnell C, Gama RE, Hughes PJ, Stanway G: Molecular relationships between 21 human rhinovirus serotypes. J Gen Virol 1995,76(Pt 10):2549-2555.PubMedView ArticleGoogle Scholar
- Laine P, Savolainen C, Blomqvist S, Hovi T: Phylogenetic analysis of human rhinovirus capsid protein VP1 and 2A protease coding sequences confirms shared genus-like relationships with human enteroviruses. J Gen Virol 2005,86(Pt 3):697-706. 10.1099/vir.0.80445-0PubMedView ArticleGoogle Scholar
- Ledford RM, Patel NR, Demenczuk TM, Watanyar A, Herbertz T, Collett MS, Pevear DC: VP1 sequencing of all human rhinovirus serotypes: insights into genus phylogeny and susceptibility to antiviral capsid-binding compounds. J Virol 2004,78(7):3663-3674. 10.1128/JVI.78.7.3663-3674.2004PubMedPubMed CentralView ArticleGoogle Scholar
- Loens K, Ieven M, Ursi D, De Laat C, Sillekens P, Oudshoorn P, Goossens H: Improved detection of rhinoviruses by nucleic acid sequence-based amplification after nucleotide sequence determination of the 5' noncoding regions of additional rhinovirus strains. J Clin Microbiol 2003,41(5):1971-1976. 10.1128/JCM.41.5.1971-1976.2003PubMedPubMed CentralView ArticleGoogle Scholar
- Savolainen C, Blomqvist S, Mulders MN, Hovi T: Genetic clustering of all 102 human rhinovirus prototype strains: serotype 87 is close to human enterovirus 70. J Gen Virol 2002,83(Pt 2):333-340.PubMedView ArticleGoogle Scholar
- Savolainen C, Laine P, Mulders MN, Hovi T: Sequence analysis of human rhinoviruses in the RNA-dependent RNA polymerase coding region reveals large within-species variation. J Gen Virol 2004,85(Pt 8):2271-2277. 10.1099/vir.0.79897-0PubMedView ArticleGoogle Scholar
- Skern T, Sommergruber W, Blaas D, Gruendler P, Fraundorfer F, Pieler C, Fogy I, Kuechler E: Human rhinovirus 2: complete nucleotide sequence and proteolytic processing signals in the capsid protein region. Nucleic Acids Res 1985,13(6):2111-2126. 10.1093/nar/13.6.2111PubMedPubMed CentralView ArticleGoogle Scholar
- Lee WM, Wang W, Rueckert RR: Complete sequence of the RNA genome of human rhinovirus 16, a clinically useful common cold virus belonging to the ICAM-1 receptor group. Virus Genes 1995,9(2):177-181. 10.1007/BF01702661PubMedView ArticleGoogle Scholar
- Hughes PJ, North C, Jellis CH, Minor PD, Stanway G: The nucleotide sequence of human rhinovirus 1B: molecular relationships within the rhinovirus genus. J Gen Virol 1988,69(Pt 1):49-58.PubMedView ArticleGoogle Scholar
- Stanway G, Hughes PJ, Mountford RC, Minor PD, Almond JW: The complete nucleotide sequence of a common cold virus: human rhinovirus 14. Nucleic Acids Res 1984,12(20):7859-7875. 10.1093/nar/12.20.7859PubMedPubMed CentralView ArticleGoogle Scholar
- Callahan PL, Mizutani S, Colonno RJ: Molecular cloning and complete sequence determination of RNA genome of human rhinovirus type 14. Proc Natl Acad Sci USA 1985,82(3):732-736. 10.1073/pnas.82.3.732PubMedPubMed CentralView ArticleGoogle Scholar
- Duechler M, Skern T, Sommergruber W, Neubauer C, Gruendler P, Fogy I, Blaas D, Kuechler E: Evolutionary relationships within the human rhinovirus genus: comparison of serotypes 89, 2, and 14. Proc Natl Acad Sci USA 1987,84(9):2605-2609. 10.1073/pnas.84.9.2605PubMedPubMed CentralView ArticleGoogle Scholar
- Harris JR, Racaniello VR: Amino acid changes in proteins 2B and 3A mediate rhinovirus type 39 growth in mouse cells. J Virol 2005,79(9):5363-5373. 10.1128/JVI.79.9.5363-5373.2005PubMedPubMed CentralView ArticleGoogle Scholar
- Louie JK, Yagi S, Nelson FA, Kiang D, Glaser CA, Rosenberg J, Cahill CK, Schnurr DP: Rhinovirus outbreak in a long term care facility for elderly persons associated with unusually high mortality. Clin Infect Dis 2005,41(2):262-265. 10.1086/430915PubMedView ArticleGoogle Scholar
- Carrillo C, Tulman ER, Delhon G, Lu Z, Carreno A, Vagnozzi A, Kutish GF, Rock DL: Comparative genomics of foot-and-mouth disease virus. J Virol 2005,79(10):6487-6504. 10.1128/JVI.79.10.6487-6504.2005PubMedPubMed CentralView ArticleGoogle Scholar
- Lukashev AN: Role of recombination in evolution of enteroviruses. Rev Med Virol 2005,15(3):157-167. 10.1002/rmv.457PubMedView ArticleGoogle Scholar
- Brown B, Oberste MS, Maher K, Pallansch MA: Complete genomic sequencing shows that polioviruses and members of human enterovirus species C are closely related in the noncapsid coding region. J Virol 2003,77(16):8973-8984. 10.1128/JVI.77.16.8973-8984.2003PubMedPubMed CentralView ArticleGoogle Scholar
- Oberste MS, Maher K, Pallansch MA: Evidence for frequent recombination within species human enterovirus B based on complete genomic sequences of all thirty-seven serotypes. J Virol 2004,78(2):855-867. 10.1128/JVI.78.2.855-867.2004PubMedPubMed CentralView ArticleGoogle Scholar
- Oberste MS, Penaranda S, Maher K, Pallansch MA: Complete genome sequences of all members of the species Human enterovirus A. J Gen Virol 2004,85(Pt 6):1597-1607. 10.1099/vir.0.79789-0PubMedView ArticleGoogle Scholar
- Simmonds P: Recombination and selection in the evolution of picornaviruses and other Mammalian positive-stranded RNA viruses. J Virol 2006,80(22):11124-11140. 10.1128/JVI.01076-06PubMedPubMed CentralView ArticleGoogle Scholar
- Oberste MS, Penaranda S, Pallansch MA: RNA recombination plays a major role in genomic change during circulation of coxsackie B viruses. J Virol 2004,78(6):2948-2955. 10.1128/JVI.78.6.2948-2955.2004PubMedPubMed CentralView ArticleGoogle Scholar
- Magden J, Kaariainen L, Ahola T: Inhibitors of virus replication: recent developments and prospects. Appl Microbiol Biotechnol 2005,66(6):612-621. 10.1007/s00253-004-1783-3PubMedView ArticleGoogle Scholar
- Appleyard G, Russell SM, Clarke BE, Speller SA, Trowbridge M, Vadolas J: Neutralization epitopes of human rhinovirus type 2. J Gen Virol 1990,71(Pt 6):1275-1282.PubMedView ArticleGoogle Scholar
- Hastings GZ, Speller SA, Francis MJ: Neutralizing antibodies to human rhinovirus produced in laboratory animals and humans that recognize a linear sequence from VP2. J Gen Virol 1990,71(Pt 12):3055-3059.PubMedView ArticleGoogle Scholar
- Hewat EA, Blaas D: Structure of a neutralizing antibody bound bivalently to human rhinovirus 2. Embo J 1996,15(7):1515-1523.PubMedPubMed CentralGoogle Scholar
- Hewat EA, Marlovits TC, Blaas D: Structure of a neutralizing antibody bound monovalently to human rhinovirus 2. J Virol 1998,72(5):4396-4402.PubMedPubMed CentralGoogle Scholar
- Speller SA, Sangar DV, Clarke BE, Rowlands DJ: The nature and spatial distribution of amino acid substitutions conferring resistance to neutralizing monoclonal antibodies in human rhinovirus type 2. J Gen Virol 1993,74(Pt 2):193-200.PubMedView ArticleGoogle Scholar
- Sherry B, Mosser AG, Colonno RJ, Rueckert RR: Use of monoclonal antibodies to identify four neutralization immunogens on a common cold picornavirus, human rhinovirus 14. J Virol 1986,57(1):246-257.PubMedPubMed CentralGoogle Scholar
- Sherry B, Rueckert R: Evidence for at least two dominant neutralization antigens on human rhinovirus 14. J Virol 1985,53(1):137-143.PubMedPubMed CentralGoogle Scholar
- Bella J, Rossmann MG: Review: rhinoviruses and their ICAM receptors. J Struct Biol 1999,128(1):69-74. 10.1006/jsbi.1999.4143PubMedView ArticleGoogle Scholar
- Hewat EA, Neumann E, Conway JF, Moser R, Ronacher B, Marlovits TC, Blaas D: The cellular receptor to human rhinovirus 2 binds around the 5-fold axis and not in the canyon: a structural view. Embo J 2000,19(23):6317-6325. 10.1093/emboj/19.23.6317PubMedPubMed CentralView ArticleGoogle Scholar
- Ledford RM, Collett MS, Pevear DC: Insights into the genetic basis for natural phenotypic resistance of human rhinoviruses to pleconaril. Antiviral Res 2005,68(3):135-138. 10.1016/j.antiviral.2005.08.003PubMedView ArticleGoogle Scholar
- Pathak HB, Ghosh SK, Roberts AW, Sharma SD, Yoder JD, Arnold JJ, Gohara DW, Barton DJ, Paul AV, Cameron CE: Structure-function relationships of the RNA-dependent RNA polymerase from poliovirus (3Dpol). A surface of the primary oligomerization domain functions in capsid precursor processing and VPg uridylylation. J Biol Chem 2002,277(35):31551-31562. 10.1074/jbc.M204408200PubMedView ArticleGoogle Scholar
- Arnold JJ, Vignuzzi M, Stone JK, Andino R, Cameron CE: Remote site control of an active site fidelity checkpoint in a viral RNA-dependent RNA polymerase. J Biol Chem 2005,280(27):25706-25716. 10.1074/jbc.M503444200PubMedPubMed CentralView ArticleGoogle Scholar
- Pfeiffer JK, Kirkegaard K: A single mutation in poliovirus RNA-dependent RNA polymerase confers resistance to mutagenic nucleotide analogs via increased fidelity. Proc Natl Acad Sci USA 2003,100(12):7289-7294. 10.1073/pnas.1232294100PubMedPubMed CentralView ArticleGoogle Scholar
- Pfeiffer JK, Kirkegaard K: Increased Fidelity Reduces Poliovirus Fitness and Virulence under Selective Pressure in Mice. PLoS Pathog 2005,1(2):e11. 10.1371/journal.ppat.0010011PubMedPubMed CentralView ArticleGoogle Scholar
- Vignuzzi M, Stone JK, Andino R: Ribavirin and lethal mutagenesis of poliovirus: molecular mechanisms, resistance and biological implications. Virus Res 2005,107(2):173-181. 10.1016/j.virusres.2004.11.007PubMedView ArticleGoogle Scholar
- Witwer C, Rauscher S, Hofacker IL, Stadler PF: Conserved RNA secondary structures in Picornaviridae genomes. Nucleic Acids Res 2001,29(24):5079-5089. 10.1093/nar/29.24.5079PubMedPubMed CentralView ArticleGoogle Scholar
- Yang Y, Rijnbrand R, McKnight KL, Wimmer E, Paul A, Martin A, Lemon SM: Sequence requirements for viral RNA replication and VPg uridylylation directed by the internal cis-acting replication element (cre) of human rhinovirus type 14. J Virol 2002,76(15):7485-7494. 10.1128/JVI.76.15.7485-7494.2002PubMedPubMed CentralView ArticleGoogle Scholar
- Gerber K, Wimmer E, Paul AV: Biochemical and genetic studies of the initiation of human rhinovirus 2 RNA replication: identification of a cis-replicating element in the coding sequence of 2A(pro). J Virol 2001,75(22):10979-10990. 10.1128/JVI.75.22.10979-10990.2001PubMedPubMed CentralView ArticleGoogle Scholar
- McKnight KL, Lemon SM: The rhinovirus type 14 genome contains an internally located RNA structure that is required for viral replication. Rna 1998,4(12):1569-1584. 10.1017/S1355838298981006PubMedPubMed CentralView ArticleGoogle Scholar
- McKnight KL: The human rhinovirus internal cis-acting replication element (cre) exhibits disparate properties among serotypes. Arch Virol 2003,148(12):2397-2418. 10.1007/s00705-003-0177-7PubMedView ArticleGoogle Scholar
- Savolainen C, Mulders MN, Hovi T: Phylogenetic analysis of rhinovirus isolates collected during successive epidemic seasons. Virus Res 2002,85(1):41-46. 10.1016/S0168-1702(02)00016-3PubMedView ArticleGoogle Scholar
- Haydon DT, Bastos AD, Knowles NJ, Samuel AR: Evidence for positive selection in foot-and-mouth disease virus capsid genes from field isolates. Genetics 2001,157(1):7-15.PubMedPubMed CentralGoogle Scholar
- Vignuzzi M, Stone JK, Arnold JJ, Cameron CE, Andino R: Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature 2006,439(7074):344-348. 10.1038/nature04388PubMedPubMed CentralView ArticleGoogle Scholar
- Yakovenko ML, Cherkasova EA, Rezapkin GV, Ivanova OE, Ivanov AP, Eremeeva TP, Baykova OY, Chumakov KM, Agol VI: Antigenic evolution of vaccine-derived polioviruses: changes in individual epitopes and relative stability of the overall immunological properties. J Virol 2006,80(6):2641-2653. 10.1128/JVI.80.6.2641-2653.2006PubMedPubMed CentralView ArticleGoogle Scholar
- Colman PM: Virus versus antibody. Structure 1997,5(5):591-593. 10.1016/S0969-2126(97)00214-1PubMedView ArticleGoogle Scholar
- Marcotte LL, Wass AB, Gohara DW, Pathak HB, Arnold JJ, Filman DJ, Cameron CE, Hogle JM: Crystal structure of poliovirus 3CD protein: virally encoded protease and precursor to the RNA-dependent RNA polymerase. J Virol 2007,81(7):3583-3596. 10.1128/JVI.02306-06PubMedPubMed CentralView ArticleGoogle Scholar
- Evans DM, Dunn G, Minor PD, Schild GC, Cann AJ, Stanway G, Almond JW, Currey K, Maizel JV Jr: Increased neurovirulence associated with a single nucleotide change in a noncoding region of the Sabin type 3 poliovaccine genome. Nature 1985,314(6011):548-550. 10.1038/314548a0PubMedView ArticleGoogle Scholar
- Kawamura N, Kohara M, Abe S, Komatsu T, Tago K, Arita M, Nomoto A: Determinants in the 5' noncoding region of poliovirus Sabin 1 RNA that influence the attenuation phenotype. J Virol 1989,63(3):1302-1309.PubMedPubMed CentralGoogle Scholar
- Minor PD, Macadam AJ, Stone DM, Almond JW: Genetic basis of attenuation of the Sabin oral poliovirus vaccines. Biologicals 1993,21(4):357-363. 10.1006/biol.1993.1096PubMedView ArticleGoogle Scholar
- Ren RB, Moss EG, Racaniello VR: Identification of two determinants that attenuate vaccine-related type 2 poliovirus. J Virol 1991,65(3):1377-1382.PubMedPubMed CentralGoogle Scholar
- Pevear DC, Hayden FG, Demenczuk TM, Barone LR, McKinlay MA, Collett MS: Relationship of pleconaril susceptibility and clinical outcomes in treatment of common colds caused by rhinoviruses. Antimicrob Agents Chemother 2005,49(11):4492-4499. 10.1128/AAC.49.11.4492-4499.2005PubMedPubMed CentralView ArticleGoogle Scholar
- Wang D, Urisman A, Liu YT, Springer M, Ksiazek TG, Erdman DD, Mardis ER, Hickenbotham M, Magrini V, Eldred J, et al.: Viral discovery and sequence recovery using DNA microarrays. PLoS Biol 2003,1(2):E2. 10.1371/journal.pbio.0000002PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990,215(3):403-410.PubMedView ArticleGoogle Scholar
- Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res 1998,8(3):195-202.PubMedView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994,22(22):4673-4680. 10.1093/nar/22.22.4673PubMedPubMed CentralView ArticleGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci 1997,13(5):555-556.PubMedGoogle Scholar
- Felsenstein J: PHYLIP-Phylogeny Inference Package (Version 3.2). Cladistics 1989, 5: 164-166.Google Scholar
- Martin DP, Williamson C, Posada D: RDP2: recombination detection and analysis from sequence alignments. Bioinformatics 2005,21(2):260-262. 10.1093/bioinformatics/bth490PubMedView ArticleGoogle Scholar
- Padidam M, Sawyer S, Fauquet CM: Possible emergence of new geminiviruses by frequent recombination. Virology 1999,265(2):218-225. 10.1006/viro.1999.0056PubMedView ArticleGoogle Scholar
- Salminen MO, Carr JK, Burke DS, McCutchan FE: Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS Res Hum Retroviruses 1995,11(11):1423-1425.PubMedView ArticleGoogle Scholar
- Smith JM: Analyzing the mosaic structure of genes. J Mol Evol 1992,34(2):126-129.PubMedGoogle Scholar
- Posada D, Crandall KA: Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci USA 2001,98(24):13757-13762. 10.1073/pnas.241370698PubMedPubMed CentralView ArticleGoogle Scholar
- Gibbs MJ, Armstrong JS, Gibbs AJ: Sister-scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics 2000,16(7):573-582. 10.1093/bioinformatics/16.7.573PubMedView ArticleGoogle Scholar
- Posada D: Evaluation of methods for detecting recombination from DNA sequences: empirical data. Mol Biol Evol 2002,19(5):708-717.PubMedView ArticleGoogle Scholar
- Nielsen R, Yang Z: Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics 1998,148(3):929-936.PubMedPubMed CentralGoogle Scholar
- Yang Z, Nielsen R, Goldman N, Pedersen AM: Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics 2000,155(1):431-449.PubMedPubMed CentralGoogle Scholar
- Pond SL, Frost SD, Muse SV: HyPhy: hypothesis testing using phylogenies. Bioinformatics 2005,21(5):676-679. 10.1093/bioinformatics/bti079PubMedView ArticleGoogle Scholar
- Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol 1994,11(5):725-736.PubMedGoogle Scholar
- Shepherd CM, Borelli IA, Lander G, Natarajan P, Siddavanahalli V, Bajaj C, Johnson JE, Brooks CL 3rd, Reddy VS: VIPERdb: a relational database for structural virology. Nucleic Acids Res 2006, (34 Database):D386-389. 10.1093/nar/gkj032
- Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE: UCSF Chimera – a visualization system for exploratory research and analysis. J Comput Chem 2004,25(13):1605-1612. 10.1002/jcc.20084PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.