Expression of L- and M-type pyruvate kinase in human tissues

Pyruvate kinase (PK) has four isozymes (L, R, M1, M2) that are encoded by two different genes of PK L and M. Differential splicing produces L- and R-type PK mRNA and M1- and M2-type PK mRNA from the PK L gene and the PK M gene, respectively. The nucleotide sequences of the 3'-noncoding region are the same between the L- and the R-type PK and between the M1- and the M2-type PK. We isolated 3'-noncoding sequences of human L- and M2-type PK cDNA to construct L-type and M-type PK specific probes. With these probes, we performed Northern blot analysis of the RNA samples extracted from human tissues. Northern blot analysis showed that both kidney and liver had mRNAs that hybridized with both the L and M probes. Small intestine, skeletal muscle, brain, testis, and lung mRNAs hybridized only with the M probe. Our probes are considered useful for the detection of the types of PK isozymes expressed in small amounts, which are very difficult to detect using the conventional PK polyacrylamide gel electrophoretic method.     

Human liver type pyruvate kinase: cDNA cloning and chromosomal assignment

Pyruvate kinase (PK) has four isozymes (L,R,M1,M2) that are encoded mainly by two different genes. We isolated a cDNA clone from a Japanese adult liver lambda gt10 cDNA library by using a rat liver(L)-type PK cDNA probe. One positively hybridizing clone, hlPK-1, which contained a 1,049-base pair cDNA insert, was subjected to DNA sequence analysis. Comparisons of the sequence data with the rat PK cDNAs indicated that the cDNA encoded information for the carboxyl terminal 105 amino acids of a human L-type PK and a 3' untranslated region of 734 nucleotides. Furthermore, the karyotype analysis of several human-mouse hybrid cells and Southern blot analysis of DNAs of the hybrids with a hlPK-1 indicated that the human L-type PK gene is located on chromosome 1.     

Exon-intron structure and evolution of the Lipocalin gene family

The Lipocalins are an ancient protein family whose expression is currently confirmed in bacteria, protoctists, plants, arthropods, and chordates. The evolution of this protein family has been assessed previously using amino acid sequence phylogenies. In this report we use an independent set of characters derived from the gene structure (exon-intron arrangement) to infer a new lipocalin phylogeny. We also present the novel gene structure of three insect lipocalins. The position and phase of introns are well preserved among lipocalin clades when mapped onto a protein sequence alignment, suggesting the homologous nature of these introns. Because of this homology, we use the intron position and phase of 23 lipocalin genes to reconstruct a phylogeny by maximum parsimony and distance methods. These phylogenies are very similar to the phylogenies derived from protein sequence. This result is confirmed by congruence analysis, and a consensus tree shows the commonalities between the two source trees. Interestingly, the intron arrangement phylogeny shows that metazoan lipocalins have more introns than other eukaryotic lipocalins, and that intron gains have occurred in the C-termini of chordate lipocalins. We also analyze the relationship of intron arrangement and protein tertiary structure, as well as the relationship of lipocalins with members of the proposed structural superfamily of calycins. Our congruence analysis validates the gene structure data as a source of phylogenetic information and helps to further refine our hypothesis on the evolutionary history of lipocalins.     

Assessment of phylogenomic and orthology approaches for phylogenetic inference

MOTIVATION: Phylogenomics integrates the vast amount of phylogenetic information contained in complete genome sequences, and is rapidly becoming the standard for reliably inferring species phylogenies. There are, however, fundamental differences between the ways in which phylogenomic approaches like gene content, superalignment, superdistance and supertree integrate the phylogenetic information from separate orthologous groups. Furthermore, they all depend on the method by which the orthologous groups are initially determined. Here, we systematically compare these four phylogenomic approaches, in parallel with three approaches for large-scale orthology determination: pairwise orthology, cluster orthology and tree-based orthology. RESULTS: Including various phylogenetic methods, we apply a total of 54 fully automated phylogenomic procedures to the fungi, the eukaryotic clade with the largest number of sequenced genomes, for which we retrieved a golden standard phylogeny from the literature. Phylogenomic trees based on gene content show, relative to the other methods, a bias in the tree topology that parallels convergence in lifestyle among the species compared, indicating convergence in gene content. CONCLUSIONS: Complete genomes are no guarantee for good or even consistent phylogenies. However, the large amounts of data in genomes enable us to carefully select the data most suitable for phylogenomic inference. In terms of performance, the superalignment approach, combined with restrictive orthology, is the most successful in recovering a fungal phylogeny that agrees with current taxonomic views, and allows us to obtain a high-resolution phylogeny. We provide solid support for what has grown to be a common practice in phylogenomics during its advance in recent years. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Structural basis for tumor pyruvate kinase M2 allosteric regulation and catalysis

Four isozymes of pyruvate kinase are differentially expressed in human tissue. Human pyruvate kinase isozyme M2 (hPKM2) is expressed in early fetal tissues and is progressively replaced by the other three isozymes, M1, R, and L, immediately after birth. In most cancer cells, hPKM2 is once again expressed to promote tumor cell proliferation. Because of its almost ubiquitous presence in cancer cells, hPKM2 has been designated as tumor specific PK-M2, and its presence in human plasma is currently being used as a molecular marker for the diagnosis of various cancers. The X-ray structure of human hPKM2 complexed with Mg(2+), K(+), the inhibitor oxalate, and the allosteric activator fructose 1,6-bisphosphate (FBP) has been determined to a resolution of 2.82 A. The active site of hPKM2 is in a partially closed conformation most likely resulting from a ligand-induced domain closure promoted by the binding of FBP. In all four subunits of the enzyme tetramer, a conserved water molecule is observed on the 2-si face of the prospective enolate and supports the hypothesis that a proton-relay system is acting as the proton donor of the reaction (1). Significant structural differences among the human M2, rabbit muscle M1, and the human R isozymes are observed, especially in the orientation of the FBP-activating loop, which is in a closed conformation when FBP is bound. The structural differences observed between the PK isozymes could potentially be exploited as unique structural templates for the design of allosteric drugs against the disease states associated with the various PK isozymes, especially cancer and nonspherocytic hemolytic anemia.     

Heterologous Expression and Biochemical Characterisation of Fourteen Esterases from Helicoverpa armigera

Esterases have recurrently been implicated in insecticide resistance in Helicoverpa armigera but little is known about the underlying molecular mechanisms. We used a baculovirus system to express 14 of 30 full-length esterase genes so far identified from midgut cDNA libraries of this species. All 14 produced esterase isozymes after native PAGE and the isozymes for seven of them migrated to two regions of the gel previously associated with both organophosphate and pyrethroid resistance in various strains. Thirteen of the enzymes obtained in sufficient yield for further analysis all showed tight binding to organophosphates and low but measurable organophosphate hydrolase activity. However there was no clear difference in activity between the isozymes from regions associated with resistance and those from elsewhere in the zymogram, or between eight of the isozymes from a phylogenetic clade previously associated with resistance in proteomic and quantitative rtPCR experiments and five others not so associated. By contrast, the enzymes differed markedly in their activities against nine pyrethroid isomers and the enzymes with highest activity for the most insecticidal isomers were from regions of the gel and, in some cases, the phylogeny that had previously been associated with pyrethroid resistance. Phospholipase treatment confirmed predictions from sequence analysis that three of the isozymes were GPI anchored. This unusual feature among carboxylesterases has previously been suggested to underpin an association that some authors have noted between esterases and resistance to the Cry1Ac toxin from Bacillus thuringiensis. However these three isozymes did not migrate to the zymogram region previously associated with Cry1Ac resistance.     

Isozymes of the glycolytic and pentose-phosphate pathways in storage tissues of different oilseeds

Isozymes of hexose-phosphate isomerase (HPI; EC 5.3.1.9), pyruvate kinase (PK; EC 2.7.1.40) and 6-phosphogluconate dehydrogenase (6PGDH; EC 1.1.1.44) have been detected in the developing cotyledons of soybean (Glycine max (L.) Merr.), safflower (Carthamnus tinctorius L.) and sunflower (Helianthus annuus L.). In each seed there are two isozymes each of PK and HPI. The isozyme patterns of 6PGDH are more complex: soybean has two forms of the enzyme, safflower three, and sunflower six. In each tissue, at least 25% of the activity of each of the three enzymes is in the plastids. This supports the proposal that the glycolytic and pentose-phosphate pathways are operating in the plastids and that the plastids are the site of long-chain fatty-acid biosynthesis in developing oilseeds.Abbreviations HPI hexose-phosphate isomerase - 6PGDH 6-phosphogluconate dehydrogenase - PK pyruvate kinase     

Purification and properties of two isozymes of pyruvate kinase from Mucor racemosus.

The dimorphic phycomycete Mucor racemosus was found to contain up to five electrophoretic forms of pyruvate kinase (ATP: pyruvate 2-O-phosphotransferase, EC 2.7.1.40) depending on growth conditions. M. racemosus hyphal cells grown on glutamic acid as the carbon source contained only the fastest electrophoretic form, designated PK1, while yeast cells grown on glucose contained only the slowest electrophoretic form, PK5. Intermediate electrophoretic forms PK2, PK3, and PK4 as well as PK1 and PK5 were found in hyphal cells grown on media containing fructose or cellibiose. All five electrophoretic forms had molecular weights of ca. 230,000 as determined from plots of log Rm versus acrylamide gel concentration. Both PK1 and PK5 were purified to homogeneity and determined to be homotetramers, with subunit molecular weights of 54,000 and 58,100, respectively. The amino acid content of PK1 and PK5 was determined and found to be similar but not identical. Analysis of limited tryptic digests and cyanogen bromide cleavage fragments of PK1 and PK5 indicate that the subunits of the two isozymes are significantly different.     

Kinetic properties of liver pyruvate kinase from the flounder (Platichthys flesus L.)

Pyruvate kinase, purified from flounder liver, in two forms, i.e. PK I and PK II, is characterized by sigmoid kinetics with phosphoenolpyruvate as substrate at pH 6.3, 6.7 and 7.7. K0.5 for PEP increases with increasing pH. PK I and PK II show hyperbolic kinetics with ADP, but are inhibited by ADP concentrations above 1-2 mM. K0.5 for ADP decreases with increasing pH. PK I and PK II differ in their K0.5 values for PEP with a factor of at least 2, showing the highest figures for the latter. K0.5 for ADP is about the same for the two enzyme forms. Other nucleotide diphosphates can replace ADP as the substrate. When the nucleoside diphosphates are arranged in a rank order showing decreasing effectiveness as substrate, different rank orders are obtained for PK I and PK II.     

Detecting non-orthology in the COGs database and other approaches grouping orthologs using genome-specific best hits

Correct orthology assignment is a critical prerequisite of numerous comparative genomics procedures, such as function prediction, construction of phylogenetic species trees and genome rearrangement analysis. We present an algorithm for the detection of non-orthologs that arise by mistake in current orthology classification methods based on genome-specific best hits, such as the COGs database. The algorithm works with pairwise distance estimates, rather than computationally expensive and error-prone tree-building methods. The accuracy of the algorithm is evaluated through verification of the distribution of predicted cases, case-by-case phylogenetic analysis and comparisons with predictions from other projects using independent methods. Our results show that a very significant fraction of the COG groups include non-orthologs: using conservative parameters, the algorithm detects non-orthology in a third of all COG groups. Consequently, sequence analysis sensitive to correct orthology assignments will greatly benefit from these findings.     

Nuclear DNA restriction site polymorphisms and the phylogeny and population structure of an intertidal snail species complex (Littorina)

Primers for amplification of four novel, unlinked nuclear DNA loci, the first reported for the rough periwinkles of the genus Littorina, are described. Patterns of restriction site polymorphism for these loci are detailed within the rough periwinkles. RFLPs are not found to be diagnostic for any of the currently accepted species within this group, nor for any of the contentious subspecies, or forms, whose taxonomic status is uncertain. However, there are important differences in allele frequencies between these taxa and certain of these mirror differences detected in a previous study of the mitochondrial DNA. These allele frequency data are used to construct a phylogeny in which groupings of the three recognised species are obvious when either Nei's genetic distances or Reynold's distances are clustered. Contentious forms (L. neglecta, L. saxatilis 'b' and L. tenebrosa) do not cluster as distinct taxa, although populations of L. neglecta have important allele frequency differences from L. saxatilis. These four loci have confirmed the consensus view of Littorina phylogeny and provided important information on population structure-however four loci is insufficient for reaching definitive conclusions. Since analysis of nuclear DNA polymorphisms such as these is invaluable for analysis of phylogeny, population structure and phylogeography, identification of additional loci is considered imperative.     

An approach of orthology detection from homologous sequences under minimum evolution

In the field of phylogenetics and comparative genomics, it is important to establish orthologous relationships when comparing homologous sequences. Due to the slight sequence dissimilarity between orthologs and paralogs, it is prone to regarding paralogs as orthologs. For this reason, several methods based on evolutionary distance, phylogeny and BLAST have tried to detect orthologs with more precision. Depending on their algorithmic implementations, each of these methods sometimes has increased false negative or false positive rates. Here, we developed a novel algorithm for orthology detection that uses a distance method based on the phylogenetic criterion of minimum evolution. Our algorithm assumes that sets of sequences exhibiting orthologous relationships are evolutionarily less costly than sets that include one or more paralogous relationships. Calculation of evolutionary cost requires the reconstruction of a neighbor-joining (NJ) tree, but calculations are unaffected by the topology of any given NJ tree. Unlike tree reconciliation, our algorithm appears free from the problem of incorrect topologies of species and gene trees. The reliability of the algorithm was tested in a comparative analysis with two other orthology detection methods using 95 manually curated KOG datasets and 21 experimentally verified EXProt datasets. Sensitivity and specificity estimates indicate that the concept of minimum evolution could be valuable for the detection of orthologs.     

Homoplasy, character function, and nemertean systematics

We question recent claims that cladistic analysis is inapplicable in nemerteans (phylum Nemertea) due to a supposedly high degree of convergence. We further argue that terms like convergence and parallelism are historical sayings and only make sense in a phylogenetic context. Therefore, an approach aiming to produce phylogenetic hypotheses cannot be rejected on the grounds of a high degree of convergence before the actual hypothesis. Convergence is not an empirical observation, but a conclusion made after an analysis. We also discuss the view that knowledge of a character's function is a prerequisite for phylogenetic analysis and conclude that this is an invalid approach. Function, like any other way of sharpening our observations, helps in formulating non-phylogenetic hypotheses of homology, but the crucial test is congruence with other characters on a phylogeny.     

Mapping adenosine cyclic 3',5'-phosphate binding sites on type I and type II adenosine cyclic 3',5'-phosphate dependent protein kinases using ribose ring and cyclic phosphate ring analogues of adenosine cyclic 3',5'-phosphate

A series of adenosine cyclic 3',5'-phosphate (cAMP) derivatives containing modifications or substitutions in either the 2',3',4', or 5' position or the phosphate were examined for their abilities to activate type I isozymes of cAMP-dependent protein kinase (PK I) from rabbit or porcine skeletal muscle and type II isozymes of cAMP-dependent protein kinase (PK II) from bovine brain and heart. The studies revealed that the activation of both PK I and PK II isozymes requires a 2'-hydroxyl group in the ribo configuration, a 3' oxygen in the ribo configuration, and a charged cyclic phosphate. The two isozymes appeared to differ in those portions of their respective cAMP-binding sites that are adjacent to the 4' position of the ribose ring and the 3' position, 5' position, and phosphate portion of the cyclic phosphate ring.     

When is enough,enough in phylogenetics? A case in point from Hordeum (Poaceae)

Direct optimization was used to reconstruct the phylogeny of the 26 diploid taxa included in the genus Hordeum. The total data set was composed of 16 nucleotide sequence regions from the nuclear as well as the plastid genome. The nine nuclear regions were from single‐copy, protein coding genes located on six of the seven chromosome pairs in the diploid H. vulgare genome. The seven plastid regions comprise protein coding genes as well as intergenic regions. Studies of character congruence between data partitions showed no correlation between chromosomal location and congruence among the nuclear sequences and a level of congruence among the plastid sequences comparable with the level among the nuclear sequences. Combined analysis of all data resolved the phylogeny completely with most clades being robust and well supported. However, due to incongruence among data partitions some relationships are still and likely to remain ambiguously inferred. Rather than adding still more genes to the phylogenetic analyses, patterns of incongruence may be better explored by adding data from multiple specimens per taxon. For some species relationships the plastid data appear positively misleading, emphasizing the need for caution if plastid data are the only or dominant type of data used for phylogenetic reconstruction and subsequent re‐classification.
© The Willi Hennig Society 2011.     

Pyruvate kinase isozymes from the green alga, Selenastrum minutum. I. Purification and physical and immunological characterization

Pyruvate kinase from the green alga Selenastrum minutum consists of two isoforms (PK1 and PK2) separable by Q-Sepharose chromatography. The two isoforms have been highly purified to respective final specific activities of 42 and 23 (mumol pyruvate produced/min)/mg protein. Purification steps included salt fractionation, anion-exchange, hydrophobic interaction, and gel filtration chromatography. The final enzyme preparations differ significantly in physical and immunological properties. PK1 is heat labile and is completely inactivated following reaction with N-ethylmaleimide. In contrast, PK2 is heat-stable and is only partially inactivated following N-ethylmaleimide treatment. PK1 appears to be homotetrameric with a native molecular mass of about 240 kDa, whereas PK2 appears to be homodecameric with a native molecular mass of approximately 590 kDa. The antigenic reaction of both final PK preparations to rabbit antiserum prepared against homogeneous germinating castor bean endosperm cytosolic pyruvate kinase was tested by immunoprecipitation and Western blotting. The two algal pyruvate kinases are immunologically unrelated as only PK2 cross-reacts with the cytosolic pyruvate kinase antibodies. These data indicate that the S. minutum pyruvate kinase isoforms, PK1 and PK2, are not interconvertible forms of the same protein, but probably represent chloroplastic and cytosolic isozymes, respectively.     

Pyruvate kinase isozymes in man

Anti human M2 type and anti human L type pyruvate kinase sera allowed us to distinguish two groups of pyruvate kinase in man. Erythrocyte and liver (L type) enzymes on the one hand were inhibited by anti L and not all by anti M2 serum; pyruvate kinase from all the other tissues on the other hand were inhibited by anti M2 and not at all by anti L serum. This latter group represent the M type pyruvate kinase isozymes. The M type isozymes have been studied by electrofocusing in thin layer acrylamide-ampholine gel. In adult tissues 4 types of isozymes were found, designated, from acid to alkaline pH, as M2 (predominant form in spleen, leukocytes, lung...), M3, M4 and M1 (predominant form in muscle and brain). In foetal tissues an extra band M2, called M2f, more anodic than M2, was added to the previously described isozymes. Except in brain (in which the isozymes M2, M3, M4 and M1 were found), the most anodic bands (M2f, M2 and M3) were predominant in all the foetal tissues. The isozymes M2f and M2 seem therefore to be the original M type pyruvate kinase forms from which the other isozymes issue. The rate of each isozyme seems to depend on tissue factors characterizing the state of differentiation of some tissues, as indicated by the ability of adult muscle extracts to change the isozymes M2 and M3 into more cathodic forms.     

Multiple copies of coding as well as pseudogene c-mos sequence exist in three lacertid species

The analysis of a 581 bp section of the nuclear gene c-mos revealed multiple copies of putative functional sequences as well as pseudogenes in three closely related lacertid species Lacerta laevis, L. kulzeri and L. cyanisparsa. A phylogenetic analysis of c-mos in comparison with a molecular phylogeny based on the mitochondrial cytochrome b gene supports our findings. The study also provides new insights into the phylogenetic relationships of L. cyanisparsa and L. laevis.Pseudogenes of the three species share 11 single-nucleotide substitutions, a 1 bp deletion and a premature stop codon but differ by group-specific mutations. This result suggests that the c-mos gene has become duplicated and subsequently silenced already in the common ancestor of the three species. Sequence divergence suggests that the duplication and the loss of function occurred in the late Miocene/early Pliocene, i.e., about 5 million years ago. Indications of gene conversion are discussed.We suggest that future studies using c-mos for phylogenetic studies should provide evidence for the orthology of the sequences compared.     

ESTimating plant phylogeny: lessons from partitioning

Background  

While Expressed Sequence Tags (ESTs) have proven a viable and efficient way to sample genomes, particularly those for which whole-genome sequencing is impractical, phylogenetic analysis using ESTs remains difficult. Sequencing errors and orthology determination are the major problems when using ESTs as a source of characters for systematics. Here we develop methods to incorporate EST sequence information in a simultaneous analysis framework to address controversial phylogenetic questions regarding the relationships among the major groups of seed plants. We use an automated, phylogenetically derived approach to orthology determination called OrthologID generate a phylogeny based on 43 process partitions, many of which are derived from ESTs, and examine several measures of support to assess the utility of EST data for phylogenies.     

The Molecular Evolution of the Vertebrate Trypsinogens

We expand the already large number of known trypsinogen nucleotide and amino acid sequences by presenting additional trypsinogen sequences from the tunicate (Boltenia villosa), the lamprey (Petromyzon marinus), the pufferfish (Fugu rubripes), and the frog (Xenopus laevis). The current array of known trypsinogen sequences now spans the entire vertebrate phylogeny. Phylogenetic analysis is made difficult by the presence of multiple isozymes within species and rates of evolution that vary highly between both species and isozymes. We nevertheless present a Fitch-Margoliash phylogeny constructed from pairwise distances. We employ this phylogeny as a vehicle for speculation on the evolution of the trypsinogen gene family as well as the general modes of evolution of multigene families. Unique attributes of the lamprey and tunicate trypsinogens are noted. Received: 12 July 1997