Integrating the universal metabolism into a phylogenetic analysis

The darwinian concept of "descent with modification" applies to metabolic pathways: pathways sharing similarities must have inherited them from an exclusive, hypothetical ancestral pathway. Comparative anatomy of biochemical pathways is performed using five criteria of homology. Primary homologies of "type I" were defined as several pathways sharing the same enzyme with high specificity for its substrate. Primary homologies of "type II" were defined as the sharing of similar enzymatic functions, cofactors, functional family, or recurrence of a set of reactions. Standard cladistic analysis is used to infer the evolutionary history of metabolic development and the relative ordering of biochemical reactions through time, from a single matrix integrating the whole basic universal metabolism. The cladogram shows that the earliest pathways to emerge are metabolism of amino acids of groups I and II (Asp, Asn, Glu, and Gln). The earliest enzymatic functions are mostly linked to amino acid catabolism: deamination, transamination, and decarboxylation. For some amino acids, catabolism and biosynthesis occur at the same time (Asp, Glu, Lys, and Met). Catabolism precedes anabolism for Asn, Gln, Arg, Trp, His, Tyr, and Phe, and anabolism precedes catabolism for Pro, Ala, Leu, Val, Ile, Cys, Gly, Ser, and Thr. The urea cycle evolves from arginine synthesis. Metabolism of fatty acids and sugars develops after the full development of metabolism of amino acids of groups I and II, and they are associated with the anabolism of amino acids of groups III and IV. Syntheses of aromatic amino acids are branched within sugar metabolism. The Krebs cycle occurs relatively late after the setting of metabolism of amino acids of groups I and II. One portion of the Krebs cycle has a catabolic origin, whereas the other portion has an anabolic origin in pathways of amino acids of groups III and IV. It is not possible to order glycolysis and gluconeogenesis with regard to the Krebs cycle, as they all belong to "period 6." Pentose-phosphate and Calvin cycles are later (periods 7 and 8, respectively). Cladistic analysis of the structure of biochemical pathways makes hypotheses in biochemical evolution explicit and parsimonious.     

Evolution of amino acid metabolism inferred through cladistic analysis

Because free amino acids were most probably available in primitive abiotic environments, their metabolism is likely to have provided some of the very first metabolic pathways of life. What were the first enzymatic reactions to emerge? A cladistic analysis of metabolic pathways of the 16 aliphatic amino acids and 2 portions of the Krebs cycle was performed using four criteria of homology. The analysis is not based on sequence comparisons but, rather, on coding similarities in enzyme properties. The properties used are shared specific enzymatic activity, shared enzymatic function without substrate specificity, shared coenzymes, and shared functional family. The tree shows that the earliest pathways to emerge are not portions of the Krebs cycle but metabolisms of aspartate, asparagine, glutamate, and glutamine. The views of Horowitz (Horowitz, N. H. (1945) Proc. Natl. Acad. Sci. U. S. A. 31, 153-157) and Cordón (Cordón, F. (1990) Tratado Evolucionista de Biologia, Aguilar, Madrid, Spain), according to which the upstream reactions in the catabolic pathways and the downstream reactions in the anabolic pathways are the earliest in evolution, are globally corroborated; however, with some exceptions. These are due to later opportunistic connections of pathways (actually already suggested by these authors). Earliest enzymatic functions are mostly catabolic; they were deaminations, transaminations, and decarboxylations. From the consensus tree we extracted four time spans for amino acid metabolism development. For some amino acids catabolism and biosynthesis occurred at the same time (Asp, Glu, Lys, Leu, Ala, Val, Ile, Pro, Arg). For others ultimate reactions that use amino acids as a substrate or as a product are distinct in time, with catabolism preceding anabolism for Asn, Gln, and Cys and anabolism preceding catabolism for Ser, Met, and Thr. Cladistic analysis of the structure of biochemical pathways makes hypotheses in biochemical evolution explicit and parsimonious.     

New prospects for deducing the evolutionary history of metabolic pathways in prokaryotes: Aromatic biosynthesis as a case-in-point

Metabolic pathways of prokaryotes are more biochemically diverse than is generally recognized. Distinctive biochemical features are shared by phylogenetic clusters. The hierarchical levels of characterstate clustering depends upon evolutionary events which fortuitously became fixed in the genome of a common ancestor. Prokaryotes can now be ordered on a phylogenetic tree. This allows the evolutionary steps that underlie the construction and regulation of appropriately complex biochemical pathways to be traced in an evolutionary progression of prokaryote types that house these pathways. Essentially the approach is to deduce ancestral character states at ever deeper phylogenetic levels, utilizing logical principles of maximum parsimony. The current perspective on the evolution of the biochemical pathway for biosynthesis of aromatic amino acids is developed as a case-in-point model for analyses that should be feasible with many major metabolic systems. Phenylalanine biosynthesis probably arose prior to the addition of branches leading to tyrosine and tryptophan. An evolutionary scenario is developed that begins with non-enzymatic reactions which may have operated in primitive systems, followed by the evolution of an enzymatic system that pre-dated the divergence of major lineages of modern eubacteria (Gram-positive bacteria, Gram-negative purple bacteria, and cyanobacteria).Florida Agricultural Experiment Station, Journal Series No. 8251.     

Structure versus time in the evolutionary diversification of avian carotenoid metabolic networks

Historical associations of genes and proteins are thought to delineate pathways available to subsequent evolution; however, the effects of past functional involvements on contemporary evolution are rarely quantified. Here, we examined the extent to which the structure of a carotenoid enzymatic network persists in avian evolution. Specifically, we tested whether the evolution of carotenoid networks was most concordant with phylogenetically structured expansion from core reactions of common ancestors or with subsampling of biochemical pathway modules from an ancestral network. We compared structural and historical associations in 467 carotenoid networks of extant and ancestral species and uncovered the overwhelming effect of pre‐existing metabolic network structure on carotenoid diversification over the last 50 million years of avian evolution. Over evolutionary time, birds repeatedly subsampled and recombined conserved biochemical modules, which likely maintained the overall structure of the carotenoid metabolic network during avian evolution. These findings explain the recurrent convergence of evolutionary distant species in carotenoid metabolism and weak phylogenetic signal in avian carotenoid evolution. Remarkable retention of an ancient metabolic structure throughout extensive and prolonged ecological diversification in avian carotenoid metabolism illustrates a fundamental requirement of organismal evolution – historical continuity of a deterministic network that links past and present functional associations of its components.     

Comprehensive proteome profiling of the Fe(III)‐reducing myxobacterium Anaeromyxobacter dehalogenans 2CP‐C during growth with fumarate and ferric citrate

Anaeromyxobacter dehalogenans is a microaerophilic member of the delta‐proteobacteria which is able to utilize a wide range of electron acceptors, including halogenated phenols, U(VI), Fe(III), nitrate, nitrite, oxygen and fumarate. To date, the knowledge regarding general metabolic activities of this ecologically relevant bacterium is limited. Here, we present a first systematic 2‐D reference map of the soluble A. dehalogenans proteome in order to provide a sound basis for further proteomic studies as well as to gain first global insights into the metabolic activities of this bacterium. Using a combination of 2‐DE and MALDI‐TOF‐MS, a total of 720 proteins spots were identified, representing 559 unique protein species. Using the proteome data, altogether 50 metabolic pathways were found to be expressed during growth with fumarate as primary electron acceptor. An analysis of the pathways revealed an extensive display of enzymes involved in the catabolism and anabolism of a variety of amino acids, including the unexpected fermentation of lysine to butyrate. Moreover, using the reference gel as basis, a semi‐quantitative analysis of protein expression changes of A. dehalogenans during growth with ferric citrate as electron acceptor was conducted. The adaptation to Fe(III) reducing conditions involved the expression changes of a total of 239 proteins. The results suggest that the adaptation to Fe(III) reductive conditions involves an increase in metabolic flux through the tricarboxylic acid cycle, which is fueled by an increased catabolism of amino acids.     

Early steps of metabolism evolution inferred by cladistic analysis of amino acid catabolic pathways

Among abiotic molecules available in primitive environments, free amino acids are good candidates as the first source of energy and molecules for early protocells. Amino acid catabolic pathways are likely to be one of the very first metabolic pathways of life. Among them, which ones were the first to emerge? A cladistic analysis of catabolic pathways of the sixteen aliphatic amino acids and two portions of the Krebs cycle is performed using four criteria of homology. The cladogram shows that the earliest pathways to emerge are not portions of the Krebs cycle but catabolisms of aspartate, asparagine, glutamate, glutamine, proline, arginine. Earliest enzymatic catabolic functions were deaminations and transaminations. Later on appeared enzymatic decarboxylations. The consensus tree allows to propose four time spans for catabolism development and corroborates the views of Cordón in 1990 about the evolution of catabolism.     

Constraining the timing of the Great Oxidation Event within the Rubisco phylogenetic tree

Ribulose 1,5‐bisphosphate (RuBP) carboxylase/oxygenase (RuBisCO, or Rubisco) catalyzes a key reaction by which inorganic carbon is converted into organic carbon in the metabolism of many aerobic and anaerobic organisms. Across the broader Rubisco protein family, homologs exhibit diverse biochemical characteristics and metabolic functions, but the evolutionary origins of this diversity are unclear. Evidence of the timing of Rubisco family emergence and diversification of its different forms has been obscured by a meager paleontological record of early Earth biota, their subcellular physiology and metabolic components. Here, we use computational models to reconstruct a Rubisco family phylogenetic tree, ancestral amino acid sequences at branching points on the tree, and protein structures for several key ancestors. Analysis of historic substitutions with respect to their structural locations shows that there were distinct periods of amino acid substitution enrichment above background levels near and within its oxygen‐sensitive active site and subunit interfaces over the divergence between Form III (associated with anoxia) and Form I (associated with oxia) groups in its evolutionary history. One possible interpretation is that these periods of substitutional enrichment are coincident with oxidative stress exerted by the rise of oxygenic photosynthesis in the Precambrian era. Our interpretation implies that the periods of Rubisco substitutional enrichment inferred near the transition from anaerobic Form III to aerobic Form I ancestral sequences predate the acquisition of Rubisco by fully derived cyanobacterial (i.e., dual photosystem‐bearing, oxygen‐evolving) clades. The partitioning of extant lineages at high clade levels within our Rubisco phylogeny indicates that horizontal transfer of Rubisco is a relatively infrequent event. Therefore, it is possible that the mutational enrichment periods between the Form III and Form I common ancestral sequences correspond to the adaptation of key oxygen‐sensitive components of Rubisco prior to, or coincident with, the Great Oxidation Event.     

Biochemical pathways in prokaryotes can be traced backward through evolutionary time

For the first time, a credible prokaryotic phylogenetic tree is being assembled by Woese and others using quantitative sequence analysis of oligonucleotides in the highly conservative rRNA. This provides an evolutionary scale against which the evolutionary steps that led to the arrangement and regulation of contemporary biochemical pathways can be measured. This paper presents an emerging evolutionary picture of aromatic amino acid biosynthesis within a large superfamily assemblage of prokaryotes that is sufficiently developed to illustrate a new perspective that will be applicable to many other biochemical pathways.      

Metabolic investigation in psychiatric disorders

A multiplicity of theories have been proposed over the years that aim to conceptualize the pathophysiology of neuropsychiatric disorders, including impaired neurotransmission, viral infections, genetic mutation, energy metabolism deficiency, excitotoxicity, oxidative stress, and others. It is likely that complex disorders such as schizophrenia, bipolar disorder, and major depression are associated with multiple etiologies and pathogenetic mechanisms. In light of the interwoven biochemistry of human organs, identifying a network of multiple interacting biochemical pathways that account for the constellation of clinical and biological features would advance our understanding of these disorders. One such approach is to evaluate simultaneously the multiple metabolites in order to uncover the dynamic relations in the relevant biochemical systems. These metabolites are a group of low-molecular-weight, redox-active compounds, such as antioxidants, amino acids, catecholamines vitamins, lipids, and nucleotides, which reflect the metabolic processes, including anabolism and catabolism as well as other related cellular processes (e.g., signal transduction, regulation, detoxification, etc.). Such an analytic approach has the potential to yield valuable insights into the likely complex pathophysiological mechanisms that affect multiple metabolic pathways and thereby offer multiple windows of therapeutic opportunities.     

The Mitochondrion: An Integration Point of Cellular Metabolism and Signalling

In addition to efficient synthesis of ATP by oxidative phosphorylation, acquisition of the mitochondrial endosymbiont brought a whole range of new metabolic capabilities to the ancestral eukaryotic cell lineage such that the mitochondrion retains an important role in numerous anabolic and catabolic processes. While respiration dominates metabolism of the mitochondrion, this organelle is also important in the catabolism of amino acids and the provision of carbon skeletons for biosynthesis of a wide range of compounds including amino acids, vitamins, lipids, and tetrapyrroles. However, mitochondrial metabolism is best understood in the context of cellular metabolism as a whole; this is particularly true in auxotrophic organisms such as plants. For this reason understanding of the integration of mitochondrial metabolism with associated metabolic pathways in distinct cellular locations is of great importance. The examples of photorespiration, proline, cysteine, branched chain amino acid, ascorbate and folate metabolism all indicate that mitochondrial steps in these pathways are critical to their function and regulation. Moreover, the central metabolic position of the mitochondrion and its key roles in bioenergetics and redox regulation, additionally mean that it is ideally placed to act as a sensor of the biochemical status of the cell. When taken together these observations suggest that the myriad nonrespiratory functions of the mitochondria are of vast importance in the coordination of plant cellular metabolism and function.     

The emergence and early evolution of biological carbon-fixation

The fixation of into living matter sustains all life on Earth, and embeds the biosphere within geochemistry. The six known chemical pathways used by extant organisms for this function are recognized to have overlaps, but their evolution is incompletely understood. Here we reconstruct the complete early evolutionary history of biological carbon-fixation, relating all modern pathways to a single ancestral form. We find that innovations in carbon-fixation were the foundation for most major early divergences in the tree of life. These findings are based on a novel method that fully integrates metabolic and phylogenetic constraints. Comparing gene-profiles across the metabolic cores of deep-branching organisms and requiring that they are capable of synthesizing all their biomass components leads to the surprising conclusion that the most common form for deep-branching autotrophic carbon-fixation combines two disconnected sub-networks, each supplying carbon to distinct biomass components. One of these is a linear folate-based pathway of reduction previously only recognized as a fixation route in the complete Wood-Ljungdahl pathway, but which more generally may exclude the final step of synthesizing acetyl-CoA. Using metabolic constraints we then reconstruct a “phylometabolic” tree with a high degree of parsimony that traces the evolution of complete carbon-fixation pathways, and has a clear structure down to the root. This tree requires few instances of lateral gene transfer or convergence, and instead suggests a simple evolutionary dynamic in which all divergences have primary environmental causes. Energy optimization and oxygen toxicity are the two strongest forces of selection. The root of this tree combines the reductive citric acid cycle and the Wood-Ljungdahl pathway into a single connected network. This linked network lacks the selective optimization of modern fixation pathways but its redundancy leads to a more robust topology, making it more plausible than any modern pathway as a primitive universal ancestral form.     

Domain loss has independently occurred multiple times in plant terpene synthase evolution

The extensive family of plant terpene synthases (TPSs) generally has a bi-domain structure, yet phylogenetic analyses consistently indicate that these synthases have evolved from larger diterpene synthases. In particular, that duplication of the diterpene synthase genes required for gibberellin phytohormone biosynthesis provided an early predecessor, whose loss of a approximately 220 amino acid 'internal sequence element' (now recognized as the γ domain) gave rise to the precursor of the modern mono- and sesqui-TPSs found in all higher plants. Intriguingly, TPSs are conserved by taxonomic relationships rather than function. This relationship demonstrates that such functional radiation has occurred both repeatedly and relatively recently, yet phylogenetic analyses assume that the 'internal/γ' domain loss represents a single evolutionary event. Here we provide evidence that such a loss was not a singular event, but rather has occurred multiple times. Specifically, we provide an example of a bi-domain diterpene synthase from Salvia miltiorrhiza, along with a sesquiterpene synthase from Triticum aestivum (wheat) that is not only closely related to diterpene synthases, but retains the ent-kaurene synthase activity relevant to the ancestral gibberellin metabolic function. Indeed, while the wheat sesquiterpene synthase clearly no longer contains the 'internal/γ' domain, it is closely related to rice diterpene synthase genes that retain the ancestral tri-domain structure. Thus, these findings provide examples of key evolutionary intermediates that underlie the bi-domain structure observed in the expansive plant TPS gene family, as well as indicating that 'internal/γ' domain loss has occurred independently multiple times, highlighting the complex evolutionary history of this important enzymatic family.     

Evolution of metabolisms: a new method for the comparison of metabolic pathways using genomics information.

The abundance of information provided by completely sequenced genomes defines a starting point for new insights in the multilevel organization of organisms and their evolution. At the lowest level enzymes and other protein complexes are formed by aggregating multiple polypeptides. At a higher level enzymes group conceptually into metabolic pathways as part of a dynamic information-processing system, and substrates are processed by enzymes yielding other substrates. A method based on a combination of sequence information with graph topology of the underlying pathway is presented. With this approach pathways of different organisms are related to each other by phylogenetic analysis, extending conventional phylogenetic analysis of individual enzymes. The new method is applied to pathways related to electron transfer and to the Krebs citric acid cycle. In addition to providing a more comprehensive understanding of similarities and differences between organisms, this method indicates different evolutionary rates between substrates and enzymes.     

The alphaD-globin gene originated via duplication of an embryonic alpha-like globin gene in the ancestor of tetrapod vertebrates

Gene duplication is thought to play an important role in the co-option of existing protein functions to new physiological pathways. The globin superfamily of genes provides an excellent example of the kind of physiological versatility that can be attained through the functional and regulatory divergence of duplicated genes that encode different subunit polypeptides of the tetrameric hemoglobin protein. In contrast to prevailing views about the evolutionary history of the alpha-globin gene family, here we present phylogenetic evidence that the alpha(A)- and alpha(D)-globin genes are not the product of a single, tandem duplication of an ancestral globin gene with adult function in the common ancestor of extant birds, reptiles, and mammals. Instead, our analysis reveals that the alpha(D)-globin gene of amniote vertebrates arose via duplication of an embryonic alpha-like globin gene that predated the radiation of tetrapods. The important evolutionary implication is that the distinct biochemical properties of alpha(D)-hemoglobin (HbD) are not exclusively derived characters that can be attributed to a post-duplication process of neofunctionalization. Rather, many of the distinct biochemical properties of HbD are retained ancestral characters that reflect the fact that the alpha(D)-globin gene arose via duplication of a gene that had a larval/embryonic function. These insights into the evolutionary origin of HbD illustrate how adaptive modifications of physiological pathways may result from the retention and opportunistic co-option of ancestral protein functions.     

The primordial metabolism: an ancestral interconnection between leucine,arginine, and lysine biosynthesis

Background

It is generally assumed that primordial cells had small genomes with simple genes coding for enzymes able to react with a wide range of chemically related substrates, interconnecting different metabolic routes. New genes coding for enzymes with a narrowed substrate specificity arose by paralogous duplication(s) of ancestral ones and evolutionary divergence. In this way new metabolic pathways were built up by primordial cells. Useful hints to disclose the origin and evolution of ancestral metabolic routes and their interconnections can be obtained by comparing sequences of enzymes involved in the same or different metabolic routes. From this viewpoint, the lysine, arginine, and leucine biosynthetic routes represent very interesting study-models. Some of the lys, arg and leu genes are paralogs; this led to the suggestion that their ancestor genes might interconnect the three pathways. The aim of this work was to trace the evolutionary pathway leading to the appearance of the extant biosynthetic routes and to try to disclose the interrelationships existing between them and other pathways in the early stages of cellular evolution.

Results

The comparative analysis of the genes involved in the biosynthesis of lysine, leucine, and arginine, their phylogenetic distribution and analysis revealed that the extant metabolic "grids" and their interrelationships might be the outcome of a cascade of duplication of ancestral genes that, according to the patchwork hypothesis, coded for unspecific enzymes able to react with a wide range of substrates. These genes belonged to a single common pathway in which the three biosynthetic routes were highly interconnected between them and also to methionine, threonine, and cell wall biosynthesis. A possible evolutionary model leading to the extant metabolic scenarios was also depicted.

Conclusion

The whole body of data obtained in this work suggests that primordial cells synthesized leucine, lysine, and arginine through a single common metabolic pathway, whose genes underwent a set of duplication events, most of which can have predated the appearance of the last common universal ancestor of the three cell domains (Archaea, Bacteria, and Eucaryotes). The model proposes a relative timing for the appearance of the three routes and also suggests a possible evolutionary pathway for the assembly of bacterial cell-wall.

     

A comparison of the composition of silk proteins produced by spiders and insects.

Proteins that are highly expressed and composed of amino acids that are costly to synthesize are likely to place a greater drain on an organism's energy resources than proteins that are composed of ingested amino acids or ones that are metabolically simple to produce. Silks are highly expressed proteins produced by all spiders and many insects. We compared the metabolic costs of silks spun by arthropods by calculating the amount of ATP required to produce their component amino acids. Although a definitive conclusion requires detailed information on the dietary pools of amino acids available to arthropods, on the basis of the central metabolic pathways, silks spun by herbivorous, Lepidoptera larvae require significantly less ATP to synthesize than the dragline silks spun by predatory spiders. While not enough data are available to draw a statistically based conclusion, comparison of homologous silks across ancestral and derived taxa of the Araneoidea seems to suggest an evolutionary trend towards reduced silk costs. However, comparison of the synthetic costs of dragline silks across all araneomorph spiders suggests a complicated evolutionary pattern that cannot be attributed to phylogenetic position alone. We propose that the diverse silk-producing systems of the araneoid spiders (including three types of protein glues and three types of silk fibroin), evolved through intra-organ competition and that taxon-specific differences in the composition of silks drawn from homologous glands may reflect limited or fluctuating amino acid availability. The different functional properties of spider silks may be a secondary result of selection acting on different polypeptide templates.     

Evolutionary model of protein secondary structure capable of revealing new biological relationships

Ancestral sequence reconstruction has had recent success in decoding the origins and the determinants of complex protein functions. However, phylogenetic analyses of remote homologues must handle extreme amino acid sequence diversity resulting from extended periods of evolutionary change. We exploited the wealth of protein structures to develop an evolutionary model based on protein secondary structure. The approach follows the differences between discrete secondary structure states observed in modern proteins and those hypothesized in their immediate ancestors. We implemented maximum likelihood-based phylogenetic inference to reconstruct ancestral secondary structure. The predictive accuracy from the use of the evolutionary model surpasses that of comparative modeling and sequence-based prediction; the reconstruction extracts information not available from modern structures or the ancestral sequences alone. Based on a phylogenetic analysis of a sequence-diverse protein family, we showed that the model can highlight relationships that are evolutionarily rooted in structure and not evident in amino acid-based analysis.     

New Methods for Detecting Positive Selection at Single Amino Acid Sites

Inferring positive selection at single amino acid sites is of particular importance for studying evolutionary mechanisms of a protein. For this purpose, Suzuki and Gojobori (1999) developed a method (SG method) for comparing the rates of synonymous and nonsynonymous substitutions at each codon site in a protein-coding nucleotide sequence, using ancestral codons at interior nodes of the phylogenetic tree as inferred by the maximum parsimony method. In the SG method, however, selective neutrality of nucleotide substitutions cannot be tested at codon sites, where only termination codons are inferred at any interior node or the number of equally parsimonious inferences of ancestral codons at all interior nodes exceeds 10,000. Here I present a modified SG method which is free from these problems. Specifically, I use the distance-based Bayesian method for inferring the single most likely ancestral codon from 61 sense codons at each interior node. In the computer simulation and real data analysis, the modified SG method showed a higher overall efficiency of detecting positive selection than the original SG method, particularly at highly polymorphic codon sites. These results indicate that the modified SG method is useful for inferring positive selection at codon sites where neutrality cannot be tested by the original SG method. I also discuss that the p-distance is preferable to the number of synonymous substitutions for inferring the phylogenetic tree in the SG method, and present a maximum likelihood method for detecting positive selection at single amino acid sites, which produced reasonable results in the real data analysis.     

Probability Steiner trees and maximum parsimony in phylogenetic analysis

The phylogenetic tree (PT) problem has been studied by a number of researchers as an application of the Steiner tree problem, a well-known network optimisation problem. Of all the methods developed for phylogenies the maximum parsimony (MP) method is a simple and commonly used method because it relies on directly observable changes in the input nucleotide or amino acid sequences. In this paper we show that the non-uniqueness of the evolutionary pathways in the MP method leads us to consider a new model of PTs. In this so-called probability representation model, for each site a node in a PT is modelled by a probability distribution of nucleotide or amino acid states, and hence the PT at a given site is a probability Steiner tree, i.e. a Steiner tree in a high-dimensional vector space. In spite of the generality of the probability representation model, in this paper we restrict our study to constructing probability phylogenetic trees (PPT) using the parsimony criterion, as well as discussing and comparing our approach with the classical MP method. We show that for a given input set although the optimal topology as well as the total tree length of the PPT is the same as the PT constructed by the classical MP method, the inferred ancestral states and branch lengths are different and the results given by our method provide a plausible alternative to the classical ones.     

Evolutionary relationships of adenylation domains in fungi

Non-ribosomal peptide synthetases (NRPSs) and NRPS-like enzymes are abundant in microbes as they are involved in the production of primary and secondary metabolites. In contrast to the well-studied NRPSs, known to produce non-ribosomal peptides, NRPS-like enzymes exhibit more diverse activities and their evolutionary relationships are unclear. Here, we present the first in-depth phylogenetic analysis of fungal NRPS-like A domains from functionally characterized pathways, and their relationships to characterized A domains found in fungal NRPSs. This study clearly differentiated amino acid reductases, including NRPSs, from CoA/AMP ligases, which could be divided into 10 distinct phylogenetic clades that reflect their conserved domain organization, substrate specificity and enzymatic activity. In particular, evolutionary relationships of adenylate forming reductases could be refined and explained the substrate specificity difference. Consistent with their phylogeny, the deduced amino acid code of A domains differentiated amino acid reductases from other enzymes. However, a diagnostic code was found for α-keto acid reductases and clade 7 CoA/AMP ligases only. Comparative genomics of loci containing these enzymes revealed that they can be independently recruited as tailoring genes in diverse secondary metabolite pathways. Based on these results, we propose a refined and clear phylogeny-based classification of A domain-containing enzymes, which will provide a robust framework for future functional analyses and engineering of these enzymes to produce new bioactive molecules.