Origin and evolution of metabolic pathways

The free energy principle induces intracellular compartmentalization

Living systems at all scales are compartmentalized into interacting subsystems. This paper reviews a mechanism that drives compartmentalization in generic systems at any scale. It first discusses three symmetries of generic physical interactions in a quantum-theoretic description. It then shows that if one of these, a permutation symmetry on the inter-system boundary, is spontaneously broken, the symmetry breaking is amplified by the Free Energy Principle (FEP). It thus shows how compartmentalization generically results from permutation symmetry breaking under the FEP. It finally notes that the FEP asymptotically restores the broken symmetry, showing that the FEP can be regarded as a theory of fluctuations away from a permutation-symmetric boundary, and hence from an entangled joint state of the interacting systems.

Chemotactic Interactions Drive Migration of Membraneless Active Droplets

In nature, chemotactic interactions are ubiquitous and play a critical role in driving the collective behavior of living organisms. Reproducing these interactions in vitro is still a paramount challenge due to the complexity of mimicking and controlling cellular features, such as tangled metabolic networks, cytosolic macromolecular crowding, and cellular migration, on a microorganism size scale. Here, we generate enzymatically active cell-sized droplets able to move freely, and by following a chemical gradient, able to interact with the surrounding droplets in a collective manner. The enzyme within the droplets generates a pH gradient that extends outside the edge of the droplets. We discovered that the external pH gradient triggers droplet migration and controls its directionality, which is selectively toward the neighboring droplets. Hence, by changing the enzyme activity inside the droplet, we tuned the droplet migration speed. Furthermore, we showed that these cellular-like features can facilitate the reconstitution of a simple and linear protometabolic pathway and increase the final reaction product generation. Our work suggests that simple and stable membraneless droplets can reproduce complex biological phenomena, opening new perspectives as bioinspired materials and synthetic biology tools.

The Sphinx and the egg: Evolutionary enigmas of the (glyco)sphingolipid biosynthetic pathway

In eukaryotes, the de novo synthesis of sphingolipids (SLs) consists of multiple sequential steps which are compartmentalized between the endoplasmic reticulum and the Golgi apparatus. Studies over many decades have identified the enzymes in the pathway, their localization, topology and an array of regulatory mechanisms. However, little is known about the evolutionary forces that underly the generation of this complex pathway or of its anteome, i.e., the metabolic pathways that converge on the SL biosynthetic pathway and are essential for its activity. After briefly describing the pathway, we discuss the mechanisms by which the enzymes of the SL biosynthetic pathway are targeted to their different subcellular locations, how the pathway per se may have evolved, including its compartmentalization, and the relationship of the pathway to eukaryogenesis. We discuss the circular interdependence of the evolution of the SL pathway, and comment on whether current Darwinian evolutionary models are able to provide genuine mechanistic insight into how the pathway came into being.

Metabolic symbiosis between oxygenated and hypoxic tumour cells: An agent-based modelling study

Deregulated metabolism is one of the hallmarks of cancer. It is well-known that tumour cells tend to metabolize glucose via glycolysis even when oxygen is available and mitochondrial respiration is functional. However, the lower energy efficiency of aerobic glycolysis with respect to mitochondrial respiration makes this behaviour, namely the Warburg effect, counter-intuitive, although it has now been recognized as source of anabolic precursors. On the other hand, there is evidence that oxygenated tumour cells could be fuelled by exogenous lactate produced from glycolysis. We employed a multi-scale approach that integrates multi-agent modelling, diffusion-reaction, stoichiometric equations, and Boolean networks to study metabolic cooperation between hypoxic and oxygenated cells exposed to varying oxygen, nutrient, and inhibitor concentrations. The results show that the cooperation reduces the depletion of environmental glucose, resulting in an overall advantage of using aerobic glycolysis. In addition, the oxygen level was found to be decreased by symbiosis, promoting a further shift towards anaerobic glycolysis. However, the oxygenated and hypoxic populations may gradually reach quasi-equilibrium. A sensitivity analysis using Latin hypercube sampling and partial rank correlation shows that the symbiotic dynamics depends on properties of the specific cell such as the minimum glucose level needed for glycolysis. Our results suggest that strategies that block glucose transporters may be more effective to reduce tumour growth than those blocking lactate intake transporters.

Evolution at the Origins of Life?

The role of evolutionary theory at the origin of life is an extensively debated topic. The origin and early development of life is usually separated into a prebiotic phase and a protocellular phase, ultimately leading to the Last Universal Common Ancestor. Most likely, the Last Universal Common Ancestor was subject to Darwinian evolution, but the question remains to what extent Darwinian evolution applies to the prebiotic and protocellular phases. In this review, we reflect on the current status of evolutionary theory in origins of life research by bringing together philosophy of science, evolutionary biology, and empirical research in the origins field. We explore the various ways in which evolutionary theory has been extended beyond biology; we look at how these extensions apply to the prebiotic development of (proto)metabolism; and we investigate how the terminology from evolutionary theory is currently being employed in state-of-the-art origins of life research. In doing so, we identify some of the current obstacles to an evolutionary account of the origins of life, as well as open up new avenues of research.

Family characteristics, phylogenetic reconstruction, and potential applications of the plant BAHD acyltransferase family

The BAHD acyltransferase family is a class of proteins in plants that can acylate a variety of primary and specialized secondary metabolites. The typically acylated products have greatly improved stability, lipid solubility, and bioavailability and thus show significant differences in their physicochemical properties and pharmacological activities. Here, we review the protein structure, catalytic mechanism, and phylogenetic reconstruction of plant BAHD acyltransferases to describe their family characteristics, acylation reactions, and the processes of potential functional differentiation. Moreover, the potential applications of the BAHD family in human activities are discussed from the perspectives of improving the quality of economic plants, enhancing the efficacy of medicinal plants, improving plant biomass for use in biofuel, and promoting stress resistance of land plants. This review provides a reference for the research and production of plant BAHD acyltransferases.

Mimicking the last step of gene elongation: hints from the bacterial hisF gene

Gene elongation consists in an in-tandem duplication of a gene and divergence and fusion of the two copies, resulting in a gene constituted by two divergent paralogous modules. Many present-day proteins show internal repeats of amino acid sequences, generated by gene elongation events; however, gene elongation is still a poorly studied evolutionary molecular mechanism. The most documented case is that of the histidine biosynthetic genes hisA and hisF, which derive from the gene elongation of an ancestral gene half the size of the extant ones. The aim of this work was to experimentally simulate the possible last step of the gene elongation event occurred during hisF gene evolution under selective pressure conditions. Azospirillum brasilense hisF gene, carrying a single nucleotide mutation that generates a stop codon between the two halves of the gene, was used to transform the histidine-auxotrophic Escherichia coli strain FB182 (hisF892). The transformed strain was subjected to selective pressure (i.e., low concentration/absence of histidine in the growth medium) and the obtained mutants were characterized. The restoration of prototrophy was strongly dependent on the time of incubation and on the strength of the selective pressure. The mutations involved the introduced stop codon with a single base substitution and none of the mutants restored the wild-type codon. Possible correlations between the different mutations and i) E. coli codon usage, ii) three-dimensional structures of the mutated HisF proteins, and iii) growth ability of the mutants were investigated. On the contrary, when the experiment was repeated by mutating a more conserved codon, only a synonymous substitution was obtained. Thus, experiments performed in this study allowed to mimic a possible gene elongation event occurred during the evolution of hisF gene, evidencing the ability of bacterial cells to modify their genome in short times under selective conditions.

On the evolution of natural product biosynthesis

Natural products are the raw material for drug discovery programmes. Bioactive natural products are used extensively in medicine and agriculture and have found utility as antibiotics, immunosuppressives, anti-cancer drugs and anthelminthics. Remarkably, the natural role and what mechanisms drive evolution of these molecules is relatively poorly understood. The exponential increase in genome and chemical data in recent years, coupled with technical advances in bioinformatics and genetics have enabled progress to be made in understanding the evolution of biosynthetic gene clusters and the products of their enzymatic machinery. Here we discuss the diversity of natural products, incorporating the mechanisms that govern evolution of metabolic pathways and how this can be applied to biosynthetic gene clusters. We build on the nomenclature of natural products in terms of primary, integrated, secondary and specialised metabolism and place this within an ecology-evolutionary-developmental biology framework. This eco-evo-devo framework we believe will help to clarify the nature and use of the term specialised metabolites in the future.

Predicting and improving the microbial removal of organic micropollutants during wastewater treatment: A review

Organic micropollutants (OMPs) consist of widely used chemicals such as pharmaceuticals and pesticides that can persist in surface and groundwaters at low concentrations (ng/L to μg/L) for a long time. The presence of OMPs in water can disrupt aquatic ecosystems and threaten the quality of drinking water sources. Wastewater treatment plants (WWTPs) rely on microorganisms to remove major nutrients from water, but their effectiveness at removing OMPs varies. Low removal efficiency might be the result of low concentrations, inherent stable chemical structures of OMPs, or suboptimal conditions in WWTPs. In this review, we discuss these factors, with special emphasis on the ongoing adaptation of microorganisms to degrade OMPs. Finally, recommendations are drawn to improve the prediction of OMP removal in WWTPs and to optimize the design of new microbial treatment strategies. OMP removal seems to be concentration-, compound-, and process-dependent, which poses a great complexity to develop accurate prediction models and effective microbial processes targeting all OMPs.

In silico analysis of sporulene biosynthesis pathway genes in the members of the class Bacilli

Sporulene, a pentacyclic triterpenoid, was discovered in Bacillus subtilis and is associated with bacterial endospores. However, the study was not further extended, leaving a trail of questions. One such question is what diversity of sporulenes exists among spore-forming members? Considering the sporulene biosynthesis pathway as a fundamental tool to survey the distribution of this molecule, a genome mining study was conducted. Mining for genes encoding putative proteins of sporulene biosynthesis pathway among the class Bacilli members revealed the presence of hepS, hepT, ytpB, and sqhC genes in the members of the family Bacillaceae, Caryophanaceae, Paenibacillaceae, and Sporolactobacillaceae. However, these genes were completely absent in the members of Staphylococcaceae, Lactobacillaceae, Aerococcaceae, Carnobacteriaceae, and Leuconostocaceae. Unlike other probable pathway related proteins, a conserved amino acid domain of putative terpenoid cyclase (YtpB) appeared deep-rooted among the genus Bacillus members. In-depth analysis showed the constant gene arrangement of hepS, hepT, ytpB, and sqhC genes in these members, there by demonstrating the conserved nature of sporulene biosynthesis pathway in the members of the genus Bacillus. Our study suggests confinement of the sporulene biosynthesis pathway to spore-forming members of the class Bacilli, majorly to the genus Bacillus.

The Operon as a Conundrum of Gene Dynamics and Biochemical Constraints: What We Have Learned from Histidine Biosynthesis

Operons represent one of the leading strategies of gene organization in prokaryotes, having a crucial influence on the regulation of gene expression and on bacterial chromosome organization. However, there is no consensus yet on why, how, and when operons are formed and conserved, and many different theories have been proposed. Histidine biosynthesis is a highly studied metabolic pathway, and many of the models suggested to explain operons origin and evolution can be applied to the histidine pathway, making this route an attractive model for the study of operon evolution. Indeed, the organization of his genes in operons can be due to a progressive clustering of biosynthetic genes during evolution, coupled with a horizontal transfer of these gene clusters. The necessity of physical interactions among the His enzymes could also have had a role in favoring gene closeness, of particular importance in extreme environmental conditions. In addition, the presence in this pathway of paralogous genes, heterodimeric enzymes and complex regulatory networks also support other operon evolution hypotheses. It is possible that histidine biosynthesis, and in general all bacterial operons, may result from a mixture of several models, being shaped by different forces and mechanisms during evolution.

Complexity and self-organization in the evolution of cell polarization

Cellular life exhibits order and complexity, which typically increase over the course of evolution. Cell polarization is a well-studied example of an ordering process that breaks the internal symmetry of a cell by establishing a preferential axis. Like many cellular processes, polarization is driven by self-organization, meaning that the macroscopic pattern emerges as a consequence of microscopic molecular interactions at the biophysical level. However, the role of self-organization in the evolution of complex protein networks remains obscure. In this Review, we provide an overview of the evolution of polarization as a self-organizing process, focusing on the model species Saccharomyces cerevisiae and its fungal relatives. Moreover, we use this model system to discuss how self-organization might relate to evolutionary change, offering a shift in perspective on evolution at the microscopic scale.

Application of Cloning-Free Genome Engineering to Escherichia coli

The propagation of foreign DNA in Escherichia coli is central to molecular biology. Recent advances have dramatically expanded the ability to engineer (bacterial) cells; however, most of these techniques remain time-consuming. The aim of the present work was to explore the possibility to use the cloning-free genome editing (CFGE) approach, proposed by Döhlemann and coworkers (2016), for E. coli genetics, and to deepen the knowledge about the homologous recombination mechanism. The E. coli auxotrophic mutant strains FB182 (hisF892) and FB181 (hisI903) were transformed with the circularized wild-type E. coli (i) hisF gene and hisF gene fragments of decreasing length, and (ii) hisIE gene, respectively. His⁺ clones were selected based on their ability to grow in the absence of histidine, and their hisF/hisIE gene sequences were characterized. CFGE method allowed the recombination of wild-type his genes (or fragments of them) within the mutated chromosomal copy, with a different recombination frequency based on the fragment length, and the generation of clones with a variable number of in tandem his genes copies. Data obtained pave the way to further evolutionary studies concerning the homologous recombination mechanism and the fate of in tandem duplicated genes.

Regulation and function of the mammalian tricarboxylic acid cycle

The tricarboxylic acid (TCA) cycle, otherwise known as the Krebs cycle, is a central metabolic pathway that performs the essential function of oxidizing nutrients to support cellular bioenergetics. More recently, it has become evident that TCA cycle behavior is dynamic, and products of the TCA cycle can be co-opted in cancer and other pathologic states. In this review, we revisit the TCA cycle, including its potential origins and the history of its discovery. We provide a detailed accounting of the requirements for sustained TCA cycle function and the critical regulatory nodes that can stimulate or constrain TCA cycle activity. We also discuss recent advances in our understanding of the flexibility of TCA cycle wiring and the increasingly appreciated heterogeneity in TCA cycle activity exhibited by mammalian cells. Deeper insight into how the TCA cycle can be differentially regulated and, consequently, configured in different contexts will shed light on how this pathway is primed to meet the requirements of distinct mammalian cell states.

On the Evolutionary History of the Twenty Encoded Amino Acids

α-Amino acids are essential molecular constituents of life, twenty of which are privileged because they are encoded by the ribosomal machinery. The question remains open as to why this number and why this 20 in particular, an almost philosophical question that cannot be conclusively resolved. They are closely related to the evolution of the genetic code and whether nucleic acids, amino acids, and peptides appeared simultaneously and were available under prebiotic conditions when the first self-sufficient complex molecular system emerged on Earth. This report focuses on prebiotic and metabolic aspects of amino acids and proteins starting with meteorites, followed by their formation, including peptides, under plausible prebiotic conditions, and the major biosynthetic pathways in the various kingdoms of life. Coenzymes play a key role in the present analysis in that amino acid metabolism is linked to glycolysis and different variants of the tricarboxylic acid cycle (TCA, rTCA, and the incomplete horseshoe version) as well as the biosynthesis of the most important coenzymes. Thus, the report opens additional perspectives and facets on the molecular evolution of primary metabolism.

The Azurin Coding Gene: Origin and Phylogenetic Distribution

Azurin is a bacterial-derived cupredoxin, which is mainly involved in electron transport reactions. Interest in azurin protein has risen in recent years due to its anticancer activity and its possible applications in anticancer therapies. Nevertheless, the attention of the scientific community only focused on the azurin protein found in Pseudomonas aeruginosa (Proteobacteria, Gammaproteobacteria). In this work, we performed the first comprehensive screening of all the bacterial genomes available in online repositories to assess azurin distribution in the three domains of life. The Azurin coding gene was not detected in the domains Archaea and Eucarya, whereas it was detected in phyla other than Proteobacteria, such as Bacteroidetes, Verrucomicrobia and Chloroflexi, and a phylogenetic analysis of the retrieved sequences was performed. Observed patchy distribution and phylogenetic data suggest that once it appeared in the bacterial domain, the azurin coding gene was lost in several bacterial phyla and/or anciently horizontally transferred between different phyla, even though a vertical inheritance appeared to be the major force driving the transmission of this gene. Interestingly, a shared conserved domain has been found among azurin members of all the investigated phyla. This domain is already known in P. aeruginosa as p28 domain and its importance for azurin anticancer activity has been widely explored. These findings may open a new and intriguing perspective in deciphering the azurin anticancer mechanisms and to develop new tools for treating cancer diseases.

Evolutionary Aspects of the Oxido-Reductive Network of Methylglyoxal

In the chemoautotrophic theory for the origin of life, offered as an alternative to broth theory, the archaic reductive citric acid cycle operating without enzymes is in the center. The non-enzymatic (methyl)glyoxalase pathway has been suggested to be the anaplerotic route for the reductive citric acid cycle. In the recent years, much has been learned about methylglyoxal, but its importance in the metabolic machinery is still uncovered. If methylglyoxal had been essential participant of the early stage of evolution, then it is a legitimate question whether it might have played a role in the early oxido-reduction network, too. Therefore, an oxido-reduction network of methylglyoxal that might have functioned under ancient circumstances without enzymes was constructed and analyzed by virtue of group contribution method. Taking methylglyoxal as input material, it turned out that the evolutionary value of reactions and biomolecules were not similar. Glycerol, glycerate, and tartonate, the output components, were conserved to different degrees. Although the tartonate route was similarly favorable from energetic point of view, its intermediates are almost not present in extant biochemistry. The presence of two carboxyl or aldehyde groups, or their combination in tricarbons of the constructed network seemed disadvantageous for selection, and the inductive effect, resulting in an asymmetry in electron cloud of chemicals, might have been important. The evolutionary role for cysteine, H₂S, and formaldehyde in the emergence of high-energy bonds in the form of thioesters and in Fe-S cluster formation as well as in imidazole synthesis was shown to bridge the gap between prebiotic chemistry and contemporary biochemistry. Overall, the ideas developed here represent an approach fitting to chemoautotrophic origin of life and implying to the role of methylglyoxal in triose formation. The proposed network is expected to have an impact upon how one may think of prebiological chemical processes on methylglyoxal, too. Finally, along the evolutionary time line, the network functioning without enzymes is situated between the formation of simple organic compounds and primeval cells, being closer to the former and well preceding the last common metabolic ancestor developed after primitive cells emerged.

The Role of Gene Duplication in the Divergence of Enzyme Function: A Comparative Approach

Gene duplication is a crucial process involved in the appearance of new genes and functions. It is thought to have played a major role in the growth of enzyme families and the expansion of metabolism at the biosphere's dawn and in recent times. Here, we analyzed paralogous enzyme content within each of the seven enzymatic classes for a representative sample of prokaryotes by a comparative approach. We found a high ratio of paralogs for three enzymatic classes: oxidoreductases, isomerases, and translocases, and within each of them, most of the paralogs belong to only a few subclasses. Our results suggest an intricate scenario for the evolution of prokaryotic enzymes, involving different fates for duplicated enzymes fixed in the genome, where around 20-40% of prokaryotic enzymes have paralogs. Intracellular organisms have a lesser ratio of duplicated enzymes, whereas free-living enzymes show the highest ratios. We also found that phylogenetically close phyla and some unrelated but with the same lifestyle share similar genomic and biochemical traits, which ultimately support the idea that gene duplication is associated with environmental adaptation.

The Histidine Biosynthetic Genes in the Superphylum Bacteroidota-Rhodothermota-Balneolota-Chlorobiota: Insights into the Evolution of Gene Structure and Organization

One of the most studied metabolic routes is the biosynthesis of histidine, especially in enterobacteria where a single compact operon composed of eight adjacent genes encodes the complete set of biosynthetic enzymes. It is still not clear how his genes were organized in the genome of the last universal common ancestor community. The aim of this work was to analyze the structure, organization, phylogenetic distribution, and degree of horizontal gene transfer (HGT) of his genes in the Bacteroidota-Rhodothermota-Balneolota-Chlorobiota superphylum, a group of phylogenetically close bacteria with different surviving strategies. The analysis of the large variety of his gene structures and organizations revealed different scenarios with genes organized in more or less compact—heterogeneous or homogeneous—operons, in suboperons, or in regulons. The organization of his genes in the extant members of the superphylum suggests that in the common ancestor of this group, genes were scattered throughout the chromosome and that different forces have driven the assembly of his genes in compact operons. Gene fusion events and/or paralog formation, HGT of single genes or entire operons between strains of the same or different taxonomic groups, and other molecular rearrangements shaped the his gene structure in this superphylum.

Recent Research Progress in Taxol Biosynthetic Pathway and Acylation Reactions Mediated by Taxus Acyltransferases

Taxol is one of the most effective anticancer drugs in the world that is widely used in the treatments of breast, lung and ovarian cancer. The elucidation of the taxol biosynthetic pathway is the key to solve the problem of taxol supply. So far, the taxol biosynthetic pathway has been reported to require an estimated 20 steps of enzymatic reactions, and sixteen enzymes involved in the taxol pathway have been well characterized, including a novel taxane-10β-hydroxylase (T10βOH) and a newly putative β-phenylalanyl-CoA ligase (PCL). Moreover, the source and formation of the taxane core and the details of the downstream synthetic pathway have been basically depicted, while the modification of the core taxane skeleton has not been fully reported, mainly concerning the developments from diol intermediates to 2-debenzoyltaxane. The acylation reaction mediated by specialized Taxus BAHD family acyltransferases (ACTs) is recognized as one of the most important steps in the modification of core taxane skeleton that contribute to the increase of taxol yield. Recently, the influence of acylation on the functional and structural diversity of taxanes has also been continuously revealed. This review summarizes the latest research advances of the taxol biosynthetic pathway and systematically discusses the acylation reactions supported by Taxus ACTs. The underlying mechanism could improve the understanding of taxol biosynthesis, and provide a theoretical basis for the mass production of taxol.

The Semi-Enzymatic Origin of Metabolic Pathways: Inferring a Very Early Stage of the Evolution of Life

The early evolution of life is a period with many important events with a lot of big and open questions. One of them is the evolution of metabolic pathways, which means the origin and assembly of enzymes that act together. The retrograde hypothesis was the first attempt to explain the origin and evolutionary history of metabolic pathways; Norman Horowitz developed this first significant hypothesis. This idea was followed by relevant proposals developed by Sam Granick, who proposed the "forward direction hypothesis," and then the successful idea of "Patchwork" assembly proposed independently by Martynas Yčas and Roy Jensen. Since then, a few new hypotheses were proposed; one of the most influential was made by Antonio Lazcano and Stanley Miller in the Journal Molecular Evolution, the "semi-enzymatic origin" of metabolic pathways. This article was cited more than 160 times, including in most papers published about the early evolution of metabolism, placing it as influential work in the field. The ideas proposed in this work and their effects on studying the origin and early evolution of life are analyzed.

Tracing molecular properties throughout evolution: A chemoinformatic approach

Evolution of metabolism is a longstanding yet unresolved question, and several hypotheses were proposed to address this complex process from a Darwinian point of view. Modern statistical bioinformatic approaches targeted to the comparative analysis of genomes are being used to detect signatures of natural selection at the gene and population level, as an attempt to understand the origin of primordial metabolism and its expansion. These studies, however, are still mainly centered on genes and the proteins they encode, somehow neglecting the small organic chemicals that support life processes. In this work, we selected steroids as an ancient family of metabolites widely distributed in all eukaryotes and applied unsupervised machine learning techniques to reveal the traits that natural selection has imprinted on molecular properties throughout the evolutionary process. Our results clearly show that sterols, the primal steroids that first appeared, have more conserved properties and that, from then on, more complex compounds with increasingly diverse properties have emerged, suggesting that chemical diversification parallels the expansion of biological complexity. In a wider context, these findings highlight the worth of chemoinformatic approaches to a better understanding the evolution of metabolism.

Ancestral zinc-finger bearing protein MucR in alpha-proteobacteria: A novel xenogeneic silencer?

The MucR/Ros family protein is conserved in alpha-proteobacteria and characterized by its zinc-finger motif that has been proposed as the ancestral domain from which the eukaryotic C2H2 zinc-finger structure evolved. In the past decades, accumulated evidences have revealed MucR as a pleiotropic transcriptional regulator that integrating multiple functions such as virulence, symbiosis, cell cycle and various physiological processes. Scattered reports indicate that MucR mainly acts as a repressor, through oligomerization and binding to multiple sites of AT-rich target promoters. The N-terminal region and zinc-finger bearing C-terminal region of MucR mediate oligomerization and DNA-binding, respectively. These features are convergent to those of xenogeneic silencers such as H-NS, MvaT, Lsr2 and Rok, which are mainly found in other lineages. Phylogenetic analysis of MucR homologs suggests an ancestral origin of MucR in alpha- and delta-proteobacteria. Multiple independent duplication and lateral gene transfer events contribute to the diversity and phyletic distribution of MucR. Finally, we posed questions which remain unexplored regarding the putative roles of MucR as a xenogeneic silencer and a general manager in balancing adaptation and regulatory integration in the pangenome context.

Evolution of linkage and genome expansion in protocells: The origin of chromosomes

Chromosomes are likely to have assembled from unlinked genes in early evolution. Genetic linkage reduces the assortment load and intragenomic conflict in reproducing protocell models to the extent that chromosomes can go to fixation even if chromosomes suffer from a replicative disadvantage, relative to unlinked genes, proportional to their length. Here we numerically show that chromosomes spread within protocells even if recurrent deleterious mutations affecting replicating genes (as ribozymes) are considered. Dosage effect selects for optimal genomic composition within protocells that carries over to the genic composition of emerging chromosomes. Lacking an accurate segregation mechanism, protocells continue to benefit from the stochastic corrector principle (group selection of early replicators), but now at the chromosome level. A remarkable feature of this process is the appearance of multigene families (in optimal genic proportions) on chromosomes. An added benefit of chromosome formation is an increase in the selectively maintainable genome size (number of different genes), primarily due to the marked reduction of the assortment load. The establishment of chromosomes is under strong positive selection in protocells harboring unlinked genes. The error threshold of replication is raised to higher genome size by linkage due to the fact that deleterious mutations affecting protocells metabolism (hence fitness) show antagonistic (diminishing return) epistasis. This result strengthens the established benefit conferred by chromosomes on protocells allowing for the fixation of highly specific and efficient enzymes.

The emergence of chromosomes harboring several genes is a crucial ingredient of the major evolutionary transition from naked replicators to cells. Linkage of replicating genes reduces conflict between them and alleviates the problem of chance loss of genes upon stochastic protocell fission. The emerging organization of protocells maintaining several segregating chromosomes with balanced gene composition also allows for an increase in the number of gene types despite recurrent deleterious mutations. We suggest that this interim genomic organization enabled protocells to evolve specific and efficient enzymes and paved the way toward an accurate mechanism for chromosome segregation later in evolution.

Exploring the role of the histidine biosynthetic hisF gene in cellular metabolism and in the evolution of (ancestral) genes: from LUCA to the extant (micro)organisms

Histidine biosynthesis is an ancestral pathway that was assembled before the appearance of the Last Universal Common Ancestor; afterwards, it remained unaltered in all the extant histidine-synthesizing (micro)organisms. It is a metabolic cross-road interconnecting histidine biosynthesis to nitrogen metabolism and the de novo synthesis of purines. This interconnection is due to the reaction catalyzed by the products of hisH and hisF genes. The latter gene is an excellent model to study which trajectories have been followed by primordial cells to build the first metabolic routes, since its evolution is the result of different molecular rearrangement events, i.e. gene duplication, gene fusion, gene elongation, and domain shuffling. Additionally, this review summarizes data concerning the involvement of hisF and its product in other different cellular processes, revealing that HisF very likely plays a role also in cell division control and involvement in virulence and nodule development in different bacteria. From the metabolic viewpoint, these results suggest that HisF plays a central role in cellular metabolism, highlighting the interconnections of different metabolic pathways.

The Role of Gene Elongation in the Evolution of Histidine Biosynthetic Genes

Gene elongation is a molecular mechanism consisting of an in-tandem duplication of a gene and divergence and fusion of the two copies, resulting in a gene constituted by two divergent paralogous modules. The aim of this work was to evaluate the importance of gene elongation in the evolution of histidine biosynthetic genes and to propose a possible evolutionary model for some of them. Concerning the genes hisA and hisF, which code for two homologous (β/α)₈-barrels, it has been proposed that the two extant genes could be the result of a cascade of gene elongation/domain shuffling events starting from an ancestor gene coding for just one (β/α) module. A gene elongation event has also been proposed for the evolution of hisB and hisD; structural analyses revealed the possibility of an early elongation event, resulting in the repetition of modules. Furthermore, it is quite possible that the gene elongations responsible for the evolution of the four proteins occurred before the earliest phylogenetic divergence. In conclusion, gene elongation events seem to have played a crucial role in the evolution of the histidine biosynthetic pathway, and they may have shaped the structures of many genes during the first steps of their evolution.

The evolution of metabolism: How to test evolutionary hypotheses at the genomic level

The origin of primordial metabolism and its expansion to form the metabolic networks extant today represent excellent systems to study the impact of natural selection and the potential adaptive role of novel compounds. Here we present the current hypotheses made on the origin of life and ancestral metabolism and present the theories and mechanisms by which the large chemical diversity of plants might have emerged along evolution. In particular, we provide a survey of statistical methods that can be used to detect signatures of selection at the gene and population level, and discuss potential and limits of these methods for investigating patterns of molecular adaptation in plant metabolism.

Conserved genomic neighborhood is a strong but no perfect indicator for a direct interaction of microbial gene products

BACKGROUND

The order of genes in bacterial genomes is not random; for example, the products of genes belonging to an operon work together in the same pathway. The cotranslational assembly of protein complexes is deemed to conserve genomic neighborhoods even stronger than a common function. This is why a conserved genomic neighborhood can be utilized to predict, whether gene products form protein complexes.

RESULTS

We were interested to assess the performance of a neighborhood-based classifier that analyzes a large number of genomes. Thus, we determined for the genes encoding the subunits of 494 experimentally verified hetero-dimers their local genomic context. In order to generate phylogenetically comprehensive genomic neighborhoods, we utilized the tools offered by the Enzyme Function Initiative. For each subunit, a sequence similarity network was generated and the corresponding genome neighborhood network was analyzed to deduce the most frequent gene product. This was predicted as interaction partner, if its abundance exceeded a threshold, which was the frequency giving rise to the maximal Matthews correlation coefficient. For the threshold of 16%, the true positive rate was 45%, the false positive rate 0.06%, and the precision 55%. For approximately 20% of the subunits, the interaction partner was not found in a neighborhood of ± 10 genes.

CONCLUSIONS

Our phylogenetically comprehensive analysis confirmed that complex formation is a strong evolutionary factor that conserves genome neighborhoods. On the other hand, for 55% of the cases analyzed here, classification failed. Either, the interaction partner was not present in a ± 10 gene window or was not the most frequent gene product.

Membranes Composed of Lipopeptides and Liponucleobases Inspired Protolife Evolution

Amino acids and peptides have been demonstrated to form lipoamino acids and lipopeptides under presumed prebiotic conditions, and readily form liposomes. Of the common nucleobases, adenine forms a liponucleobase even below 100 °C. Adenine as well as other nucleobases can also be derivatized with ethylene carbonate (and likely other similar compounds) onto which fatty acids can be attached. The fatty acid tails along with appropriately functionalized nucleobases provide some solubility of liponucleobases in membranes. Such membranes would provide a structure in which three of biology's major components are closely associated and available for chemical interactions. Nucleobase-to-nucleobase interactions would ensure that the liponucleobases would have a uniquely different head-group relationship than other amphiphiles within a membrane, likely forming rafts due their π-π interactions and providing surface discontinuities that could serve as catalytic sites. The π-π bond distance in aromatic compounds is typically 0.34 nm, commensurate with that of the amine to carboxylate distance in alpha amino acids. This would have provided opportunity for hydrogen bonding between amino acids and the distal primary amines or tautomeric carbonyl/hydroxyl groups of two π-bonded nucleobases. Such bonding would weaken the covalent linkages within the amino acids, making them susceptible to forming peptide bonds with an adjacent amino acid, likely a lipoamino acid or lipopeptide. Were this second lipoamino acid bound to a third π-bonded nucleobase, it could result in orientation, destabilization and peptide formation. The stacked triplet of nucleobases might constitute the primordial codon triplet from which peptides were synthesized: primordial translation.

Evolutionary Developments in Plant Specialized Metabolism, Exemplified by Two Transferase Families

Plant specialized metabolism emerged from the land colonization by ancient plants, becoming diversified along with plant evolution. To date, more than 1 million metabolites have been predicted to exist in the plant kingdom, and their metabolic processes have been revealed on the molecular level. Previous studies have reported that rates of evolution are greater for genes involved in plant specialized metabolism than in primary metabolism. This perspective introduces topics on the enigmatic molecular evolution of some plant specialized metabolic processes. Two transferase families, BAHD acyltransferases and aromatic prenyltransferases, which are involved in the biosynthesis of paclitaxel and meroterpenes, respectively, have shown apparent expansion. The latter family has been shown to beinvolved in the biosynthesis of a variety of aromatic substances, including prenylated coumarins in citrus plants and shikonin in Lithospermum erythrorhizon. These genes have evolved in the development of each special subfamily within the plant lineage. The broadness of substrate specificity and the exon-intron structure of their genes may provide hints to explain the evolutionary process underlying chemodiversity in plants.

Sensing between reactions - how the metabolic microenvironment shapes immunity

Perception of potential threat is key for survival. The immune system constantly patrols the organism scanning for potential pathogenic or malignant danger. Recent evidence suggests that immunosurveillance not only relies on classic receptors [e.g. Toll-like receptors (TLRs) or antibodies] but is also based on sensing of the metabolic environment. Metabolites interact in numerous ways with immune cells, and are therefore more than just reaction intermediates. This new perspective opens the door for potential, future therapeutic strategies. Here we describe how immune functionality during infections, cancer or autoimmunity, as exemplified by short-chain fatty acids, lactate and reactive oxygen species (ROS), can be shaped by metabolic intermediates.

Selection for novel metabolic capabilities in Salmonella enterica

Bacteria are known to display extensive metabolic diversity and many studies have shown that they can use an extensive repertoire of small molecules as carbon‐ and energy sources. However, it is less clear to what extent a bacterium can expand its existing metabolic capabilities by acquiring mutations that, for example, rewire its metabolic pathways. To investigate this capability and potential for evolution of novel phenotypes, we sampled large populations of mutagenized Salmonella enterica to select very rare mutants that can grow on minimal media containing 124 low molecular weight compounds as sole carbon sources. We found mutants growing on 18 of these novel carbon sources, and identified the causal mutations that allowed growth for four of them. Mutations that relieve physiological constraints or increase expression of existing pathways were found to be important contributors to the novel phenotypes. For the remaining 14 novel phenotypes, whole genome sequencing of independent mutants and genetic analysis suggested that these novel metabolic phenotypes result from a combination of multiple mutations. This work, by virtue of identifying the genetic and mechanistic basis for new metabolic capabilities, sheds light on the properties of adaptive landscapes underlying the evolution of novel phenotypes.

The Diversity of Nutritional Metabolites: Origin, Dissection, and Application in Crop Breeding

The chemical diversity of plants is very high, and plant-based foods provide almost all the nutrients necessary for human health, either directly or indirectly. With advancements in plant metabolomics studies, the concept of nutritional metabolites has been expanded and updated. Because the concentration of many nutrients is usually low in plant-based foods, especially those from crops, metabolome-assisted breeding techniques using molecular markers associated with the synthesis of nutritional metabolites have been developed and used to improve nutritional quality of crops. Here, we review the origins of the diversity of nutrient metabolites from a genomic perspective and the role of gene duplication and divergence. In addition, we systematically review recent advances in the metabolomic and genetic basis of metabolite production in major crops. With the development of genome sequencing and metabolic detection technologies, multi-omic integrative analysis of genomes, transcriptomes, and metabolomes has greatly facilitated the deciphering of the genetic basis of metabolic pathways and the diversity of nutrient metabolites. Finally, we summarize the application of nutrient diversity in crop breeding and discuss the future development of a viable alternative to metabolome-assisted breeding techniques that can be used to improve crop nutrient quality.

Influence of pathway topology and functional class on the molecular evolution of human metabolic genes

Metabolic networks comprise thousands of enzymatic reactions functioning in a controlled manner and have been shaped by natural selection. Thanks to the genome data, the footprints of adaptive (positive) selection are detectable, and the strength of purifying selection can be measured. This has made possible to know where, in the metabolic network, adaptive selection has acted and where purifying selection is more or less strong and efficient. We have carried out a comprehensive molecular evolutionary study of all the genes involved in the human metabolism. We investigated the type and strength of the selective pressures that acted on the enzyme-coding genes belonging to metabolic pathways during the divergence of primates and rodents. Then, we related those selective pressures to the functional and topological characteristics of the pathways. We have used DNA sequences of all enzymes (956) of the metabolic pathways comprised in the HumanCyc database, using genome data for humans and five other mammalian species. We have found that the evolution of metabolic genes is primarily constrained by the layer of the metabolism in which the genes participate: while genes encoding enzymes of the inner core of metabolism are much conserved, those encoding enzymes participating in the outer layer, mediating the interaction with the environment, are evolutionarily less constrained and more plastic, having experienced faster functional evolution. Genes that have been targeted by adaptive selection are endowed by higher out-degree centralities than non-adaptive genes, while genes with high in-degree centralities are under stronger purifying selection. When the position along the pathway is considered, a funnel-like distribution of the strength of the purifying selection is found. Genes at bottom positions are highly preserved by purifying selection, whereas genes at top positions, catalyzing the first steps, are open to evolutionary changes. These results show how functional and topological characteristics of metabolic pathways contribute to shape the patterns of evolutionary pressures driven by natural selection and how pathway network structure matters in the evolutionary process that shapes the evolution of the system.

Metabolic models and gene essentiality data reveal essential and conserved metabolism in prokaryotes

Essential metabolic reactions are shaping constituents of metabolic networks, enabling viable and distinct phenotypes across diverse life forms. Here we analyse and compare modelling predictions of essential metabolic functions with experimental data and thereby identify core metabolic pathways in prokaryotes. Simulations of 15 manually curated genome-scale metabolic models were integrated with 36 large-scale gene essentiality datasets encompassing a wide variety of species of bacteria and archaea. Conservation of metabolic genes was estimated by analysing 79 representative genomes from all the branches of the prokaryotic tree of life. We find that essentiality patterns reflect phylogenetic relations both for modelling and experimental data, which correlate highly at the pathway level. Genes that are essential for several species tend to be highly conserved as opposed to non-essential genes which may be conserved or not. The tRNA-charging module is highlighted as ancestral and with high centrality in the networks, followed closely by cofactor metabolism, pointing to an early information processing system supplied by organic cofactors. The results, which point to model improvements and also indicate faults in the experimental data, should be relevant to the study of centrality in metabolic networks and ancient metabolism but also to metabolic engineering with prokaryotes.

If we tried to list every known chemical reaction within an organism–human, plant or even bacteria–we would get quite a long and confusing read. But when this information is represented in so-called genome-scale metabolic networks, we have the means to access computationally each of those reactions and their interconnections. Some parts of the network have alternatives, while others are unique and therefore can be essential for growth. Here, we simulate growth and compare essential reactions and genes for the simplest type of unicellular species–prokaryotes–to understand which parts of their metabolism are universally essential and potentially ancestral. We show that similar patterns of essential reactions echo phylogenetic relationships (this makes sense, as the genome provides the building plan for the enzymes that perform those reactions). Our computational predictions correlate strongly with experimental essentiality data. Finally, we show that a crucial step of protein synthesis (tRNA charging) and the synthesis and transformation of small molecules that enzymes require (cofactors) are the most essential and conserved parts of metabolism in prokaryotes. Our results are a step further in understanding the biology and evolution of prokaryotes but can also be relevant in applied studies including metabolic engineering and antibiotic design.

Scots pine aminopropyltransferases shed new light on evolution of the polyamine biosynthesis pathway in seed plants

Background and Aims

Polyamines are small metabolites present in all living cells and play fundamental roles in numerous physiological events in plants. The aminopropyltransferases (APTs), spermidine synthase (SPDS), spermine synthase (SPMS) and thermospermine synthase (ACL5), are essential enzymes in the polyamine biosynthesis pathway. In angiosperms, SPMS has evolved from SPDS via gene duplication, whereas in gymnosperms APTs are mostly unexplored and no SPMS gene has been reported. The present study aimed to investigate the functional properties of the SPDS and ACL5 proteins of Scots pine (Pinus sylvestris L.) in order to elucidate the role and evolution of APTs in higher plants.

Methods

Germinating Scots pine seeds and seedlings were analysed for polyamines by high-performance liquid chromatography (HPLC) and the expression of PsSPDS and PsACL5 genes by in situ hybridization. Recombinant proteins of PsSPDS and PsACL5 were produced and investigated for functional properties. Also gene structures, promoter regions and phylogenetic relationships of PsSPDS and PsACL5 genes were analysed.

Key Results

Scots pine tissues were found to contain spermidine, spermine and thermospermine. PsSPDS enzyme catalysed synthesis of both spermidine and spermine. PsACL5 was found to produce thermospermine, and PsACL5 gene expression was localized in the developing procambium in embryos and tracheary elements in seedlings.

Conclusions

Contrary to previous views, our results demonstrate that SPMS activity is not a novel feature developed solely in the angiosperm lineage of seed plants but also exists as a secondary property in the Scots pine SPDS enzyme. The discovery of bifunctional SPDS from an evolutionarily old conifer reveals the missing link in the evolution of the polyamine biosynthesis pathway. The finding emphasizes the importance of pre-existing secondary functions in the evolution of new enzyme activities via gene duplication. Our results also associate PsACL5 with the development of vascular structures in Scots pine.

BioPAXViz: a cytoscape application for the visual exploration of metabolic pathway evolution

Summary

BioPAXViz is a Cytoscape (version 3) application, providing a comprehensive framework for metabolic pathway visualization. Beyond the basic parsing, viewing and browsing roles, the main novel function that BioPAXViz provides is a visual comparative analysis of metabolic pathway topologies across pre-computed pathway phylogenomic profiles given a species phylogeny. Furthermore, BioPAXViz supports the display of hierarchical trees that allow efficient navigation through sets of variants of a single reference pathway. Thus, BioPAXViz can significantly facilitate, and contribute to, the study of metabolic pathway evolution and engineering.

Availability and Implementation

BioPAXViz has been developed as a Cytoscape app and is available at: https://github.com/CGU-CERTH/BioPAX.Viz. The software is distributed under the MIT License and is accompanied by example files and data. Additional documentation is available at the aforementioned GitHub repository.

Contact

ouzounis@certh.gr.

Whither life? Conjectures on the future evolution of biochemistry

Life has existed on the Earth for approximately four billion years. The sheer depth of evolutionary time, and the diversity of extant species, makes it tempting to assume that all the key biochemical innovations underpinning life have already happened. But we are only a little over halfway through the trajectory of life on our planet. In this Opinion piece, we argue: (i) that sufficient time remains for the evolution of new processes at the heart of metabolic biochemistry and (ii) that synthetic biology is providing predictive insights into the nature of these innovations. By way of example, we focus on engineered solutions to existing inefficiencies in energy generation, and on the complex, synthetic regulatory circuits that are currently being implemented.

A two-step evolutionary process establishes a non-native vitamin B6 pathway in Bacillus subtilis

Pyridoxal 5'-phosphate (PLP), the most important form of vitamin B6 serves as a cofactor for many proteins. Two alternative pathways for de novo PLP biosynthesis are known: the short deoxy-xylulose-5-phosphate (DXP)-independent pathway, which is present in the Gram-positive model bacterium Bacillus subtilis and the longer DXP-dependent pathway, which has been intensively studied in the Gram-negative model bacterium Escherichia coli. Previous studies revealed that bacteria contain many promiscuous enzymes causing a so-called 'underground metabolism', which can be important for the evolution of novel pathways. Here, we evaluated the potential of B. subtilis to use a truncated non-native DXP-dependent PLP pathway from E. coli for PLP synthesis. Adaptive laboratory evolution experiments revealed that two non-native enzymes catalysing the last steps of the DXP-dependent PLP pathway and two genomic alterations are sufficient to allow growth of vitamin B6 auxotrophic bacteria as rapid as the wild type. Thus, the existence of an underground metabolism in B. subtilis facilitates the generation of a pathway for synthesis of PLP using parts of a non-native vitamin B6 pathway. The introduction of non-native enzymes into a metabolic network and rewiring of native metabolism could be helpful to generate pathways that might be optimized for producing valuable substances.

Molecular evolution of acetohydroxyacid synthase in bacteria

Acetohydroxyacid synthase (AHAS) is the key enzyme in the biosynthetic pathways of branched chain amino acids in bacteria. Since it does not exist in animal and plant cells, AHAS is an attractive target for developing antimicrobials and herbicides. In some bacteria, there is a single copy of AHAS, while in others there are multiple copies. Therefore, it is necessary to investigate the origin and evolutionary pathway of various AHASs in bacteria. In this study, all the available protein sequences of AHAS in bacteria were investigated, and an evolutionary model of AHAS in bacteria is proposed, according to gene structure, organization and phylogeny. Multiple copies of AHAS in some bacteria might be evolved from the single copy of AHAS, the ancestor. Gene duplication, domain deletion and horizontal gene transfer might occur during the evolution of this enzyme. The results show the biological significance of AHAS, help to understand the functions of various AHASs in bacteria, and would be useful for developing industrial production strains of branched chain amino acids or novel antimicrobials.

Metabolomics analysis: Finding out metabolic building blocks

In this paper we propose a new methodology for the analysis of metabolic networks. We use the notion of strongly connected components of a graph, called in this context metabolic building blocks. Every strongly connected component is contracted to a single node in such a way that the resulting graph is a directed acyclic graph, called a metabolic DAG, with a considerably reduced number of nodes. The property of being a directed acyclic graph brings out a background graph topology that reveals the connectivity of the metabolic network, as well as bridges, isolated nodes and cut nodes. Altogether, it becomes a key information for the discovery of functional metabolic relations. Our methodology has been applied to the glycolysis and the purine metabolic pathways for all organisms in the KEGG database, although it is general enough to work on any database. As expected, using the metabolic DAGs formalism, a considerable reduction on the size of the metabolic networks has been obtained, specially in the case of the purine pathway due to its relative larger size. As a proof of concept, from the information captured by a metabolic DAG and its corresponding metabolic building blocks, we obtain the core of the glycolysis pathway and the core of the purine metabolism pathway and detect some essential metabolic building blocks that reveal the key reactions in both pathways. Finally, the application of our methodology to the glycolysis pathway and the purine metabolism pathway reproduce the tree of life for the whole set of the organisms represented in the KEGG database which supports the utility of this research.

Origin and evolution of transporter substrate specificity within the NPF family

Despite vast diversity in metabolites and the matching substrate specificity of their transporters, little is known about how evolution of transporter substrate specificities is linked to emergence of substrates via evolution of biosynthetic pathways. Transporter specificity towards the recently evolved glucosinolates characteristic of Brassicales is shown to evolve prior to emergence of glucosinolate biosynthesis. Furthermore, we show that glucosinolate transporters belonging to the ubiquitous NRT1/PTR FAMILY (NPF) likely evolved from transporters of the ancestral cyanogenic glucosides found across more than 2500 species outside of the Brassicales. Biochemical characterization of orthologs along the phylogenetic lineage from cassava to A. thaliana, suggests that alterations in the electrogenicity of the transporters accompanied changes in substrate specificity. Linking the evolutionary path of transporter substrate specificities to that of the biosynthetic pathways, exemplify how transporter substrate specificities originate and evolve as new biosynthesis pathways emerge.

Investigating Evolutionary Dynamics of RHA1 Operons

Grouping genes as operons is an important genomic feature of prokaryotic organisms. The comprehensive understanding of the operon organizations would be helpful to decipher transcriptional mechanisms, cellular pathways, and the evolutionary landscape of prokaryotic genomes. Although thousands of prokaryotes have been sequenced, genome-wide investigation of the evolutionary dynamics (division and recombination) of operons among these genomes remains unexplored. Here, we systematically analyzed the operon dynamics of Rhodococcus jostii RHA1 (RHA1), an oleaginous bacterium with high potential applications in biofuel, by comparing 340 prokaryotic genomes that were carefully selected from different genera. Interestingly, 99% of RHA1 operons were observed to exhibit evolutionary events of division and recombination among the 340 compared genomes. An operon that encodes all enzymes related to histidine biosynthesis in RHA1 (His-operon) was found to be segmented into smaller gene groups (sub-operons) in diverse genomes. These sub-operons were further reorganized with different functional genes as novel operons that are related to different biochemical processes. Comparatively, the operons involved in the functional categories of lipid transport and metabolism are relatively conserved among the 340 compared genomes. At the pathway level, RHA1 operons found to be significantly conserved were involved in ribosome synthesis, oxidative phosphorylation, and fatty acid synthesis. These analyses provide evolutionary insights of operon organization and the dynamic associations of various biochemical pathways in different prokaryotes.

Evolution of initiator tRNAs and selection of methionine as the initiating amino acid

Transfer RNAs (tRNAs) have been important in shaping biomolecular evolution. Initiator tRNAs (tRNAi), a special class of tRNAs, carry methionine (or its derivative, formyl-methionine) to ribosomes to start an enormously energy consuming but a highly regulated process of protein synthesis. The processes of tRNAi evolution, and selection of methionine as the universal initiating amino acid remain an enigmatic problem. We constructed phylogenetic trees using the whole sequence, the acceptor-TψC arm ('minihelix'), and the anticodon-dihydrouridine arm regions of tRNAi from 158 species belonging to all 3 domains of life. All the trees distinctly assembled into 3 domains of life. Large trees, generated using data for all the tRNAs of a vast number of species, fail to reveal the major evolutionary events and identity of the probable elongator tRNA sequences that could be ancestor of tRNAi. Therefore, we constructed trees using the minihelix or the whole sequence of species specific tRNAs, and iterated our analysis on 50 eubacterial species. We identified tRNA(Pro), tRNA(Glu), or tRNA(Thr) (but surprisingly not elongator tRNA(Met)) as probable ancestors of tRNAi. We then determined the factors imposing selection of methionine as the initiating amino acid. Overall frequency of occurrence of methionine, whose metabolic cost of synthesis is the highest among all amino acids, remains almost unchanged across the 3 domains of life. Our correlation analysis shows that its high metabolic cost is independent of many physicochemical properties of the side chain. Our results indicate that selection of methionine, as the initiating amino acid was possibly a consequence of the evolution of one-carbon metabolism, which plays an important role in regulating translation initiation.

Xenogeneic Silencing and Its Impact on Bacterial Genomes

The H-NS (heat-stable nucleoid structuring) protein affects both nucleoid compaction and global gene regulation. H-NS appears to act primarily as a silencer of AT-rich genetic material acquired by horizontal gene transfer. As such, it is key in the regulation of most genes involved in virulence and in adaptation to new environmental niches. Here we review recent progress in understanding the biochemistry of H-NS and how xenogeneic silencing affects bacterial evolution. We highlight the strengths and weaknesses of some of the models proposed in H-NS-mediated nucleoprotein complex formation. Based on recent single-molecule studies, we also propose a novel mode of DNA compaction by H-NS termed intrabridging to explain over two decades of observations of the H-NS molecule.

Genome-Wide Dosage-Dependent and -Independent Regulation Contributes to Gene Expression and Evolutionary Novelty in Plant Polyploids

Polyploidy provides evolutionary and morphological novelties in many plants and some animals. However, the role of genome dosage and composition in gene expression changes remains poorly understood. Here, we generated a series of resynthesized Arabidopsis tetraploids that contain 0-4 copies of Arabidopsis thaliana and Arabidopsis arenosa genomes and investigated ploidy and hybridity effects on gene expression. Allelic expression can be defined as dosage dependent (expression levels correlate with genome dosages) or otherwise as dosage independent. Here, we show that many dosage-dependent genes contribute to cell cycle, photosynthesis, and metabolism, whereas dosage-independent genes are enriched in biotic and abiotic stress responses. Interestingly, dosage-dependent genes tend to be preserved in ancient biochemical pathways present in both plant and nonplant species, whereas many dosage-independent genes belong to plant-specific pathways. This is confirmed by an independent analysis using Arabidopsis phylostratigraphic map. For A. thaliana loci, the dosage-dependent alleles are devoid of TEs and tend to correlate with H3K9ac, H3K4me3, and CG methylation, whereas the majority of dosage-independent alleles are enriched with TEs and correspond to H3K27me1, H3K27me3, and CHG (H = A, T, or C) methylation. Furthermore, there is a parent-of-origin effect on nonadditively expressed genes in the reciprocal allotetraploids especially when A. arenosa is used as the pollen donor, leading to metabolic and morphological changes. Thus, ploidy, epigenetic modifications, and cytoplasmic-nuclear interactions shape gene expression diversity in polyploids. Dosage-dependent expression can maintain growth and developmental stability, whereas dosage-independent expression can facilitate functional divergence between homeologs (subfunctionalization and/or neofunctionalization) during polyploid evolution.