The Molecular Basis of Hydrolethalus Syndrome

Hydrolethalus syndrome (HLS) is a severe fetal malformation syndrome that is inherited by an autosomal recessive manner. HLS belongs to the Finnish disease heritage, an entity of rare diseases that are more prevalent in Finland than in other parts of the world. The phenotypic spectrum of the syndrome is wide and it is characterized by several developmental abnormalities, including hydrocephalus and absent midline structures in the brain, abnormal lobation of the lungs, polydactyly as well as micrognathia and other craniofacial anomalies. Polyhydramnios are relatively frequent during pregnancy. HLS can nowadays be effectively identified by ultrasound scan already at the end of the first trimester of pregnancy. One of the main goals in this study was to identify and characterize the gene defect underlying HLS. The defect was found from a previously unknown gene that was named HYLS1. Identification of the gene defect made it possible to confirm the HLS diagnosis genetically, an aspect that provides valuable information for the families in which a fetus is suspected to have HLS. Neuropathological findings of mutation confirmed HLS cases were described for the first time in detail in this study. Also, detailed general pathological findings were described. Since HYLS1 was an unknown gene with no relatives in the known gene families, many functional studies were performed in order to unravel the function of the gene and of the protein it codes for. Studies showed, for example, that the subcellular localization of the HYLS1 protein was different when the normal and the defective forms were compared. In addition, HYLS1 was shown to possess transactivation potential which was significantly diminished in the defective form. According to the results of this study it can be stated that HYLS1 most likely participates in transcriptional regulation and also in the regulation of cholesterol metabolism and that the function of HYLS1 is critical for normal fetal development.

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.

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.

Phylogeography and hybrid swarms: history of brackish water bivalve diversity in North European marginal seas

This study addressed the large-scale molecular zoogeography in two brackish water bivalve molluscs, Macoma balthica and Cerastoderma glaucum, and genetic signatures of the postglacial colonization of Northern Europe by them. The traditional view poses that M. balthica in the Baltic, White and Barents seas (i.e. marginal seas) represent direct postglacial descendants of the adjacent Northeast Atlantic populations, but this has recently been challenged by observations of close genetic affinities between these marginal populations and those of the Northeast Pacific. The primary aim of the thesis was to verify, quantify and characterize the Pacific genetic contribution across North European populations of M. balthica and to resolve the phylogeographic histories of the two bivalve taxa in range-wide studies using information from mitochondrial DNA (mtDNA) and nuclear allozyme polymorphisms. The presence of recent Pacific genetic influence in M. balthica of the Baltic, White and Barents seas, along with an Atlantic element, was confirmed by mtDNA sequence data. On a broader temporal and geographical scale, altogether four independent trans-Arctic invasions of Macoma from the Pacific since the Miocene seem to have been involved in generating the current North Atlantic lineage diversity. The latest trans-Arctic invasion that affected the current Baltic, White and Barents Sea populations probably took place in the early post-glacial. The nuclear genetic compositions of these marginal sea populations are intermediate between those of pure Pacific and Atlantic subspecies. In the marginal sea populations of mixed ancestry (Barents, White and Northern Baltic seas), the Pacific and Atlantic components are now randomly associated in the genomes of individual clams, which indicates both pervasive historical interbreeding between the previously long-isolated lineages (subspecies), and current isolation of these populations from the adjacent pure Atlantic populations. These mixed populations can be characterized as self-supporting hybrid swarms, and they arguably represent the most extensive marine animal hybrid swarms so far documented. Each of the three swarms still has a distinct genetic composition, and the relative Pacific contributions vary from 30 to 90 % in local populations. This diversity highlights the potential of introgressive hybridization to rapidly give rise to new evolutionarily and ecologically significant units in the marine realm. In the south of the Danish straits and in the Southern Baltic Sea, a broad genetic transition zone links the pure North Sea subspecies M. balthica rubra to the inner Baltic hybrid swarm, which has about 60 % of Pacific contribution in its genome. This transition zone has no regular smooth clinal structure, but its populations show strong genotypic disequilibria typical of a hybrid zone maintained by the interplay of selection and gene flow by dispersing pelagic larvae. The structure of the genetic transition is partly in line with features of Baltic water circulation and salinity stratification, with greater penetration of Atlantic genes on the Baltic south coast and in deeper water populations. In all, the scenarios of historical isolation and secondary contact that arise from the phylogeographic studies of both Macoma and Cerastoderma shed light to the more general but enigmatic patterns seen in marine phylogeography, where deep genetic breaks are often seen in species with high dispersal potential.

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.

Functional genes and gene array analysis as tools for monitoring hydrocarbon biodegradation

Bioremediation, which is the exploitation of the intrinsic ability of environmental microbes to degrade and remove harmful compounds from nature, is considered to be an environmentally sustainable and cost-effective means for environmental clean-up. However, a comprehensive understanding of the biodegradation potential of microbial communities and their response to decontamination measures is required for the effective management of bioremediation processes. In this thesis, the potential to use hydrocarbon-degradative genes as indicators of aerobic hydrocarbon biodegradation was investigated. Small-scale functional gene macro- and microarrays targeting aliphatic, monoaromatic and low molecular weight polyaromatic hydrocarbon biodegradation were developed in order to simultaneously monitor the biodegradation of mixtures of hydrocarbons. The validity of the array analysis in monitoring hydrocarbon biodegradation was evaluated in microcosm studies and field-scale bioremediation processes by comparing the hybridization signal intensities to hydrocarbon mineralization, real-time polymerase chain reaction (PCR), dot blot hybridization and both chemical and microbiological monitoring data. The results obtained by real-time PCR, dot blot hybridization and gene array analysis were in good agreement with hydrocarbon biodegradation in laboratory-scale microcosms. Mineralization of several hydrocarbons could be monitored simultaneously using gene array analysis. In the field-scale bioremediation processes, the detection and enumeration of hydrocarbon-degradative genes provided important additional information for process optimization and design. In creosote-contaminated groundwater, gene array analysis demonstrated that the aerobic biodegradation potential that was present at the site, but restrained under the oxygen-limited conditions, could be successfully stimulated with aeration and nutrient infiltration. During ex situ bioremediation of diesel oil- and lubrication oil-contaminated soil, the functional gene array analysis revealed inefficient hydrocarbon biodegradation, caused by poor aeration during composting. The functional gene array specifically detected upper and lower biodegradation pathways required for complete mineralization of hydrocarbons. Bacteria representing 1 % of the microbial community could be detected without prior PCR amplification. Molecular biological monitoring methods based on functional genes provide powerful tools for the development of more efficient remediation processes. The parallel detection of several functional genes using functional gene array analysis is an especially promising tool for monitoring the biodegradation of mixtures of hydrocarbons.

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.

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.