Rampant horizontal gene transfer and phospho-donor change in the evolution of the phosphofructokinase

Genome-wide identification of the PFK gene family and their expression analysis in Quercus rubra

The glycolytic pathway involves phosphofructokinase (PFK), a rate-limiting enzyme that catalyzes the phosphorylation of fructose-6-phosphate. In plants, the two PFK members are ATP-dependent phosphofructokinase (PFK) and pyrophosphate-fructose-6-phosphate phosphotransferase (PFP). However, the functions of the PFK family members in Quercus rubra are not well understood. The purpose of this study was to investigate the genome-wide distribution of the PFK family members and their roles in Q. rubra by performing a systematic study of the phylogenetic relationships, molecular characteristics, motifs, chromosomal and subcellular locations, and cis-elements of QrPFKs. We identified 14 QrPFK genes in the genome of Q. rubra, followed by examining their expression in different tissues, including the roots, stems, and leaves. The phylogenetic tree divided the 14 QrPFK genes into two groups: 11 belonging to PFK and three belonging to PFP. The expression profiles of all 14 proteins were relatively the same in leaves but differed between stems and roots. Four genes (Qurub.02G189400.1, Qurub.02G189400.2, Qurub.09G134300.1, and Qurub.09G134300.2) were expressed at very low levels in both stems and roots, while two (Qurub.05G235500.1 and Qurub.05G235500.1) were expressed at low levels and the others showed relatively high expression in all tissues.

Delineating transitions during the evolution of specialised peroxisomes: Glycosome formation in kinetoplastid and diplonemid protists

One peculiarity of protists belonging to classes Kinetoplastea and Diplonemea within the phylum Euglenozoa is compartmentalisation of most glycolytic enzymes within peroxisomes that are hence called glycosomes. This pathway is not sequestered in peroxisomes of the third Euglenozoan class, Euglenida. Previous analysis of well-studied kinetoplastids, the ‘TriTryps’ parasites Trypanosoma brucei, Trypanosoma cruzi and Leishmania spp., identified within glycosomes other metabolic processes usually not present in peroxisomes. In addition, trypanosomatid peroxins, i.e. proteins involved in biogenesis of these organelles, are divergent from human and yeast orthologues. In recent years, genomes, transcriptomes and proteomes for a variety of euglenozoans have become available. Here, we track the possible evolution of glycosomes by querying these databases, as well as the genome of Naegleria gruberi, a non-euglenozoan, which belongs to the same protist supergroup Discoba. We searched for orthologues of TriTryps proteins involved in glycosomal metabolism and biogenesis. Predicted cellular location(s) of each metabolic enzyme identified was inferred from presence or absence of peroxisomal-targeting signals. Combined with a survey of relevant literature, we refine extensively our previously postulated hypothesis about glycosome evolution. The data agree glycolysis was compartmentalised in a common ancestor of the kinetoplastids and diplonemids, yet additionally indicates most other processes found in glycosomes of extant trypanosomatids, but not in peroxisomes of other eukaryotes were either sequestered in this ancestor or shortly after separation of the two lineages. In contrast, peroxin divergence is evident in all euglenozoans. Following their gain of pathway complexity, subsequent evolution of peroxisome/glycosome function is complex. We hypothesize compartmentalisation in glycosomes of glycolytic enzymes, their cofactors and subsequently other metabolic enzymes provided selective advantage to kinetoplastids and diplonemids during their evolution in changing marine environments. We contend two specific properties derived from the ancestral peroxisomes were key: existence of nonselective pores for small solutes and the possibility of high turnover by pexophagy. Critically, such pores and pexophagy are characterised in extant trypanosomatids. Increasing amenability of free-living kinetoplastids and recently isolated diplonemids to experimental study means our hypothesis and interpretation of bioinformatic data are suited to experimental interrogation.

Detection of horizontal gene transfer in the genome of the choanoflagellate Salpingoeca rosetta

Horizontal gene transfer (HGT), the movement of heritable materials between distantly related organisms, is crucial in eukaryotic evolution. However, the scale of HGT in choanoflagellates, the closest unicellular relatives of metazoans, and its possible roles in the evolution of animal multicellularity remains unexplored. We identified at least 175 candidate HGTs in the genome of the colonial choanoflagellate Salpingoeca rosetta using sequence-based tests. The majority of these were orthologous to genes in bacterial and microalgal lineages, yet displayed genomic features consistent with the rest of the S. rosetta genome-evidence of ancient acquisition events. Putative functions include enzymes involved in amino acid and carbohydrate metabolism, cell signaling, and the synthesis of extracellular matrix components. Functions of candidate HGTs may have contributed to the ability of choanoflagellates to assimilate novel metabolites, thereby supporting adaptation, survival in diverse ecological niches, and response to external cues that are possibly critical in the evolution of multicellularity in choanoflagellates.

Phosphofructokinases A and B from Mycobacterium tuberculosis Display Different Catalytic Properties and Allosteric Regulation

Tuberculosis (TB) remains one of the major health concerns worldwide. Mycobacterium tuberculosis (Mtb), the causative agent of TB, can flexibly change its metabolic processes during different life stages. Regulation of key metabolic enzyme activities by intracellular conditions, allosteric inhibition or feedback control can effectively contribute to Mtb survival under different conditions. Phosphofructokinase (Pfk) is one of the key enzymes regulating glycolysis. Mtb encodes two Pfk isoenzymes, Pfk A/Rv3010c and Pfk B/Rv2029c, which are differently expressed upon transition to the hypoxia-induced non-replicating state of the bacteria. While pfkB gene and protein expression are upregulated under hypoxic conditions, Pfk A levels decrease. Here, we present biochemical characterization of both Pfk isoenzymes, revealing that Pfk A and Pfk B display different kinetic properties. Although the glycolytic activity of Pfk A is higher than that of Pfk B, it is markedly inhibited by an excess of both substrates (fructose-6-phosphate and ATP), reaction products (fructose-1,6-bisphosphate and ADP) and common metabolic allosteric regulators. In contrast, synthesis of fructose-1,6-bisphosphatase catalyzed by Pfk B is not regulated by higher levels of substrates, and metabolites. Importantly, we found that only Pfk B can catalyze the reverse gluconeogenic reaction. Pfk B thus can support glycolysis under conditions inhibiting Pfk A function.

Biochemical and transcript level differences between the three human phosphofructokinases show optimisation of each isoform for specific metabolic niches

6-Phosphofructokinase-1-kinase (PFK) tetramers catalyse the phosphorylation of fructose 6-phosphate (F6P) to fructose 1,6-bisphosphate (F16BP). Vertebrates have three PFK isoforms (PFK-M, PFK-L, and PFK-P). This study is the first to compare the kinetics, structures, and transcript levels of recombinant human PFK isoforms. Under the conditions tested PFK-M has the highest affinities for F6P and ATP (K0.5ATP 152 µM; K0.5F6P 147 µM), PFK-P the lowest affinities (K0.5ATP 276 µM; K0.5F6P 1333 µM), and PFK-L demonstrates a mixed picture of high ATP affinity and low F6P affinity (K0.5ATP 160 µM; K0.5F6P 1360 µM). PFK-M is more resistant to ATP inhibition compared with PFK-L and PFK-P (respectively, 23%, 31%, 50% decreases in specificity constants). GTP is an alternate phospho donor. Interface 2, which regulates the inactive dimer to active tetramer equilibrium, differs between isoforms, resulting in varying tetrameric stability. Under the conditions tested PFK-M is less sensitive to fructose 2,6-bisphosphate (F26BP) allosteric modulation than PFK-L or PFK-P (allosteric constants [K0.5ATP+F26BP/K0.5ATP] 1.10, 0.92, 0.54, respectively). Structural analysis of two allosteric sites reveals one may be specialised for AMP/ADP and the other for smaller/flexible regulators (citrate or phosphoenolpyruvate). Correlations between PFK-L and PFK-P transcript levels indicate that simultaneous expression may expand metabolic capacity for F16BP production whilst preserving regulatory capabilities. Analysis of cancer samples reveals intriguing parallels between PFK-P and PKM2 (pyruvate kinase M2), and simultaneous increases in PFK-P and PFKFB3 (responsible for F26BP production) transcript levels, suggesting prioritisation of metabolic flexibility in cancers. Our results describe the kinetic and transcript level differences between the three PFK isoforms, explaining how each isoform may be optimised for distinct roles.

The pentose phosphate pathway of cellulolytic clostridia relies on 6-phosphofructokinase instead of transaldolase

The genomes of most cellulolytic clostridia do not contain genes annotated as transaldolase. Therefore, for assimilating pentose sugars or for generating C5 precursors (such as ribose) during growth on other (non-C5) substrates, they must possess a pathway that connects pentose metabolism with the rest of metabolism. Here we provide evidence that for this connection cellulolytic clostridia rely on the sedoheptulose 1,7-bisphosphate (SBP) pathway, using pyrophosphate-dependent phosphofructokinase (PPi-PFK) instead of transaldolase. In this reversible pathway, PFK converts sedoheptulose 7-phosphate (S7P) to SBP, after which fructose-bisphosphate aldolase cleaves SBP into dihydroxyacetone phosphate and erythrose 4-phosphate. We show that PPi-PFKs of Clostridium thermosuccinogenes and Clostridium thermocellum indeed can convert S7P to SBP, and have similar affinities for S7P and the canonical substrate fructose 6-phosphate (F6P). By contrast, (ATP-dependent) PfkA of Escherichia coli, which does rely on transaldolase, had a very poor affinity for S7P. This indicates that the PPi-PFK of cellulolytic clostridia has evolved the use of S7P. We further show that C. thermosuccinogenes contains a significant SBP pool, an unusual metabolite that is elevated during growth on xylose, demonstrating its relevance for pentose assimilation. Last, we demonstrate that a second PFK of C. thermosuccinogenes that operates with ATP and GTP exhibits unusual kinetics toward F6P, as it appears to have an extremely high degree of cooperative binding, resulting in a virtual on/off switch for substrate concentrations near its K½ value. In summary, our results confirm the existence of an SBP pathway for pentose assimilation in cellulolytic clostridia.

Structural comparisons of phosphoenolpyruvate carboxykinases reveal the evolutionary trajectories of these phosphodiester energy conversion enzymes

Inorganic pyrophosphate (PP_i) consists of two phosphate molecules and can act as an energy and phosphate donor in cellular reactions, similar to ATP. Several kinases use PP_i as a substrate, and these kinases have recently been suggested to have evolved from ATP-dependent functional homologs, which have significant amino acid sequence similarity to PP_i-utilizing enzymes. In contrast, phosphoenolpyruvate carboxykinase (PEPCK) can be divided into three types according to the phosphate donor (ATP, GTP, or PPi), and the amino acid sequence similarity of these PEPCKs is too low to confirm that they share a common ancestor. Here we solved the crystal structure of a PP_i-PEPCK homolog from the bacterium Actinomyces israelii at 2.6 Å resolution and compared it with previously reported structures from ATP- and GTP-specific PEPCKs to assess the degrees of similarities and divergences among these PEPCKs. These comparisons revealed that they share a tertiary structure with significant value and that amino acid residues directly contributing to substrate recognition, except for those that recognize purine moieties, are conserved. Furthermore, the order of secondary structural elements between PP_i-, ATP-, and GTP-specific PEPCKs was strictly conserved. The structure-based comparisons of the three PEPCK types provide key insights into the structural basis of PP_i specificity and suggest that all of these PEPCKs are derived from a common ancestor.

Assessing Cofactor Usage in Pseudoclostridium thermosuccinogenes via Heterologous Expression of Central Metabolic Enzymes

Pseudoclostridium thermosuccinogenes and Hungateiclostridium thermocellum are being studied for their potential to contribute to a more sustainable bio-based economy. Both species were shown previously to rely on GTP or pyrophosphate instead of ATP as cofactors in specific reactions of central energy metabolism for reasons that are not well understood yet. Since it is often impossible to predict cofactor specificity from the primary protein structure, thirteen enzymes from P. thermosuccinogenes were cloned and heterologous expressed in Escherichia coli to assess the cofactor usage in vitro and paint a more complete picture of the cofactor usage in the central metabolism of P. thermosuccinogenes. The assays were conducted with heat-treated E. coli cell-free extract devoid of background activity to allow the quick assessment of a relatively large number of (thermophilic) enzymes. Selected enzymes were also purified to allow the determination of the enzyme kinetics for competing cofactors. Following the results of the glucokinase (GK), galactokinase, xylulokinase (XK), and ribokinase assays, it seems that phosphorylation of monosaccharides by and large is mainly GTP-dependent. Some possible implications of this relating to the adenylate/guanylate energy charge are discussed here. Besides the highly expressed pyrophosphate-dependent 6-phosphofructokinase, another 6-phosphofructokinase was found to be equally dependent on ATP and GTP, while no 6-phosphofructokinase activity could be demonstrated for a third. Both type I glyceraldehyde 3-phosphate dehydrogenases were found to be NAD+-dependent, and further, acetate kinase, isocitrate dehydrogenase, and three enzymes predicted to be responsible for the interconversion of phosphoenolpyruvate and pyruvate (i.e., pyruvate kinase; pyruvate, phosphate dikinase; phosphoenolpyruvate synthase), were also assessed.

Genome-wide identification, expression and functional analysis of the phosphofructokinase gene family in Chinese white pear (Pyrus bretschneideri)

Phosphofructokinase plays an essential role in sugar metabolism in plants. Plants possess two types of phosphofructokinase proteins for phosphorylation of fructose-6-phosphate, the pyrophosphate-dependent fructose-6-phosphate phosphotransferase (PFP), and the ATP-dependent phosphofructokinase (PFK). Until now, the gene evolution, expression patterns, and functions of phosphofructokinase proteins were unknown in pear. In this report, 14 phosphofructokinase genes were identified in pear. The phylogenetic tree indicated that the phosphofructokinase gene family could be grouped into two subfamilies, with 10 genes belonging to the PbPFK subfamily, and 4 genes belonging to the PbPFP subfamily. Conserved motifs and exon numbers of the phosphofructokinase were found in pear and other six species. The evolution analysis indicated that WGD/Segmental and dispersed duplications were the main duplication models for the phosphofructokinase genes expansion in pear and other six species. Analysis of cis-regulatory element sequences of all phosphofructokinase genes identified light regulation and the MYB binding site in the promoter of all pear phosphofructokinase genes, suggesting that phosphofructokinase might could be regulated by light and MYB transcription factors (TFs). Gene expression patterns revealed that PbPFP1 showed similar pattern with sorbitol contents, suggesting important contributions to sugar accumulation during fruit development. Further functional analysis indicated that the phosphofructokinase gene PbPFP1 was localized on plasma membrane compartment, indicating that PbPFP1 had function in plasma membrane. Transient transformation of PbPFP1 in pear fruits led to significant increases of fructose and sorbitol compared to controls. Overall, our study provides important insights into the gene expression patterns and important potential functions of phosphofructokinase for sugar accumulation in pear fruits, which will help to enrich understanding of sugar-related bio-pathways and lay the molecular basis for fruit quality improvement.

The kinetic characteristics of human and trypanosomatid phosphofructokinases for the reverse reaction

Eukaryotic ATP-dependent phosphofructokinases (PFKs) are often considered unidirectional enzymes catalysing the transfer of a phospho moiety from ATP to fructose 6-phosphate to produce ADP and fructose 1,6-bisphosphate. The reverse reaction is not generally considered to occur under normal conditions and has never been demonstrated for any eukaryotic ATP-dependent PFKs, though it does occur in inorganic pyrophosphate-dependent PFKs and has been experimentally shown for bacterial ATP-dependent PFKs. The evidence is provided via two orthogonal assays that all three human PFK isoforms can catalyse the reverse reaction in vitro, allowing determination of kinetic properties. Additionally, the reverse reaction was shown possible for PFKs from three clinically important trypanosomatids; these enzymes are contained within glycosomes in vivo. This compartmentalisation may facilitate reversal, given the potential for trypanosomatids to have an altered ATP/ADP ratio in glycosomes compared with the cytosol. The kinetic properties of each trypanosomatid PFK were determined, including the response to natural and artificial modulators of enzyme activity. The possible physiological relevance of the reverse reaction in trypanosomatid and human PFKs is discussed.

Primary Metabolism in Citrus Fruit as Affected by Its Unique Structure

Citrus is one of the world's most important fruit crops, contributing essential nutrients, such as vitamin C and minerals, to the human diet. It is characterized by two important traits: first, its major edible part is composed of juice sacs, a unique structure among fruit, and second, relatively high levels of citric acid are accumulated in the vacuole of the juice sac cell. Although the major routes of primary metabolism are generally the same in citrus fruit and other plant systems, the fruit's unique structural features challenge our understanding of carbon flow into the fruit and its movement through all of its parts. In fact, acid metabolism and accumulation have only been summarized in a few reviews. Here we present a comprehensive view of sugar, acid and amino acid metabolism and their connections within the fruit, all in relation to the fruit's unique structure.

Horizontally transferred genes in the ctenophore Mnemiopsis leidyi

Horizontal gene transfer (HGT) has had major impacts on the biology of a wide range of organisms from antibiotic resistance in bacteria to adaptations to herbivory in arthropods. A growing body of literature shows that HGT between non-animals and animals is more commonplace than previously thought. In this study, we present a thorough investigation of HGT in the ctenophore Mnemiopsis leidyi. We applied tests of phylogenetic incongruence to identify nine genes that were likely transferred horizontally early in ctenophore evolution from bacteria and non-metazoan eukaryotes. All but one of these HGTs (an uncharacterized protein) are homologous to characterized enzymes, supporting previous observations that genes encoding enzymes are more likely to be retained after HGT events. We found that the majority of these nine horizontally transferred genes were expressed during development, suggesting that they are active and play a role in the biology of M. leidyi. This is the first report of HGT in ctenophores, and contributes to an ever-growing literature on the prevalence of genetic information flowing between non-animals and animals.

Identification of a pyrophosphate-dependent kinase and its donor selectivity determinants

Almost all kinases utilize ATP as their phosphate donor, while a few kinases utilize pyrophosphate (PPi) instead. PPi-dependent kinases are often homologous to their ATP-dependent counterparts, but determinants of their different donor specificities remain unclear. We identify a PPi-dependent member of the ribokinase family, which differs from known PPi-dependent kinases, and elucidate its PPi-binding mode based on the crystal structures. Structural comparison and sequence alignment reveal five important residues: three basic residues specifically recognizing PPi and two large hydrophobic residues occluding a part of the ATP-binding pocket. Two of the three basic residues adapt a conserved motif of the ribokinase family for the PPi binding. Using these five key residues as a signature pattern, we discover additional PPi-specific members of the ribokinase family, and thus conclude that these residues are the determinants of PPi-specific binding. Introduction of these residues may enable transformation of ATP-dependent ribokinase family members into PPi-dependent enzymes.

While most kinases are ATP-dependent some utilize pyrophosphate (PPi) instead. Here the authors structurally characterize a PPi-dependent kinase, identify its key recognition residues and find further PPi-dependent ribokinase family members with this signature pattern.

Genomes of ubiquitous marine and hypersaline Hydrogenovibrio, Thiomicrorhabdus and Thiomicrospira spp. encode a diversity of mechanisms to sustain chemolithoautotrophy in heterogeneous environments

Chemolithoautotrophic bacteria from the genera Hydrogenovibrio, Thiomicrorhabdus and Thiomicrospira are common, sometimes dominant, isolates from sulfidic habitats including hydrothermal vents, soda and salt lakes and marine sediments. Their genome sequences confirm their membership in a deeply branching clade of the Gammaproteobacteria. Several adaptations to heterogeneous habitats are apparent. Their genomes include large numbers of genes for sensing and responding to their environment (EAL- and GGDEF-domain proteins and methyl-accepting chemotaxis proteins) despite their small sizes (2.1-3.1 Mbp). An array of sulfur-oxidizing complexes are encoded, likely to facilitate these organisms' use of multiple forms of reduced sulfur as electron donors. Hydrogenase genes are present in some taxa, including group 1d and 2b hydrogenases in Hydrogenovibrio marinus and H. thermophilus MA2-6, acquired via horizontal gene transfer. In addition to high-affinity cbb₃ cytochrome c oxidase, some also encode cytochrome bd-type quinol oxidase or ba₃ -type cytochrome c oxidase, which could facilitate growth under different oxygen tensions, or maintain redox balance. Carboxysome operons are present in most, with genes downstream encoding transporters from four evolutionarily distinct families, which may act with the carboxysomes to form CO₂ concentrating mechanisms. These adaptations to habitat variability likely contribute to the cosmopolitan distribution of these organisms.

Rethinking the evolution of eukaryotic metabolism: novel cellular partitioning of enzymes in stramenopiles links serine biosynthesis to glycolysis in mitochondria

Background

An important feature of eukaryotic evolution is metabolic compartmentalization, in which certain pathways are restricted to the cytosol or specific organelles. Glycolysis in eukaryotes is described as a cytosolic process. The universality of this canon has been challenged by recent genome data that suggest that some glycolytic enzymes made by stramenopiles bear mitochondrial targeting peptides.

Results

Mining of oomycete, diatom, and brown algal genomes indicates that stramenopiles encode two forms of enzymes for the second half of glycolysis, one with and the other without mitochondrial targeting peptides. The predicted mitochondrial targeting was confirmed by using fluorescent tags to localize phosphoglycerate kinase, phosphoglycerate mutase, and pyruvate kinase in Phytophthora infestans, the oomycete that causes potato blight. A genome-wide search for other enzymes with atypical mitochondrial locations identified phosphoglycerate dehydrogenase, phosphoserine aminotransferase, and phosphoserine phosphatase, which form a pathway for generating serine from the glycolytic intermediate 3-phosphoglycerate. Fluorescent tags confirmed the delivery of these serine biosynthetic enzymes to P. infestans mitochondria. A cytosolic form of this serine biosynthetic pathway, which occurs in most eukaryotes, is missing from oomycetes and most other stramenopiles. The glycolysis and serine metabolism pathways of oomycetes appear to be mosaics of enzymes with different ancestries. While some of the noncanonical oomycete mitochondrial enzymes have the closest affinity in phylogenetic analyses with proteins from other stramenopiles, others cluster with bacterial, plant, or animal proteins. The genes encoding the mitochondrial phosphoglycerate kinase and serine-forming enzymes are physically linked on oomycete chromosomes, which suggests a shared origin.

Conclusions

Stramenopile metabolism appears to have been shaped through the acquisition of genes by descent and lateral or endosymbiotic gene transfer, along with the targeting of the proteins to locations that are novel compared to other eukaryotes. Colocalization of the glycolytic and serine biosynthesis enzymes in mitochondria is apparently necessary since they share a common intermediate. The results indicate that descriptions of metabolism in textbooks do not cover the full diversity of eukaryotic biology.

Electronic supplementary material

The online version of this article (10.1186/s12862-017-1087-8) contains supplementary material, which is available to authorized users.

Glycolysis as the Central Core of Fermentation

The increasing concerns of greenhouse gas emissions have increased the interest in dark fermentation as a means of productions for industrial chemicals, especially from renewable cellulosic biomass. However, the metabolism, including glycolysis, of many candidate organisms for cellulosic biomass conversion through consolidated bioprocessing is still poorly understood and the genomes have only recently been sequenced. Because a variety of industrial chemicals are produced directly from sugar metabolism, the careful understanding of glycolysis from a genomic and biochemical point of view is essential in the development of strategies for increasing product yields and therefore increasing industrial potential. The current review discusses the different pathways available for glycolysis along with unexpected variations from traditional models, especially in the utilization of alternate energy intermediates (GTP, pyrophosphate). This reinforces the need for a careful description of interactions between energy metabolites and glycolysis enzymes for understanding carbon and electron flux regulation.

N-Terminal Presequence-Independent Import of Phosphofructokinase into Hydrogenosomes of Trichomonas vaginalis

Mitochondrial evolution entailed the origin of protein import machinery that allows nuclear-encoded proteins to be targeted to the organelle, as well as the origin of cleavable N-terminal targeting sequences (NTS) that allow efficient sorting and import of matrix proteins. In hydrogenosomes and mitosomes, reduced forms of mitochondria with reduced proteomes, NTS-independent targeting of matrix proteins is known. Here, we studied the cellular localization of two glycolytic enzymes in the anaerobic pathogen Trichomonas vaginalis: PPi-dependent phosphofructokinase (TvPPi-PFK), which is the main glycolytic PFK activity of the protist, and ATP-dependent PFK (TvATP-PFK), the function of which is less clear. TvPPi-PFK was detected predominantly in the cytosol, as expected, while all four TvATP-PFK paralogues were imported into T. vaginalis hydrogenosomes, although none of them possesses an NTS. The heterologous expression of TvATP-PFK in Saccharomyces cerevisiae revealed an intrinsic capability of the protein to be recognized and imported into yeast mitochondria, whereas yeast ATP-PFK resides in the cytosol. TvATP-PFK consists of only a catalytic domain, similarly to "short" bacterial enzymes, while ScATP-PFK includes an N-terminal extension, a catalytic domain, and a C-terminal regulatory domain. Expression of the catalytic domain of ScATP-PFK and short Escherichia coli ATP-PFK in T. vaginalis resulted in their partial delivery to hydrogenosomes. These results indicate that TvATP-PFK and the homologous ATP-PFKs possess internal structural targeting information that is recognized by the hydrogenosomal import machinery. From an evolutionary perspective, the predisposition of ancient ATP-PFK to be recognized and imported into hydrogenosomes might be a relict from the early phases of organelle evolution.

Discovery of PPi-type Phosphoenolpyruvate Carboxykinase Genes in Eukaryotes and Bacteria

Phosphoenolpyruvate carboxykinase (PEPCK) is one of the pivotal enzymes that regulates the carbon flow of the central metabolism by fixing CO2 to phosphoenolpyruvate (PEP) to produce oxaloacetate or vice versa. Whereas ATP- and GTP-type PEPCKs have been well studied, and their protein identities are established, inorganic pyrophosphate (PPi)-type PEPCK (PPi-PEPCK) is poorly characterized. Despite extensive enzymological studies, its protein identity and encoding gene remain unknown. In this study, PPi-PEPCK has been identified for the first time from a eukaryotic human parasite, Entamoeba histolytica, by conventional purification and mass spectrometric identification of the native enzyme, followed by demonstration of its enzymatic activity. A homolog of the amebic PPi-PEPCK from an anaerobic bacterium Propionibacterium freudenreichii subsp. shermanii also exhibited PPi-PEPCK activity. The primary structure of PPi-PEPCK has no similarity to the functional homologs ATP/GTP-PEPCKs and PEP carboxylase, strongly suggesting that PPi-PEPCK arose independently from the other functional homologues and very likely has unique catalytic sites. PPi-PEPCK homologs were found in a variety of bacteria and some eukaryotes but not in archaea. The molecular identification of this long forgotten enzyme shows us the diversity and functional redundancy of enzymes involved in the central metabolism and can help us to understand the central metabolism more deeply.

HAMAP in 2015: updates to the protein family classification and annotation system

HAMAP (High-quality Automated and Manual Annotation of Proteins—available at http://hamap.expasy.org/) is a system for the automatic classification and annotation of protein sequences. HAMAP provides annotation of the same quality and detail as UniProtKB/Swiss-Prot, using manually curated profiles for protein sequence family classification and expert curated rules for functional annotation of family members. HAMAP data and tools are made available through our website and as part of the UniRule pipeline of UniProt, providing annotation for millions of unreviewed sequences of UniProtKB/TrEMBL. Here we report on the growth of HAMAP and updates to the HAMAP system since our last report in the NAR Database Issue of 2013. We continue to augment HAMAP with new family profiles and annotation rules as new protein families are characterized and annotated in UniProtKB/Swiss-Prot; the latest version of HAMAP (as of 3 September 2014) contains 1983 family classification profiles and 1998 annotation rules (up from 1780 and 1720). We demonstrate how the complex logic of HAMAP rules allows for precise annotation of individual functional variants within large homologous protein families. We also describe improvements to our web-based tool HAMAP-Scan which simplify the classification and annotation of sequences, and the incorporation of an improved sequence-profile search algorithm.

Bacterial influences on animal origins

Animals evolved in seas teeming with bacteria, yet the influences of bacteria on animal origins are poorly understood. Comparisons among modern animals and their closest living relatives, the choanoflagellates, suggest that the first animals used flagellated collar cells to capture bacterial prey. The cell biology of prey capture, such as cell adhesion between predator and prey, involves mechanisms that may have been co-opted to mediate intercellular interactions during the evolution of animal multicellularity. Moreover, a history of bacterivory may have influenced the evolution of animal genomes by driving the evolution of genetic pathways for immunity and facilitating lateral gene transfer. Understanding the interactions between bacteria and the progenitors of animals may help to explain the myriad ways in which bacteria shape the biology of modern animals, including ourselves.

Enhanced whole genome sequence and annotation of Clostridium stercorarium DSM8532T using RNA-seq transcriptomics and high-throughput proteomics

BACKGROUND

Growing interest in cellulolytic clostridia with potential for consolidated biofuels production is mitigated by low conversion of raw substrates to desired end products. Strategies to improve conversion are likely to benefit from emerging techniques to define molecular systems biology of these organisms. Clostridium stercorarium DSM8532T is an anaerobic thermophile with demonstrated high ethanol production on cellulose and hemicellulose. Although several lignocellulolytic enzymes in this organism have been well-characterized, details concerning carbohydrate transporters and central metabolism have not been described. Therefore, the goal of this study is to define an improved whole genome sequence (WGS) for this organism using in-depth molecular profiling by RNA-seq transcriptomics and tandem mass spectrometry-based proteomics.

RESULTS

A paired-end Roche/454 WGS assembly was closed through application of an in silico algorithm designed to resolve repetitive sequence regions, resulting in a circular replicon with one gap and a region of 2 kilobases with 10 ambiguous bases. RNA-seq transcriptomics resulted in nearly complete coverage of the genome, identifying errors in homopolymer length attributable to 454 sequencing. Peptide sequences resulting from high-throughput tandem mass spectrometry of trypsin-digested protein extracts were mapped to 1,755 annotated proteins (68% of all protein-coding regions). Proteogenomic analysis confirmed the quality of annotation and improvement pipelines, identifying a missing gene and an alternative reading frame. Peptide coverage of genes hypothetically involved in substrate hydrolysis, transport and utilization confirmed multiple pathways for glycolysis, pyruvate conversion and recycling of intermediates. No sequences homologous to transaldolase, a central enzyme in the pentose phosphate pathway, were observed by any method, despite demonstrated growth of this organism on xylose and xylan hemicellulose.

CONCLUSIONS

Complementary omics techniques confirm the quality of genome sequence assembly, annotation and error-reporting. Nearly complete genome coverage by RNA-seq likely indicates background DNA in RNA extracts, however these preps resulted in WGS enhancement and transcriptome profiling in a single Illumina run. No detection of transaldolase by any method despite xylose utilization by this organism indicates an alternative pathway for sedoheptulose-7-phosphate degradation. This report combines next-generation omics techniques to elucidate previously undefined features of substrate transport and central metabolism for this organism and its potential for consolidated biofuels production from lignocellulose.

Genomic evaluation of Thermoanaerobacter spp. for the construction of designer co-cultures to improve lignocellulosic biofuel production

The microbial production of ethanol from lignocellulosic biomass is a multi-component process that involves biomass hydrolysis, carbohydrate transport and utilization, and finally, the production of ethanol. Strains of the genus Thermoanaerobacter have been studied for decades due to their innate abilities to produce comparatively high ethanol yields from hemicellulose constituent sugars. However, their inability to hydrolyze cellulose, limits their usefulness in lignocellulosic biofuel production. As such, co-culturing Thermoanaerobacter spp. with cellulolytic organisms is a plausible approach to improving lignocellulose conversion efficiencies and yields of biofuels. To evaluate native lignocellulosic ethanol production capacities relative to competing fermentative end-products, comparative genomic analysis of 11 sequenced Thermoanaerobacter strains, including a de novo genome, Thermoanaerobacter thermohydrosulfuricus WC1, was conducted. Analysis was specifically focused on the genomic potential for each strain to address all aspects of ethanol production mentioned through a consolidated bioprocessing approach. Whole genome functional annotation analysis identified three distinct clades within the genus. The genomes of Clade 1 strains encode the fewest extracellular carbohydrate active enzymes and also show the least diversity in terms of lignocellulose relevant carbohydrate utilization pathways. However, these same strains reportedly are capable of directing a higher proportion of their total carbon flux towards ethanol, rather than non-biofuel end-products, than other Thermoanaerobacter strains. Strains in Clade 2 show the greatest diversity in terms of lignocellulose hydrolysis and utilization, but proportionately produce more non-ethanol end-products than Clade 1 strains. Strains in Clade 3, in which T. thermohydrosulfuricus WC1 is included, show mid-range potential for lignocellulose hydrolysis and utilization, but also exhibit extensive divergence from both Clade 1 and Clade 2 strains in terms of cellular energetics. The potential implications regarding strain selection and suitability for industrial ethanol production through a consolidated bioprocessing co-culturing approach are examined throughout the manuscript.

Characterization of recombinant PPi-dependent 6-phosphofructokinases from Methylosinus trichosporium OB3b and Methylobacterium nodulans ORS 2060

The properties of the purified recombinant PPi-dependent 6-phosphofructokinases (PPi-PFKs) from the methanotroph Methylosinus trichosporium OB3b and rhizospheric phytosymbiont Methylobacterium nodulans ORS 2060 were determined. The dependence of activities of PPi-PFK-His(6)-tag from Ms. trichosporium OB3b (6 × 45 kDa) and PPi-PFK from Mb. nodulans ORS 2060 (4 × 43 kDa) on the concentrations of substrates of forward and reverse reactions conformed to Michaelis-Menten kinetics. Besides fructose-6-phosphate, the enzymes also phosphorylated sedoheptulose-7-phosphate. ADP or AMP (1 mM each) inhibited activity of the Ms. trichosporium PPi-PFK but did not affect the activity of the Mb. nodulans enzyme. Preference of PPi-PFKs to fructose-1,6-bisphosphate implied a predominant function of the enzymes in hexose phosphate synthesis in these bacteria. PPi-PFKs from the methylotrophs have low similarity of translated amino acid sequences (17% identity) and belong to different phylogenetic subgroups of type II 6-phosphofructokinases. The relationship of PPi-PFKs with microaerophilic character of Ms. trichosporium OB3b and adaptation of Mb. nodulans ORS 2060 to anaerobic phase of phytosymbiosis are discussed.

Linking genome content to biofuel production yields: a meta-analysis of major catabolic pathways among select H2 and ethanol-producing bacteria

Background

Fermentative bacteria offer the potential to convert lignocellulosic waste-streams into biofuels such as hydrogen (H₂) and ethanol. Current fermentative H₂ and ethanol yields, however, are below theoretical maxima, vary greatly among organisms, and depend on the extent of metabolic pathways utilized. For fermentative H₂ and/or ethanol production to become practical, biofuel yields must be increased. We performed a comparative meta-analysis of (i) reported end-product yields, and (ii) genes encoding pyruvate metabolism and end-product synthesis pathways to identify suitable biomarkers for screening a microorganism’s potential of H₂ and/or ethanol production, and to identify targets for metabolic engineering to improve biofuel yields. Our interest in H₂ and/or ethanol optimization restricted our meta-analysis to organisms with sequenced genomes and limited branched end-product pathways. These included members of the Firmicutes, Euryarchaeota, and Thermotogae.

Results

Bioinformatic analysis revealed that the absence of genes encoding acetaldehyde dehydrogenase and bifunctional acetaldehyde/alcohol dehydrogenase (AdhE) in Caldicellulosiruptor, Thermococcus, Pyrococcus, and Thermotoga species coincide with high H₂ yields and low ethanol production. Organisms containing genes (or activities) for both ethanol and H₂ synthesis pathways (i.e. Caldanaerobacter subterraneus subsp. tengcongensis, Ethanoligenens harbinense, and Clostridium species) had relatively uniform mixed product patterns. The absence of hydrogenases in Geobacillus and Bacillus species did not confer high ethanol production, but rather high lactate production. Only Thermoanaerobacter pseudethanolicus produced relatively high ethanol and low H₂ yields. This may be attributed to the presence of genes encoding proteins that promote NADH production. Lactate dehydrogenase and pyruvate:formate lyase are not conducive for ethanol and/or H₂ production. While the type(s) of encoded hydrogenases appear to have little impact on H₂ production in organisms that do not encode ethanol producing pathways, they do influence reduced end-product yields in those that do.

Conclusions

Here we show that composition of genes encoding pathways involved in pyruvate catabolism and end-product synthesis pathways can be used to approximate potential end-product distribution patterns. We have identified a number of genetic biomarkers for streamlining ethanol and H₂ producing capabilities. By linking genome content, reaction thermodynamics, and end-product yields, we offer potential targets for optimization of either ethanol or H₂ yields through metabolic engineering.

Horizontal gene transfer in choanoflagellates

Horizontal gene transfer (HGT), also known as lateral gene transfer, results in the rapid acquisition of genes from another organism. HGT has long been known to be a driving force in speciation in prokaryotes, and there is evidence for HGT from symbiotic and infectious bacteria to metazoans, as well as from protists to bacteria. Recently, it has become clear that as many as a 1,000 genes in the genome of the choanoflagellate Monosiga brevicollis may have been acquired by HGT. Interestingly, these genes reportedly come from algae, bacteria, and other choanoflagellate prey. Some of these genes appear to have allowed an ancestral choanoflagellate to exploit nutrient-poor environments and were not passed on to metazoan descendents. However, some of these genes are also found in animal genomes, suggesting that HGT into a common ancestor of choanozoans and animals may have contributed to metazoan evolution.

Characterization of the recombinant pyrophosphate-dependent 6-phosphofructokinases from Methylomicrobium alcaliphilum 20Z and Methylococcus capsulatus Bath

The Embden-Meyerhof-Parnas (EMP) glycolysis is the starting point of the core carbon metabolism. Aerobic methanotrophs possessing activity of the pyrophosphate-dependent 6-phosphofructokinase (PPi-PFK) instead of the classical glycolytic enzyme ATP-dependent 6-phosphofructokinase (ATP-PFK) are promising model bacteria for elucidation of the role of inorganic pyrophosphate (PPi) and PPi-dependent glycolysis in microorganisms. Characterization of the His(6)-tagged PPi-PFKs from two methanotrophs, halotolerant alkaliphilic Methylomicrobium alcaliphilum 20Z and thermotolerant Methylococcus capsulatus Bath, showed differential capabilities of PPi-PFKs to phosphorylate sedoheptulose-7-phosphate and this property correlated well with the metabolic patterns of these bacteria assimilating C(1) substrate either via the ribulosemonophosphate (RuMP) pathway (Mm. alcaliphilum 20Z) or simultaneously via the RuMP and serine pathways and the Calvin cycle (Mc. capsulatus Bath). Analysis of the genomic draft of Mm. alcaliphilum 20Z (https://www.genoscope.cns.fr/agc/mage) has provided in silico evidence for the existence of a PPi-dependent pyruvate-phosphate dikinase (PPDK). Expression of the ppdk gene at oxygen limitation along with the presence of PPi-PFK in Mm. alcaliphilum 20Z implied functioning of PPi-dependent glycolysis and PPi recycling under conditions when oxidative phosphorylation is hampered.

Characterization of recombinant pyrophosphate-dependent 6-phosphofructokinase from halotolerant methanotroph Methylomicrobium alcaliphilum 20Z

Pyrophosphate-dependent 6-phosphofructokinase (PPi-PFK) was obtained as His₆-tagged protein by cloning of the pfp gene from the aerobic obligate methanotroph Methylomicrobium alcaliphilum 20Z and characterized. The recombinant PPi-PFK (4×45 kDa) was highly active, non-allosteric and stringently specific to pyrophosphate as the phosphoryl donor. The enzyme was more specific for the reverse reaction substrate fructose-1,6-bisphosphate (K(m) 0.095 mM, V(max) 805 U/mg of protein) than for the forward reaction substrate fructose-6-phosphate (K(m) 0.64 mM, V(max) 577 U/mg of protein). It also phosphorylated sedoheptulose-7-phosphate with much lower efficiency (K(m) 1.01 mM, V(max) 0.118 U/mg of protein). The kinetic properties of the M. alcaliphilum PP(i)-PFK were analyzed and compared with those of PP(i)-PFKs from other methanotrophs. The PP(i)-PFK from M. alcaliphilum shows highest sequence identity to PPi-PFK from obligate mesophilic methanotroph Methylomonas methanica (89%), and only low identity to the enzyme from thermotolerant Methylococcus capsulatus Bath (16%). This extensive sequence divergence of PPi-PFKs correlated with differential ability to phosphorylate sedoheptulose-7-phosphate and with the metabolic patterns of these bacteria assimilating C₁ substrate either via the ribulose monophoshate (RuMP) cycle or simultaneously via the RuMP and the Calvin cycles. Based on enzymic and genomic data, the involvement of PPi-PFK in pyrophosphate-dependent glycolysis in M. alcaliphilum 20Z was fist proposed.

Pyrophosphate as a central energy carrier in the hydrogen-producing extremely thermophilic Caldicellulosiruptor saccharolyticus

The role of inorganic pyrophosphate (PPi) as an energy carrier in the central metabolism of the extremely thermophilic bacterium Caldicellulosiruptor saccharolyticus was investigated. In agreement with its annotated genome sequence, cell extracts were shown to exhibit PPi-dependent phosphofructokinase and pyruvate phosphate dikinase activity. In addition, membrane-bound pyrophosphatase activity was demonstrated, while no significant cytosolic pyrophosphatase activity was detected. During the exponential growth phase, high PPi levels (approximately 4 +/- 2 mM) and relatively low ATP levels (0.43 +/- 0.07 mM) were found, and the PPi/ATP ratio decreased 13-fold when the cells entered the stationary phase. Pyruvate kinase activity appeared to be allosterically affected by PPi. Altogether, these findings suggest an important role for PPi in the central energy metabolism of C. saccharolyticus.

Acetate and succinate production in amoebae, helminths, diplomonads, trichomonads and trypanosomatids: common and diverse metabolic strategies used by parasitic lower eukaryotes

Parasites that often grow anaerobically in their hosts have adopted a fermentative strategy relying on the production of partially oxidized end products, including lactate, glycerol, ethanol, succinate and acetate. This review focuses on recent progress in understanding acetate production in protist parasites, such as amoebae, diplomonads, trichomonads, trypanosomatids and in the metazoan parasites helminths, as well as the succinate production pathway(s) present in some of them. We also describe the unconventional organisation of the tricarboxylic acid cycle associated with the fermentative strategy adopted by the procyclic trypanosomes, which may resemble the probable structure of the primordial TCA cycle in prokaryotes.

Expression of phosphofructokinase in Neisseria meningitidis

Neisseria meningitidis serogroup B is a pathogen that can infect diverse sites within the human host. According to the N. meningitidis genomic information and experimental observations, glucose can be completely catabolized through the Entner-Doudoroff pathway and the pentose phosphate pathway. The Embden-Meyerhof-Parnas pathway is not functional, because the gene for phosphofructokinase (PFK) is not present. The phylogenetic distribution of PFK indicates that in most obligate aerobic organisms, PFK is lacking. We conclude that this is because of the limited contribution of PFK to the energy supply in aerobically grown organisms in comparison with the energy generated through oxidative phosphorylation. Under anaerobic or microaerobic conditions, the available energy is limiting, and PFK provides an advantage, which explains the presence of PFK in many (facultatively) anaerobic organisms. In accordance with this, in silico flux balance analysis predicted an increase of biomass yield as a result of PFK expression. However, analysis of a genetically engineered N. meningitidis strain that expressed a heterologous PFK showed that the yield of biomass on substrate decreased in comparison with a pfkA-deficient control strain, which was associated mainly with an increase in CO(2) production, whereas production of by-products was similar in the two strains. This might explain why the pfkA gene has not been obtained by horizontal gene transfer, since it is initially unfavourable for biomass yield. No large effects related to heterologous expression of pfkA were observed in the transcriptome. Although our results suggest that introduction of PFK does not contribute to a more efficient strain in terms of biomass yield, achievement of a robust, optimal metabolic network that enables a higher growth rate or a higher biomass yield might be possible after adaptive evolution of the strain, which remains to be investigated.

Cross-species analysis of the glycolytic pathway by comparison of molecular interaction fields

The electrostatic potential of an enzyme is a key determinant of its substrate interactions and catalytic turnover. Here we invoke comparative analysis of protein electrostatic potentials, along with sequence and structural analysis, to classify and characterize all the enzymes in an entire pathway across a set of different organisms. The electrostatic potentials of the enzymes from the glycolytic pathway of 11 eukaryotes were analyzed by qPIPSA (quantitative protein interaction property similarity analysis). The comparison allows the functional assignment of neuron-specific isoforms of triosephosphate isomerase from zebrafish, the identification of unusual protein surface interaction properties of the mosquito glucose-6-phosphate isomerase and the functional annotation of ATP-dependent phosphofructokinases and cofactor-dependent phosphoglycerate mutases from plants. We here show that plants possess two parallel pathways to convert glucose. One is similar to glycolysis in humans, the other is specialized to let plants adapt to their environmental conditions. We use differences in electrostatic potentials to estimate kinetic parameters for the triosephosphate isomerases from nine species for which published parameters are not available. Along the core glycolytic pathway, phosphoglycerate mutase displays the most conserved electrostatic potential. The largest cross-species variations are found for glucose-6-phosphate isomerase, enolase and fructose-1,6-bisphosphate aldolase. The extent of conservation of electrostatic potentials along the pathway is consistent with the absence of a single rate-limiting step in glycolysis.

Plant-like phosphofructokinase from Plasmodium falciparum belongs to a novel class of ATP-dependent enzymes

Malaria parasite-infected erythrocytes exhibit enhanced glucose utilisation and 6-phospho-1-fructokinase (PFK) is a key enzyme in glycolysis. Here we present the characterisation of PFK from the human malaria parasite Plasmodium falciparum. Of the two putative PFK genes on chromosome 9 (PfPFK9) and 11 (PfPFK11), only the PfPFK9 gene appeared to possess all the catalytic features appropriate for PFK activity. The deduced PfPFK proteins contain domains homologous to the plant-like pyrophosphate (PPi)-dependent PFK beta and alpha subunits, which are quite different from the human erythrocyte PFK protein. The PfPFK9 gene beta and alpha regions were cloned and expressed as His(6)- and GST-tagged proteins in Escherichia coli. Complementation of PFK-deficient E. coli and activity analysis of purified recombinant proteins confirmed that PfPFK9beta possessed catalytic activity. Monoclonal antibodies against the recombinant beta protein confirmed that the PfPFK9 protein has beta and alpha domains fused into a 200 kDa protein, as opposed to the independent subunits found in plants. Despite an overall structural similarity to plant PPi-PFKs, the recombinant protein and the parasite extract exhibited only ATP-dependent enzyme activity, and none with PPi. Unlike host PFK, the Plasmodium PFK was insensitive to fructose-2,6-bisphosphate (F-2,6-bP), phosphoenolpyruvate (PEP) and citrate. A comparison of the deduced PFK proteins from several protozoan PFK genome databases implicates a unique class of ATP-dependent PFK present amongst the apicomplexan protozoans.

Characterization of the pyrophosphate-dependent 6-phosphofructokinase from Methylococcus capsulatus Bath

An active pyrophosphate-dependent 6-phosphofructokinase (PPi-PFK) from the thermotolerant methanotroph Methylococcus capsulatus Bath, containing a six-residue polyhistidine tag, was characterized. The enzyme was homodimeric (2 x 45 kDa), nonallosteric and most active at pH 7.0. PPi-PFK catalyzed reactions of PPi-dependent phosphorylation of fructose-6-phosphate (F-6-P) (K(m) 2.27 mM and V(max) 7.6 U mg(-1) of protein), sedoheptulose-7-phosphate (K(m) 0.027 mM and V(max) 31 U mg(-1)) and ribulose-5-phosphate. In the reaction with F-6-P, the apparent K(m) for PPi was 0.027 mM, while in the reverse reaction, K(m) for orthophosphate was 8.69 mM and that for fructose-1,6-bisphosphate 0.328 mM (V(max) 9.0 U mg(-1)). Phylogenetically, M. capsulatus PPi-PFK was most similar to PPi-PFKs from the lithoautotrophic ammonia oxidizers Nitrosomonas europaea (74.0%), Nitrosospira multiformis (73.6%) and Betaproteobacterial methylotroph Methylibium petroleiphilum PM1 (71.6% identity). Genes coding PPi-PFK and a putative V-type H(+)-translocating pyrophosphatase (H(+)-PPi-ase) were cotranscribed as an operon. The potential significance of the PPi-PFK for regulation of carbon and energy fluxes in M. capsulatus Bath is discussed.

The crystal structure of ATP-bound phosphofructokinase from Trypanosoma brucei reveals conformational transitions different from those of other phosphofructokinases

The crystal structure of the ATP-bound form of the tetrameric phosphofructokinase (PFK) from Trypanosoma brucei enables detailed comparisons to be made with the structures of the apoenzyme form of the same enzyme, as well as with those of bacterial ATP-dependent and PP(i)-dependent PFKs. The active site of T. brucei PFK (which is strictly ATP-dependent but belongs to the PP(i)-dependent family by sequence similarities) is a chimera of the two types of PFK. In particular, the active site of T. brucei PFK possesses amino acid residues and structural features characteristic of both types of PFK. Conformational changes upon ATP binding are observed that include the opening of the active site to accommodate the two substrates, MgATP and fructose 6-phosphate, and a dramatic ordering of the C-terminal helices, which act like reaching arms to hold the tetramer together. These conformational transitions are fundamentally different from those of other ATP-dependent PFKs. The substantial differences in structure and mechanism of T. brucei PFK compared with bacterial and mammalian PFKs give optimism for the discovery of species-specific drugs for the treatment of diseases caused by protist parasites of the trypanosomatid family.

Immunoproteomic assay of extracellular proteins in Streptococcus equi ssp. zooepidemicus

A proteomic approach combining two-dimensional electrophoresis, Western blot and matrix-assisted laser desorption tandem time-of-flight mass spectrometry has been used to map the extracellular proteins of Streptococcus equi ssp. zooepidemicus (S. zooepidemicus) strain ATCC 35246. These bioinformatic technologies facilitated the identification of novel S. zooepidemicus vaccine candidate antigens and therapeutic agents. Despite the limitations posed by the unavailability of complete genome and proteome data for S. zooepidemicus, seven of 15 chosen immunogenic spots were successfully identified as streptococcal proteins (AE1 and AE4 c. 10) from homologous Streptococcus species. Among these, AE6 and AE7 were identified as S. zooepidemicus UDP-N-acetyl-glucosamine pyrophosphorylase and UDP-glucose pyrophosphorylase proteins. In addition, AE4 was determined to be glyceraldehyde-3-phosphate dehydrogenase from Enterococcus faecalis. Following SIGNALIP 3.0 (http://www.cbs.dtu.dk/servicess/SignalIP) prediction, data suggested that AE5, AE7 and AE9 contained signal peptides. BLAST (http://www.sanger.ac.uk) results found that nucleotide sequences of all identified proteins shared high homology (> or = 65%) with S. zooepidemicus. The majority of proteins identified in our study remain formally unreported in S. zooepidemicus. However, these proteins serve a vital role in the immune system and reproduction of host species. Therefore, we further evaluated the proteins as vaccine candidates in this study.

Phylogeny vs genome reshuffling: horizontal gene transfer

The evolutionary events in organisms can be tracked to the transfer of genetic material. The inheritance of genetic material among closely related organisms is a slow evolutionary process. On the other hand, the movement of genes among distantly related species can account for rapid evolution. The later process has been quite evident in the appearance of antibiotic resistance genes among human and animal pathogens. Phylogenetic trees based on such genes and those involved in metabolic activities reflect the incongruencies in comparison to the 16S rDNA gene, generally used for taxonomic relationships. Such discrepancies in gene inheritance have been termed as horizontal gene transfer (HGT) events. In the post-genomic era, the explosion of known sequences through large-scale sequencing projects has unraveled the weakness of traditional 16S rDNA gene tree based evolutionary model. Various methods to scrutinize HGT events include atypical composition, abnormal sequence similarity, anomalous phylogenetic distribution, unusual phyletic patterns, etc. Since HGT generates greater genetic diversity, it is likely to increase resource use and ecosystem resilience.

Does anoxia tolerance involve altering the energy currency towards PPi?

Oxygen deficit is an important abiotic stress influencing plants, because this condition results in an 'energy crisis'. Most species only survive short periods of anoxia, but several wetland species tolerate prolonged anoxia. Transcriptomic and proteomic studies, using anoxia-tolerant rice and anoxia-intolerant Arabidopsis, have provided evidence for the selective adoption of pyrophosphate (PPi) over ATP as high-energy donor molecules, which may contribute to anoxia tolerance. The use of PPi in some tolerant plant species is similar to that observed in many anaerobic prokaryotes. Investigations are being performed to better understand the origin and regulation of reversible PPi-dependent glycolytic enzymes such as cytosolic pyruvate phosphate dikinase, as well as PPi-consuming enzymes, which are engaged during the anoxic energy crisis. This will be crucial in unraveling this currency switch and its contribution to anoxia tolerance.

Characterisation of the ATP-dependent phosphofructokinase gene family from Arabidopsis thaliana

Plants possess two different types of phosphofructokinases, an ATP-dependent (PFK) and a pyrophosphate-dependent form (PFP). While plant PFPs have been investigated in detail, cDNA clones coding for PFK have not been identified in Arabidopsis thaliana. Searching the A. thaliana genome revealed 11 putative members of a phosphofructokinase gene family. Among those, four sequences showed high homology to the alpha- or beta-subunits of plant PFPs. Seven cDNAs resulted in elevated PFK, but not PFP activity after transient expression in tobacco leaves suggesting that they encode Arabidopsis PFKs. RT-PCR revealed different tissue-specific expression of the individual forms. Furthermore, analysis of GFP fusion proteins indicated their presence in different sub-cellular compartments.

A novel form of 6-phosphofructokinase. Identification and functional relevance of a third type of subunit in Pichia pastoris

Classically, 6-phosphofructokinases are homo- and hetero-oligomeric enzymes consisting of alpha subunits and alpha/beta subunits, respectively. Herein, we describe a new form of 6-phosphofructokinase (Pfk) present in several Pichia species, which is composed of three different types of subunit, alpha, beta, and gamma. The sequence of the gamma subunit shows no similarity to classic Pfk subunits or to other known protein sequences. In-depth structural and functional studies revealed that the gamma subunit is a constitutive component of Pfk from Pichia pastoris (PpPfk). Analyses of the purified PpPfk suggest a heterododecameric assembly from the three different subunits. Accordingly, it is the largest and most complex Pfk identified yet. Although, the gamma subunit is not required for enzymatic activity, the gamma subunit-deficient mutant displays a decreased growth on nutrient limitation and reduced cell flocculation when compared with the P. pastoris wild-type strain. Subsequent characterization of purified Pfks from wild-type and gamma subunit-deficient strains revealed that the allosteric regulation of the PpPfk by ATP, fructose 2,6-bisphosphate, and AMP is fine-tuned by the gamma subunit. Therefore, we suggest that the gamma subunit contributes to adaptation of P. pastoris to energy resources.

Purification, microsequencing and cloning of spinach ATP-dependent phosphofructokinase link sequence and function for the plant enzyme

Despite its importance in plant metabolism, no sequences of higher plant ATP-dependent phosphofructokinase (EC 2.7.1.11) are annotated in the databases. We have purified the enzyme from spinach leaves 309-fold to electrophoretic homogeneity. The purified enzyme was a homotetramer of approximately 52 kDa subunits with a specific activity of 600 mU x mg(-1) and a Km value for ATP of 81 microm. The purified enzyme was not activated by phosphate, but slightly inhibited instead, suggesting that it was the chloroplast isoform. The inclusion of adenosine 5'-(beta,gamma-imido)triphosphate was conducive to enzyme activity during the purification protocol. The sequences of eight tryptic peptides from the final protein preparation, which did not utilize pyrophosphate as a phosphoryl donor, were determined and an exactly corresponding cDNA was cloned. The sequence of enzymatically active spinach ATP-dependent phosphofructokinase suggests that a large family of genomics-derived higher plant sequences currently annotated in the databases as putative pyrophosphate-dependent phosphofructokinases according to sequence similarity is misannotated with respect to the cosubstrate.

The glycolytic pathway of Trimastix pyriformis is an evolutionary mosaic

Background

Glycolysis and subsequent fermentation is the main energy source for many anaerobic organisms. The glycolytic pathway consists of ten enzymatic steps which appear to be universal amongst eukaryotes. However, it has been shown that the origins of these enzymes in specific eukaryote lineages can differ, and sometimes involve lateral gene transfer events. We have conducted an expressed sequence tag (EST) survey of the anaerobic flagellate Trimastix pyriformis to investigate the nature of the evolutionary origins of the glycolytic enzymes in this relatively unstudied organism.

Results

We have found genes in the Trimastix EST data that encode enzymes potentially catalyzing nine of the ten steps of the glycolytic conversion of glucose to pyruvate. Furthermore, we have found two different enzymes that in principle could catalyze the conversion of phosphoenol pyruvate (PEP) to pyruvate (or the reverse reaction) as part of the last step in glycolysis. Our phylogenetic analyses of all of these enzymes revealed at least four cases where the relationship of the Trimastix genes to homologs from other species is at odds with accepted organismal relationships. Although lateral gene transfer events likely account for these anomalies, with the data at hand we were not able to establish with confidence the bacterial donor lineage that gave rise to the respective Trimastix enzymes.

Conclusion

A number of the glycolytic enzymes of Trimastix have been transferred laterally from bacteria instead of being inherited from the last common eukaryotic ancestor. Thus, despite widespread conservation of the glycolytic biochemical pathway across eukaryote diversity, in a number of protist lineages the enzymatic components of the pathway have been replaced by lateral gene transfer from disparate evolutionary sources. It remains unclear if these replacements result from selectively advantageous properties of the introduced enzymes or if they are neutral outcomes of a gene transfer 'ratchet' from food or endosymbiotic organisms or a combination of both processes.

Reconstructing the mosaic glycolytic pathway of the anaerobic eukaryote Monocercomonoides

All eukaryotes carry out glycolysis, interestingly, not all using the same enzymes. Anaerobic eukaryotes face the challenge of fewer molecules of ATP extracted per molecule of glucose due to their lack of a complete tricarboxylic acid cycle. This may have pressured anaerobic eukaryotes to acquire the more ATP-efficient alternative glycolytic enzymes, such as pyrophosphate-fructose 6-phosphate phosphotransferase and pyruvate orthophosphate dikinase, through lateral gene transfers from bacteria and other eukaryotes. Most studies of these enzymes in eukaryotes involve pathogenic anaerobes; Monocercomonoides, an oxymonad belonging to the eukaryotic supergroup Excavata, is a nonpathogenic anaerobe representing an evolutionarily and ecologically distinct sampling of an anaerobic glycolytic pathway. We sequenced cDNA encoding glycolytic enzymes from a previously established cDNA library of Monocercomonoides and analyzed the relationships of these enzymes to those from other organisms spanning the major groups of Eukaryota, Bacteria, and Archaea. We established that, firstly, Monocercomonoides possesses alternative versions of glycolytic enzymes: fructose-6-phosphate phosphotransferase, both pyruvate kinase and pyruvate orthophosphate dikinase, cofactor-independent phosphoglycerate mutase, and fructose-bisphosphate aldolase (class II, type B). Secondly, we found evidence for the monophyly of oxymonads, kinetoplastids, diplomonads, and parabasalids, the major representatives of the Excavata. We also found several prokaryote-to-eukaryote as well as eukaryote-to-eukaryote lateral gene transfers involving glycolytic enzymes from anaerobic eukaryotes, further suggesting that lateral gene transfer was an important factor in the evolution of this pathway for denizens of this environment.

The first crystal structure of phosphofructokinase from a eukaryote: Trypanosoma brucei

The crystal structure of the ATP-dependent phosphofructokinase (PFK) from Trypanosoma brucei provides the first detailed description of a eukaryotic PFK, and enables comparisons to be made with the crystal structures of bacterial ATP-dependent and PPi-dependent PFKs. The structure reveals that two insertions (the 17-20 and 329-348 loops) that are characteristic of trypanosomatid PFKs, but absent from bacterial and mammalian ATP-dependent PFKs, are located within and adjacent to the active site, and are in positions to play important roles in the enzyme's mechanism. The 90 residue N-terminal extension forms a novel domain that includes an "embracing arm" across the subunit boundary to the symmetry-related subunit in the tetrameric enzyme. Comparisons with the PPi-dependent PFK from Borrelia burgdorferi show that several features thought to be characteristic of PPi-dependent PFKs are present in the trypanosome ATP-dependent PFK. These two enzymes are generally more similar to each other than to the bacterial or mammalian ATP-dependent PFKs. However, there are critical differences at the active site of PPi-dependent PFKs that are sufficient to prevent the binding of ATP. This crystal structure of a eukaryotic PFK has enabled us to propose a detailed model of human muscle PFK that shows active site and other differences that offer opportunities for structure-based drug discovery for the treatment of sleeping sickness and other diseases caused by the trypanosomatid family of protozoan parasites.

Pyruvate-phosphate dikinase of oxymonads and parabasalia and the evolution of pyrophosphate-dependent glycolysis in anaerobic eukaryotes

In pyrophosphate-dependent glycolysis, the ATP/ADP-dependent enzymes phosphofructokinase (PFK) and pyruvate kinase are replaced by the pyrophosphate-dependent PFK and pyruvate phosphate dikinase (PPDK), respectively. This variant of glycolysis is widespread among bacteria, but it also occurs in a few parasitic anaerobic eukaryotes such as Giardia and Entamoeba spp. We sequenced two genes for PPDK from the amitochondriate oxymonad Streblomastix strix and found evidence for PPDK in Trichomonas vaginalis and other parabasalia, where this enzyme was thought to be absent. The Streblomastix and Giardia genes may be related to one another, but those of Entamoeba and perhaps Trichomonas are distinct and more closely related to bacterial homologues. These findings suggest that pyrophosphate-dependent glycolysis is more widespread in eukaryotes than previously thought, enzymes from the pathway coexists with ATP-dependent more often than previously thought and may be spread by lateral transfer of genes for pyrophosphate-dependent enzymes from bacteria.

Horizontal gene transfer from Bacteria to rumen Ciliates indicates adaptation to their anaerobic, carbohydrates-rich environment

Background

The horizontal transfer of expressed genes from Bacteria into Ciliates which live in close contact with each other in the rumen (the foregut of ruminants) was studied using ciliate Expressed Sequence Tags (ESTs). More than 4000 ESTs were sequenced from representatives of the two major groups of rumen Cilates: the order Entodiniomorphida (Entodinium simplex, Entodinium caudatum, Eudiplodinium maggii, Metadinium medium, Diploplastron affine, Polyplastron multivesiculatum and Epidinium ecaudatum) and the order Vestibuliferida, previously called Holotricha (Isotricha prostoma, Isotricha intestinalis and Dasytricha ruminantium).

Results

A comparison of the sequences with the completely sequenced genomes of Eukaryotes and Prokaryotes, followed by large-scale construction and analysis of phylogenies, identified 148 ciliate genes that specifically cluster with genes from the Bacteria and Archaea. The phylogenetic clustering with bacterial genes, coupled with the absence of close relatives of these genes in the Ciliate Tetrahymena thermophila, indicates that they have been acquired via Horizontal Gene Transfer (HGT) after the colonization of the gut by the rumen Ciliates.

Conclusion

Among the HGT candidates, we found an over-representation (>75%) of genes involved in metabolism, specifically in the catabolism of complex carbohydrates, a rich food source in the rumen. We propose that the acquisition of these genes has greatly facilitated the Ciliates' colonization of the rumen providing evidence for the role of HGT in the adaptation to new niches.