The evolution of molecular genetic pathways and networks

There is growing interest in the evolutionary dynamics of molecular genetic pathways and networks, and the extent to which the molecular evolution of a gene depends on its position within a pathway or network, as well as over-all network topology. Investigations on the relationships between network organization, topological architecture and evolutionary dynamics provide intriguing hints as to how networks evolve. Recent studies also suggest that genetic pathway and network structures may influence the action of evolutionary forces, and may play a role in maintaining phenotypic robustness in organisms.     

Patterns of variation during adaptation in functionally linked loci

An understanding of the distribution of natural patterns of genetic variation is relevant to such fundamental biological fields as evolution and development. One recent approach to understanding such patterns has been to focus on the constraints that may arise as a function of the network or pathway context in which genes are embedded. Despite theoretical expectations of higher evolutionary constraint for genes encoding upstream versus downstream enzymes in metabolic pathways, empirical results have varied. Here we combine two complementary models from population genetics and enzyme kinetics to explore genetic variation as a function of pathway position when selection acts on whole-pathway flux. We are able to qualitatively reproduce empirically observed patterns of polymorphism and divergence and suggest that expectations should vary depending on the evolutionary trajectory of a population. Upstream genes are initially more polymorphic and diverge faster after an environmental change, while we see the opposite trend as the population approaches its fitness optimum.     

METAMODELS AND PHYLOGENETIC REPLICATION: A SYSTEMATIC APPROACH TO THE EVOLUTION OF DEVELOPMENTAL PATHWAYS

Molecular genetic analysis of phenotypic variation has revealed many examples of evolutionary change in the developmental pathways that control plant and animal morphology. A major challenge is to integrate the information from diverse organisms and traits to understand the general patterns of developmental evolution. This integration can be facilitated by evolutionary metamodels—traits that have undergone multiple independent changes in different species and whose development is controlled by well-studied regulatory pathways. The metamodel approach provides the comparative equivalent of experimental replication, allowing us to test whether the evolution of each developmental pathway follows a consistent pattern, and whether different pathways are predisposed to different modes of evolution by their intrinsic organization. A review of several metamodels suggests that the structure of developmental pathways may bias the genetic basis of phenotypic evolution, and highlights phylogenetic replication as a value-added approach that produces deeper insights into the mechanisms of evolution than single-species analyses.     

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.     

Targets of selection in a disease resistance network in wild tomatoes

Studies combining comparative genomics and information on biochemical pathways have revealed that protein evolution can be affected by the amount of pleiotropy associated with a particular gene. The amount of pleiotropy, in turn, can be a function of the position at which a gene operates in a pathway and the pathway structure. Genes that serve as convergence points and have several partners (so-called hubs) often show the greatest constraint and hence the slowest rate of protein evolution. In this article, we have studied five genes (Pto, Fen, Rin4, Prf and Pfi) in a defence signalling network in a wild tomato species, Solanum peruvianum. These proteins operate together and contribute to bacterial resistance in tomato. We predicted that Prf (and possibly Pfi), which serves as a convergence point for upstream signals, should show greater evolutionary constraint. However, we found instead that two of the genes which potentially interact with pathogen ligands, Rin4 and Fen, have evolved under strong evolutionary constraint, whereas Prf and Pfi, which probably function further downstream in the network, show evidence of balancing selection. This counterintuitive observation may be probable in pathogen defence networks, because pathogens may target positions throughout resistance networks to manipulate or nullify host resistance, thereby leaving a molecular signature of host-parasite co-evolution throughout a single network.     

The Evolution and Origin of Animal Toll-Like Receptor Signaling Pathway Revealed by Network-Level Molecular Evolutionary Analyses

Genes carry out their biological functions through pathways in complex networks consisting of many interacting molecules. Studies on the effect of network architecture on the evolution of individual proteins will provide valuable information for understanding the origin and evolution as well as functional conservation of signaling pathways. However, the relationship between the network architecture and the individual protein sequence evolution is yet little known. In current study, we carried out network-level molecular evolution analysis on TLR (Toll-like receptor ) signaling pathway, which plays an important role in innate immunity in insects and mammals, and we found that: 1) The selection constraint of genes was negatively correlated with its position along TLR signaling pathway; 2) all genes in TLR signaling pathway were highly conserved and underwent strong purifying selection; 3) the distribution of selective pressure along the pathway was driven by differential nonsynonymous substitution levels; 4) The TLR signaling pathway might present in a common ancestor of sponges and eumetazoa, and evolve via the TLR, IKK, IκB and NF-κB genes underwent duplication events as well as adaptor molecular enlargement, and gene structure and conservation motif of NF-κB genes shifted in their evolutionary history. Our results will improve our understanding on the evolutionary history of animal TLR signaling pathway as well as the relationship between the network architecture and the sequences evolution of individual protein.     

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.     

Rapid Pathway Evolution Facilitated by Horizontal Gene Transfers across Prokaryotic Lineages

The evolutionary history of biological pathways is of general interest, especially in this post-genomic era, because it may provide clues for understanding how complex systems encoded on genomes have been organized. To explain how pathways can evolve de novo, some noteworthy models have been proposed. However, direct reconstruction of pathway evolutionary history both on a genomic scale and at the depth of the tree of life has suffered from artificial effects in estimating the gene content of ancestral species. Recently, we developed an algorithm that effectively reconstructs gene-content evolution without these artificial effects, and we applied it to this problem. The carefully reconstructed history, which was based on the metabolic pathways of 160 prokaryotic species, confirmed that pathways have grown beyond the random acquisition of individual genes. Pathway acquisition took place quickly, probably eliminating the difficulty in holding genes during the course of the pathway evolution. This rapid evolution was due to massive horizontal gene transfers as gene groups, some of which were possibly operon transfers, which would convey existing pathways but not be able to generate novel pathways. To this end, we analyzed how these pathways originally appeared and found that the original acquisition of pathways occurred more contemporaneously than expected across different phylogenetic clades. As a possible model to explain this observation, we propose that novel pathway evolution may be facilitated by bidirectional horizontal gene transfers in prokaryotic communities. Such a model would complement existing pathway evolution models.     

Making new molecules - evolution of pathways for novel metabolites in plants

Plants have adapted to their environments by diversifying in various ways. This diversification is reflected at the phytochemical level in their production of numerous specialized secondary metabolites that provide protection against biotic and abiotic stresses. Plant speciation is therefore intimately linked to metabolic diversification, yet we do not currently have a deep understanding of how new metabolic pathways evolve. Recent evidence indicates that genes for individual secondary metabolic pathways can be either distributed throughout the genome or clustered, but the relative frequencies of these two pathway organizations remain to be established. While it is possible that clustering is a feature of pathways that have evolved in recent evolutionary time, the answer to this and how dispersed and clustered pathways may be related remain to be addressed. Recent advances enabled by genomics and systems biology are beginning to yield the first insights into network evolution in plant metabolism. This review focuses on recent progress in understanding the evolution of clustered and dispersed pathways for new secondary metabolites in plants.     

Metabolic systems cost-benefit analysis for interpreting network structure and regulation

MOTIVATION: Interpretation of bioinformatics data in terms of cellular function is a major challenge facing systems biology. This question is complicated by robust metabolic networks filled with structural features like parallel pathways and isozymes. Under conditions of nutrient sufficiency, metabolic networks are well known to be regulated for thermodynamic efficiency however; efficient biochemical pathways are anabolically expensive to construct. While parameters like thermodynamic efficiency have been extensively studied, a systems-based analysis of anabolic proteome synthesis 'costs' and the cellular function implications of these costs has not been reported. RESULTS: A cost-benefit analysis of an in silico Escherichia coli network revealed the relationship between metabolic pathway proteome synthesis requirements, DNA-coding sequence length, thermodynamic efficiency and substrate affinity. The results highlight basic metabolic network design principles. Pathway proteome synthesis requirements appear to have shaped biochemical network structure and regulation. Under conditions of nutrient scarcity and other general stresses, E. coli expresses pathways with relatively inexpensive proteome synthesis requirements instead of more efficient but also anabolically more expensive pathways. This evolutionary strategy provides a cellular function-based explanation for common network motifs like isozymes and parallel pathways and possibly explains 'overflow' metabolisms observed during nutrient scarcity. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Selection and the Cell Cycle: Positive Darwinian Selection in a Well-Known DNA Damage Response Pathway

Cancer is a common occurrence in multi-cellular organisms and is not strictly limited to the elderly in a population. It is therefore possible that individuals with genotypes that protect against early onset cancers have a selective advantage. In this study the patterns of mutation in the proteins of a well-studied DNA damage response pathway have been examined for evidence of adaptive evolutionary change. Using a maximum likelihood framework and the mammalian species phylogeny, together with codon models of evolution, selective pressure variation across the interacting network of proteins has been detected. The presence of signatures of adaptive evolution in BRCA1 and BRCA2 has already been documented but the effect on the entire network of interacting proteins in this damage response pathway has, until now, been unknown. Positive selection is evident throughout the network with a total of 11 proteins out of 15 examined displaying patterns of substitution characteristic of positive selection. It is also shown here that modern human populations display evidence of an ongoing selective sweep in 9 of these DNA damage repair proteins. The results presented here provide the community with new residues that may be relevant to cancer susceptibility while also highlighting those proteins where human and mouse have undergone lineage-specific functional shift. An understanding of this damage response pathway from an evolutionary perspective will undoubtedly contribute to future cancer treatment approaches.     

Optimal stoichiometric designs of ATP-producing systems as determined by an evolutionary algorithm.

The design of metabolic pathways is thought to be the result of an optimization process such that the structure of contemporary metabolic routes maximizes a particular objective function. Recently, it has been shown that some essential stoichiometric properties of glycolysis can be explained on the basis of the requirement for a high ATP production rate. Because the number of stoichiometrically feasible designs increases strongly with the number of reactions involved, a systematic analysis of all the possibilities turns out to be inaccessible beyond a certain system size. We present, therefore, an alternative approach to compute in a more efficient way the optimal design of glycolysis interacting with an external ATP-consuming reaction. The algorithm is based on the laws of evolution by natural selection, and may be viewed as a particular version of evolutionary algorithms. The following conclusions are derived: (a) evolutionary algorithms are very useful search strategies in determining optimal stoichiometries of metabolic pathways. (b) Essential topological features of the glycolytic network may be explained on the basis of flux optimization. (c) There is a strong interrelation between the optimal stoichiometries and the thermodynamic and kinetic properties of the participating reactions. (d) Some subsequences of reactions in optimal pathways are strongly conserved at variation of system parameters, which may be understood by applying principles of metabolic control analysis.     

Cysteine and iron accelerate the formation of ribose-5-phosphate,providing insights into the evolutionary origins of the metabolic network structure

The structure of the metabolic network is highly conserved, but we know little about its evolutionary origins. Key for explaining the early evolution of metabolism is solving a chicken–egg dilemma, which describes that enzymes are made from the very same molecules they produce. The recent discovery of several nonenzymatic reaction sequences that topologically resemble central metabolism has provided experimental support for a “metabolism first” theory, in which at least part of the extant metabolic network emerged on the basis of nonenzymatic reactions. But how could evolution kick-start on the basis of a metal catalyzed reaction sequence, and how could the structure of nonenzymatic reaction sequences be imprinted on the metabolic network to remain conserved for billions of years? We performed an in vitro screening where we add the simplest components of metabolic enzymes, proteinogenic amino acids, to a nonenzymatic, iron-driven reaction network that resembles glycolysis and the pentose phosphate pathway (PPP). We observe that the presence of the amino acids enhanced several of the nonenzymatic reactions. Particular attention was triggered by a reaction that resembles a rate-limiting step in the oxidative PPP. A prebiotically available, proteinogenic amino acid cysteine accelerated the formation of RNA nucleoside precursor ribose-5-phosphate from 6-phosphogluconate. We report that iron and cysteine interact and have additive effects on the reaction rate so that ribose-5-phosphate forms at high specificity under mild, metabolism typical temperature and environmental conditions. We speculate that accelerating effects of amino acids on rate-limiting nonenzymatic reactions could have facilitated a stepwise enzymatization of nonenzymatic reaction sequences, imprinting their structure on the evolving metabolic network.

The evolutionary origins of metabolism are largely unknown. This study shows that the prebiotically available proteinogenic amino acid cysteine can promote the metabolism-like rate-limiting formation of ribose-5-phosphate, suggesting that early metabolic pathways could have emerged thought the stepwise enzymatization of non-enzymatic reaction sequences.