Detection of HTLV-I in peripheral blood lymphocytes from patients with chronic HTLV-I-associated myelopathy/tropical spastic paraparesis and asymptomatic carriers by PCR-in situ hybridization

Less than 5% of people infected with human T-lymphotropic virus type I (HTLV-I) develop HTLV-I-associated myelopathy/tropical spastic paraparesis (HAM/TSP), a chronic progressive neurologic disease. A number of factors have been implicated in the development of HAM/TSP including heterogeneity of viral sequences, host-genetic background, viral-specific cellular immune responses and viral load. This study examined the presence of HTLV-Itax DNA in peripheral blood lymphocytes (PBL) from 2 chronic HAM/TSP patients and 2 asymptomatic HTLV-I carriers by using PCR-in situ hybridization (PCR-ISH) for the in situ presence of proviral HTLV-Itax DNA. By this technique, rare PBL from these HTLV-I-infected individuals contained HTLV-I DNA. PCR-ISH did not detect any difference in the number of infected cells between HAM/TSP patients and asymptomatic carriers.     

Lack of tax diversity for tropical spastic paraparesis/human T-cell lymphotropic virus type 1 (HTLV-I) associated myelopathy development in HTLV-I-infected subjects in São Paulo, Brazil

The product of human T-cell lymphotropic virus type 1 (HTLV-1) tax gene has a transactivating effect of the viral and cellular gene expression. Genetic variations in this gene have been correlated with differences in clinical outcomes. Based upon its diversity, two closely related substrains, namely tax A and tax B, have been described. The tax A substrain has been found at a higher frequency among human T-cell leukemia virus type 1 (TSP/HAM) patients than among healthy HTLV-I-infected asymptomatic subjects in Japan. In this study, we determined the distribution of tax substrains in HTLV-I-infected subjects in the city of S?o Paulo, Brazil. Using the ACCII restriction enzyme site, we detected only tax A substrain from 48 TSP/HAM patients and 28 healthy HTLV-I carriers. The sequenced tax genes from nine TSP/HAM patients and five asymptomatic HTLV-I carriers showed a similar pattern of mutation, which characterizes tax A. Our results indicate that HTLV-I tax subtypes have no significant influences on TSP/HAM disease progression. Furthermore, monophyletic introduction of HTLV-I to Brazil probably occurred during the African slave trade many years ago.     

Ferritin gene organization: Differences between plants and animals suggest possible kingdom-specific selective constraints

Ferritin, a protein widespread in nature, concentrates iron ∼10¹¹–10¹²-fold above the solubility within a spherical shell of 24 subunits; it derives in plants and animals from a common ancestor (based on sequence) but displays a cytoplasmic location in animals compared to the plastid in contemporary plants. Ferritin gene regulation in plants and animals is altered by development, hormones, and excess iron; iron signals target DNA in plants but mRNA in animals. Evolution has thus conserved the two end points of ferritin gene expression, the physiological signals and the protein structure, while allowing some divergence of the genetic mechanisms. Comparison of ferritin gene organization in plants and animals, made possible by the cloning of a dicot (soybean) ferritin gene presented here and the recent cloning of two monocot (maize) ferritin genes, shows evolutionary divergence in ferritin gene organization between plants and animals but conservation among plants or among animals; divergence in the genetic mechanism for iron regulation is reflected by the absence in all three plant genes of the IRE, a highly conserved, noncoding sequence in vertebrate animal ferritin mRNA. In plant ferritin genes, the number of introns (n= 7) is higher than in animals (n= 3). Second, no intron positions are conserved when ferritin genes of plants and animals are compared, although all ferritin gene introns are in the coding region; within kingdoms, the intron positions in ferritin genes are conserved. Finally, secondary protein structure has no apparent relationship to intron/exon boundaries in plant ferritin genes, whereas in animal ferritin genes the correspondence is high. The structural differences in introns/exons among phylogenetically related ferritin coding sequences and the high conservation of the gene structure within plant or animal kingdoms suggest that kingdom-specific functional constraints may exist to maintain a particular intron/exon pattern within ferritin genes. In the case of plants, where ferritin gene intron placement is unrelated to triplet codons or protein structure, and where ferritin is targeted to the plastid, the selection pressure on gene organization may relate to RNA function and plastid/nuclear signaling. Received: 25 July 1995 / Accepted: 3 October 1995     

The transactivator gene of human T-cell leukemia virus type I is more variable within and between healthy carriers than patients with tropical spastic paraparesis.

Human T-cell leukemia virus type I (HTLV-1) causes T-cell leukemia and tropical spastic paraparesis (TSP) in a minority of infected people, whereas the majority remain healthy. No association between a particular HTLV-I sequence and disease manifestation has been found in previous studies. We studied here the sequence variability of the gene for the HTLV-I Tax protein, which is the dominant target antigen of the very strong cytotoxic T-lymphocyte response to the virus. In HTLV-I infection, the intraisolate nucleotide variability is much greater than the variability between isolates. The predicted protein sequence of Tax was significantly more variable in the healthy seropositive individuals' provirus than in those of the patients with TSP. Thus, tax sequence heterogeneity, rather than the presence of particular sequences, distinguishes healthy HTLV-I-seropositive individuals from patients with TSP.     

Multistability in a Model for CTL Response to HTLV-I Infection and Its Implications to HAM/TSP Development and Prevention

Human T-cell leukaemia/lymphoma virus type I (HTLV-I) is a retrovirus that has been identified as the causative agent of HTLV-I-associated myelopathy/tropical spastic paraparesis (HAM/TSP) and other illnesses. HTLV-I infects primarily CD4⁺ T cells and the transmission occurs through direct cell-to-cell contact. HAM/TSP patients harbor higher proviral loads in peripheral blood lymphocytes than asymptomatic carriers. Also, HAM/TSP patients exhibit a remarkably high number of circulating HTLV-I-specific CD8⁺ cytotoxic T lymphocytes (CTLs) in the peripheral blood. While CTLs have a protective role by killing the infected cells and lowering the proviral load, a high level of CTLs and their cytotoxicity are believed to be a main cause of the development of HAM/TSP. A mathematical model for HTLV-I infection of CD4⁺ T cells that incorporates the CD8⁺ cytotoxic T-cell (CTL) response is investigated. Our mathematical analysis reveals that the system can stabilize at a carrier steady-state with persistent viral infection but no CTL response, or at a HAM/TSP steady-state at which both the viral infection and CTL response are persistent. We also establish two threshold parameters R ₀ and R ₁, the basic reproduction numbers for viral persistence and for CTL response, respectively. We show that the parameter R ₁ can be used to distinguish asymptomatic carriers from HAM/TSP patients, and as an important control parameter for preventing the development of HAM/TSP.     

Heterogeneity of the glutathione transferase genes encoding enzymes responsible for insecticide degradation in the housefly

One of the four glutathione-S-transferases (GST) that is overproduced in the insecticide-resistant Cornell-R strain of the housefly (Musca domestica) produces an activity that degrades the insecticide dimethyl parathion and conjugates glutathione to lindane. In earlier work, it was shown that the resistant Cornell-R carries an amplification, probably a duplication, of one or more of its GST loci and that this amplification is directly related to resistance. Using polymerase chain reaction (PCR) amplification with genomic DNA, multiple copies of the gene encoding the parathion-degrading activity (called MdGst-3) were subcloned from both the ancestral, insecticide-susceptible strain BPM and from the insecticide-resistant Cornell-R. In BPM, three different MdGst-3 genes were identified while in Cornell-R, 12 different MdGst-3 sequences were found that, though closely related to ancestral genes, had diverged by a few nucleotides. This diversity in MdGst-3 genomic sequences in Cornell-R is reflected in the expressed sequences, as sampled through a cDNA bank. Population heterozygosity cannot account for these multiple GST genes. We suggest that selection for resistance to insecticides has resulted in not only amplification of the MdGst-3 genes but also in the divergence of sequence between the amplified copies. Received: 22 November 1995 / Accepted: 23 February 1996     

Evolutionary dynamics of the human endogenous retrovirus family HERV-K inferred from full-length proviral genomes

Several distinct families of endogenous retroviruses exist in the genomes of primates. Most of them are remnants of ancient germ-line infections. The human endogenous retrovirus family HERV-K represents the unique known case of endogenous retrovirus that amplified in the human genome after the divergence of human and chimpanzee lineages. There are two types of HERV-K proviral genomes differing by the presence or absence of 292 bp in the pol-env boundary. Human-specific insertions exist for both types. The analyses shown in the present work reveal that several lineages of type 1 and type 2 HERV-K proviruses remained transpositionally active after the human/chimpanzee split. The data also reflect the important role of mosaic evolution (either by recombination or gene conversion) during the evolutionary history of HERV-K. Received: 5 February 2001 / Accepted: 22 March 2001     

Of Worms and Men: An Evolutionary Perspective on the Fibroblast Growth Factor (FGF) and FGF Receptor Families

FGFs (fibroblast growth factors) play major roles in a number of developmental processes. Recent studies of several human disorders, and concurrent analysis of gene knock-out and properties of the corresponding recombinant proteins have shown that FGFs and their receptors are prominently involved in the development of the skeletal system in mammals. We have compared the sequences of the nine known mammalian FGFs, FGFs from other vertebrates, and three additional sequences that we extracted from existing databases: two human FGF sequences that we tentatively designated FGF10 and FGF11, and an FGF sequence from C?norhabditis elegans. Similarly, we have compared the sequences of the four FGF receptor paralogs found in chordates with four non-chordate FGF receptors, including one recently identified in C. elegans. The comparison of FGF and FGF receptor sequences in vertebrates and nonvertebrates shows that the FGF and FGF receptor families have evolved through phases of gene duplications, one of which may have coincided with the emergence of vertebrates, in relation with their new system of body scaffold. Received: 6 April 1996 / Accepted: 5 July 1996     

Molecular detection and isolation of human T-cell lymphotropic virus type I (HTLV-I) from patients with HAM/TSP in São Paulo,Brazil

Background: Infection with HTLV-I is etiologically linked with HTLV-I-associated myelopathy/tropical spastic paraparesis (HAM/TSP). However some patients with chronic progressive paraparesis resembling HAM/TSP have been shown to be infected with HTLV-II.Objective: To clarify the role of each of these human retroviruses in the etiology of HAM/TSP in São Paulo, Brazil.Study design: A detailed serological and molecular analysis of HTLV-I/II infection was performed in a cohort of 19 patients with HAM/TSP attending a neurological clinic.Results: Plasma samples analyzed for anti-HTLV-I/II antibodies using a Western blot assay, comprising HTLV-I (rgp46^I)- and HTLV-II (rgp46^II)-specific recombinant env epitopes, demonstrated reactivity to rgp46^I and hence were typed as seropositive for HTLV-I. Presence of HTLV genomic sequences in peripheral blood mononuclear cells (PBMC) was sought after by PCR using consensus primers SK 110 and SK 111 for the pol region of HTLV proviral DNA followed by hybridization with type-specific probes—SK 112 (HTLV-I) and SK 188 (HTLV-II). Southern blots from all individuals hybridized with SK 112 but not with SK 188, further confirming HTLV-I infection. Cocultivation of PBMC from eight of these patients with activated lymphocytes from normal individuals resulted in active viral production, detected as presence of soluble p24^gag antigen in culture supernatants. Investigation of risk factors for HTLV-I infection in these individuals revealed that five out of 19 patients studied (26.3%) had received blood transfusions previous to disease onset.Conclusions: We demonstrate HTLV-I as the only viral type involved in the etiology of HAM/TSP in a cohort from São Paulo, Brazil, and emphasize that prevention measures, including widespread routine screening of blood donations for HTLV should be conducted in Brazil.     

The Pyrophosphate-Dependent Phosphofructokinase of the Protist, Trichomonas vaginalis, and the Evolutionary Relationships of Protist Phosphofructokinases

The pyrophosphate-dependent phosphofructokinase (PP_i-PFK) of the amitochondriate protist Trichomonas vaginalis has been purified. The enzyme is a homotetramer of about 50 kDa subunits and is not subject to allosteric regulation. The protein was fragmented and a number of peptides were sequenced. Based on this information a PCR product was obtained from T. vaginalis gDNA and used to isolate corresponding cDNA and gDNA clones. Southern analysis indicated the presence of five genes. One open reading frame (ORF) was completely sequenced and for two others the 5′ half of the gene was determined. The sequences were highly similar. The complete ORF corresponded to a polypeptide of about 46 kDa. All the peptide sequences obtained were present in the derived sequences. The complete ORF was highly similar to that of other PFKs, primarily in its amino-terminal half. The T. vaginalis enzyme was most similar to PP_i-PFK of the mitochondriate heterolobosean, Naegleria fowleri. Most of the residues shown or assumed to be involved in substrate binding in other PP_i-PFKs were conserved in the T. vaginalis enzyme. Direct comparison and phylogenetic reconstruction revealed a significant divergence among PP_i-PFKs and related enzymes, which can be assigned to at least four distantly related groups, three of which contain enzymes of protists. The separation of these groups is supported with a high percentage of bootstrap proportions. The short T. vaginalis PFK shares a most recent common ancestor with the enzyme from N. fowleri. This pair is clearly separated from a group comprising the long (>60-kDa) enzymes from Giardia lamblia, Entamoeba histolytica pfk2, the spirochaetes Borrelia burgdorferi and Trepomena pallidum, as well as the α- and β-subunits of plant PP_i-PFKs. The third group (``X') containing protist sequences includes the glycosomal ATP-PFK of Trypanosoma brucei, E. histolytica pfk1, and a second sequence from B. burgdorferi. The fourth group (``Y') comprises cyanobacterial and high-G + C, Gram-positive eubacterial sequences. The well-studied PP_i-PFK of Propionibacterium freudenreichii is highly divergent and cannot be assigned to any of these groups. These four groups are well separated from typical ATP-PFKs, the phylogenetic analysis of which confirmed relationships established earlier. These findings indicate a complex history of a key step of glycolysis in protists with several early gene duplications and possible horizontal gene transfers. Received: 5 December 1997 / Accepted: 18 March 1998     

The Heat-Shock Gene hsp83 of Drosophila auraria: Genomic Organization, Nucleotide Sequence, and Long Antiparallel Coupled ORFs (LAC ORFs)

The genomic organization of the hsp83 gene of Drosophila auraria, a far-eastern endemic species belonging to the montium subgroup of the melanogaster species group, is presented here. Based on in situ hybridization on polytene chromosomes, cDNA and genomic clone mapping, nucleotide sequencing, and genomic Southern analysis, hsp83 is shown to be present as a single-copy gene at locus 64B on the 3L chromosome arm in D. auraria. This gene is organized into two exons separated by a 929-bp intron. The first exon represents the mRNA leader sequence and is not translated, while the coding region, having a length of 2,151 bp, is solely included in the second exon. Nucleotide sequence comparisons of D. auraria hsp83 with homologous sequences from other organisms show high conservation of the coding region (88–92% identity) in the genus Drosophila, in addition to the conserved genomic organization of two-exons–one-intron, of comparable size and arrangement. A phylogenetic tree based on the protein sequences of homologous genes from representative organisms is in accord with the accredited phylogenetic position of D. auraria. In the hsp83 gene region, a second case of long antiparallel coupled open reading frames (LAC ORFs) for this species was found. The antiparallel to the hsp83 gene ORF is 1,554 bases long, while the two ORFs overlap has a size of 1,548 bp. The anti-hsp83 ORF does not show significant homology to any known gene sequences. In addition, no similar LAC ORF structures were found in homologous gene regions of other organisms. Received: 18 April 1997 / Accepted: 1 August 1997     

Evolution of the Larval Cuticle Proteins Coded by the Secondary Sex Chromosome Pair: X2 and Neo-Y of Drosophila miranda: II. Comparison at the Amino Acid Sequence Level

The larval cuticle proteins (LCPs) are encoded by a multigene family, Lcp1–4, located at the right arm of the metacentric autosome 2 (2R) in Drosophila melanogaster. Due to a chromosome fusion the Lcp locus of Drosophila miranda is situated on a pair of secondary sex chromosomes, the X2 and neo-Y chromosomes. Comparing the deduced amino acid sequences of the autosomal D. melanogaster loci with the sex-chromosomal loci of D. miranda, we were able to trace the evolution of the Lcp loci with respect to their different chromosomal inheritance. The length of the signal peptide is conserved in all four LCPs, while the size of the mature LCPs varies. Conserved protein motifs became obvious from the alignment, indicating regions of structural and functional importance. Analyzing intra- and interspecific sequence similarities of the Lcp gene families allowed us to reconstruct the phylogeny of the gene cluster. Alignment with cuticle amino acid sequences originating from divergent insect species reveals motifs already present in the primordial insect LCPs. These motifs indicate different levels of constraint acting during the evolution of the LCPs. Received: 27 December 1995 / Accepted: 30 April 1996     

Evolutionary Dynamics of HTLV-I

Using mathematical models to describe the in vivo dynamics of HTLV-I infection, an explanation is offered for the slow rate of evolution of HTLV-I relative to HIV-1. In agreement with experimental findings, it is assumed that cell activation is required for successful replication in T helper cells and that HTLV-I induces a significant degree of bystander activation. It is found that the rate of evolution of HTLV-I is limited by the restricted availability of activated uninfected T cells, both at high and low proviral loads. This limits the within-host sequence diversity of HTLV-I and may therefore account for the slow rate of evolution of the virus in the population. Specific differences in the in vivo dynamics of HTLV-I and HIV-1 are identified which may account for the discrepancy in the rate of evolution of these two retroviruses. Received: 7 September 1999 / Accepted: 6 December 1999     

Accelerated Evolution of Cytochrome b in Simian Primates: Adaptive Evolution in Concert with Other Mitochondrial Proteins?

We have sequenced the cytochrome b gene of Horsfield's tarsier, Tarsius bancanus, to complete a data set of sequences for this gene from representatives of each primate infraorder. These primate cytochrome b sequences were combined with those from representatives of three other mammalian orders (cat, whale, and rat) in an analysis of relative evolutionary rates. The nonsynonymous nucleotide substitution rate of the cytochrome b gene has increased approximately twofold along lineages leading to simian primates compared to that of the tarsier and other primate and nonprimate mammalian species. However, the rate of transversional substitutions at fourfold degenerate sites has remained uniform among all lineages. This increase in the evolutionary rate of cytochrome b is similar in character and magnitude to that described previously for the cytochrome c oxidase subunit II gene. We propose that the evolutionary rate increase observed for cytochrome b and cytochrome c oxidase subunit II may underlie an episode of coadaptive evolution of these two proteins in the mitochondria of simian primates. Received: 15 December 1997 / Accepted: 24 February 1998     

Genetic stability of human T lymphotropic virus type I despite antiviral pressures by CTLs

Human T lymphotropic virus type I (HTLV-I)-associated myelopathy/tropical spastic paraparesis (HAM/TSP) is an inflammatory neurological disease. Patients with HAM/TSP show high proviral load despite increased HTLV-I Tax-specific CTL. It is still unknown whether the CTL efficiently eliminate the virus in vivo and/or whether a naturally occurring variant virus becomes predominant by escaping from the CTL. To address these issues, we sequenced a large number of HTLV-I tax genes from HLA-A*02 HAM/TSP patients and estimated synonymous and nonsynonymous changes of the genes to detect positive selection pressure on the virus. We found the pressures in three of six CTL epitopes in HTLV-I Tax, where amino acid substitutions preferentially occurred. Although some of variant viruses were not recognized by the CTL, no variant viruses accumulated within 3-8 years, indicating genetic stability of HTLV-I tax gene. These results suggest that CTL eliminate the infected cells in vivo and naturally occurring variant viruses do not predominate. As Tax is a regulatory protein which controls viral replication, the amino acid substitutions in Tax may reduce viral fitness for replication. Viral fitness and host immune response may contribute to the viral evolution within the infected individuals. Furthermore, the genetic stability in the epitopes despite the antiviral pressures suggests that the three epitopes can be the candidate targets for HTLV-I vaccine development.     

Intron Conservation in a UV-Specific DNA Repair Gene Encoded by Chlorella Viruses

Large dsDNA-containing chlorella viruses encode a pyrimidine dimer-specific glycosylase (PDG) that initiates repair of UV-induced pyrimidine dimers. The PDG enzyme is a homologue of the bacteriophage T4-encoded endonuclease V. The pdg gene was cloned and sequenced from 42 chlorella viruses isolated over a 12-year period from diverse geographic regions. Surprisingly, the pdg gene from 15 of these 42 viruses contain a 98-nucleotide intron that is 100% conserved among the viruses and another 4 viruses contain an 81-nucleotide intron, in the same position, that is nearly 100% identical (one virus differed by one base). In contrast, the nucleotides in the pdg coding regions (exons) from the intron-containing viruses are 84 to 100% identical. The introns in the pdg gene have 5′-AG/GTATGT and 3′-TTGCAG/AA splice site sequences which are characteristic of nuclear-located, spliceosomal processed pre-mRNA introns. The 100% identity of the 98-nucleotide intron sequence in the 15 viruses and the near-perfect identity of an 81-nucleotide intron sequence in another 4 viruses imply strong selective pressure to maintain the DNA sequence of the intron when it is in the pdg gene. However, the ability of intron-plus and intron-minus viruses to repair UV-damaged DNA in the dark was nearly identical. These findings contradict the widely accepted dogma that intron sequences are more variable than exon sequences. Received: 13 May 1999 / Accepted: 20 August 1999     

The tryptophan biosynthetic pathway of aphid endosymbionts (Buchnera): Genetics and evolution of plasmid-associated anthranilate synthase (trpEG) within the aphididae

The bacterial endosymbionts (Buchnera) from the aphids Rhopalosiphum padi, R. maidis, Schizaphis graminum, and Acyrthosiphon pisum contain the genes for anthranilate synthase (trpEG) on plasmids made up of one or more 3.6-kb units. Anthranilate synthase is the first as well as the rate-limiting enzyme in the tryptophan biosynthetic pathway. The amplification of trpEG on plasmids may result in an increase of enzyme protein and overproduction of this essential amino acid, which is required by the aphid host. The nucleotide sequence of trpEG from endosymbionts of different species of aphids is highly conserved, as is an approximately 500-bp upstream DNA segment which has the characteristics of an origin of replication. Phylogenetic analyses were performed using trpE and trpG from the endosymbionts of these four aphids as well as from the endosymbiont of Schlechtendalia chinensis, in which trpEG occurs on the chromosome. The resulting phylogeny was congruent with trees derived from sequences of two chromosome-located bacterial genes (part of trpB and 16S ribosomal DNA). In turn, trees obtained from plasmid-borne and bacterial chromosome-borne sequences were congruent with the tree resulting from phylogenetic analysis of three aphid mitochondrial regions (portions of the small and large ribosomal DNA subunits, as well as cytochrome oxidase II). Congruence of trees based on genes from host mitochondria and from bacteria adds to previous support for exclusively vertical transmission of the endosymbionts within aphid lineages. Congruence with trees based on plasmid-borne genes supports the origin of the plasmid-borne trpEG from the chromosomal genes of the same lineage and the absence of subsequent plasmid exchange among endosymbionts of different species of aphids. Received: 22 August 1995 / Accepted: 6 September 1995     

The Dentin Matrix Protein 1 Gene of Prototherian and Metatherian Mammals

Mineralization of tooth dentin (the deposition of hydroxyapatite crystals in and around collagen type I fibers of the extracellular matrix) requires the involvement of several genes, among them the gene coding for the dentin matrix protein 1, DMP1. We determined the exon–intron organization of the cattle DMP1 gene and used this information to amplify by the polymerase chain reaction homologous gene fragments from the genomic DNA of two species of metatherian (marsupial) mammals and one prototherian (monotreme) species. The translated proto- and metatherian protein sequences are highly divergent from the eutherian sequences but retain the general characteristics of the DMP1 (high acidity, serine-richness, multiple glycosylation sites, and the presence of the RGD cell attachment tripeptide). They therefore appear to be functional even though, evolutionarily, teeth are in a regression phase in prototherians. It is possible, therefore, that DMP1 is also involved in other functions besides dentinogenesis. The DMP1 gene appears to evolve rapidly and apparently tolerates non-frame-shifting insertions/deletions throughout the coding sequence. Received: 22 February 1998 / Accepted: 11 July 1998     

Multiple Substitutions Affect the Phylogenetic Utility of Cytochrome b and 12S rDNA Data: Examining a Rapid Radiation in Leporid (Lagomorpha) Evolution

Partial sequences of two mitochondrial genes, the 12S ribosomal gene (739 bp) and the cytochrome b gene (672 bp), were analyzed in hopes of reconstructing the evolutionary relationships of 11 leporid species, representative of seven genera. However, partial cytochrome b sequences were of little phylogenetic value in this study. A suite of pairwise comparisons between taxa revealed that at the intergeneric level, the cytochrome b gene is saturated at synonymous coding positions due to multiple substitution events. Furthermore, variation at the nonsynonymous positions is limited, rendering the cytochrome b gene of little phylogenetic value for assessing the relationships between leporid genera. If the cytochrome b data are analyzed without accounting for these two classes of nucleotides (i.e., synonymous and nonsynonymous sites), one may incorrectly conclude that signal exists in the cytochrome b data. The mitochondrial 12S rRNA gene, on the other hand, has not experienced excessive saturation at either stem or loop positions. Phylogenies reconstructed from the 12S rDNA data support hypotheses based on fossil evidence that African rock rabbits (Pronolagus) are outside of the main leporid stock and that leporids experienced a rapid radiation. However, the molecular data suggest that this radiation event occurred in the mid-Miocene several millions of years earlier than the Pleistocene dates suggested by paleontological evidence. Received: 23 April 1998 / Accepted: 14 May 1998     

Molecular Evolutionary Analysis of the Thiamine-Diphosphate-Dependent Enzyme,Transketolase

Members of the transketolase group of thiamine-diphosphate-dependent enzymes from 17 different organisms including mammals, yeast, bacteria, and plants have been used for phylogenetic reconstruction. Alignment of the amino acid and DNA sequences for 21 transketolase enzymes and one putative transketolase reveals a number of highly conserved regions and invariant residues that are of predicted importance for enzyme activity, based on the crystal structure of yeast transketolase. One particular sequence of 36 residues has some similarities to the nucleotide-binding motif and we designate it as the transketolase motif. We report further evidence that the recP protein from Streptococcus pneumoniae might be a transketolase and we list a number of invariant residues which might be involved in substrate binding. Phylogenies derived from the nucleotide and the amino acid sequences by various methods show a conventional clustering for mammalian, plant, and gram-negative bacterial transketolases. The branching order of the gram-positive bacteria could not be inferred reliably. The formaldehyde transketolase (sometimes known as dihydroxyacetone synthase) of the yeast Hansenula polymorpha appears to be orthologous to the mammalian enzymes but paralogous to the other yeast transketolases. The occurrence of more than one transketolase gene in some organisms is consistent with several gene duplications. The high degree of similarity in functionally important residues and the fact that the same kinetic mechanism is applicable to all characterized transketolase enzymes is consistent with the proposition that they are all derived from one common ancestral gene. Transketolase appears to be an ancient enzyme that has evolved slowly and might serve as a model for a molecular clock, at least within the mammalian clade. Received: 13 September 1995 / Accepted: 14 November 1996