Evolution of variation in presence and absence of genes in bacterial pathways

ABSTRACT: BACKGROUND: Bacterial genomes exhibit a remarkable degree of variation in the presence and absence of genes, which probably extends to the level of individual pathways. This variation may be a consequence of the significant evolutionary role played by horizontal gene transfer, but might also be explained by the loss of genes through mutation. A challenge is to understand why there would be variation in gene presence within pathways if they confer a benefit only when complete. RESULTS: Here, we develop a mathematical model to study how variation in pathway content is produced by horizontal transfer, gene loss and partial exposure of a population to a novel environment. CONCLUSIONS: We discuss the possibility that variation in gene presence acts as cryptic genetic variation on which selection acts when the appropriate environment occurs. We find that a high level of variation in gene presence can be readily explained by decay of the pathway through mutation when there is no longer exposure to the selective environment, or when selection becomes too weak to maintain the genes. In the context of pathway variation the role of horizontal gene transfer is probably the initial introduction of a complete novel pathway rather than in building up the variation in a genome without the pathway.     

Yeast metabolic innovations emerged via expanded metabolic network and gene positive selection

Yeasts are known to have versatile metabolic traits, while how these metabolic traits have evolved has not been elucidated systematically. We performed integrative evolution analysis to investigate how genomic evolution determines trait generation by reconstructing genome‐scale metabolic models (GEMs) for 332 yeasts. These GEMs could comprehensively characterize trait diversity and predict enzyme functionality, thereby signifying that sequence‐level evolution has shaped reaction networks towards new metabolic functions. Strikingly, using GEMs, we can mechanistically map different evolutionary events, e.g. horizontal gene transfer and gene duplication, onto relevant subpathways to explain metabolic plasticity. This demonstrates that gene family expansion and enzyme promiscuity are prominent mechanisms for metabolic trait gains, while GEM simulations reveal that additional factors, such as gene loss from distant pathways, contribute to trait losses. Furthermore, our analysis could pinpoint to specific genes and pathways that have been under positive selection and relevant for the formulation of complex metabolic traits, i.e. thermotolerance and the Crabtree effect. Our findings illustrate how multidimensional evolution in both metabolic network structure and individual enzymes drives phenotypic variations.     

Repeated horizontal gene transfer of GALactose metabolism genes violates Dollo’s law of irreversible loss

Dollo’s law posits that evolutionary losses are irreversible, thereby narrowing the potential paths of evolutionary change. While phenotypic reversals to ancestral states have been observed, little is known about their underlying genetic causes. The genomes of budding yeasts have been shaped by extensive reductive evolution, such as reduced genome sizes and the losses of metabolic capabilities. However, the extent and mechanisms of trait reacquisition after gene loss in yeasts have not been thoroughly studied. Here, through phylogenomic analyses, we reconstructed the evolutionary history of the yeast galactose utilization pathway and observed widespread and repeated losses of the ability to utilize galactose, which occurred concurrently with the losses of GALactose (GAL) utilization genes. Unexpectedly, we detected multiple galactose-utilizing lineages that were deeply embedded within clades that underwent ancient losses of galactose utilization. We show that at least two, and possibly three, lineages reacquired the GAL pathway via yeast-to-yeast horizontal gene transfer. Our results show how trait reacquisition can occur tens of millions of years after an initial loss via horizontal gene transfer from distant relatives. These findings demonstrate that the losses of complex traits and even whole pathways are not always evolutionary dead-ends, highlighting how reversals to ancestral states can occur.     

Evolution of tryptophan biosynthetic pathway in microbial genomes: a comparative genetic study

Biosynthetic pathway evolution needs to consider the evolution of a group of genes that code for enzymes catalysing the multiple chemical reaction steps leading to the final end product. Tryptophan biosynthetic pathway has five chemical reaction steps that are highly conserved in diverse microbial genomes, though the genes of the pathway enzymes show considerable variations in arrangements, operon structure (gene fusion and splitting) and regulation. We use a combined bioinformatic and statistical analyses approach to address the question if the pathway genes from different microbial genomes, belonging to a wide range of groups, show similar evolutionary relationships within and between them. Our analyses involved detailed study of gene organization (fusion/splitting events), base composition, relative synonymous codon usage pattern of the genes, gene expressivity, amino acid usage, etc. to assess inter- and intra-genic variations, between and within the pathway genes, in diverse group of microorganisms. We describe these genetic and genomic variations in the tryptophan pathway genes in different microorganisms to show the similarities across organisms, and compare the same genes across different organisms to find the possible variability arising possibly due to horizontal gene transfers. Such studies form the basis for moving from single gene evolution to pathway evolutionary studies that are important steps towards understanding the systems biology of intracellular pathways.     

Evidences of lateral gene transfer between archaea and pathogenic bacteria

Acquisition of new genetic material through horizontal gene transfer has been shown to be an important feature in the evolution of many pathogenic bacteria. Changes in the genetic repertoire, occurring through gene acquisition and deletion, are the major events underlying the emergence and evolution of bacterial pathogens. However, horizontal gene transfer across the domains i.e. archaea and bacteria is not so common. In this context, we explore events of horizontal gene transfer between archaea and bacteria. In order to determine whether the acquisition of archaeal genes by lateral gene transfer is an important feature in the evolutionary history of the pathogenic bacteria, we have developed a scheme of stepwise eliminations that identifies archaeal-like genes in various bacterial genomes. We report the presence of 9 genes of archaeal origin in the genomes of various bacteria, a subset of which is also unique to the pathogenic members and are not found in respective non-pathogenic counterparts. We believe that these genes, having been retained in the respective genomes through selective advantage, have key functions in the organism’s biology and may play a role in pathogenesis.     

Horizontal transfer of a pathway for coumarate catabolism unexpectedly inhibits purine nucleotide biosynthesis

A microbe’s ecological niche and biotechnological utility are determined by its specific set of co‐evolved metabolic pathways. The acquisition of new pathways, through horizontal gene transfer or genetic engineering, can have unpredictable consequences. Here we show that two different pathways for coumarate catabolism failed to function when initially transferred into Escherichia coli. Using laboratory evolution, we elucidated the factors limiting activity of the newly acquired pathways and the modifications required to overcome these limitations. Both pathways required host mutations to enable effective growth with coumarate, but the necessary mutations differed. In one case, a pathway intermediate inhibited purine nucleotide biosynthesis, and this inhibition was relieved by single amino acid replacements in IMP dehydrogenase. A strain that natively contains this coumarate catabolism pathway, Acinetobacter baumannii, is resistant to inhibition by the relevant intermediate, suggesting that natural pathway transfers have faced and overcome similar challenges. Molecular dynamics simulation of the wild type and a representative single‐residue mutant provide insight into the structural and dynamic changes that relieve inhibition. These results demonstrate how deleterious interactions can limit pathway transfer, that these interactions can be traced to specific molecular interactions between host and pathway, and how evolution or engineering can alleviate these limitations.     

Evolution of vitamin B6 (pyridoxine) metabolism by gain and loss of genes

Vitamin B(6) (VB6) functions as a cofactor of many diverse enzymes in amino acid metabolism. Three metabolic pathways for pyridoxal 5'-phosphate (PLP; the active form of VB6) are known: the de novo pathway, the salvage pathway, and the fungal type pathway. Most unicellular organisms and plants biosynthesize VB6 using one or two of these three biosynthetic pathways. However, animals such as insects and mammals do not possess any of the pathways and, thus, need to intake VB6 in their diet to survive. It is conceivable that breakdowns of these pathways occurred in the evolutionary lineages of insects and mammals, and one of the major reasons for this would be the loss of pertinent genes. We studied the evolution of VB6 biosynthesis from the view of the gain and loss of 10 pertinent genes in 122 species whose genome sequences were completely determined. The results revealed that each gene in the pathways was lost more than once in the entire evolutionary lineages of the 122 species. We also found the following three points regarding the evolution of PLP biosynthesis: (1) the breakdown of the PLP biosynthetic pathways occurred independently at least three times in animal lineages, (2) the de novo pathway was formed by the generation of pdxB in gamma-proteobacteria, and (3) the order of the gene loss in VB6 metabolism was conserved among different evolutionary lineages. These results suggest that the evolution of VB6 metabolism was subject to gains and frequent losses of related genes in the 122 species examined. This dynamic nature of the evolutionary changes must have been responsible for the breakdowns of the pathways, resulting in profound differentiation of heterotrophy among the species.     

Eukaryotic origin of a metabolic pathway in virus by horizontal gene transfer

Horizontal gene transfer, the movement of genetic materials across the normal mating barriers between organisms occurs frequently and contributes significantly to the evolution of both eukaryotic and prokaryotic genomes. However, few concurrent transfers of functionally related genes implemented in a pathway from eukaryotes to prokaryotes are observed. Here, we did phylogenetic analyses to support that the genes, i.e. dihydrofolate reductase, glycine hydroxymethyltransferase, and thymidylate synthase involved in thymidylate metabolism, in Hz-1 virus were obtained from insect genome recently by independent horizontal gene transfers. In addition, five other related genes in nucleotide metabolism show evidences of horizontal gene transfers. These genes demonstrate similar expression pattern, and they may have formatted a functionally related pathway (e.g. thymidylate synthesis, and DNA replication) in Hz-1 virus. In conclusion, we provide an example of horizontal gene transfer of functionally related genes in a pathway to prokaryote from eukaryote.     

Caught in the act: Rapid,symbiont‐driven evolution

Facultative bacterial endosymbionts can transfer horizontally among lineages of their arthropod hosts, providing the recipient with a suite of traits that can lead to rapid evolutionary response, as has been recently demonstrated. But how common is symbiont‐driven evolution? Evidence suggests that successful symbiont transfers are most likely within a species or among closely related species, although more distant transfers have occurred over evolutionary history. Symbiont‐driven evolution need not be a function of a recent horizontal transfer, however. Many endosymbionts infect only a small proportion of a host population, but could quickly increase in frequency under favorable selection regimes. Some host species appear to accumulate a diversity of facultative endosymbionts, and it is among these species that symbiont‐driven evolution should be most prevalent. It remains to be determined how frequently symbionts enable rapid evolutionary response by their hosts, but substantial ecological effects are a likely consequence whenever it does occur.     

Molecular clocks and evolutionary relationships: Possible distortions due to horizontal gene flow

Summary This paper discusses recent evidence suggesting that genetic information from one species occasionally transfers to another remotely related species. Besides addressing the issue of whether or not the molecular data are consistent with a wide-spread influence of horizontal gene transfer, the paper shows that horizontal gene flow would not necessarily preclude a linear molecular clock or change the rate of molecular evolution (assuming the neutral allele theory). A pervasive influence of horizontal gene transfer is more than just consistent with the data of molecular evolution, it also provides a unique explanation for a number of possibly conflicting phylogenies and contradictory clocks. This phenomenon might explain why some protein clocks are linear while the superoxide dismutase clock is not, how the molecular data on the phylogeny of apes and Australian song birds are not necessarily in conflict with those based on morphology, and, finally, why the mycoplasmas have an accelerated molecular clock.     

Phylogenetic analysis based on genome-scale metabolic pathway reaction content

Phylogenetic classifications based on single genes such as rRNA genes do not provide a complete and accurate picture of evolution because they do not account for evolutionary leaps caused by gene transfer, duplication, deletion and functional replacement. Here, we present a whole-genome-scale phylogeny based on metabolic pathway reaction content. From the genome sequences of 42 microorganisms, we deduced the metabolic pathway reactions and used the relatedness of these contents to construct a phylogenetic tree that represents the similarity of metabolic profiles (relatedness) as well as the extent of metabolic pathway similarity (evolutionary distance). This method accounts for horizontal gene transfer and specific gene loss by comparison of whole metabolic subpathways, and allows evaluation of evolutionary relatedness and changes in metabolic pathways. Thus, a tree based on metabolic pathway content represents both the evolutionary time scale (changes in genetic content) and the evolutionary process (changes in metabolism).     

Unintentional selection,unanticipated insights: introductions,stocking and the evolutionary ecology of fishes

Natural environmental change has produced countless opportunities for species to disperse into and persist in habitats where they previously did not exist. Introduction and stocking programmes have facilitated similar sorts of colonization opportunities across considerably greater geographical scales and often in much shorter periods of time. Even though the mechanism of colonization differs, the result can be the same: evolutionary change in the colonizing population in response to novel selection pressures. As a consequence, some human‐mediated fish transfers have unintentionally yielded novel research opportunities to study how phenotypes and genes interact with their environment and affect ecological and evolutionary change. The primary purpose here is to explore how work, directly or indirectly involved with human‐mediated transfers, has unintentionally yielded novel research and research opportunities in fish ecology and evolution. Insights have produced new knowledge or altered previously held perceptions on topics such as local adaptation, rate of evolutionary change, phenotypic plasticity, alternative reproductive strategies, population structure and colonization probability. Well‐documented stocking programmes, especially in terms of history, numbers and original population sources, can provide highly fertile ground for generating further insights on the ecology and evolution of fishes and of the factors likely to influence the success of conservation‐based, restoration programmes.     

Inverse relationship between evolutionary rate and age of mammalian genes

A large number of genes is shared by all living organisms, whereas many others are unique to some specific lineages, indicating their different times of origin. The availability of a growing number of eukaryotic genomes allows us to estimate which mammalian genes are novel genes and, approximately, when they arose. In this article, we classify human genes into four different age groups and estimate evolutionary rates in human and mouse orthologs. We show that older genes tend to evolve more slowly than newer ones; that is, proteins that arose earlier in evolution currently have a larger proportion of sites subjected to negative selection. Interestingly, this property is maintained when a fraction of the fastest-evolving genes is excluded or when only genes belonging to a given functional class are considered. One way to explain this relationship is by assuming that genes maintain their functional constraints along all their evolutionary history, but the nature of more recent evolutionary innovations is such that the functional constraints operating on them are increasingly weaker. Alternatively, our results would also be consistent with a scenario in which the functional constraints acting on a gene would not need to be constant through evolution. Instead, starting from weak functional constraints near the time of origin of a gene-as supported by mechanisms proposed for the origin of orphan genes-there would be a gradual increase in selective pressures with time, resulting in fewer accepted mutations in older versus more novel genes.     

Comparative genomics and bioenergetics.

Bacterial and archaeal complete genome sequences have been obtained from a wide range of evolutionary lines, which allows some general conclusions about the phylogenetic distribution and evolution of bioenergetic pathways to be drawn. In particular, I searched in the complete genomes for key enzymes involved in aerobic and anaerobic respiratory pathways and in photosynthesis, and mapped them into an rRNA tree of sequenced species. The phylogenetic distribution of these enzymes is very irregular, and clearly shows the diverse strategies of energy conservation used by prokaryotes. In addition, a thorough phylogenetic analysis of other bioenergetic protein families of wide distribution reveals a complex evolutionary history for the respective genes. A parsimonious explanation for these complex phylogenetic patterns and for the irregular distribution of metabolic pathways is that the last common ancestor of Bacteria and Archaea contained several members of every gene family as a consequence of previous gene or genome duplications, while different patterns of gene loss occurred during the evolution of every gene family. This would imply that the last universal ancestor was a bioenergetically sophisticated organism. Finally, important steps that occurred during the evolution of energetic machineries, such as the early evolution of aerobic respiration and the acquisition of eukaryotic mitochondria from a proteobacterium ancestor, are supported by the analysis of the complete genome sequences.     

Frequent,independent transfers of a catabolic gene from bacteria to contrasted filamentous eukaryotes

Even genetically distant prokaryotes can exchange genes between them, and these horizontal gene transfer events play a central role in adaptation and evolution. While this was long thought to be restricted to prokaryotes, certain eukaryotes have acquired genes of bacterial origin. However, gene acquisitions in eukaryotes are thought to be much less important in magnitude than in prokaryotes. Here, we describe the complex evolutionary history of a bacterial catabolic gene that has been transferred repeatedly from different bacterial phyla to stramenopiles and fungi. Indeed, phylogenomic analysis pointed to multiple acquisitions of the gene in these filamentous eukaryotes—as many as 15 different events for 65 microeukaryotes. Furthermore, once transferred, this gene acquired introns and was found expressed in mRNA databases for most recipients. Our results show that effective inter-domain transfers and subsequent adaptation of a prokaryotic gene in eukaryotic cells can happen at an unprecedented magnitude.     

Evolution of arginine deiminase (ADI) pathway genes

We have analyzed the evolution of the three genes encoding structural enzymes of the arginine deiminase (ADI) pathway, arginine deiminase (ADI), ornithine transcarbamoylase (OTC), and carbamate kinase (CK) in a wide range of organisms, including Archaea, Bacteria, and Eukarya. This catabolic route was probably present in the last common ancestor to all the domains of life. The results obtained indicate that these genes have undergone a complex evolutionary history, including horizontal transfer events, duplications, and losses. Therefore, these genes are not adequate to infer organismal relationships at deep branching levels, but they provide an insight into how catabolic genes evolved and were assembled into metabolic pathways. Our results suggest that the three genes evolved independently and were later assembled into a single cluster with functional interdependence, thus, providing support for the gene recruitment hypothesis. Furthermore, the molecular phylogenetic analysis of OTC suggests a new classification of these genes into three subfamilies.     

Inferring phylogenetic networks by the maximum parsimony criterion: a case study

Horizontal gene transfer (HGT) may result in genes whose evolutionary histories disagree with each other, as well as with the species tree. In this case, reconciling the species and gene trees results in a network of relationships, known as the "phylogenetic network" of the set of species. A phylogenetic network that incorporates HGT consists of an underlying species tree that captures vertical inheritance and a set of edges which model the "horizontal" transfer of genetic material. In a series of papers, Nakhleh and colleagues have recently formulated a maximum parsimony (MP) criterion for phylogenetic networks, provided an array of computationally efficient algorithms and heuristics for computing it, and demonstrated its plausibility on simulated data. In this article, we study the performance and robustness of this criterion on biological data. Our findings indicate that MP is very promising when its application is extended to the domain of phylogenetic network reconstruction and HGT detection. In all cases we investigated, the MP criterion detected the correct number of HGT events required to map the evolutionary history of a gene data set onto the species phylogeny. Furthermore, our results indicate that the criterion is robust with respect to both incomplete taxon sampling and the use of different site substitution matrices. Finally, our results show that the MP criterion is very promising in detecting HGT in chimeric genes, whose evolutionary histories are a mix of vertical and horizontal evolution. Besides the performance analysis of MP, our findings offer new insights into the evolution of 4 biological data sets and new possible explanations of HGT scenarios in their evolutionary history.