Unravelling Molecular and Cellular Disease Mechanisms in Infantile Neuronal Ceroid Lipofuscinosis (INCL)

Studying neurodegeneration provides an opportunity to gain insights into normal cell physiology, and not just pathophysiology. In this thesis work the focus is on Infantile Neuronal Ceroid Lipofuscinosis (INCL). It is a recessively inherited lysosomal storage disorder. The disease belongs to the neuronal ceroid lipofuscinoses (NCLs), a group of common progressive neurodegenerative diseases of the childhood. Characteristic accumulation of autofluorescent storage material is seen in most tissues but only neurons of the central nervous system are damaged and eventually lost during the course of the disease leaving most other cell types unaffected. The disease is caused by mutations in the CLN1 gene, but the physiological function of the corresponding protein the palmitoyl protein thioesterase (PPT1) has remained elusive. The aim of this thesis work was to shed light on the molecular and cell biological mechanisms behind INCL. This study pinpointed the localization of PPT1 in axonal presynapses of neurons. It also established the role of PPT1 in early neuronal maturation as well as importance in mature neuronal synapses. This study revealed an endocytic defect in INCL patient cells manifesting itself as delayed trafficking of receptor and non-receptor mediated endocytic markers. Furthermore, this study was the first to connect the INCL storage proteins the sphingolipid activator proteins (SAPs) A and D to pathological events on the cellular level. Abnormal endocytic processing and intracellular re-localization was demonstrated in patient cells and disease model knock-out mouse neurons. To identify early affected cellular and metabolic pathways in INCL, knock-out mouse neurons were studied by global transcript profiling and functional analysis. The gene expression analysis revealed changes in neuronal maturation and cell communication strongly associated with the regulated secretory system. Furthermore, cholesterol metabolic pathways were found to be affected. Functional studies with the knock-out mouse model revealed abnormalities in neuronal maturation as well as key neuronal functions including abnormalities in intracellular calcium homeostasis and cholesterol metabolism. Together the findings, introduced in this thesis work, support the essential role of PPT1 in developing neurons as well as synaptic sites of mature neurons. Results of this thesis also elucidate early events in INCL pathogenesis revealing defective pathways ultimately leading to the neurodegenerative process. These results contribute to the understanding of the vital physiological function of PPT1 and broader knowledge of common cellular mechanisms behind neurodegeneration. These results add to the knowledge of these severe diseases offering basis for new approaches in treatment strategies.

Mechanisms and molecular regulation of mammalian tooth replacement

In most non-mammalian vertebrates, such as fish and reptiles, teeth are replaced continuously. However, tooth replacement in most mammals, including human, takes place only once and further renewal is apparently inhibited. It is not known how tooth replacement is genetically regulated, and little is known on the physiological mechanism and evolutionary reduction of tooth replacement in mammals. In this study I have attempted to address these questions. In a rare human condition cleidocranial dysplasia, caused by a mutation in a Runt domain transcription factor Runx2, tooth replacement is continued. Runx2 mutant mice were used to investigate the molecular mechanisms of Runx2 function. Microarray analysis from dissected embryonic day 14 Runx2 mutant and wild type dental mesenchymes revealed many downstream targets of Runx2, which were validated using in situ hybridization and tissue culture methods. Wnt signaling inhibitor Dkk1 was identified as a candidate target, and in tissue culture conditions it was shown that Dkk1 is induced by FGF4 and this induction is Runx2 dependent. These experiments demonstrated a connection between Runx2, FGF and Wnt signaling in tooth development and possibly also in tooth replacement. The role of Wnt signaling in tooth replacement was further investigated by using a transgenic mouse model where Wnt signaling mediator β-catenin is continuously stabilized in dental epithelium. This stabilization led to activated Wnt signaling and to the formation of multiple enamel knots. In vitro and transplantation experiments were performed to examine the process of extra tooth formation. We showed that new teeth were continuously generated and that new teeth form from pre-existing teeth. A morphodynamic activator-inhibitor model was used to simulate enamel knot formation. By increasing the intrinsic production rate of the activator (β-catenin), the multiple enamel knot phenotype was reproduced by computer simulations. It was thus concluded that β-catenin acts as an upstream activator of enamel knots, closely linking Wnt signaling to the regulation of tooth renewal. As mice do not normally replace teeth, we used other model animals to investigate the physiological and genetic mechanisms of tooth replacement. Sorex araneus, the common shrew was earlier reported to have non-functional tooth replacement in all antemolar tooth positions. We showed by histological and gene expression studies that there is tooth replacement only in one position, the premolar 4 and that the deciduous tooth is diminished in size and disappears during embryogenesis without becoming functional. The growth rates of deciduous and permanent premolar 4 were measured and it was shown by competence inference that the early initiation of the replacement tooth in relation to the developmental stage of the deciduous tooth led to the inhibition of deciduous tooth morphogenesis. It was concluded that the evolutionary loss of deciduous teeth may involve the early activation of replacement teeth, which in turn suppress their predecessors. Mustela putorius furo, the ferret, has a dentition that resembles that of the human as ferrets have teeth that belong to all four tooth families, and all the antemolar teeth are replaced once. To investigate the replacement mechanism, histological serial sections from different embryonic stages were analyzed. It was noticed that tooth replacement is a process which involves the growth and detachment of the dental lamina from the lingual cervical loop of the deciduous tooth. Detachment of the deciduous tooth leads to a free successional dental lamina, which grows deeper into the mesenchyme, and later buds the replacement tooth. A careful 3D analysis of serial histological sections was performed and it was shown that replacement teeth are initiated from the successional dental lamina and not from the epithelium of the deciduous tooth. The molecular regulation of tooth replacement was studied and it was shown by examination of expression patterns of candidate regulatory genes that BMP/Wnt inhibitor Sostdc1 was strongly expressed in the buccal aspect of the dental lamina, and in the intersection between the detaching deciduous tooth and the successional dental lamina, suggesting a role for Sostdc1 in the process of detachment. Shh was expressed in the enamel knot and in the inner enamel epithelium in both generations of teeth supporting the view that the morphogenesis of both generations of teeth is regulated by similar mechanisms. In summary, histological and molecular studies on different model animals and transgenic mouse models were used to investigate tooth replacement. This thesis work has significantly contributed to the knowledge on the physiological mechanisms and molecular regulation of tooth replacement and its evolutionary suppression in mammals.

Substrate Selectivity and Molecular Adaptation in the Outer Membrane Proteases of Yersinia pestis and Salmonella enterica

Proteolysis is important in bacterial pathogenesis and colonization of animal and plant hosts. In this work I have investigated the functions of the bacterial outer membrane proteases, omptins, of Yersinia pestis and Salmonella enterica. Y. pestis is a zoonotic pathogen that causes plague and has evolved from gastroenteritis-causing Yersinia pseudotuberculosis about 13 000 years ago. S. enterica causes gastroenteritis and typhoid fever in humans. Omptins are transmembrane β-barrels with ten antiparallel β-strands and five surface-exposed loops. The loops are important in substrate recognition, and variation in the loop sequences leads to different substrate selectivities between omptins, which makes omptins an ideal platform to investigate functional adaptation and to alter their polypeptide substrate preferences. The omptins Pla of Y. pestis and PgtE of S. enterica are 75% identical in their amino acid sequences. Pla is a multifunctional protein with proteolytic and non-proteolytic functions, and it increases bacterial penetration and proliferation in the host. Functions of PgtE increase migration of S. enterica in vivo and bacterial survival in mouse macrophages, thus enhancing bacterial spread within the host. Mammalian plasminogen/fibrinolytic system maintains the balance between coagulation and fibrinolysis and participates in several cellular processes, e.g., cell migration and degradation of extracellular matrix proteins. This system consists of activation cascades, which are strictly controlled by several regulators, such as plasminogen activator inhibitor 1 (PAI-1), α2-antiplasmin (α2AP), and thrombin-activatable fibrinolysis inhibitor (TAFI). This work reveals novel interactions of the omptins of Y. pestis and S. enterica with the regulators of the plasminogen/fibrinolytic system: Pla and PgtE inactivate PAI-1 by cleavage at the reactive site peptide bond, and degrade TAFI, preventing its activation to TAFIa. Structure-function relationship studies with Pla showed that threonine 259 of Pla is crucial in plasminogen activation, as it prevents degradation of the plasmin catalytic domain by the omptin and thus maintains plasmin stability. In this work I constructed chimeric proteins between Pla and Epo of Erwinia pyrifoliae that share 78% sequence identity to find out which amino acids and regions in Pla are important for its functions. Epo is neither a plasminogen activator nor an invasin, but it degrades α2AP and PAI-1. Cumulative substitutions towards Pla sequence turned Epo into a Pla-like protein. In addition to threonine 259, loops 3 and 5 are critical in plasminogen activation by Pla. Turning Epo into an invasin required substitution of 31 residues located at the extracellular side of the Epo protein above the lipid bilayer, and also of the β1-strand in the N-terminal transmembrane region of the protein. These studies give an example of how omptins adapt to novel functions that advantage their host bacteria in different ecological niches.

Molecular Alterations in Asbestos-Related Lung Cancer

Background: Asbestos is a well known cancer-causing mineral fibre, which has a synergistic effect on lung cancer risk in combination with tobacco smoking. Several in vitro and in vivo experiments have demonstrated that asbestos can evoke chromosomal damage and cause alterations as well as gene expression changes. Lung tumours, in general, have very complex karyotypes with several recurrently gained and lost chromosomal regions and this has made it difficult to identify specific molecular changes related primarily to asbestos exposure. The main aim of these studies has been to characterize asbestos-related lung cancer at a molecular level. Methods: Samples from asbestos-exposed and non-exposed lung cancer patients were studied using array comparative genomic hybridization (aCGH) and fluorescent in situ hybridization (FISH) to detect copy number alterations (CNA) as well as microsatellite analysis to detect allelic imbalance (AI). In addition, asbestos-exposed cell lines were studied using gene expression microarrays. Results: Eighteen chromosomal regions showing differential copy number in the lung tumours of asbestos-exposed patients compared to those of non-exposed patients were identified. The most significant differences were detected at 2p21-p16.3, 5q35.3, 9q33.3-q34.11, 9q34.13-q34.3, 11p15.5, 14q11.2 and 19p13.1-p13.3 (p<0.005). The alterations at 2p and 9q were validated and characterized in detail using AI and FISH analysis in a larger study population. Furthermore, in vitro studies were performed to examine the early gene expression changes induced by asbestos in three different lung cell lines. The results revealed specific asbestos-associated gene expression profiles and biological processes as well as chromosomal regions enriched with genes believed to contribute to the common asbestos-related responses in the cell lines. Interestingly, the most significant region enriched with asbestos-response genes was identified at 2p22, close to the previously identified region showing asbestos-related CNA in lung tumours. Additionally, in this thesis, the dysregulated biological processes (Gene Ontology terms) detected in the cell line experiment were compared to dysregulated processes identified in patient samples in a later study (Ruosaari et al., 2008a). Commonly affected processes such as those related to protein ubiquitination, ion transport and surprisingly sensory perception of smell were identified. Conclusions: The identification of specific CNA and dysregulated biological processes shed some light on the underlying genes acting as mediators in asbestos-related lung carcinogenesis. It is postulated that the combination of several asbestos-specific molecular alterations could be used to develop a diagnostic method for the identification of asbestos-related lung cancer.

Exposure-related human cancer: Molecular changes in sinonasal cancer and lung cancer, with focus on TP53 mutations

Cancer is a leading cause of death worldwide and the total number of cancer cases continues to increase. Many cancers, for example sinonasal cancer and lung cancer, have clear external risk factors and so are potentially preventable. The occurrence of sinonasal cancer is strongly associated with wood dust exposure and the main risk factor for lung cancer is tobacco smoking. Although the molecular mechanisms involved in lung carcinogenesis have been widely studied, very little is known about the molecular changes leading to sinonasal cancer. In this work, mutations in the tumour suppressor TP53 gene in cases of sinonasal cancer and lung cancer and the associations of these mutations with exposure factors were studied. In addition, another important mechanism in many cancers, inflammation, was explored by analyzing the expression of the inflammation related enzyme, COX-2, in sinonasal cancer. The results demonstrate that TP53 mutations are frequent in sinonasal cancer and lung cancer and in both cancers they are associated with exposure. In sinonasal cancer, the occurrence of TP53 mutation significantly increased in relation to long duration and high level of exposure to wood dust. Smoking was not associated with the overall occurrence of the TP53 mutation in sinonasal cancer, but was associated with multiple TP53 mutations. Furthermore, inflammation appears to play a part in sinonasal carcinogenesis as indicated by our results showing that the expression of COX-2 was associated with adenocarcinoma type of tumours, wood dust exposure and non-smoking. In lung cancer, we detected statistically significant associations between TP53 mutations and duration of smoking, gender and histology. We also found that patients with a tumour carrying a G to T transversion, a mutation commonly found in association with tobacco smoking, had a high level of smoking-related bulky DNA adducts in their non-tumorous lung tissue. Altogether, the information on molecular changes in exposure induced cancers adds to the observations from epidemiological studies and helps to understand the role and impact of different etiological factors, which in turn can be beneficial for risk assessment and prevention.

Molecular details of phage phi6 RNA-dependent RNA synthesis

For most RNA viruses RNA-dependent RNA polymerases (RdRPs) encoded by the virus are responsible for the entire RNA metabolism. Thus, RdRPs are critical components in the viral life cycle. However, it is not fully understood how these important enzymes function during viral replication. Double-stranded RNA (dsRNA) viruses perform the synthesis of their RNA genome within a proteinacous viral particle containing an RdRP as a minor constituent. The phi6 bacteriophage is the best-studied dsRNA virus, providing an excellent background for studies of its RNA synthesis. The purified recombinant phi6 RdRP is highly active in vitro and it possesses both RNA replication and transcription activities. The crystal structure of the phi6 polymerase, solved in complex with a number of ligands, provides a working model for detailed in vitro studies of RNA-dependent RNA polymerization. In this thesis, the primer-independent initiation of the phi6 RdRP was studied in vitro using biochemical and structural methods. A C-terminal, four-amino-acid-long loop protruding into the central cavity of the phi6 RdRP has been suggested to stabilize the incoming nucleotides of the initiation complex formation through stacking interactions. A similar structural element has been found from several other viral RdRPs. In this thesis, this so-called initiation platform loop was subjected to site-directed mutagenesis to address its role in the initiation. It was found that the initiation mode of the mutants is primer-dependent, requiring either an oligonucleotide primer or a back-priming initiation mechanism for the RNA synthesis. The crystal structure of a mutant RdRP with altered initiation platform revealed a set of contacts important for primer-independent initiation. Since phi6 RdRP is structurally and functionally homologous to several viral RdRPs, among them the hepatitis C virus RdRP, these results provide further general insight to understand primer-independent initiation. In this study it is demonstrated that manganese phasing could be used as a practical tool for solving structures of large proteins with a bound manganese ion. The phi6 RdRP was used as a case study to obtain phases for crystallographic analysis. Manganese ions are naturally bound to the phi6 RdRP at the palm domain of the enzyme. In a crystallographic experiment, X-ray diffraction data from a phi6 RdRP crystal were collected at a wavelength of 1.89 Å, which is the K edge of manganese. With this data an automatically built model of the core region of the protein could be obtained. Finally, in this work terminal nucleotidyl transferase (TNTase) activity of the phi6 RdRP was documented in the isolated polymerase as well as in the viral particle. This is the first time that such an activity has been reported in a polymerase of a dsRNA virus. The phi6 RdRP used uridine triphosphates as the sole substrate in a TNTase reaction but could accept several heterologous templates. The RdRP was able to add one or a few non-templated nucleotides to the 3' end of the single- or double-stranded RNA substrate. Based on the results on particle-mediated TNTase activity and previous structural information of the polymerase, a model for termination of the RNA-dependent RNA synthesis is suggested in this thesis.

Molecular Tuning of Rhodopsins for High Visual Sensitivity in Different Light Environments : Variation in Absorbance Spectrum and Opsin Sequence within and between Species

Visual pigments of different animal species must have evolved at some stage to match the prevailing light environments, since all visual functions depend on their ability to absorb available photons and transduce the event into a reliable neural signal. There is a large literature on correlation between the light environment and spectral sensitivity between different fish species. However, little work has been done on evolutionary adaptation between separated populations within species. More generally, little is known about the rate of evolutionary adaptation to changing spectral environments. The objective of this thesis is to illuminate the constraints under which the evolutionary tuning of visual pigments works as evident in: scope, tempo, available molecular routes, and signal/noise trade-offs. Aquatic environments offer Nature s own laboratories for research on visual pigment properties, as naturally occurring light environments offer an enormous range of variation in both spectral composition and intensity. The present thesis focuses on the visual pigments that serve dim-light vision in two groups of model species, teleost fishes and mysid crustaceans. The geographical emphasis is in the brackish Baltic Sea area with its well-known postglacial isolation history and its aquatic fauna of both marine and fresh-water origin. The absorbance spectrum of the (single) dim-light visual pigment were recorded by microspectrophotometry (MSP) in single rods of 26 fish species and single rhabdoms of 8 opossum shrimp populations of the genus Mysis inhabiting marine, brackish or freshwater environments. Additionally, spectral sensitivity was determined from six Mysis populations by electroretinogram (ERG) recording. The rod opsin gene was sequenced in individuals of four allopatric populations of the sand goby (Pomatoschistus minutus). Rod opsins of two other goby species were investigated as outgroups for comparison. Rod absorbance spectra of the Baltic subspecies or populations of the primarily marine species herring (Clupea harengus membras), sand goby (P. minutus), and flounder (Platichthys flesus) were long-wavelength-shifted compared to their marine populations. The spectral shifts are consistent with adaptation for improved quantum catch (QC) as well as improved signal-to-noise ratio (SNR) of vision in the Baltic light environment. Since the chromophore of the pigment was pure A1 in all cases, this has apparently been achieved by evolutionary tuning of the opsin visual pigment. By contrast, no opsin-based differences were evident between lake and sea populations of species of fresh-water origin, which can tune their pigment by varying chromophore ratios. A more detailed analysis of differences in absorbance spectra and opsin sequence between and within populations was conducted using the sand goby as model species. Four allopatric populations from the Baltic Sea (B), Swedish west coast (S), English Channel (E), and Adriatic Sea (A) were examined. Rod absorbance spectra, characterized by the wavelength of maximum absorbance (λmax), differed between populations and correlated with differences in the spectral light transmission of the respective water bodies. The greatest λmax shift as well as the greatest opsin sequence difference was between the Baltic and the Adriatic populations. The significant within-population variation of the Baltic λmax values (506-511 nm) was analyzed on the level of individuals and was shown to correlate well with opsin sequence substitutions. The sequences of individuals with λmax at shorter wavelengths were identical to that of the Swedish population, whereas those with λmax at longer wavelengths additionally had substitution F261F/Y in the sixth transmembrane helix of the protein. This substitution (Y261) was also present in the Baltic common gobies and is known to redshift spectra. The tuning mechanism of the long-wavelength type Baltic sand gobies is assumed to be the co-expression of F261 and Y261 in all rods to produce ≈ 5 nm redshift. The polymorphism of the Baltic sand goby population possibly indicates ambiguous selection pressures in the Baltic Sea. The visual pigments of all lake populations of the opossum shrimp (Mysis relicta) were red-shifted by 25 nm compared with all Baltic Sea populations. This is calculated to confer a significant advantage in both QC and SNR in many humus-rich lakes with reddish water. Since only A2 chromophore was present, the differences obviously reflect evolutionary tuning of the visual protein, the opsin. The changes have occurred within the ca. 9000 years that the lakes have been isolated from the Sea after the most recent glaciation. At present, it seems that the mechanism explaining the spectral differences between lake and sea populations is not an amino acid substitution at any other conventional tuning site, but the mechanism is yet to be found.

Molecular genetics of RECQL4 syndromes

RAPADILINO syndrome is an autosomally resessively inherited condition that belongs to a group of rare syndromes more common in Finland than in other parts of the world. RAPADILINO is characterized by pre- and postnatal growth retardation, radial ray defects, diarrhoea of unknown aetiology during chilhood, a facial resemblance with other patients and normal intelligence. In Finland, 15 patients with this condition have been found which compares with only five patients in other parts of the world. We found RECQL4 gene mutations in RAPADILINO patients and proved this syndrome to be allelic with a subgroup of Rothmund-Thomson syndrome (RTS). Later we found RECQL4 mutations in patients with Baller-Gerold syndrome (BGS). These three syndromes share clinical findings and differential diagnostics rely on poikiloderma and craniosynostosis not seen in RAPADILINO syndrome. We found five different mutations in the Finnish RAPADILINO patients. The g.2545delT mutation is the founder mutation in the Finnish population as all the patients are either homozygotes or compound heterozygotes for it. This mutation leads to the inframe skipping of exon seven from mRNA. The protein encoded by this mutant mRNA lacks the nuclear retention signal and thus leads to the mislocalization of the mutant protein. The genotype-phenotype correlation is not straightforward but it seems that RAPADILINO could be due to alteration in protein function and truncating mutations in both alleles are more common among RTS patients. RTS patients with RECQL4 mutations have an elevated risk for osteosarcoma, but their risk to develop other types of malignancies is not increased.Two Finnish RAPADILINO patients have been diagnosed with osteosarcoma, but in addition to this we have found an excess of lymphoma cases among the Finnish RAPADILINO patients. This difference between cancer types could be due to different mutations found in these syndromes. The mutation screening of the patients will help to differentiate patients who have RECQL4 mutations and thus the elevated cancer risk. Patients will benefit from the follow up since early detection of malignancies is important for the treatment.

Phospholipids of lipid-containing bacteriophages and their transbilayer distribution

In this study we used electro-spray ionization mass-spectrometry to determine phospholipid class and molecular species compositions in bacteriophages PM2, PRD1, Bam35 and phi6 as well as their hosts. To obtain compositional data of the individual leaflets, phospholipid transbilayer distribution in the viral membranes was studied. We found that 1) the membranes of all studied bacteriophage are enriched in PG as compared to the host membranes, 2) molecular species compositions in the phage and host membranes are similar, and 3) phospholipids in the viral membranes are distributed asymmetrically with phosphatidylglycerol enriched in the outer leaflet and phosphatidylethanolamine in the inner one (except Bam35). Alternative models for selective incorporation of phospholipids to phages and for the origins of the asymmetric phospholipid transbilayer distribution are discussed. Notably, the present data are also useful when constructing high resolution structural models of bacteriophages, since diffraction methods cannot provide a detailed structure of the membrane due to high motility of the lipids and lack of symmetric organization of membrane proteins.

Molecular epidemiology of human rhinoviruses

The first part of this work investigates the molecular epidemiology of a human enterovirus (HEV), echovirus 30 (E-30). This project is part of a series of studies performed in our research team analyzing the molecular epidemiology of HEV-B viruses. A total of 129 virus strains had been isolated in different parts of Europe. The sequence analysis was performed in three different genomic regions: 420 nucleotides (nt) in the VP4/VP2 capsid protein coding region, the entire VP1 capsid protein coding gene of 876 nt, and 150 nt in the VP1/2A junction region. The analysis revealed a succession of dominant sublineages within a major genotype. The temporally earlier genotypes had been replaced by a genetically homogenous lineage that has been circulating in Europe since the late 1970s. The same genotype was found by other research groups in North America and Australia. Globally, other cocirculating genetic lineages also exist. The prevalence of a dominant genotype makes E-30 different from other previously studied HEVs, such as polioviruses and coxsackieviruses B4 and B5, for which several coexisting genetic lineages have been reported. The second part of this work deals with molecular epidemiology of human rhinoviruses (HRVs). A total of 61 field isolates were studied in the 420-nt stretch in the capsid coding region of VP4/VP2. The isolates were collected from children under two years of age in Tampere, Finland. Sequences from the clinical isolates clustered in the two previously known phylogenetic clades. Seasonal clustering was found. Also, several distinct serotype-like clusters were found to co-circulate during the same epidemic season. Reappearance of a cluster after disappearing for a season was observed. The molecular epidemiology of the analyzed strains turned out to be complex, and we decided to continue our studies of HRV. Only five previously published complete genome sequences of HRV prototype strains were available for analysis. Therefore, all designated HRV prototype strains (n=102) were sequenced in the VP4/VP2 region, and the possibility of genetic typing of HRV was evaluated. Seventy-six of the 102 prototype strains clustered in HRV genetic group A (HRV-A) and 25 in group B (HRV-B). Serotype 87 clustered separately from other HRVs with HEV species D. The field strains of HRV represented as many as 19 different genotypes, as judged with an approximate demarcation of a 20% nt difference in the VP4/VP2 region. The interserotypic differences of HRV were generally similar to those reported between different HEV serotypes (i.e. about 20%), but smaller differences, less than 10%, were also observed. Because some HRV serotypes are genetically so closely related, we suggest that the genetic typing be performed using the criterion "the closest prototype strain". This study is the first systematic genetic characterization of all known HRV prototype strains, providing a further taxonomic proposal for classification of HRV. We proposed to divide the genus Human rhinoviruses into HRV-A and HRV-B. The final part of the work comprises a phylogenetic analysis of a subset (48) of HRV prototype strains and field isolates (12) in the nonstructural part of the genome coding for the RNA-dependent RNA polymerase (3D). The proposed division of the HRV strains in the species HRV-A and HRV-B was also supported by 3D region. HRV-B clustered closer to HEV species B, C, and also to polioviruses than to HRV-A. Intraspecies variation within both HRV-A and HRV-B was greater in the 3D coding region than in the VP4/VP2 coding region, in contrast to HEV. Moreover, the diversity of HRV in 3D exceeded that of HEV. One group of HRV-A, designated HRV-A', formed a separate cluster outside other HRV-A in the 3D region. It formed a cluster also in the capsid region, but located within HRV-A. This may reflect a different evolutionary history of distinct genomic regions among HRV-A. Furthermore, the tree topology within HRV-A in the 3D region differed from that in the VP4/VP2, suggesting possible recombination events in the evolution of the strains. No conflicting phylogenies were observed in any of the 12 field isolates. Possible recombination was further studied using the Similarity and Bootscanning analyses of the complete genome sequences of HRV available in public databases. Evidence for recombination among HRV-A was found, as HRV2 and HRV39 showed higher similarity in the nonstructural part of the genome. Whether HRV2 and HRV39 strains - and perhaps also some other HRV-A strains not yet completely sequenced - are recombinants remains to be determined.