Evolution of the enzymes of the citric acid cycle and the glyoxylate cycle of higher plants. A case study of endosymbiotic gene transfer.

The citric acid or tricarboxylic acid cycle is a central element of higher-plant carbon metabolism which provides, among other things, electrons for oxidative phosphorylation in the inner mitochondrial membrane, intermediates for amino-acid biosynthesis, and oxaloacetate for gluconeogenesis from succinate derived from fatty acids via the glyoxylate cycle in glyoxysomes. The tricarboxylic acid cycle is a typical mitochondrial pathway and is widespread among alpha-proteobacteria, the group of eubacteria as defined under rRNA systematics from which mitochondria arose. Most of the enzymes of the tricarboxylic acid cycle are encoded in the nucleus in higher eukaryotes, and several have been previously shown to branch with their homologues from alpha-proteobacteria, indicating that the eukaryotic nuclear genes were acquired from the mitochondrial genome during the course of evolution. Here, we investigate the individual evolutionary histories of all of the enzymes of the tricarboxylic acid cycle and the glyoxylate cycle using protein maximum likelihood phylogenies, focusing on the evolutionary origin of the nuclear-encoded proteins in higher plants. The results indicate that about half of the proteins involved in this eukaryotic pathway are most similar to their alpha-proteobacterial homologues, whereas the remainder are most similar to eubacterial, but not specifically alpha-proteobacterial, homologues. A consideration of (a) the process of lateral gene transfer among free-living prokaryotes and (b) the mechanistics of endosymbiotic (symbiont-to-host) gene transfer reveals that it is unrealistic to expect all nuclear genes that were acquired from the alpha-proteobacterial ancestor of mitochondria to branch specifically with their homologues encoded in the genomes of contemporary alpha-proteobacteria. Rather, even if molecular phylogenetics were to work perfectly (which it does not), then some nuclear-encoded proteins that were acquired from the alpha-proteobacterial ancestor of mitochondria should, in phylogenetic trees, branch with homologues that are no longer found in most alpha-proteobacterial genomes, and some should reside on long branches that reveal affinity to eubacterial rather than archaebacterial homologues, but no particular affinity for any specific eubacterial donor.     

Cloning, expression, and substrate specificity of a fungal chymotrypsin. Evidence for lateral gene transfer from an actinomycete bacterium

Unlike trypsins, chymotrypsins have not until now been found in fungi. Expressed sequence tag analysis of the deuteromycete Metarhizium anisopliae identified two trypsins (family S1) and a novel chymotrypsin (CHY1). CHY1 resembles actinomycete (bacterial) chymotrypsins (family S2) rather than other eukaryote enzymes (family S1) in being synthesized as a precursor species (374 amino acids, pI/MW: 5.07/38,279) containing a large N-terminal fragment (186 amino acids). Chy1 was expressed in Pichia pastoris yielding an enzyme with a chymotrypsin specificity for branched aliphatic and aromatic C-terminal amino acids. This is predictable as key catalytic residues determining the specificity of Streptomyces griseus chymotrypsins are conserved with CHY1. Mature (secreted) CHY1 (pI/MW: 8.29/18,499) shows closest overall amino acid identity to S. griseus protease C (55%) and clustered with other secreted bacterial S2 chymotrypsins that diverged widely from animal and endocellular bacterial enzymes in phylogenetic trees of the chymotrypsin superfamily. Conversely, actinomycete chymotrypsins are much more closely related to fungal proteases than to other eubacterial sequences. Complete genomes of yeast, gram eubacteria, archaebacteria, and mitochondria do not contain paralogous genes. Expressed sequence tag data bases from other fungi also lack chymotrypsin homologs. In light of this patchy distribution, we conclude that chy1 probably arose by lateral gene transfer from an actinomycete bacterium.     

Spirochete contributions to the eukaryotic genome

Expanded genome/proteome databases and effective use of sequence alignment tools make it possible to trace the phylogeny of individual eukaryotic proteins and ultimately to identify the prokaryotes that contributed to the last eukaryotic common ancestor (LECA). I developed an application of reciprocal BLASTp that identifies (1) the prokaryotic lineages that have contributed to the nuclear genome and (2) the specific proteins acquired from prokaryotic ancestors. Eight complete eubacterial proteomes were analyzed: two free-living spirochetes, two clostridia, two actinobacteria, and two proteobacteria (one alpha and one gamma). The data reveal a spirochete genetic contribution to the eukaryotic genome including essential proteins involved in DNA binding and repair, cyclic nucleotide metabolism, acyltransferase, and signal transduction. My results, consistent with the sulfur syntrophy hypothesis that posits LECA evolved from a merger of spirochetes (eubacteria) with sulfidogenic eocytes (archaebacteria), confirm the contribution of mitochondrial genes from alpha-proteobacteria. A contribution from clostridia to eukaryote genomes was also detected whereas none was seen from either actinobacterium or Escherichia coli. The complete spirochete and clostridial genetic contributions to eukaryotes and those of other eu-and archaebacteria can be identified by this method.     

A novel group of type I polyketide synthases (PKS) in animals and the complex phylogenomics of PKSs

Type I polyketide synthases (PKSs), and related fatty acid synthases (FASs), represent a large group of proteins encoded by a diverse gene family that occurs in eubacteria and eukaryotes (mainly in fungi). Collectively, enzymes encoded by this gene family produce a wide array of polyketide compounds that encompass a broad spectrum of biological activity including antibiotic, antitumor, antifungal, immunosuppressive, and predator defense functional roles. We employed a phylogenomics approach to estimate relationships among members of this gene family from eubacterial and eukaryotic genomes. Our results suggest that some animal genomes (sea urchins, birds, and fish) possess a previously unidentified group of pks genes, in addition to possessing fas genes used in fatty acid metabolism. These pks genes in the chicken, fish, and sea urchin genomes do not appear to be closely related to any other animal or fungal genes, and instead are closely related to pks genes from the slime mold Dictyostelium and eubacteria. Continued accumulation of genome sequence data from diverse animal lineages is required to clarify whether the presence of these (non-fas) pks genes in animal genomes owes their origins to horizontal gene transfer (from eubacterial or Dictostelium genomes) or to more conventional patterns of vertical inheritance coupled with massive gene loss in several animal lineages. Additionally, results of our broad-scale phylogenetic analyses bolster the support for previous hypotheses of horizontal gene transfer of pks genes from bacterial to fungal and protozoan lineages.     

Nucleotide bias causes a genomewide bias in the amino acid composition of proteins

We analyzed the nucleotide contents of several completely sequenced genomes, and we show that nucleotide bias can have a dramatic effect on the amino acid composition of the encoded proteins. By surveying the genes in 21 completely sequenced eubacterial and archaeal genomes, along with the entire Saccharomyces cerevisiae genome and two Plasmodium falciparum chromosomes, we show that biased DNA encodes biased proteins on a genomewide scale. The predicted bias affects virtually all genes within the genome, and it could be clearly seen even when we limited the analysis to sets of homologous gene sequences. Parallel patterns of compositional bias were found within the archaea and the eubacteria. We also found a positive correlation between the degree of amino acid bias and the magnitude of protein sequence divergence. We conclude that mutational bias can have a major effect on the molecular evolution of proteins. These results could have important implications for the interpretation of protein-based molecular phylogenies and for the inference of functional protein adaptation from comparative sequence data.     

Molecular phylogenies based on ribosomal protein L11, L1, L10, and L12 sequences

Summary Available sequences that correspond to the E. coli ribosomal proteins L11, L1, L10, and L12 from eubacteria, archaebacteria, and eukaryotes have been aligned. The alignments were analyzed qualitatively for shared structural features and for conservation of deletions or insertions. The alignments were further subjected to quantitative phylogenetic analysis, and the amino acid identity between selected pairs of sequences was calculated. In general, eubacteria, archaebacteria, and eukaryotes each form coherent and well-resolved nonoverlapping phylogenetic domains. The degree of diversity of the four proteins between the three groups is not uniform. For L11, the eubacterial and archaebacterial proteins are very similar whereas the eukaryotic L11 is clearly less similar. In contrast, in the case of the L12 proteins and to a lesser extent the L10 proteins, the archaebacterial and eukaryotic proteins are similar whereas the eubacterial proteins are different. The eukaryotic L1 equivalent protein has yet to be identified. If the root of the universal tree is near or within the eubacterial domain, our ribosomal protein-based phylogenies indicate that archaebacteria are monophyletic. The eukaryotic lineage appears to originate either near or within the archaebacterial domain. Correspondence to: P. Dennis     

Origin and evolution of the mitochondrial proteome.

The endosymbiotic theory for the origin of mitochondria requires substantial modification. The three identifiable ancestral sources to the proteome of mitochondria are proteins descended from the ancestral alpha-proteobacteria symbiont, proteins with no homology to bacterial orthologs, and diverse proteins with bacterial affinities not derived from alpha-proteobacteria. Random mutations in the form of deletions large and small seem to have eliminated nonessential genes from the endosymbiont-mitochondrial genome lineages. This process, together with the transfer of genes from the endosymbiont-mitochondrial genome to nuclei, has led to a marked reduction in the size of mitochondrial genomes. All proteins of bacterial descent that are encoded by nuclear genes were probably transferred by the same mechanism, involving the disintegration of mitochondria or bacteria by the intracellular membranous vacuoles of cells to release nucleic acid fragments that transform the nuclear genome. This ongoing process has intermittently introduced bacterial genes to nuclear genomes. The genomes of the last common ancestor of all organisms, in particular of mitochondria, encoded cytochrome oxidase homologues. There are no phylogenetic indications either in the mitochondrial proteome or in the nuclear genomes that the initial or subsequent function of the ancestor to the mitochondria was anaerobic. In contrast, there are indications that relatively advanced eukaryotes adapted to anaerobiosis by dismantling their mitochondria and refitting them as hydrogenosomes. Accordingly, a continuous history of aerobic respiration seems to have been the fate of most mitochondrial lineages. The initial phases of this history may have involved aerobic respiration by the symbiont functioning as a scavenger of toxic oxygen. The transition to mitochondria capable of active ATP export to the host cell seems to have required recruitment of eukaryotic ATP transport proteins from the nucleus. The identity of the ancestral host of the alpha-proteobacterial endosymbiont is unclear, but there is no indication that it was an autotroph. There are no indications of a specific alpha-proteobacterial origin to genes for glycolysis. In the absence of data to the contrary, it is assumed that the ancestral host cell was a heterotroph.     

Comprehensive analysis of the origin of eukaryotic genomes

There is currently no consensus on the evolutionary origin of eukaryotes. In the search of the ancestors of eukaryotes, we analyzed the phylogeny of 46 genomes, including those of 2 eukaryotes, 8 archaea, and 36 eubacteria. To avoid the effects of gene duplications, we used inparalog pairs of genes with orthologous relationships. First, we grouped these inparalogs into the functional categories of the nucleus, cytoplasm, and mitochondria. Next, we counted the sister groups of eukaryotes in prokaryotic phyla and plotted them on a standard phylogenetic tree. Finally, we used Pearson's chi-square test to estimate the origin of the genomes from specific prokaryotic ancestors. The results suggest the eukaryotic nuclear genome descends from an archaea that was neither euryarchaeota nor crenarchaeota and that the mitochondrial genome descends from alpha-proteobacteria. In contrast, genes related to the cytoplasm do not appear to originate from a specific group of prokaryotes.     

Origin and Evolution of the Mitochondrial Proteome

The endosymbiotic theory for the origin of mitochondria requires substantial modification. The three identifiable ancestral sources to the proteome of mitochondria are proteins descended from the ancestral α-proteobacteria symbiont, proteins with no homology to bacterial orthologs, and diverse proteins with bacterial affinities not derived from α-proteobacteria. Random mutations in the form of deletions large and small seem to have eliminated nonessential genes from the endosymbiont-mitochondrial genome lineages. This process, together with the transfer of genes from the endosymbiont-mitochondrial genome to nuclei, has led to a marked reduction in the size of mitochondrial genomes. All proteins of bacterial descent that are encoded by nuclear genes were probably transferred by the same mechanism, involving the disintegration of mitochondria or bacteria by the intracellular membranous vacuoles of cells to release nucleic acid fragments that transform the nuclear genome. This ongoing process has intermittently introduced bacterial genes to nuclear genomes. The genomes of the last common ancestor of all organisms, in particular of mitochondria, encoded cytochrome oxidase homologues. There are no phylogenetic indications either in the mitochondrial proteome or in the nuclear genomes that the initial or subsequent function of the ancestor to the mitochondria was anaerobic. In contrast, there are indications that relatively advanced eukaryotes adapted to anaerobiosis by dismantling their mitochondria and refitting them as hydrogenosomes. Accordingly, a continuous history of aerobic respiration seems to have been the fate of most mitochondrial lineages. The initial phases of this history may have involved aerobic respiration by the symbiont functioning as a scavenger of toxic oxygen. The transition to mitochondria capable of active ATP export to the host cell seems to have required recruitment of eukaryotic ATP transport proteins from the nucleus. The identity of the ancestral host of the α-proteobacterial endosymbiont is unclear, but there is no indication that it was an autotroph. There are no indications of a specific α-proteobacterial origin to genes for glycolysis. In the absence of data to the contrary, it is assumed that the ancestral host cell was a heterotroph.     

NADH: ubiquinone oxidoreductase from bovine heart mitochondria. A fourth nuclear encoded subunit with a homologue encoded in chloroplast genomes.

The amino acid sequence has been determined of the precursor of a nuclear encoded 20 kDa subunit of complex I from bovine heart mitochondria. The sequence of the mature protein is related to a protein of uncertain function, hitherto known as psbG, encoded in the chloroplast genomes of higher plants. Open reading frames encoding homologues of psbG have also been detected in bacteria and in the mitochondrial genome of Paramecium tetraurelia. The chloroplast psbG gene is found between ndhC and ndhJ, which encode homologues of ND3, a hydrophobic subunit of complex I encoded in the bovine mitochondrial genome, and of the nuclear encoded 30 kDa subunit of complex I. This 20 kDa protein is the eleventh out of the forty or more subunits of bovine complex I with a chloroplast encoded homologue, and its sequence provides further support for the presence in chloroplasts of a multisubunit enzyme related to complex I that could be involved in chlororespiration. The strict conservation of three cysteines suggests that the subunit might be an iron-sulphur protein.     

Cloning of the HSP70 gene from Halobacterium marismortui: relatedness of archaebacterial HSP70 to its eubacterial homologs and a model for the evolution of the HSP70 gene.

Heat shock induces the synthesis of a set of proteins in Halobacterium marismortui whose molecular sizes correspond to the known major heat shock proteins. By using the polymerase chain reaction and degenerate oligonucleotide primers for conserved regions of the 70-kDa heat shock protein (HSP70) family, we have successfully cloned and sequenced a gene fragment containing the entire coding sequence for HSP70 from H. marismortui. HSP70 from H. marismortui shows between 44 and 47% amino acid identity with various eukaryotic HSP70s and between 51 and 58% identity with its eubacterial and archaebacterial homologs. On the basis of a comparison of all available HSP70 sequences, we have identified a number of unique sequence signatures in this protein family that provide a clear distinction between eukaryotic organisms and prokaryotic organisms (archaebacteria and eubacteria). The archaebacterial (viz., H. marismortui and Methanosarcina mazei) HSP70s have been found to contain all of the signature sequences characteristic of eubacteria (particularly the gram-positive bacteria), which suggests a close evolutionary relationship between these groups. In addition, detailed analyses of HSP70 sequences that we have carried out have revealed a number of additional novel features of the HSP70 protein family. These include (i) the presence of an insertion of about 25 to 27 amino acids in the N-terminal quadrants of all known eukaryotic and prokaryotic HSP70s except those from archaebacteria and the gram-positive group of bacteria, (ii) significant sequence similarity in HSP70 regions comprising its first and second quadrants from organisms lacking the above insertion, (iii) highly significant similarity between a protein, MreB, of Escherichia coli and the N-terminal half of HSP70s, (iv) significant sequence similarity between the N-terminal quadrant of HSP70 (from gram-positive bacteria and archaebacteria) and the m-type thioredoxin of plant chloroplasts. To account for these and other observations, a model for the evolution of HSP70 proteins involving gene duplication is proposed. The model proposes that HSP70 from archaebacteria (H. marismortui and M. mazei) and the gram-positive group of bacteria constitutes the ancestral form of the protein and that all other HSP70s (viz., other eubacteria as well as eukaryotes) containing the insert have evolved from this ancient protein.     

The ubiquitous potential Z-forming sequence of eucaryotes, (dT-dG)n . (dC-dA)n, is not detectable in the genomes of eubacteria, archaebacteria, or mitochondria.

The potential Z-forming sequence (dT-dG)n . (dC-dA)n is an abundant, interspersed repeat element that is ubiquitous in eucaryotic nuclear genomes. We report that in contrast to eucaryotic nuclear DNA, the genomes of eubacteria, archaebacteria, and mitochondria lack this sequence, since even a single tract of greater than or equal to 14 base pairs in length is not detectable through either hybridization or sequence analysis. Interestingly, the phylogenetic distribution of the (dT-dG)n . (dC-dA)n repeat exhibits a striking parallel to that of (dT-dC)n . (dG-dA)n, but not to other homocopolymeric sequences such as (dC-dG)n . (dC-dG)n or (dT-dA)n . (dT-dA)n.     

Conservation of ribosomal protein gene ordering in 16 complete genomes

The organization of ribosomal proteins in 16 prokaryotic genomes was studied as an example of comparative genome analyses of gene systems. Hypothetical ribosomal protein-containing operons were constructed. These operons also contained putative genes and other non-ribosomal genes. The correspondences among these genes across different organisms were clarified by sequence homology computations. In this way a cross tabulation of 70 ribosomal proteins genes was constructed. On average, these were organized into 9-14 operons in each genome. There were also 25 non-ribosomal or putative genes in these mainly ribosomal protein operons. Hence the table contains 95 genes in total. It was found that: (i) the conservation of the block of about 20 r-proteins in the L3 and L4 operons across almost the entire eubacteria and archaebacteria is remarkable; (ii) some operons only belong to eubacteria or archaebacteria; (iii) although the ribosomal protein operons are highly conserved within domain, there are fine variations in some operons across different organisms within each domain, and these variations are informative on the evolutionary relations among the organisms. This method provides a new potential for studying the origin and evolution of old species.     

Sequence comparison of glyceraldehyde-3-phosphate dehydrogenases from the three urkingdoms: evolutionary implication

The primary structure of the glyceraldehyde-3-phosphate dehydrogenase from the archaebacteria shows striking deviation from the known sequences of eubacterial and eukaryotic sequences, despite unequivocal homologies in functionally important regions. Thus, the structural similarity between the eubacterial and eukaryotic enzymes is significantly higher than that between the archaebacterial enzymes and the eubacterial and eukaryotic enzymes. This preferred similarity of eubacterial and eukaryotic glyceraldehyde-3-phosphate dehydrogenase structures does not correspond to the phylogenetic distances among the three urkingdoms as deduced from comparisons of ribosomal ribonucleic acid sequences. Indications will be presented that the closer relationship of the eubacterial and eukaryotic glyceraldehyde-3-phosphate dehydrogenase resulted from a gene transfer from eubacteria to eukaryotes after the segregation of the three urkingdoms.     

Bifunctional aldehyde/alcohol dehydrogenase (ADHE) in chlorophyte algal mitochondria

Protein profiles of mitochondria isolated from the heterotrophic chlorophyte Polytomella sp. grown on ethanol at pH 6.0 and pH 3.7 were analyzed by Blue Native and denaturing polyacrylamide gel electrophoresis. Steady-state levels of oxidative phosphorylation complexes were influenced by external pH. Levels of an abundant, soluble, mitochondrial protein of 85 kDa and its corresponding mRNA increased at pH 6.0 relative to pH 3.7. N-terminal and internal sequencing of the 85 kDa mitochondrial protein together with the corresponding cDNA identified it as a bifunctional aldehyde/alcohol dehydrogenase (ADHE) with strong similarity to homologues from eubacteria and amitochondriate protists. A mitochondrial targeting sequence of 27 amino acids precedes the N-terminus of the mature mitochondrial protein. A gene encoding an ADHE homologue was also identified in the genome of Chlamydomonas reinhardtii, a photosynthetic relative of Polytomella. ADHE reveals a complex picture of sequence similarity among homologues. The lack of ADHE from archaebacteria indicates a eubacterial origin for the eukaryotic enzyme. Among eukaryotes, ADHE has hitherto been characteristic of anaerobes since it is essential to cytosolic energy metabolism of amitochondriate protists such as Giardia intestinalis and Entamoeba histolytica. Its abundance and expression pattern suggest an important role for ADHE in mitochondrial metabolism of Polytomella under the conditions studied. The current data are compatible with the view that Polytomella ADHE could be involved either in ethanol production or assimilation, or both, depending upon environmental conditions. Presence of ADHE in an oxygen-respiring algal mitochondrion and co-expression at ambient oxygen levels with respiratory chain components is unexpected with respect to the view that eukaryotes acquired ADHE genes specifically as an adaptation to an anaerobic lifestyle.     

Phylogenomic evidence for a myxococcal contribution to the mitochondrial fatty acid beta-oxidation

Background

The origin of eukaryotes remains a fundamental question in evolutionary biology. Although it is clear that eukaryotic genomes are a chimeric combination of genes of eubacterial and archaebacterial ancestry, the specific ancestry of most eubacterial genes is still unknown. The growing availability of microbial genomes offers the possibility of analyzing the ancestry of eukaryotic genomes and testing previous hypotheses on their origins.

Methodology/Principal Findings

Here, we have applied a phylogenomic analysis to investigate a possible contribution of the Myxococcales to the first eukaryotes. We conducted a conservative pipeline with homologous sequence searches against a genomic sampling of 40 eukaryotic and 357 prokaryotic genomes. The phylogenetic reconstruction showed that several eukaryotic proteins traced to Myxococcales. Most of these proteins were associated with mitochondrial lipid intermediate pathways, particularly enzymes generating reducing equivalents with pivotal roles in fatty acid β-oxidation metabolism. Our data suggest that myxococcal species with the ability to oxidize fatty acids transferred several genes to eubacteria that eventually gave rise to the mitochondrial ancestor. Later, the eukaryotic nucleocytoplasmic lineage acquired those metabolic genes through endosymbiotic gene transfer.

Conclusions/Significance

Our results support a prokaryotic origin, different from α-proteobacteria, for several mitochondrial genes. Our data reinforce a fluid prokaryotic chromosome model in which the mitochondrion appears to be an important entry point for myxococcal genes to enter eukaryotes.     

Coverage of whole proteome by structural genomics observed through protein homology modeling database

We have been developing FAMSBASE, a protein homology-modeling database of whole ORFs predicted from genome sequences. The latest update of FAMSBASE (), which is based on the protein three-dimensional (3D) structures released by November 2003, contains modeled 3D structures for 368,724 open reading frames (ORFs) derived from genomes of 276 species, namely 17 archaebacterial, 130 eubacterial, 18 eukaryotic and 111 phage genomes. Those 276 genomes are predicted to have 734,193 ORFs in total and the current FAMSBASE contains protein 3D structure of approximately 50% of the ORF products. However, cases that a modeled 3D structure covers the whole part of an ORF product are rare. When portion of an ORF with 3D structure is compared in three kingdoms of life, in archaebacteria and eubacteria, approximately 60% of the ORFs have modeled 3D structures covering almost the entire amino acid sequences, however, the percentage falls to about 30% in eukaryotes. When annual differences in the number of ORFs with modeled 3D structure are calculated, the fraction of modeled 3D structures of soluble protein for archaebacteria is increased by 5%, and that for eubacteria by 7% in the last 3 years. Assuming that this rate would be maintained and that determination of 3D structures for predicted disordered regions is unattainable, whole soluble protein model structures of prokaryotes without the putative disordered regions will be in hand within 15 years. For eukaryotic proteins, they will be in hand within 25 years. The 3D structures we will have at those times are not the 3D structure of the entire proteins encoded in single ORFs, but the 3D structures of separate structural domains. Measuring or predicting spatial arrangements of structural domains in an ORF will then be a coming issue of structural genomics.     

Organelle origins and ribosomal RNA

As the detailed molecular biology of organelle genomes has unfolded, there has been a general acceptance of the view that plastids and mitochondria are of endosymbiotic, eubacterial origin. Plastid genes are strikingly similar to their eubacterial (particularly cyanobacterial) counterparts in sequence, organization, and mode of expression, and such features strongly support the hypothesis that the plastid and its genome were derived in evolution from a blue-green alga-like endosymbiont. Mitochondria, on the other hand, are problematic: mitochondrial genes are organized and expressed in remarkably diverse ways in the different major groups of eukaryotes, and in no case are these features particularly characteristic of either bacterial or nuclear genomes. There is, however, clear evidence derived from gene sequence supporting the eubacterial ancestry of mitochondria, and some of the most compelling data have come from analyses of mitochondrial ribosomal RNA (rRNA). Plant mitochondrial rRNA genes diverge in sequence at a particularly slow rate, and these genes have proven to be especially supportive of the endosymbiont hypothesis, pointing to an origin of mitochondria from within the alpha subdivision of the purple bacteria. Ribosomal RNA sequences provide a basis for the construction of global phylogenetic trees that probe the evolutionary history of organelles, and that address the question of whether mitochondria and plastids are monophyletic or polyphyletic in origin. Such studies raise the possibility that the rRNA genes of plant mitochondria originated separately from the mitochondrial rRNA genes of other eukaryotes.     

ApA Dinucleotide Periodicity in Prokaryote, Eukaryote, and Organelle Genomes

Computer analyses of various genome sequences revealed the existence of certain periodical patterns of adenine–adenine dinucleotides (ApA). For each genome sequence of 13 eubacteria, 3 archaebacteria, 10 eukaryotes, 60 mitochondria, and 9 chloroplasts, we counted frequencies of ApA dinucleotides at each downstream position within 50 bp from every ApA. We found that the complete genomes of all three archaebacteria have clear ApA periodicities of about 10 bps. On the other hand, all of the 13 eubacteria we analyzed were found to have an ApA periodicity of about 11 bp. Similar periodicities exist in the 10 eukaryotes, although higher organisms such as primates tend to have weaker periodic patterns. None of the mitochondria and chroloplasts we analyzed showed an evident periodic pattern. Received: 3 November 1998 / Accepted: 24 March 1999