Identification of Sequences Encoding the Detoxification Metalloisomerase Glyoxalase I in Microbial Genomes from Several Pathogenic Organisms

The ubiquitous glyoxalase system, which is composed of two enzymes, removes cellular cytotoxic methylglyoxal (MG). In an effort to identify critical residues conserved in the evolution of the first enzyme in this system, glyoxalase I (GlxI), as well as the structural implications of sequence alterations in this enzyme, a search of the National Center for Biotechnology Information (NCBI) database of unfinished genomes was undertaken. Eleven putative GlxI sequences from pathogenic organisms were identified and analyses of these sequences in relation to the known and previously identified GlxI enzymes were performed. Several of these sequences show a very high similarity to the Escherichia coli GlxI sequence, most notably the 79% identity of the sequence identified from Yersinia pestis, the causative agent of bubonic plague. In addition to the conservation of residues critical to binding the catalytic metal in all of the proposed GlxI enzymes, four regions in the Homo sapiens GlxI enzyme are absent in all of the bacterial GlxI sequences, with the exception of Pseudomonas putida. Removal of these regions may alter the active-site conformation of the bacterial enzymes in relation to that of the H. sapiens. These differences may be targeted for the development of inhibitors selective to the bacterial enzymes. Received: 13 October 1999 / Accepted: 17 January 2000     

Phylogenetic Position of the Hexactinellida Within the Phylum Porifera Based on the Amino Acid Sequence of the Protein Kinase C from Rhabdocalyptus dawsoni

Recent analyses of genes encoding proteins typical for multicellularity, especially adhesion molecules and receptors, favor the conclusion that all metazoan phyla, including the phylum Porifera (sponges), are of monophyletic origin. However, none of these data includes cDNA encoding a protein from the sponge class Hexactinellida. We have now isolated and characterized the cDNA encoding a protein kinase C, belonging to the C subfamily (cPKC), from the hexactinellid sponge Rhabdocalyptus dawsoni. The two conserved regions, the regulatory part with the pseudosubstrate site, the two zinc fingers, and the C2 domain, as well as the catalytic domain were used for phylogenetic analyses. Sequence alignment and construction of a phylogenetic tree from the catalytic domains revealed that the yeast Saccharomyces cerevisiae and the protozoan Trypanosoma brucei are at the base of the tree, while the hexactinellid R. dawsoni branches off first among the metazoan sequences; the other two classes of the Porifera, the Calcarea (the sequence from Sycon raphanus was used) and the Demospongiae (sequences from Geodia cydonium and Suberites domuncula were used), branch off later. The statistically robust tree also shows that the two cPKC sequences from the higher invertebrates Drosophila melanogaster and Lytechinus pictus are most closely related to the calcareous sponge. This finding was also confirmed by comparing the regulatory part of the kinase gene. We suggest, that (i) within the phylum Porifera, the class Hexactinellida diverged first from a common ancestor to the Calcarea and the Demospongiae, which both appeared later, and (ii) the higher invertebrates are more closely related to the calcareous sponges. Received: 6 August 1997 / Accepted: 24 October 1997     

Signature Sequences in Diverse Proteins Provide Evidence of a Close Evolutionary Relationship Between the Deinococcus-Thermus Group and Cyanobacteria

A number of proteins have been identified that contain prominent sequence signatures that are uniquely shared by the members of the Deinococcus-Thermus genera and the cyanobacterial species but which are not found in any of the other eubacterial or archaebacterial homologs. The proteins containing such sequence signatures include (1) the DnaJ/Hsp40 family of proteins, (2) DNA polymerase I, (3) the protein synthesis elongation factor EF-Tu, and (4) the elongation factor EF-Ts. A strong affinity of the Deinococcus-Thermus species to cyanobacteria is also seen in the phylogenetic trees based on Hsp70 and DnaJ sequences. These results provide strong evidence of a close and specific evolutionary relationship between species belonging to these two eubacterial divisions. Received: 10 September 1997 / Accepted: 15 December 1997     

Did Archaeal and Bacterial Cells Arise Independently from Noncellular Precursors? A Hypothesis Stating That the Advent of Membrane Phospholipid with Enantiomeric Glycerophosphate Backbones Caused the Separation of the Two Lines of Descent

One of the most remarkable biochemical differences between the members of two domains Archaea and Bacteria is the stereochemistry of the glycerophosphate backbone of phospholipids, which are exclusively opposite. The enzyme responsible to the formation of Archaea-specific glycerophosphate was found to be NAD(P)-linked sn-glycerol-1-phosphate (G-1-P) dehydrogenase and it was first purified from Methanobacterium thermoautotrophicum cells and its gene was cloned. This structure gene named egsA (enantiomeric glycerophosphate synthase) consisted of 1,041 bp and coded the enzyme with 347 amino acid residues. The amino acid sequence deduced from the base sequence of the cloned gene (egsA) did not share any sequence similarity except for NAD-binding region with that of NAD(P)-linked sn-glycerol-3-phosphate (G-3-P) dehydrogenase of Escherichia coli which catalyzes the formation of G-3-P backbone of bacterial phospholipids, while the deduced protein sequence of the enzyme revealed some similarity with bacterial glycerol dehydrogenases. Because G-1-P dehydrogenase and G-3-P dehydrogenase would originate from different ancestor enzymes and it would be almost impossible to interchange stereospecificity of the enzymes, it seems likely that the stereostructure of membrane phospholipids of a cell must be maintained from the time of birth of the first cell. We propose here the hypothesis that Archaea and Bacteria were differentiated by the occurrence of cells enclosed by membranes of phospholipids with G-1-P and G-3-P as a backbone, respectively. Received: 24 March 1997 / Accepted: 21 May 1997     

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     

Unidentified Open Reading Frames in the Genome of Methanococcus jannaschii Are Similar in Sequence to an Archaebacterial Gene for tRNA Nucleotidyltransferase

The protein sequence of ATP/CTP:tRNA nucleotidyltransferase (cca) from Sulfolobus shibatae was used to search open reading frames in the genome of Methanococcus jannaschii. Translations of two unidentified open reading frames showed significant sequence similarity to portions of the Sulfolobus cca protein. When the two open reading frames were joined together, the expanded open reading frame was similar in sequence to the entire Sulfolobus cca protein and displayed features of the active site signature sequence proposed for members of class I enzymes within the superfamily of nucleotidyltransferases (Yue et al., 1996, RNA 2, 895–908). A possible UUG start codon was identified based on significant sequence similarity of the resulting amino-terminal region to that of Sulfolobus, and on a six-base complementarity between an adjacent upstream sequence and Methanococcus 16S rRNA. Received: 10 February 1997     

The sequences of heat shock protein 40 (DnaJ) homologs provide evidence for a close evolutionary relationship between the Deinococcus- Thermus group and cyanobacteria

The genes encoding for heat shock protein 40 (Hsp40 or DnaJ) homologs were cloned and sequenced from the archaebacterium Halobacterium cutirubrum and the eubacterium Deinococcus proteolyticus to add to sequences from the gene banks. These genes were identified downstream of the Hsp70 (or DnaK) genes in genomic fragments spanning this region and, as in other prokaryotic species, Hsp70-Hsp40 genes are likely part of the same operon. The Hsp40 homolog from D. proteolyticus was found to be lacking a central 204 base pair region present in H. cutirubrum that encodes for the four cysteine-rich domains of the repeat consensus sequence CxxCxGxG (where x is any amino acid), present in most Hsp40 homologs. The available sequences from various archaebacteria, eubacteria, and eukaryotes show that the same deletion is also present in the homologs from Thermus aquaticus and two cyanobacteria, but in no other species tested. This unique deletion and the clustering of homologs from the Deinococcus–Thermus group and cyanobacterial species in the Hsp40 phylogenetic trees suggest a close evolutionary relationship between these groups as was also shown recently for Hsp70 sequences (R.S. Gupta et al., J Bacteriol 179:345–357, 1997). Sequence comparisons indicate that the Hsp40 homologs are not as conserved as the Hsp70 sequences. Phylogenetic analysis provides no reliable information concerning evolutionary relationship between prokaryotes and eukaryotes and their usefulness in this regard is limited. However, in phylogenetic trees based on Hsp40 sequences, the two archaebacterial homologs showed a polyphyletic branching within Gram-positive bacteria, similar to that seen with Hsp70 sequences. Received: 30 January 1997 / Accepted: 22 March 1997     

Heat Shock Protein 70 Family: Multiple Sequence Comparisons, Function, and Evolution

The heat shock protein 70 kDa sequences (HSP70) are of great importance as molecular chaperones in protein folding and transport. They are abundant under conditions of cellular stress. They are highly conserved in all domains of life: Archaea, eubacteria, eukaryotes, and organelles (mitochondria, chloroplasts). A multiple alignment of a large collection of these sequences was obtained employing our symmetric-iterative ITERALIGN program (Brocchieri and Karlin 1998). Assessments of conservation are interpreted in evolutionary terms and with respect to functional implications. Many archaeal sequences (methanogens and halophiles) tend to align best with the Gram-positive sequences. These two groups also miss a signature segment [about 25 amino acids (aa) long] present in all other HSP70 species (Gupta and Golding 1993). We observed a second signature sequence of about 4 aa absent from all eukaryotic homologues, significantly aligned in all prokaryotic sequences. Consensus sequences were developed for eight groups [Archaea, Gram-positive, proteobacterial Gram-negative, singular bacteria, mitochondria, plastids, eukaryotic endoplasmic reticulum (ER) isoforms, eukaryotic cytoplasmic isoforms]. All group consensus comparisons tend to summarize better the alignments than do the individual sequence comparisons. The global individual consensus ``matches' 87% with the consensus of consensuses sequence. A functional analysis of the global consensus identifies a (new) highly significant mixed charge cluster proximal to the carboxyl terminus of the sequence highlighting the hypercharge run EEDKKRRER (one-letter aa code used). The individual Archaea and Gram-positive sequences contain a corresponding significant mixed charge cluster in the location of the charge cluster of the consensus sequence. In contrast, the four Gram-negative proteobacterial sequences of the alignment do not have a charge cluster (even at the 5% significance level). All eukaryotic HSP70 sequences have the analogous charge cluster. Strikingly, several of the eukaryotic isoforms show multiple mixed charged clusters. These clusters were interpreted with supporting data related to HSP70 activity in facilitating chaperone, transport, and secretion function. We observed that the consensus contains only a single tryptophan residue and a single conserved cysteine. This is interpreted with respect to the target rule for disaggregating misfolded proteins. The mitochondrial HSP70 connections to bacterial HSP70 are analyzed, suggesting a polyphyletic split of Trypanosoma and Leishmania protist mitochondrial (Mt) homologues separated from Mt-animal/fungal/plant homologues. Moreover, the HSP70 sequences from the amitochondrial Entamoeba histolytica and Trichomonas vaginalis species were analyzed. The E. histolytica HSP70 is most similar to the higher eukaryotic cytoplasmic sequences, with significantly weaker alignments to ER sequences and much diminished matching to all eubacterial, mitochondrial, and chloroplast sequences. This appears to be at variance with the hypothesis that E. histolytica rather recently lost its mitochondrial organelle. T. vaginalis contains two HSP70 sequences, one Mt-like and the second similar to eukaryotic cytoplasmic sequences suggesting two diverse origins. Received: 29 January 1998 / Accepted: 14 May 1998     

Similarity of Structural Features and Evolution of Satellite DNAs from Palorus subdepressus (Coleoptera) and Related Species

A novel highly abundant satellite DNA comprising 20% of the genome has been characterized in Palorus subdepressus (Insecta, Coleoptera). The 72-bp-long monomer sequence is composed of two copies of T₂A₅T octanucleotide alternating with 22-nucleotide-long elements of an inverted repeat. Phylogenetic analysis revealed clustering of monomer sequence variants into two clades. Two types of variants are prevalently organized in an alternating pattern, thus showing a tendency to generate a new complex repeating unit 144 bp in length. Fluorescent in situ hybridization revealed even distribution of the satellite in the region of pericentric heterochromatin of all 20 chromosomes. P. subdepressus satellite sequence is clearly species specific, lacking similarity even with the satellite from congeneric species P. ratzeburgii. However, on the basis of similarity in predicted tertiary structure induced by intrinsic DNA curvature and in repeat length, P. subdepressus satellite can be classified into the same group with satellites from related tenebrionid species P. ratzeburgii, Tenebrio molitor, and T. obscurus. It can be reasonably inferred that repetitive sequences of different origin evolve under constraints to adopt and conserve particular features. Obtained results suggest that the higher-order structure and repeat length, but not the nucleotide sequence itself, are maintained through evolution of these species. Received: 23 April 1997 / Accepted: 11 July 1997     

Comparative Analysis of Evolution in a Rodent Histone H2a Pseudogene

Sequences were obtained from five species of rodents that are orthologous to an H2a histone pseudogene from Mus musculus. The pseudogene is part of the cluster of replication-dependent histone genes found on Mus musculus chromosome 13. Comparative analysis of these five sequences together with the previously published sequence from M. musculus shows that this gene has likely been a pseudogene throughout the evolution of the genus Mus, while the gene from Rattus norvegicus is likely functional. Three large (>20 bp) deletions were found among the Mus pseudogenes, a feature that is very unusual compared to surveys of processed pseudogenes. In addition, there are two single-base deletions and one 4-bp insertion among the Mus pseudogenes. The species distributions of one of the large deletions and the 4-bp insertion require either independent insertions of an identical sequence, independent deletions with identical boundaries, or a deletion followed by precise reintegration of the original sequence. The evidence favors the hypothesis of multiple deletions with identical boundaries. The ``coding' regions of the Mus pseudogenes show a much reduced level of among-species variability in the 3′ half of the pseudogene, compared both to the 5′ half and to flanking sequences. This supports a hypothesis that the 3′ end of the pseudogene is the target of frequent gene conversion by functional H2a genes. Received: 1 April 1997 / Accepted: 12 June 1997     

Phylogenetic Relationships of the Glycolytic Enzyme, Glyceraldehyde-3-Phosphate Dehydrogenase, from Parabasalid Flagellates

Over 90% of the open reading frame of gap genes for glycolytic glyceraldehyde-3-phosphate dehydrogenase (GAPDH; EC 1.2.1.12) was obtained with PCR from five species of Parabasala. With gap1 from Trichomonas vaginalis obtained earlier, the data include two sequences each for three species. All sequences were colinear with T. vaginalis gap1 and shared with it as a synapomorphy a 10- to 11-residue insertion not found in any other gap and an S-loop with characteristic features of eubacterial GAPDH. All residues known to be highly conserved in this enzyme were present. The parabasalid sequences formed a robust monophyletic group in phylogenetic reconstructions with distance-based, maximum-parsimony, and maximum-likelihood methods. The two genes of the amphibian commensal, Trichomitus batrachorum, shared a common ancestor with the rest, which separate into two well-supported lineages. T. vaginalis and Tetratrichomonas gallinarum (both representatives of Trichomonadinae) formed one, while Monocercomonas sp. and Tritrichomonas foetus formed the other. These data agreed with and/or were close to published reconstructions based on other macromolecules. They did not support the ancestral position of Monocercomonas sp. proposed on the basis of morphological characteristics but confirmed an early emergence of Trichomitus batrachorum. The sequence pairs obtained from three species indicated either gene duplications subsequent to the divergence of the corresponding lineages or a strong gene conversion later in these lineages. The parabasalid clade was a robust part of the eubacterial radiation of GAPDH and showed no relationships to the clade that contained all other eukaryotic gap genes. The data clearly reveal that the members of this lineage use in their glycolytic pathway a GAPDH species with properties and an evolutionary history that are unique among all eukaryotes studied so far. Received: 28 April 1997 / Accepted: 14 October 1997     

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     

A Reexamination of the Phylogenetic Position of Callimico (Primates) Incorporating New Mitochondrial DNA Sequence Data

The New World monkeys are divided into two main groups, Callitrichidae and Cebidae. Callimico goeldii shares traits with both the Cebidae and the Callitrichidae. Recent morphological phyletic studies generally place Callimico as the most basal member of the Callitrichidae. In contrast, genetic studies (immunological, restriction fragment, and sequence data) have consistently placed Callimico somewhere within the Callitrichidae, not basal to this clade. A DNA sequence data set from the terminal 236 codons of the mitochondrial ND4 gene and the tRNA^His, tRNA^Ser, and tRNA^Leu genes was generated to clarify the position of Callimico. The sequences of 887 base pairs were analyzed by maximum-parsimony, neighbor-joining, and maximum-likelihood methods. The results of these various methods are generally congruent and place Callimico within the Callitrichidae between the marmosets (Callithrix and Cebuella) and the tamarins (Saguinus and Leontopithecus). Combined analyses of all suitable nuclear and mitochondrial gene sequences confirm the position of Callimico between the marmosets and the tamarins. As available molecular evidence indicates that Callimico is more closely related to the marmosets than to the tamarins, a reconsideration of the morphological evidence in light of the consensus tree from DNA sequence analyses is warranted. The marmosets and tamarins share four morphological characters (loss of the third molar, loss of the hypocone, reduced body size, reproductive twinning). Dwarfism may have evolved repeatedly among the Callitrichidae. It is well-known that the loss of a character can occur many times independently. The reproduction of marmosets and tamarins is extremely specialized and it is difficult to imagine that this complex and unique twinning system evolved separately in marmosets and tamarins. However, it is possible that a secondary reversal to single offspring took place in Callimico. Received: 20 March 1997 / Accepted: 17 December 1997     

Molecular Phylogenetic Inference of the Woolly Mammoth Mammuthus primigenius, Based on Complete Sequences of Mitochondrial Cytochrome b and 12S Ribosomal RNA Genes

Complete sequences of cytochrome b (1,137 bases) and 12S ribosomal RNA (961 bases) genes in mitochondrial DNA were successfully determined from the woolly mammoth (Mammuthus primigenius), African elephant (Loxodonta africana), and Asian elephant (Elephas maximus). From these sequence data, phylogenetic relationships among three genera were examined. Molecular phylogenetic trees reconstructed by the neighbor-joining and the maximum parsimony methods provided an identical topology both for cytochrome b and 12S rRNA genes. These results support the ``Mammuthus-Loxodonta' clade, which is contrary to some previous morphological reports that Mammuthus is more closely related to Elephas than to Loxodonta. Received: 8 April 1997 / Accepted: 23 July 1997     

Equally Parsimonious Pathways Through an RNA Sequence Space Are Not Equally Likely

An experimental system for determining the potential ability of sequences resembling 5S ribosomal RNA (rRNA) to perform as functional 5S rRNAs in vivo in the Escherichia coli cellular environment was devised previously. Presumably, the only 5S rRNA sequences that would have been fixed by ancestral populations are ones that were functionally valid, and hence the actual historical paths taken through RNA sequence space during 5S rRNA evolution would have most likely utilized valid sequences. Herein, we examine the potential validity of all sequence intermediates along alternative equally parsimonious trajectories through RNA sequence space which connect two pairs of sequences that had previously been shown to behave as valid 5S rRNAs in E. coli. The first trajectory requires a total of four changes. The 14 sequence intermediates provide 24 apparently equally parsimonious paths by which the transition could occur. The second trajectory involves three changes, six intermediate sequences, and six potentially equally parsimonious paths. In total, only eight of the 20 sequence intermediates were found to be clearly invalid. As a consequence of the position of these invalid intermediates in the sequence space, seven of the 30 possible paths consisted of exclusively valid sequences. In several cases, the apparent validity/invalidity of the intermediate sequences could not be anticipated on the basis of current knowledge of the 5S rRNA structure. This suggests that the interdependencies in RNA sequence space may be more complex than currently appreciated. If ancestral sequences predicted by parsimony are to be regarded as actual historical sequences, then the present results would suggest that they should also satisfy a validity requirement and that, in at least limited cases, this conjecture can be tested experimentally. Received: 27 August 1996 / Accepted: 14 April 1997     

Molecular Evolution of the Aldo-keto Reductase Gene Superfamily

The aldo-keto reductase enzymes comprise a functionally diverse gene family which catalyze the NADPH-dependant reduction of a variety of carbonyl compounds. The protein sequences of 45 members of this family were aligned and phylogenetic trees were deduced from this alignment using the neighbor-joining and Fitch algorithms. The branching order of these trees indicates that the vertebrate enzymes cluster in three groups, which have a monophyletic origin distinct from the bacterial, plant, and invertebrate enzymes. A high level of conservation was observed between the vertebrate hydroxysteroid dehydrogenase enzymes, prostaglandin F synthase, and ρ-crystallin of Xenopus laevis. We infer from the phylogenetic analysis that prostaglandin F synthase may represent a recent recruit to the eicosanoid biosynthetic pathway from the hydroxysteroid dehydrogenase pathway and furthermore that, in the context of gene recruitment, Xenopus laevisρ-crystallin may represent a shared gene. Received: 26 August 1996 / Accepted: 5 June 1997     

Analysis of the Cluster of Ribosomal Protein Genes in the Plastid Genome of a Unicellular Red Alga Cyanidioschyzon merolae: Translocation of the str Cluster as an Early Event in the Rhodophyte-Chromophyte Lineage of Plastid Evolution

The nucleotide sequence of a cluster of ribosomal protein genes in the plastid genome of a unicellular red alga, Cyanidioschyzon merolae, which has been supposed to be the most primitive alga, was determined. The phylogenetic tree inferred from the amino acid sequence of ribosomal proteins of two rhodophytes, a chromophyte, a glaucophyte, two chlorophytes (land plants), a cyanobacterium, and three eubacteria suggested a close relationship between the cyanobacterium Synechocystis PCC6803 and the plastids of various species in the kingdom Plantae, which is consistent with the hypothesis of the endosymbiotic origin of plastids. In this tree, the two species of rhodophytes were grouped with the chromophyte, and the glaucophyte was grouped with the chlorophytes. Analysis of the organization of the genes encoding the ribosomal proteins suggested that the translocation of the str cluster occurred early in the lineage of rhodophytes and chromophytes after these groups had been separated from chlorophytes and glaucophytes. Received: 2 June 1997 / Accepted: 15 July 1997     

Extensive Amplification and Transposition of a Novel Repetitive Element, Xstir, Together with Its Terminal Inverted Repeat in the Evolution of Xenopus

A DNA fragment containing short tandem repeat sequences (approximately 86-bp repeat) was isolated from a Xenopus laevis cDNA library. Southern blot and in situ hybridization analyses revealed that the repeat was highly dispersed in the genome and was present at approximately 1 million copies per haploid genome. We named this element Xstir (Xenopus short tandemly and invertedly repeating element) after its arrangement in the genome. The majority of the genomic Xstir sequences were digested to monomer and dimer sizes with several restriction enzymes. Their sequences were found to be highly homogeneous and organized into tandem arrays in the genome. Alignment analyses of several known sequences showed that some of the Xstir-like sequences were also organized into interspersed inverted repeats. The inverted repeats consisted of an inverted pair of two differently modified Xstirs separated by a short insert. In addition, these were framed by another novel inverted repeat (Xstir-TIR). The Xstir-TIR sequence was also found at the ends of tandem Xstir arrays. Furthermore, we found that Xstir-TIR was linked to a motif characterizing the T2 family which belonged to a vertebrate MITE (miniature inverted-repeat transposable element) family, suggesting the importance of Xstir-TIR for their amplification and transposition. The present study of 11 anuran and 2 urodele species revealed that Xstir or Xstir-like sequences were extensively amplified in the three Xenopus species. Genomic Xstir populations of X. borealis and X. laevis were mutually indistinguishable but significantly different from that of X. tropicalis. Received: 5 April 2000 / Accepted: 3 August 2000     

Molecular Characterization of the Principal Symbiotic Bacteria of the Weevil Sitophilus oryzae: A Peculiar G + C Content of an Endocytobiotic DNA

The principal intracellular symbiotic bacteria of the cereal weevil Sitophilus oryzae were characterized using the sequence of the 16S rDNA gene (rrs gene) and G + C content analysis. Polymerase chain reaction amplification with universal eubacterial primers of the rrs gene showed a single expected sequence of 1,501 bp. Comparison of this sequence with the available database sequences placed the intracellular bacteria of S. oryzae as members of the Enterobacteriaceae family, closely related to the free-living bacteria, Erwinia herbicola and Escherichia coli, and the endocytobiotic bacteria of the tsetse fly and aphids. Moreover, by high-performance liquid chromatography, we measured the genomic G + C content of the S. oryzae principal endocytobiotes (SOPE) as 54%, while the known genomic G + C content of most intracellular bacteria is about 39.5%. Furthermore, based on the third codon position G + C content and the rrs gene G + C content, we demonstrated that most intracellular bacteria except SOPE are A + T biased irrespective of their phylogenetic position. Finally, using the hsp60 gene sequence, the codon usage of SOPE was compared with that of two phylogenetically closely related bacteria: E. coli, a free-living bacterium, and Buchnera aphidicola, the intracellular symbiotic bacteria of aphids. Taken together, these results show a peculiar and distinctly different DNA composition of SOPE with respect to the other obligate intracellular bacteria, and, combined with biological and biochemical data, they elucidate the evolution of symbiosis in S. oryzae. Received: 8 September 1997 / Accepted: 24 October 1997