Higher-plant chloroplast and cytosolic 3-phosphoglycerate kinases: a case of endosymbiotic gene replacement

Exploring New Routes for Genetic Resistances to Potyviruses: The Case of the Arabidopsis thaliana Phosphoglycerates Kinases (PGK) Metabolic Enzymes

The development of recessive resistance by loss of susceptibility is a consistent strategy to combat and limit damages caused by plant viruses. Susceptibility genes can be turned into resistances, a feat that can either be selected among the plant’s natural diversity or engineered by biotechnology. Here, we summarize the current knowledge on the phosphoglycerate kinases (PGK), which have emerged as a new class of susceptibility factors to single-stranded positive RNA viruses, including potyviruses. PGKs are metabolic enzymes involved in glycolysis and the carbon reduction cycle, encoded by small multigene families in plants. To fulfil their role in the chloroplast and in the cytosol, PGKs genes encode differentially addressed proteins. Here, we assess the diversity and homology of chloroplastic and cytosolic PGKs sequences in several crops and review the current knowledge on their redundancies during plant development, taking Arabidopsis as a model. We also show how PGKs have been shown to be involved in susceptibility—and resistance—to viruses. Based on this knowledge, and drawing from the experience with the well-characterized translation initiation factors eIF4E, we discuss how PGKs genes, in light of their subcellular localization, function in metabolism, and susceptibility to viruses, could be turned into efficient genetic resistances using genome editing techniques.

The Phosphoglycerate Kinase (PGK) Gene Family of Maize (Zea mays var. B73)

Phosphoglycerate kinase (PGK, E.C. 2.7.2.3) interconverts ADP + 1,3-bisphospho-glycerate (1,3-bPGA) to ATP + 3-phosphoglycerate (3PGA). While most bacteria have a single pgk gene and mammals possess two copies, plant genomes contain three or more PGK genes. In this study, we identified five Pgk genes in the Zea mays var. B73 genome, predicted to encode proteins targeted to different subcellular compartments: ZmPgk1, ZmPgk2, and ZmPgk4 (chloroplast), ZmPgk3 (cytosol), and ZmPgk5 (nucleus). The expression of ZmPgk3 was highest in non-photosynthetic tissues (roots and cobs), where PGK activity was also greatest, consistent with a function in glycolysis. Green tissues (leaf blade and husk leaf) showed intermediate levels of PGK activity, and predominantly expressed ZmPgk1 and ZmPgk2, suggesting involvement in photosynthetic metabolism. ZmPgk5 was weakly expressed and ZmPgk4 was not detected in any tissue. Phylogenetic analysis showed that the photosynthetic and glycolytic isozymes of plants clustered together, but were distinct from PGKs of animals, fungi, protozoa, and bacteria, indicating that photosynthetic and glycolytic isozymes of plants diversified after the divergence of the plant lineage from other groups. These results show the distinct role of each PGK in maize and provide the basis for future studies into the regulation and function of this key enzyme.

Pooled CRISPR/Cas9 reveals redundant roles of plastidial phosphoglycerate kinases in carbon fixation and metabolism

Phosphoglycerate kinase (PGK) is a highly conserved reversible enzyme that participates in both glycolysis and photosynthesis. In Arabidopsis thaliana, one cytosolic PGK (PGKc) and two plastidial PGKs (PGKp) are known. It remains debatable whether the two PGKp isozymes are functionally redundant or specialized in plastidial carbon metabolism and fixation. Here, using a pooled clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) strategy, we found that plants with single mutations in pgkp1 or pgkp2 were not significantly affected, whereas a pgkp1pgkp2 double mutation was lethal due to retarded carbon fixation, suggesting that PGKp isozymes play redundant functional roles. Metabolomic analysis demonstrated that the sugar-deficient pgkp1pgkp2 double mutation was partially complemented by exogenous sugar, although respiration intermediates were not rescued. Chloroplast development was defective in pgkp1pgkp2, due to a deficiency in glycolysis-dependent galactoglycerolipid biosynthesis. Ectopic expression of a plastid targeting PGKc did not reverse the pgkp1pgkp2 double-mutant phenotypes. Therefore, PGKp1 and PGKp2 play redundant roles in carbon fixation and metabolism, whereas the molecular function of PGKc is more divergent. Our study demonstrated the functional conservation and divergence of glycolytic enzymes.

DNA synthesis pattern, proteome, and ABA and GA signalling in developing seeds of Norway maple (Acer platanoides)

Mature seeds of Norway maple exhibit desiccation tolerance and deep physiological dormancy. Flow cytometry, proteomics, and immunodetection have been combined to investigate seed development of this species. DNA content analysis revealed that cell cycle/endoreduplication activity differs between seed organs and developmental stages. In the embryo axis, the proportion of the nuclei with the highest DNA content (4C) increases at the beginning of maturation (17 weeks after flowering; WAF), and then is stable until the end of maturation, to increase again after drying. In cotyledons, during maturation endopolyploid nuclei (8C) occur and the intensity of endoreduplication increases up to 21 WAF, and then is stable until development is completed. In dry mature seeds, the proportion of 4C nuclei is high, and reaches 36% in the embryo axis and 52% in cotyledons. Proteomic studies revealed that energy and carbon metabolism, fatty acid biosynthesis, storage and antioxidant proteins are associated with seed development. Study of the ABI5 protein, a transcription factor involved in ABA signalling, and the RGL2 protein, a repressor of the GA signalling indicates that the highest accumulation of these proteins occurs in fully-matured and dried seeds. It is suggested that this increase in accumulation can be associated with completion of maturation, mainly with desiccation and dormancy acquisition.

Grain width 2 (GW2) and its interacting proteins regulate seed development in rice (Oryza sativa L.)

BACKGROUND

Seed size has been extensively studied in crop plants, as it determines crop yield. However, the mechanism of seed development remains elusive. In this study, we explored the mechanism of seed development in rice (Oryza sativa L.), and identified proteins affecting seed size.

RESULTS

Proteomic analysis showed that glyceraldehyde 3-phosphate dehydrogenase, chitinase 14 (CHT14), and phosphoglycerate kinase (PGK) accumulated to high levels in the seeds of the natural japonica rice mutant Oochikara, which carries a loss-of-function mutation in the grain width 2 (GW2) gene; GW2 encodes a RING-type E3 ubiquitin ligase. In vitro pull-down and ubiquitination assays showed that CHT14 and PGK directly interacted with GW2 but were not ubiquitinated by GW2. Immunoblot analysis revealed that protein disulfide isomerase-like 1-1 accumulated to high levels in young developing seeds of the gw2 mutant compared with the wild type. Histochemical β-glucuronidase staining showed strong expression of GW2 in leaf and root tissues but weak expression in leaf sheaths and internodes. In addition, transformation of the green fluorescent protein (GFP) gene under the control of the GW2 promoter in rice revealed GFP expression in the aleurone layer of seeds.

CONCLUSIONS

Collectively, these results suggest that GW2 regulates seed size through direct interactions with proteins involved in carbohydrate metabolism by modulating their activity or stability and controlling disulfide bond formation in various proteins during seed development. Additionally, GW2 participates in vegetative as well as reproductive growth, and protects the seed from pathogen attack.

Regulation of the Calvin-Benson-Bassham cycle in the enigmatic diatoms: biochemical and evolutionary variations on an original theme

In Plantae, the Calvin-Benson-Bassham (CBB) cycle is highly regulated and most of its enzymes have been thoroughly studied. Since diatoms arose as a result of secondary endosymbiosis with one or more Plantae ancestors, their precise evolutionary history is enigmatic and complex resulting in biochemical variations on the original CBB cycle theme. The Rubisco Michaelis constant for CO2 is higher in diatoms than land plants and the nuclear-encoded Rubisco activase in Plantae is replaced by an analogous chloroplast-encoded CbbX (Calvin-Benson-Bassham protein X) in diatoms. In the CBB cycle reduction phase, phosphoglycerate kinase in diatoms is redox-regulated and similar to that in red algae; however, glyceraldehyde phosphate dehydrogenase (GAPDH) is not redox-regulated, unlike in Plantae. The phosphoribulokinase (PRK)-GAPDH-CP12 complex found in many photosynthetic organisms has not yet been found in diatoms, but a ferredoxin-NADP reductase (FNR)-GAPDH-CP12 complex has been found in one species. In the CBB cycle regeneration phase, sedoheptulose 1,7-bisphosphatase and PRK are not redox-regulated in diatoms, unlike in Plantae. Regulation at the transcriptional level seems to be important in diatoms. CBB cycle enzyme properties appear to be variable among diatoms, but this view relies on results from a few model species: a greater range of diatoms need to be studied to test this.This article is part of the themed issue 'The peculiar carbon metabolism in diatoms'.

Metabolic Plasticity and Inter-Compartmental Interactions in Rice Metabolism: An Analysis from Reaction Deletion Study

More than 20% of the total caloric intake of human population comes from rice. The expression of rice genes and hence, the concentration of enzymatic proteins might vary due to several biotic and abiotic stresses. It in turn, can influence the overall metabolism and survivability of rice plant. Thus, understanding the rice cellular metabolism, its plasticity and potential readjustments under different perturbations can help rice biotechnologists to design efficient rice cultivars. Here, using the flux balance analysis (FBA) method, with the help of in-silico reaction deletion strategy, we study the metabolic plasticity of genome-scale metabolic model of rice leaf. A set of 131 reactions, essential for the production of primary biomass precursors is identified; deletion of any of them can inhibit the overall biomass production. Usability Index (I_U) for the rest of the reactions are estimated and based on this parameter, they are classified into three categories—maximally-favourable, quasi-favourable and unfavourable for the primary biomass production. The lower value of 1 − I_U of a reaction suggests that the cell cannot easily bypass it for biomass production. While some of the alternative paths are energetically equally efficient, others demand for higher photon. The variations in (i) ATP/NADPH ratio, (ii) exchange of metabolites through chloroplastic transporters and (iii) total biomass production are also presented here. Mutual metabolic dependencies of different cellular compartments are also demonstrated.

Endosymbiotic gene transfer from prokaryotic pangenomes: Inherited chimerism in eukaryotes

Endosymbiotic theory in eukaryotic-cell evolution rests upon a foundation of three cornerstone partners--the plastid (a cyanobacterium), the mitochondrion (a proteobacterium), and its host (an archaeon)--and carries a corollary that, over time, the majority of genes once present in the organelle genomes were relinquished to the chromosomes of the host (endosymbiotic gene transfer). However, notwithstanding eukaryote-specific gene inventions, single-gene phylogenies have never traced eukaryotic genes to three single prokaryotic sources, an issue that hinges crucially upon factors influencing phylogenetic inference. In the age of genomes, single-gene trees, once used to test the predictions of endosymbiotic theory, now spawn new theories that stand to eventually replace endosymbiotic theory with descriptive, gene tree-based variants featuring supernumerary symbionts: prokaryotic partners distinct from the cornerstone trio and whose existence is inferred solely from single-gene trees. We reason that the endosymbiotic ancestors of mitochondria and chloroplasts brought into the eukaryotic--and plant and algal--lineage a genome-sized sample of genes from the proteobacterial and cyanobacterial pangenomes of their respective day and that, even if molecular phylogeny were artifact-free, sampling prokaryotic pangenomes through endosymbiotic gene transfer would lead to inherited chimerism. Recombination in prokaryotes (transduction, conjugation, transformation) differs from recombination in eukaryotes (sex). Prokaryotic recombination leads to pangenomes, and eukaryotic recombination leads to vertical inheritance. Viewed from the perspective of endosymbiotic theory, the critical transition at the eukaryote origin that allowed escape from Muller's ratchet--the origin of eukaryotic recombination, or sex--might have required surprisingly little evolutionary innovation.

Photosynthetic gene expression in higher plants

Within the chloroplasts of higher plants and algae, photosynthesis converts light into biological energy, fueling the assimilation of atmospheric carbon dioxide into biologically useful molecules. Two major steps, photosynthetic electron transport and the Calvin-Benson cycle, require many gene products encoded from chloroplast as well as nuclear genomes. The expression of genes in both cellular compartments is highly dynamic and influenced by a diverse range of factors. Light is the primary environmental determinant of photosynthetic gene expression. Working through photoreceptors such as phytochrome, light regulates photosynthetic genes at transcriptional and posttranscriptional levels. Other processes that affect photosynthetic gene expression include photosynthetic activity, development, and biotic and abiotic stress. Anterograde (from nucleus to chloroplast) and retrograde (from chloroplast to nucleus) signaling insures the highly coordinated expression of the many photosynthetic genes between these different compartments. Anterograde signaling incorporates nuclear-encoded transcriptional and posttranscriptional regulators, such as sigma factors and RNA-binding proteins, respectively. Retrograde signaling utilizes photosynthetic processes such as photosynthetic electron transport and redox signaling to influence the expression of photosynthetic genes in the nucleus. The basic C3 photosynthetic pathway serves as the default form used by most of the plant species on earth. High temperature and water stress associated with arid environments have led to the development of specialized C4 and CAM photosynthesis, which evolved as modifications of the basic default expression program. The goal of this article is to explain and summarize the many gene expression and regulatory processes that work together to support photosynthetic function in plants.

How protein targeting to primary plastids via the endomembrane system could have evolved? A new hypothesis based on phylogenetic studies

Background

It is commonly assumed that a heterotrophic ancestor of the supergroup Archaeplastida/Plantae engulfed a cyanobacterium that was transformed into a primary plastid; however, it is still unclear how nuclear-encoded proteins initially were imported into the new organelle. Most proteins targeted to primary plastids carry a transit peptide and are transported post-translationally using Toc and Tic translocons. There are, however, several proteins with N-terminal signal peptides that are directed to higher plant plastids in vesicles derived from the endomembrane system (ES). The existence of these proteins inspired a hypothesis that all nuclear-encoded, plastid-targeted proteins initially carried signal peptides and were targeted to the ancestral primary plastid via the host ES.

Results

We present the first phylogenetic analyses of Arabidopsis thaliana α-carbonic anhydrase (CAH1), Oryza sativa nucleotide pyrophosphatase/phosphodiesterase (NPP1), and two O. sativa α-amylases (αAmy3, αAmy7), proteins that are directed to higher plant primary plastids via the ES. We also investigated protein disulfide isomerase (RB60) from the green alga Chlamydomonas reinhardtii because of its peculiar dual post- and co-translational targeting to both the plastid and ES. Our analyses show that these proteins all are of eukaryotic rather than cyanobacterial origin, and that their non-plastid homologs are equipped with signal peptides responsible for co-translational import into the host ES. Our results indicate that vesicular trafficking of proteins to primary plastids evolved long after the cyanobacterial endosymbiosis (possibly only in higher plants) to permit their glycosylation and/or transport to more than one cellular compartment.

Conclusions

The proteins we analyzed are not relics of ES-mediated protein targeting to the ancestral primary plastid. Available data indicate that Toc- and Tic-based translocation dominated protein import into primary plastids from the beginning. Only a handful of host proteins, which already were targeted through the ES, later were adapted to reach the plastid via the vesicular trafficking. They represent a derived class of higher plant plastid-targeted proteins with an unusual evolutionary history.

Reviewers

This article was reviewed by Prof. William Martin, Dr. Philippe Deschamps (nominated by Dr. Purificacion Lopez-Garcia) and Dr Simonetta Gribaldo.

An evolutionary network of genes present in the eukaryote common ancestor polls genomes on eukaryotic and mitochondrial origin

To test the predictions of competing and mutually exclusive hypotheses for the origin of eukaryotes, we identified from a sample of 27 sequenced eukaryotic and 994 sequenced prokaryotic genomes 571 genes that were present in the eukaryote common ancestor and that have homologues among eubacterial and archaebacterial genomes. Maximum-likelihood trees identified the prokaryotic genomes that most frequently contained genes branching as the sister to the eukaryotic nuclear homologues. Among the archaebacteria, euryarchaeote genomes most frequently harbored the sister to the eukaryotic nuclear gene, whereas among eubacteria, the α-proteobacteria were most frequently represented within the sister group. Only 3 genes out of 571 gave a 3-domain tree. Homologues from α-proteobacterial genomes that branched as the sister to nuclear genes were found more frequently in genomes of facultatively anaerobic members of the rhiozobiales and rhodospirilliales than in obligate intracellular ricketttsial parasites. Following α-proteobacteria, the most frequent eubacterial sister lineages were γ-proteobacteria, δ-proteobacteria, and firmicutes, which were also the prokaryote genomes least frequently found as monophyletic groups in our trees. Although all 22 higher prokaryotic taxa sampled (crenarchaeotes, γ-proteobacteria, spirochaetes, chlamydias, etc.) harbor genes that branch as the sister to homologues present in the eukaryotic common ancestor, that is not evidence of 22 different prokaryotic cells participating at eukaryote origins because prokaryotic “lineages” have laterally acquired genes for more than 1.5 billion years since eukaryote origins. The data underscore the archaebacterial (host) nature of the eukaryotic informational genes and the eubacterial (mitochondrial) nature of eukaryotic energy metabolism. The network linking genes of the eukaryote ancestor to contemporary homologues distributed across prokaryotic genomes elucidates eukaryote gene origins in a dialect cognizant of gene transfer in nature.

The eukaryotic tree of life: endosymbiosis takes its TOL

Resolving the structure of the eukaryotic tree of life remains one of the most important and challenging tasks facing biologists. The notion of six eukaryotic 'supergroups' has recently gained some acceptance, and several papers in 2007 suggest that resolution of higher taxonomic levels is possible. However, in organisms that acquired photosynthesis via secondary (i.e. eukaryote-eukaryote) endosymbiosis, the host nuclear genome is a mosaic of genes derived from two (or more) nuclei, a fact that is often overlooked in studies attempting to reconstruct the deep evolutionary history of eukaryotes. Accurate identification of gene transfers and replacements involving eukaryotic donor and recipient genomes represents a potentially formidable challenge for the phylogenomics community as more protist genomes are sequenced and concatenated data sets grow.

Horizontal Gene Transfer in Aminoacyl-tRNA Synthetases Including Leucine-Specific Subtypes

Aminoacyl-tRNA synthetases catalyze a fundamental reaction for the flow of genetic information from RNA to protein. Their presence in all organisms known today highlights their important role in the early evolution of life. We investigated the evolutionary history of aminoacyl-tRNA synthetases on the basis of sequence data from more than 200 Archaea, Bacteria, and Eukaryota. Phylogenetic profiles are in agreement with previous observations that many genes for aminoacyl-tRNA synthetases were transferred horizontally between species from all domains of life. We extended these findings by a detailed analysis of the history of leucyl-tRNA synthetases. Thereby, we identified a previously undetected case of horizontal gene transfer from Bacteria to Archaea based on phylogenetic profiles, trees, and networks. This means that, finally, the last subfamily of aminoacyl-tRNA synthetases has lost its exceptional position as the sole subfamily that is devoid of horizontal gene transfer. Furthermore, the leucyl-tRNA synthetase phylogenetic tree suggests a dichotomy of the archaeal/eukaryotic-cytosolic and bacterial/eukaryotic-mitochondrial proteins. We argue that the traditional division of life into Prokaryota (non-chimeric) and Eukaryota (chimeric) is favorable compared to Woese's trichotomy into Archaea/Bacteria/Eukaryota.

Jumping Genes and Shrinking Genomes – Probing the Evolution of Eukaryotic Photosynthesis with Genomics

The advent of comparative genomics has revolutionized the study of the origin and evolution of eukaryotes and their organelles. Genomic analysis has revealed that the endosymbiosis that gave rise to plastids--the light-harvesting apparatus of photosynthetic eukaryotes--had a profound impact on the genetic composition of the host, far beyond the contribution of cyanobacterial genes for plastid-specific functions. Here I discuss recent advances in our appreciation of the mosaic nature of the eukaryotic nuclear genome, and the ongoing role endosymbiosis plays in shaping its content.

Distribution and phylogenies of enzymes of the Embden-Meyerhof-Parnas pathway from archaea and hyperthermophilic bacteria support a gluconeogenic origin of metabolism

Enzymes of the gluconeogenic/glycolytic pathway (the Embden-Meyerhof-Parnas (EMP) pathway), the reductive tricarboxylic acid cycle, the reductive pentose phosphate cycle and the Entner-Doudoroff pathway are widely distributed and are often considered to be central to the origins of metabolism. In particular, several enzymes of the lower portion of the EMP pathway (the so-called trunk pathway), including triosephosphate isomerase (TPI; EC 5.3.1.1), glyceraldehyde-3-phosphate dehydrogenase (GAPDH; EC 1.2.1.12/13), phosphoglycerate kinase (PGK; EC 2.7.2.3) and enolase (EC 4.2.1.11), are extremely well conserved and universally distributed among the three domains of life. In this paper, the distribution of enzymes of gluconeogenesis/glycolysis in hyperthermophiles--microorganisms that many believe represent the least evolved organisms on the planet--is reviewed. In addition, the phylogenies of the trunk pathway enzymes (TPIs, GAPDHs, PGKs and enolases) are examined. The enzymes catalyzing each of the six-carbon transformations in the upper portion of the EMP pathway, with the possible exception of aldolase, are all derived from multiple gene sequence families. In contrast, single sequence families can account for the archaeal and hyperthermophilic bacterial enzyme activities of the lower portion of the EMP pathway. The universal distribution of the trunk pathway enzymes, in combination with their phylogenies, supports the notion that the EMP pathway evolved in the direction of gluconeogenesis, i.e., from the bottom up.

Nucleus-encoded genes for plastid-targeted proteins in Helicosporidium: functional diversity of a cryptic plastid in a parasitic alga

Plastids are the organelles of plants and algae that house photosynthesis and many other biochemical pathways. Plastids contain a small genome, but most of their proteins are encoded in the nucleus and posttranslationally targeted to the organelle. When plants and algae lose photosynthesis, they virtually always retain a highly reduced "cryptic" plastid. Cryptic plastids are known to exist in many organisms, although their metabolic functions are seldom understood. The best-studied example of a cryptic plastid is from the intracellular malaria parasite, Plasmodium, which has retained a plastid for the biosynthesis of fatty acids, isoprenoids, and heme by the use of plastid-targeted enzymes. To study a completely independent transformation of a photosynthetic plastid to a cryptic plastid in another alga-turned-parasite, we conducted an expressed sequence tag (EST) survey of Helicosporidium. This parasite has recently been recognized as a highly derived green alga. Based on phylogenetic relationships to other plastid homologues and the presence of N-terminal transit peptides, we have identified 20 putatively plastid-targeted enzymes that are involved in a wide variety of metabolic pathways. Overall, the metabolic diversity of the Helicosporidium cryptic plastid exceeds that of the Plasmodium plastid, as it includes representatives of most of the pathways known to operate in the Plasmodium plastid as well as many others. In particular, several amino acid biosynthetic pathways have been retained, including the leucine biosynthesis pathway, which was only recently recognized in plant plastids. These two parasites represent different evolutionary trajectories in plastid metabolic adaptation.

Chloroplast phosphoglycerate kinase from Euglena gracilis

Two chloroplast phosphoglycerate kinase isoforms from the photosynthetic flagellate Euglena gracilis were purified to homogeneity, partially sequenced, and subsequently cDNAs encoding phosphoglycerate kinase isoenzymes from both the chloroplast and cytosol of E. gracilis were cloned and sequenced. Chloroplast phosphoglycerate kinase, a monomeric enzyme, was encoded as a polyprotein precursor of at least four mature subunits that were separated by conserved tetrapeptides. In a Neighbor-Net analysis of sequence similarity with homologues from numerous prokaryotes and eukaryotes, cytosolic phosphoglycerate kinase of E. gracilis showed the highest similarity to cytosolic and glycosomal homologues from the Kinetoplastida. The chloroplast isoenzyme of E. gracilis did not show a close relationship to sequences from other photosynthetic organisms but was most closely related to cytosolic homologues from animals and fungi.

Temporal changes in phosphoglycerate kinase coding sequences: a quantitative measure

The ratio of the average of the square of the number of the nucleotides to that of the random sequence of the same strand bias is proposed as a quantitative measure of evolution in some coding DNA sequences. Applying this measure to the phosphoglycerate kinase gene we observe a monotonic rise of the ratio with evolution. We present an interpretation of this data on some bacteria.

GAPDH gene diversity in spirochetes: a paradigm for genetic promiscuity

In this study we have determined gap sequences from nine different spirochetes. Phylogenetic analyses of these sequences in the context of all other available eubacterial and a selection of eukaryotic Gap sequences demonstrated that the eubacterial glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene diversity encompasses at least five highly distinct gene families. Within these gene families, spirochetes show an extreme degree of sequence divergence that is probably the result of several lateral gene transfer events between spirochetes and other eubacterial phyla, and early gene duplications in the eubacterial ancestor. A Gap1 sequence from the syphilis spirochete Treponema pallidum has recently been shown to be closely related to GapC sequences from Euglenozoa. Here we demonstrate that several other spirochetal species are part of this cluster, supporting the conclusion that an interkingdom gene transfer from spirochetes to Euglenozoa must have occurred. Furthermore, we provide evidence that the GAPDH genes present in the protists Parabasalia may also be of spirochetal descent.

An overview of endosymbiotic models for the origins of eukaryotes, their ATP-producing organelles (mitochondria and hydrogenosomes), and their heterotrophic lifestyle

The evolutionary processes underlying the differentness of prokaryotic and eukaryotic cells and the origin of the latter's organelles are still poorly understood. For about 100 years, the principle of endosymbiosis has figured into thoughts as to how these processes might have occurred. A number of models that have been discussed in the literature and that are designed to explain this difference are summarized. The evolutionary histories of the enzymes of anaerobic energy metabolism (oxygen-independent ATP synthesis) in the three basic types of heterotrophic eukaryotes those that lack organelles of ATP synthesis, those that possess mitochondria and those that possess hydrogenosomes--play an important role in this issue. Traditional endosymbiotic models generally do not address the origin of the heterotrophic lifestyle and anaerobic energy metabolism in eukaryotes. Rather they take it as a given, a direct inheritance from the host that acquired mitochondria. Traditional models are contrasted to an alternative endosymbiotic model (the hydrogen hypothesis), which addresses the origin of heterotrophy and the origin of compartmentalized energy metabolism in eukaryotes.

Purification and cloning of chloroplast 6-phosphogluconate dehydrogenase from spinach. Cyanobacterial genes for chloroplast and cytosolic isoenzymes encoded in eukaryotic chromosomes

Previous attempts to purify chloroplast 6-phosphogluconate dehydrogenase (cp6PGDH), a key enzyme of the oxidative pentose phosphate pathway, have been unsuccessful due to rapid activity loss. An efficient purification protocol was developed and the enzyme from spinach leaves was purified 1000-fold to apparent homogeneity with a specific activity of 60 U.mg-1. The enzyme is a homodimer with subunits of 50 kDa. Antibodies raised against the purified cp6PGDH detected a 53-kDa protein from a crude extract, indicating alterations during purification. Purified cp6PGDH was microsequenced and the corresponding spinach cDNA was cloned using PCR techniques and degenerate primers. The cDNA for cytosolic 6PGDH from spinach was cloned for comparison. Phylogenetic analysis in the context of available homologues from eukaryotes and eubacteria revealed that animal and fungal cytosolic 6PGDH sequences are more similar to their homologues from gamma-proteobacteria, whereas plant 6PGDH is more similar to its cyanobacterial homologues. The ancestral gene for higher plant 6PGDH was acquired from the antecedent of plastids through endosymbiosis and gene transfer to the nucleus. A subsequent gene duplication gave rise to higher plant cytosolic 6PGDH, which assumed the function of its pre-existing cytosolic homologue through endosymbiotic gene replacement. The protein phylogeny of both 6PGDH and of the first enzyme of the oxidative pentose phosphate pathway, glucose-6-phosphate dehydrogenase, indicate a surprisingly close relationship between the plant and Trypanosoma brucei lineages, suggesting that T. brucei (a relative of Euglena gracilis) may be secondarily nonphotosynthetic.

Enolase from Trypanosoma brucei, from the amitochondriate protist Mastigamoeba balamuthi, and from the chloroplast and cytosol of Euglena gracilis: pieces in the evolutionary puzzle of the eukaryotic glycolytic pathway

Genomic or cDNA clones for the glycolytic enzyme enolase were isolated from the amitochondriate pelobiont Mastigamoeba balamuthi, from the kinetoplastid Trypanosoma brucei, and from the euglenid Euglena gracilis. Clones for the cytosolic enzyme were found in all three organisms, whereas Euglena was found to also express mRNA for a second isoenzyme that possesses a putative N-terminal plastid-targeting peptide and is probably targeted to the chloroplast. Database searching revealed that Arabidopsis also possesses a second enolase gene that encodes an N-terminal extension and is likely targeted to the chloroplast. A phylogeny of enolase amino acid sequences from 6 archaebacteria, 24 eubacteria, and 32 eukaryotes showed that the Mastigamoeba enolase tended to branch with its homologs from Trypanosoma and from the amitochondriate protist Entamoeba histolytica. The compartment-specific isoenzymes in Euglena arose through a gene duplication independent of that which gave rise to the compartment-specific isoenzymes in Arabidopsis, as evidenced by the finding that the Euglena enolases are more similar to the homolog from the eubacterium Treponema pallidum than they are to homologs from any other organism sampled. In marked contrast to all other glycolytic enzymes studied to date, enolases from all eukaryotes surveyed here (except Euglena) are not markedly more similar to eubacterial than to archaebacterial homologs. An intriguing indel shared by enolase from eukaryotes, from the archaebacterium Methanococcus jannaschii, and from the eubacterium Campylobacter jejuni maps to the surface of the three-dimensional structure of the enzyme and appears to have occurred at the same position in parallel in independent lineages.

Pyruvate kinase of the hyperthermophilic crenarchaeote Thermoproteus tenax: physiological role and phylogenetic aspects

Pyruvate kinase (PK; EC 2.7.1.40) of Thermoproteus tenax was purified to homogeneity, and its coding gene was cloned and expressed in Escherichia coli. It represents a homomeric tetramer with a molecular mass of 49 kDa per subunit. PK exhibits positive binding cooperativity with respect to phosphoenolpyruvate and metal ions such as Mg(2+) and Mn(2+). Heterotropic effects, as commonly found for PKs from bacterial and eucaryal sources, could not be detected. The enzyme does not depend on K(+) ions. Heterotrophically grown cells exhibit specific activity of PK four times higher than autotrophically grown cells. Since the mRNA level of the PK coding gene is also accordingly higher in heterotrophic cells, we conclude that the PK activity is adjusted to growth conditions mainly on the transcript level. The enzymic properties of the PK and the regulation of its expression are discussed with respect to the physiological framework given by the T. tenax-specific variant of the Embden-Meyerhof-Parnas pathway. T. tenax PK shows moderate overall sequence similarity (25 to 40% identity) to its bacterial and eucaryal pendants. Phylogenetic analyses of the known PK sequences result in a dichotomic tree topology that divides the enzymes into two major PK clusters, probably diverged by an early gene duplication event. The phylogenetic divergence is paralleled by a striking phenotypic differentiation of PKs: PKs of cluster I, which occur in eucaryal cytoplasm, some gamma proteobacteria, and low-GC gram-positive bacteria, are only active in the presence of fructose-1,6-bisphosphate or other phosphorylated sugars, whereas PKs of cluster II, found in various bacterial phyla, plastids, and in Archaea, show activity without effectors but are commonly regulated by the energy charge of the cell.

The crystal structure of d-glyceraldehyde-3-phosphate dehydrogenase from the hyperthermophilic archaeon Methanothermus fervidus in the presence of NADP(+) at 2.1 A resolution

The crystal structure of the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from the archaeon Methanothermus fervidus has been solved in the holo form at 2.1 A resolution by molecular replacement. Unlike bacterial and eukaryotic homologous enzymes which are strictly NAD(+)-dependent, GAPDH from this organism exhibits a dual-cofactor specificity, with a marked preference for NADP(+) over NAD(+). The present structure is the first archaeal GAPDH crystallized with NADP(+). GAPDH from M. fervidus adopts a homotetrameric quaternary structure which is topologically similar to that observed for its bacterial and eukaryotic counterparts. Within the cofactor-binding site, the positively charged side-chain of Lys33 decisively contributes to NADP(+) recognition through a tight electrostatic interaction with the adenosine 2'-phosphate group. Like other GAPDHs, GAPDH from archaeal sources binds the nicotinamide moiety of NADP(+) in a syn conformation with respect to the adjacent ribose and so belongs to the B-stereospecific class of oxidoreductases. Stabilization of the syn conformation is principally achieved through hydrogen bonding of the carboxamide group with the side-chain of Asp171, a structural feature clearly different from what is observed in all presently known GAPDHs from bacteria and eukaryotes. Within the catalytic site, the reported crystal structure definitively confirms the essential role previously assigned to Cys140 by site-directed mutagenesis studies. In conjunction with new mutation results reported in this paper, inspection of the crystal structure gives reliable evidence for the direct implication of the side-chain of His219 in the catalytic mechanism. M. fervidus grows optimally at 84 degrees C with a maximal growth temperature of 97 degrees C. The paper includes a detailed comparison of the present structure with four other homologous enzymes extracted from mesophilic as well as thermophilic organisms. Among the various phenomena related to protein thermostabilization, reinforcement of electrostatic and hydrophobic interactions as well as a more efficient molecular packing appear to be essentially promoted by the occurrence of two additional alpha-helices in the archaeal GAPDHs. The first one, named alpha4, is located in the catalytic domain and participates in the enzyme architecture at the quaternary structural level. The second one, named alphaJ, occurs at the C terminus and contributes to the molecular packing within each monomer by filling a peripherical pocket in the tetrameric assembly.

Functional characterization of the phosphorylating D-glyceraldehyde 3-phosphate dehydrogenase from the archaeon Methanothermus fervidus by comparative molecular modelling and site-directed mutagenesis

Phosphorylating archaeal D-glyceraldehyde 3-phosphate dehydrogenases (GraP-DHs) share only 15-20% identity with their glycolytic bacterial and eukaryotic counterparts. Unlike the latter which are NAD-specific, archaeal GraP-DHs exhibit a dual-cofactor specificity with a marked preference for NADP. In the present study, we have constructed a three-dimensional model of the Methanothermus fervidus GraP-DH based upon the X-ray structures of the Bacillus stearothermophilus and Escherichia coli GraP-DHs. The overall structure of the archaeal enzyme is globally similar to homology modelling-derived structures, in particular for the cofactor binding domain, which might adopt a classical Rossmann fold. M. fervidus GraP-DH can be considered as a dimer of dimers which exhibits negative and positive cooperativity in binding the coenzymes NAD and NADP, respectively. As expected, the differences between the model and the templates are located mainly within the loops. Based on the predictions derived from molecular modelling, site-directed mutagenesis was performed to characterize better the cofactor binding pocket and the catalytic domain. The Lys32Ala, Lys32Glu and Lys32Asp mutants led to a drastic increase in the Km value for NADP (i.e. 165-, 500- and 1000-fold, respectively), thus demonstrating that the invariant Lys32 residue is one of the most important determinants favouring the adenosine 2'-PO42- binding of NADP. The involvement of the side chain of Asn281, which was postulated to play a role equivalent to that of the Asn313 of bacterial and eukaryotic GraP-DHs in fixing the position of the nicotinamide ring in a syn orientation [Fabry, S. & Hensel, R. (1988) Gene 64, 189-197], was ruled out. Most of the amino acids involved in catalysis and in substrate recognition in bacterial and eukaryotic GraP-DHs are not conserved in the archaeal enzyme except for the essential Cys149. Inspection of our model suggests that side chains of invariant residues Asn150, Arg176, Arg177 and His210 are located in or near the active site pocket. The Arg177Asn mutation induced strong allosteric properties with the Pi, indicating that this residue should be located near to the intersubunit interfaces. The Arg176Asn mutation led to a 10-fold decrease in the kcat, a 35-fold increase in the Km value for D-glyceraldehyde 3-phosphate and a 1000-fold decrease in the acylation rate. These results strongly suggest that Arg176 is involved in the Ps site. The His210Asn mutation increased the pKapp of the catalytic Cys149 from 6.3 to 7.6, although no Cys-/His+ ion pair was detectable [Talfournier, F., Colloc'h, N., Mornon, J.P. & Branlant, G. (1998) Eur. J. Biochem. 252, 447-457]. No other invariant amino acid which can play a role as a base catalyst to favour the hydride transfer is located in the active site. The fact that the efficiency of phosphorolysis is 1000-fold lower when compared to the B. stearothermophilus GraP-DH suggests significant differences in the nature of the Pi site. Despite these differences, it is likely that the archaeal GraP-DHs and their bacterial and eukaryotic counterparts have evolved from a common ancestor.

RPE, a plant gene involved in early developmental steps of nematode feeding cells

Sedentary plant-parasitic nematodes are able to induce the redifferentiation of root cells into multinucleate nematode feeding sites (NFSs). We have isolated by promoter trapping an Arabidopsis thaliana gene that is essential for the early steps of NFS formation induced by the root-knot nematode Meloidogyne incognita. Its pattern of expression is similar to that of key regulators of the cell cycle, but it is not observed with the cyst nematode. Later in NFS development, this gene is induced by both root-knot and cyst nematodes. It encodes a protein similar to the D-ribulose-5-phosphate 3-epimerase (RPE) (EC 5.1.3.1), a key enzyme in the reductive Calvin cycle and the oxidative pentose phosphate pathway (OPPP). Quantitative RT-PCR showed the accumulation of RPE transcripts in potato, as in Arabidopsis NFS. Homozygous rpe plants have a germination mutant phenotype that can be rescued in dwarf plants on sucrose-supplemented medium. During root development, this gene is expressed in the meristems and initiation sites of lateral roots. These results suggest that the genetic control of NFSs and the first stages of meristem formation share common steps and confirms the previous cytological observations which indicate that root cells undergo metabolic reprogramming when they turn into NFSs.

NAD+-dependent glyceraldehyde-3-phosphate dehydrogenase from Thermoproteus tenax. The first identified archaeal member of the aldehyde dehydrogenase superfamily is a glycolytic enzyme with unusual regulatory properties

The hyperthermophilic archaeum Thermoproteus tenax possesses two glyceraldehyde-3-phosphate dehydrogenases differing in cosubstrate specificity and phosphate dependence of the catalyzed reaction. NAD+-dependent glyceraldehyde-3-phosphate dehydrogenase catalyzes the phosphate-independent irreversible oxidation of D-glyceraldehyde 3-phosphate to 3-phosphoglycerate. The coding gene was cloned, sequenced, and expressed in Escherichia coli. Sequence comparisons showed no similarity to phosphorylating glyceraldehyde-3-phosphate dehydrogenases but revealed a relationship to aldehyde dehydrogenases, with the highest similarity to the subgroup of nonphosphorylating glyceraldehyde-3-phosphate dehydrogenases. The activity of the enzyme is affected by a series of metabolites. All effectors tested influence the affinity of the enzyme for its cosubstrate NAD+. Whereas NADP(H), NADH, and ATP reduce the affinity for the cosubstrate, AMP, ADP, glucose 1-phosphate, and fructose 6-phosphate increase the affinity for NAD+. Additionally, most of the effectors investigated induce cooperativity of NAD+ binding. The irreversible catabolic oxidation of glyceraldehyde 3-phosphate, the control of the enzyme by energy charge of the cell, and the regulation by intermediates of glycolysis and glucan degradation identify the NAD+-dependent glyceraldehyde-3-phosphate dehydrogenase as an integral constituent of glycolysis in T. tenax. Its regulatory properties substitute for those lacking in the reversible nonregulated pyrophosphate-dependent phosphofructokinase in this variant of the Embden-Meyerhof-Parnas pathway.

The hydrogen hypothesis for the first eukaryote

A new hypothesis for the origin of eukaryotic cells is proposed, based on the comparative biochemistry of energy metabolism. Eukaryotes are suggested to have arisen through symbiotic association of an anaerobic, strictly hydrogen-dependent, strictly autotrophic archaebacterium (the host) with a eubacterium (the symbiont) that was able to respire, but generated molecular hydrogen as a waste product of anaerobic heterotrophic metabolism. The host's dependence upon molecular hydrogen produced by the symbiont is put forward as the selective principle that forged the common ancestor of eukaryotic cells.

Archaea and the prokaryote-to-eukaryote transition

Since the late 1970s, determining the phylogenetic relationships among the contemporary domains of life, the Archaea (archaebacteria), Bacteria (eubacteria), and Eucarya (eukaryotes), has been central to the study of early cellular evolution. The two salient issues surrounding the universal tree of life are whether all three domains are monophyletic (i.e., all equivalent in taxanomic rank) and where the root of the universal tree lies. Evaluation of the status of the Archaea has become key to answering these questions. This review considers our cumulative knowledge about the Archaea in relationship to the Bacteria and Eucarya. Particular attention is paid to the recent use of molecular phylogenetic approaches to reconstructing the tree of life. In this regard, the phylogenetic analyses of more than 60 proteins are reviewed and presented in the context of their participation in major biochemical pathways. Although many gene trees are incongruent, the majority do suggest a sisterhood between Archaea and Eucarya. Altering this general pattern of gene evolution are two kinds of potential interdomain gene transferrals. One horizontal gene exchange might have involved the gram-positive Bacteria and the Archaea, while the other might have occurred between proteobacteria and eukaryotes and might have been mediated by endosymbiosis.

The tigA gene is a transcriptional fusion of glycolytic genes encoding triose-phosphate isomerase and glyceraldehyde-3-phosphate dehydrogenase in oomycota

Genes encoding triose-phosphate isomerase (TPI) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) are fused and form a single transcriptional unit (tigA) in Phytophthora species, members of the order Pythiales in the phylum Oomycota. This is the first demonstration of glycolytic gene fusion in eukaryotes and the first case of a TPI-GAPDH fusion in any organism. The tigA gene from Phytophthora infestans has a typical Oomycota transcriptional start point consensus sequence and, in common with most Phytophthora genes, has no introns. Furthermore, Southern and PCR analyses suggest that the same organization exists in other closely related genera, such as Pythium, from the same order (Oomycota), as well as more distantly related genera, Saprolegnia and Achlya, in the order Saprolegniales. Evidence is provided that in P. infestans, there is at least one other discrete copy of a GAPDH-encoding gene but not of a TPI-encoding gene. Finally, a phylogenetic analysis of TPI does not place Phytophthora within the assemblage of crown eukaryotes and suggests TPI may not be particularly useful for resolving relationships among major eukaryotic groups.

Choline monooxygenase, an unusual iron-sulfur enzyme catalyzing the first step of glycine betaine synthesis in plants: prosthetic group characterization and cDNA cloning

Plants synthesize the osmoprotectant glycine betaine via the route choline --> betaine aldehyde --> glycine betaine. In spinach, the first step is catalyzed by choline monooxygenase (CMO), a ferredoxin-dependent stromal enzyme that has been hypothesized to be an oligomer of identical subunits and to be an Fe-S protein. Analysis by HPLC and matrix-assisted laser desorption ionization MS confirmed that native CMO contains only one type of subunit (Mr 42,864). Determination of acid-labile sulfur and nonheme iron demonstrated that there is one [2Fe-2S] cluster per subunit, and EPR spectral data indicated that this cluster is of the Rieske type--i.e., coordinated by two Cys and two His ligands. A full-length CMO cDNA (1,622 bp) was cloned from spinach using a probe generated by PCR amplification for which the primers were based on internal peptide sequences. The ORF encoded a 440-amino acid polypeptide that included a 60-residue transit peptide. The deduced amino acid sequence included two Cys-His pairs spaced 16 residues apart, a motif characteristic of Rieske-type Fe-S proteins. Larger regions that included this motif also showed some sequence similarity (approximately 40%) to Rieske-type proteins, particularly bacterial oxygenases. Otherwise there was very little similarity between CMO and proteins from plants or other organisms. RNA and immunoblot analyses showed that the expression of CMO in leaves increased several-fold during salinization. We conclude that CMO is a stress-inducible representative of a new class of plant oxygenases.

Evidence that eukaryotic triosephosphate isomerase is of alpha-proteobacterial origin

We have cloned and sequenced genes for triosephosphate isomerase (TPI) from the gamma-proteobacterium Francisella tularensis, the green non-sulfur bacterium Chloroflexus aurantiacus, and the alpha-proteobacterium Rhizobium etli and used these in phylogenetic analysis with TPI sequences from other members of the Bacteria, Archaea, and Eukarya. These analyses show that eukaryotic TPI genes are most closely related to the homologue from the alpha-proteobacterium and most distantly related to archaebacterial homologues. This relationship suggests that the TPI genes present in modern eukaryotic genomes were derived from an alpha-proteobacterial genome (possibly that of the protomitochondrial endosymbiont) after the divergence of Archaea and Eukarya. Among these eukaryotic genes are some from deeply branching, amitochondrial eukaryotes (namely Giardia), which further suggests that this event took place quite early in eukaryotic evolution.

Evidence for the Heterolobosea from phylogenetic analysis of genes encoding glyceraldehyde-3-phosphate dehydrogenase

The phylogenetic relationships between major slime mould groups and the identification of their unicellular relatives has been a subject of controversy for many years. Traditionally, it has been assumed that two slime mould groups, the acrasids and the dictyostelids were related by virtue of their cellular slime mould habit; a view still endorsed by at least one current classification scheme. However, a decade ago, on the basis of detailed ultrastructural resemblances it was proposed that acrasids of the family Acrasidae were not relatives of other slime moulds but instead related to a group of mostly free-living unicellular amoebae, the Schizopyrenida. The class Heterolobosea was created to contain these organisms and has since figured in many discussions of protist evolution. We sought to test the validity of Heterolobosea by characterizing homologs of the highly conserved glycolytic enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH) from an acrasid, Acrasis rosea; a dictyostelid, Dictyostelium discoideum; and the schizopyrenid Naegleria andersoni. Phylogenetic analysis of these and other GAPDH sequences, using maximum parsimony, neighbour-joining distance and maximum likelihood methods strongly supports the Heterolobosea hypothesis and discredits the concept of a cellular slime mould grouping. Moreover, all of our analyses place Dictyostelium discoideum as a relatively recently originating lineage, most closely related to the Metazoa, similar to other recently published phylogenies of protein-coding genes. However, GAPDH phylogenies do not show robust branching orders for most of the relationships between major groups. We propose that several of the incongruencies observed between GAPDH and other molecular phylogenies are artifacts resulting from substitutional saturation of this enzyme.

Molecular characterization of transketolase (EC 2.2.1.1) active in the Calvin cycle of spinach chloroplasts

A cDNA encoding the Calvin cycle enzyme transketolase (TKL; EC 2.2.1.1) was isolated from Sorghum bicolor via subtractive differential hybridization, and used to isolate several full-length cDNA clones for this enzyme from spinach. Functional identity of the encoded mature subunit was shown by an 8.6-fold increase of TKL activity upon induction of Escherichia coli cells that overexpress the spinach TKL subunit under the control of the bacteriophage T7 promoter. Chloroplast localization of the cloned enzyme is shown by processing of the in vitro synthesized precursor upon uptake by isolated chloroplasts. Southern blot-analysis suggests that TKL is encoded by a single gene in the spinach genome. TKL proteins of both higher-plant chloroplasts and the cytosol of non-photosynthetic eukaryotes are found to be unexpectedly similar to eubacterial homologues, suggesting a possible eubacterial origin of these nuclear genes. Chloroplast TKL is the last of the demonstrably chloroplast-localized Calvin cycle enzymes to have been cloned and thus completes the isolation of gene probes for all enzymes of the pathway in higher plants.

Higher-plant chloroplast and cytosolic fructose-1,6-bisphosphatase isoenzymes: origins via duplication rather than prokaryote-eukaryote divergence

Full-size cDNAs encoding the precursors of chloroplast fructose-1,6-bisphosphatase (FBP), sedoheptulose-1,7-bisphosphatase (SBP), and the small subunit of Rubisco (RbcS) from spinach were cloned. These cDNAs complete the set of homologous probes for all nuclear-encoded enzymes of the Calvin cycle from spinach (Spinacia oleracea L.). FBP enzymes not only of higher plants but also of non-photosynthetic eukaryotes are found to be unexpectedly similar to eubacterial homologues, suggesting a eubacterial origin of these eukaryotic nuclear genes. Chloroplast and cytosolic FBP isoenzymes of higher plants arose through a gene duplication event which occurred early in eukaryotic evolution. Both FBP and SBP of higher plant chloroplasts have acquired substrate specificity, i.e. have undergone functional specialization since their divergence from bifunctional FBP/SBP enzymes of free-living eubacteria.

Microsequencing and cDNA cloning of the Calvin cycle/OPPP enzyme ribose-5-phosphate isomerase (EC 5.3.1.6) from spinach chloroplasts

Ribose-5-phosphate isomerase (RPI) catalyses the interconversion of ribose-5-phosphate and ribulose-5-phosphate in the reductive and oxidative pentose phosphate pathways in plants. RPI from spinach chloroplasts was purified and microsequenced. Via PCR with degenerate primers designed against microsequenced peptides, a hybridisation probe was obtained and used to isolate several cDNA clones which encode RPI. The nuclear-encoded 239 amino acid mature RPI subunit has a predicted size of 25.3 kDa and is translated as a cytosolic precursor possessing a 50 amino acid transit peptide. The processing site of the transit peptide was identified from protein sequence data. Spinach leaves possess only one type of homodimeric RPI enzyme which is localized in chloroplasts and is encoded by a single nuclear gene. Molecular characterization of RPI supports the view that a single amphibolic RPI enzyme functions in the oxidative and reductive pentose phosphate pathways of spinach plastids.

Cloning of the amphibolic Calvin cycle/OPPP enzyme D-ribulose-5-phosphate 3-epimerase (EC 5.1.3.1) from spinach chloroplasts: functional and evolutionary aspects

Exploiting the differential expression of genes for Calvin cycle enzymes in bundle-sheath and mesophyll cells of the C4 plant Sorghum bicolor L., we isolated via subtractive hybridization a molecular probe for the Calvin cycle enzyme D-ribulose-5-phosphate 3-epimerase (R5P3E)(EC 5.1.3.1), with the help of which several full-size cDNAs were isolated from spinach. Functional identity of the encoded mature subunit was shown by R5P3E activity found in affinity-purified glutatione S-transferase fusions expressed in Escherichia coli and by three-fold increase of R5P3E activity upon induction of E. coli overexpressing the spinach subunit under the control of the bacteriophage T7 promoter, demonstrating that we have cloned the first functional ribulose-5-phosphate 3-epimerase from any eukaryotic source. The chloroplast enzyme from spinach shares about 50% amino acid identity with its homologues from the Calvin cycle operons of the autotrophic purple bacteria Alcaligenes eutrophus and Rhodospirillum rubrum. A R5P3E-related eubacterial gene family was identified which arose through ancient duplications in prokaryotic chromosomes, three R5P3E-related genes of yet unknown function have persisted to the present within the E. coli genome. A gene phylogeny reveals that spinach R5P3E is more similar to eubacterial homologues than to the yeast sequence, suggesting a eubacterial origin for this plant nuclear gene.