Comparative phage genomics and the evolution of Siphoviridae: insights from dairy phages

Comparative phage genomics can retrace part of the evolutionary history of phage modules encoding phage-specific functions such as capsid building or establishment of the lysogenic state. The diagnosis of relatedness is not based exclusively on sequence similarity, but includes topological considerations of genome organization. The gene maps from the lambda-, psiM2-, L5-, Sfi21-, Sfi11-, phiC31-, sk1- and TM4-like phages showed a remarkable synteny of their structural genes defining a lambda supergroup within Siphoviridae (Caudovirales with long non-contractile tails). A hierarchy of relatedness within the lambda supergroup suggested elements of vertical evolution in the capsid module of Siphoviridae. Links to P22-like Podoviridae and P2-like Myoviridae were also detected. Numerous cases of horizontal gene transfer were observed, but recent transfers were limited to interbreeding phage populations. We suggest that tailed phages are the result of both vertical and horizontal evolution and are thus a good model system for web-like phylogenies.     

Evolution of metabolic pathways by chance assembly of enzyme proteins generated from sense and antisense strands of pre-existing genes.

In order to get an insight into the evolutionary aspect of metabolic pathways, especially of the ubiquitous glycolytic pathway, we have carried out an extensive search of sense-sense and sense-antisense similarities for enzyme proteins in the glycolytic pathway, the pentose phosphate cycle, alcohol and lactate fermentation pathways and the TCA cycle. This investigation of amino acid sequences reveals a curious pattern of similarity relations; no similarity can be found between the enzyme proteins in a section of the glycolytic pathway where the glyceraldehyde-3-phosphate or even glycerol-3-phosphate is converted into the pyruvate while many examples of sense-sense and sense-antisense similarities are found even between enzyme proteins in distant blocks, e.g. between the proteins in the TCA cycle and those in the pentose phosphate cycle, as well as between the functionally associated proteins in each of these blocks. Complementary to this characteristic pattern of amino acid sequence similarity, the search for similarities of nucleotide sequences also finds that the similarities of glycolytic enzyme genes, some sense-sense and others sense-antisense similarities, are concentrated on the nucleotide sequences of prokaryotic 16S or eukaryotic 18S ribosomal RNA gene with its flanks, although some of the copy sequences are also found in transfer RNA genes as well as in 23S or 26S ribosomal RNA gene. These results strongly suggest that the metabolic pathways have been developed by the chance assembly of enzyme proteins generated from the sense and antisense strands of pre-existing genes, e.g. the fermentation pathways and pentose phosphate cycle by the proteins from the genes of enzymes in the glycolytic pathway and the TCA cycle from all these successively increased genes, ascribing the origin of metabolic enzyme genes to the close relation between the glycolytic enzyme protein genes and the RNA gene cluster.     

Replacement of the Arginine Biosynthesis Operon in Xanthomonadales by Lateral Gene Transfer

The role of lateral gene transfer (LGT) in prokaryotes has been shown to rapidly change the genome content, providing new gene tools for environmental adaptation. Features related to pathogenesis and resistance to strong selective conditions have been widely shown to be products of gene transfer between bacteria. The genomes of the γ-proteobacteria from the genus Xanthomonas, composed mainly of phytopathogens, have potential genomic islands that may represent imprints of such evolutionary processes. In this work, the evolution of genes involved in the pathway responsible for arginine biosynthesis in Xanthomonadales was investigated, and several lines of evidence point to the foreign origin of the arg genes clustered within a potential operon. Their presence inside a potential genomic island, bordered by a tRNA gene, the unusual ranking of sequence similarity, and the atypical phylogenies indicate that the metabolic pathway for arginine biosynthesis was acquired through LGT in the Xanthomonadales group. Moreover, although homologues were also found in Bacteroidetes (Flavobacteria group), for many of the genes analyzed close homologues are detected in different life domains (Eukarya and Archaea), indicating that the source of these arg genes may have been outside the Bacteria clade. The possibility of replacement of a complete primary metabolic pathway by LGT events supports the selfish operon hypothesis and may occur only under very special environmental conditions. Such rare events reveal part of the history of these interesting mosaic Xanthomonadales genomes, disclosing the importance of gene transfer modifying primary metabolism pathways and extending the scenario for bacterial genome evolution.     

Phylogeny vs genome reshuffling: horizontal gene transfer

The evolutionary events in organisms can be tracked to the transfer of genetic material. The inheritance of genetic material among closely related organisms is a slow evolutionary process. On the other hand, the movement of genes among distantly related species can account for rapid evolution. The later process has been quite evident in the appearance of antibiotic resistance genes among human and animal pathogens. Phylogenetic trees based on such genes and those involved in metabolic activities reflect the incongruencies in comparison to the 16S rDNA gene, generally used for taxonomic relationships. Such discrepancies in gene inheritance have been termed as horizontal gene transfer (HGT) events. In the post-genomic era, the explosion of known sequences through large-scale sequencing projects has unraveled the weakness of traditional 16S rDNA gene tree based evolutionary model. Various methods to scrutinize HGT events include atypical composition, abnormal sequence similarity, anomalous phylogenetic distribution, unusual phyletic patterns, etc. Since HGT generates greater genetic diversity, it is likely to increase resource use and ecosystem resilience.     

Inter-phylum HGT has shaped the metabolism of many mesophilic and anaerobic bacteria

Genome sequencing has revealed that horizontal gene transfer (HGT) is a major evolutionary process in bacteria. Although it is generally assumed that closely related organisms engage in genetic exchange more frequently than distantly related ones, the frequency of HGT among distantly related organisms and the effect of ecological relatedness on the frequency has not been rigorously assessed. Here, we devised a novel bioinformatic pipeline, which minimized the effect of over-representation of specific taxa in the available databases and other limitations of homology-based approaches by analyzing genomes in standardized triplets, to quantify gene exchange between bacterial genomes representing different phyla. Our analysis revealed the existence of networks of genetic exchange between organisms with overlapping ecological niches, with mesophilic anaerobic organisms showing the highest frequency of exchange and engaging in HGT twice as frequently as their aerobic counterparts. Examination of individual cases suggested that inter-phylum HGT is more pronounced than previously thought, affecting up to ∼16% of the total genes and ∼35% of the metabolic genes in some genomes (conservative estimation). In contrast, ribosomal and other universal protein-coding genes were subjected to HGT at least 150 times less frequently than genes encoding the most promiscuous metabolic functions (for example, various dehydrogenases and ABC transport systems), suggesting that the species tree based on the former genes may be reliable. These results indicated that the metabolic diversity of microbial communities within most habitats has been largely assembled from preexisting genetic diversity through HGT and that HGT accounts for the functional redundancy among phyla.     

Optimal gene partition into operons correlates with gene functional order

Gene arrangement into operons varies between bacterial species. Genes in a given system can be on one operon in some organisms and on several operons in other organisms. Existing theories explain why genes that work together should be on the same operon, since this allows for advantageous lateral gene transfer and accurate stoichiometry. But what causes the frequent separation into multiple operons of co-regulated genes that act together in a pathway? Here we suggest that separation is due to benefits made possible by differential regulation of each operon. We present a simple mathematical model for the optimal distribution of genes into operons based on a balance of the cost of operons and the benefit of regulation that provides 'just-when-needed' temporal order. The analysis predicts that genes are arranged such that genes on the same operon do not skip functional steps in the pathway. This prediction is supported by genomic data from 137 bacterial genomes. Our work suggests that gene arrangement is not only the result of random historical drift, genome re-arrangement and gene transfer, but has elements that are solutions of an evolutionary optimization problem. Thus gene functional order may be inferred by analyzing the operon structure across different genomes.     

A SINGLE ORIGIN OF PLASTIDS REVISITED: CONVERGENT EVOLUTION IN ORGANELLAR GENOME CONTENT1

In recent years a consensus has emerged from molecular phylogenetic investigations favoring a common endosymbiotic ancestor for all chloroplasts. It is within this conceptual framework that most comparative analyses of eukaryotic biochemistry and genetics now are interpreted. One of the first and most influential sources of data leading to this consensus is the remarkable similarity in genome content among all major plastid lineages. Here we report statistical analyses of two sequence data sets, genes encoding ribosomal proteins and transfer RNAs, from representatives of the three primary plastid lineages and a mitochondrion. The latter almost certainly originated in an independent endosymbiotic association and serves as a control for similarity due to convergent evolution. When genes related to organelle‐specific function are factored out, plastid genomes appear to be no more similar to each other than they are to the mitochondrion. Total similarities in gene content, measured as deviations from the expectation from a process of random gene loss, are correlated with the extent of reduction in the two genomes compared. They do not appear to reflect putative evolutionary relationships among plastids. These analyses indicate that similarities in plastid genome content are better explained by convergent evolution due to constraint on gene loss than by a shared evolutionary history. A review of other data cited as support for a single plastid origin suggests that the alternative hypothesis of multiple origins is at least equally consistent in most cases.     

Higher Gene Duplicabilities for Metabolic Proteins Than for Nonmetabolic Proteins in Yeast and <Emphasis Type="Italic">E. coli</Emphasis>

Although the evolutionary significance of gene duplication has long been appreciated, it remains unclear what factors determine gene duplicability. In this study we investigated whether metabolism is an important determinant of gene duplicability because cellular metabolism is crucial for the survival and reproduction of an organism. Using genomic data and metabolic pathway data from the yeast (Saccharomyces cerevisiae) and Escherichia coli, we found that metabolic proteins indeed tend to have higher gene duplicability than nonmetabolic proteins. Moreover, a detailed analysis of metabolic pathways in these two organisms revealed that genes in the central metabolic pathways and the catabolic pathways have, on average, higher gene duplicability than do other genes and that most genes in anabolic pathways are single-copy genes.Reviewing Editor: Dr. Rüdiger Cerff     

Prediction of operon-like gene clusters in the Arabidopsis thaliana genome based on co-expression analysis of neighboring genes

Operon-like arrangements of genes occur in eukaryotes ranging from yeasts and filamentous fungi to nematodes, plants, and mammals. In plants, several examples of operon-like gene clusters involved in metabolic pathways have recently been characterized, e.g. the cyclic hydroxamic acid pathways in maize, the avenacin biosynthesis gene clusters in oat, the thalianol pathway in Arabidopsis thaliana, and the diterpenoid momilactone cluster in rice. Such operon-like gene clusters are defined by their co-regulation or neighboring positions within immediate vicinity of chromosomal regions. A comprehensive analysis of the expression of neighboring genes therefore accounts a crucial step to reveal the complete set of operon-like gene clusters within a genome. Genome-wide prediction of operon-like gene clusters should contribute to functional annotation efforts and provide novel insight into evolutionary aspects acquiring certain biological functions as well. We predicted co-expressed gene clusters by comparing the Pearson correlation coefficient of neighboring genes and randomly selected gene pairs, based on a statistical method that takes false discovery rate (FDR) into consideration for 1469 microarray gene expression datasets of A. thaliana. We estimated that A. thaliana contains 100 operon-like gene clusters in total. We predicted 34 statistically significant gene clusters consisting of 3 to 22 genes each, based on a stringent FDR threshold of 0.1. Functional relationships among genes in individual clusters were estimated by sequence similarity and functional annotation of genes. Duplicated gene pairs (determined based on BLAST with a cutoff of E<10(-5)) are included in 27 clusters. Five clusters are associated with metabolism, containing P450 genes restricted to the Brassica family and predicted to be involved in secondary metabolism. Operon-like clusters tend to include genes encoding bio-machinery associated with ribosomes, the ubiquitin/proteasome system, secondary metabolic pathways, lipid and fatty-acid metabolism, and the lipid transfer system.     

METABOLIC FLUX IS A DETERMINANT OF THE EVOLUTIONARY RATES OF ENZYME‐ENCODING GENES

Relationships between evolutionary rates and gene properties on a genomic, functional, pathway, or system level are being explored to unravel the principles of the evolutionary process. In particular, functional network properties have been analyzed to recognize the constraints they may impose on the evolutionary fate of genes. Here we took as a case study the core metabolic network in human erythrocytes and we analyzed the relationship between the evolutionary rates of its genes and the metabolic flux distribution throughout it. We found that metabolic flux correlates with the ratio of nonsynonymous to synonymous substitution rates. Genes encoding enzymes that carry high fluxes have been more constrained in their evolution, while purifying selection is more relaxed in genes encoding enzymes carrying low metabolic fluxes. These results demonstrate the importance of considering the dynamical functioning of gene networks when assessing the action of selection on system‐level properties.     

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.     

GO4genome: A Prokaryotic Phylogeny Based on Genome Organization

Determining the phylogeny of closely related prokaryotes may fail in an analysis of rRNA or a small set of sequences. Whole-genome phylogeny utilizes the maximally available sample space. For a precise determination of genome similarity, two aspects have to be considered when developing an algorithm of whole-genome phylogeny: (1) gene order conservation is a more precise signal than gene content; and (2) when using sequence similarity, failures in identifying orthologues or the in situ replacement of genes via horizontal gene transfer may give misleading results. GO4genome is a new paradigm, which is based on a detailed analysis of gene function and the location of the respective genes. For characterization of genes, the algorithm uses gene ontology enabling a comparison of function independent of evolutionary relationship. After the identification of locally optimal series of gene functions, their length distribution is utilized to compute a phylogenetic distance. The outcome is a classification of genomes based on metabolic capabilities and their organization. Thus, the impact of effects on genome organization that are not covered by methods of molecular phylogeny can be studied. Genomes of strains belonging to Escherichia coli, Shigella, Streptococcus, Methanosarcina, and Yersinia were analyzed. Differences from the findings of classical methods are discussed. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Measuring genome conservation across taxa: divided strains and united kingdoms

Species evolutionary relationships have traditionally been defined by sequence similarities of phylogenetic marker molecules, recently followed by whole-genome phylogenies based on gene order, average ortholog similarity or gene content. Here, we introduce genome conservation--a novel metric of evolutionary distances between species that simultaneously takes into account, both gene content and sequence similarity at the whole-genome level. Genome conservation represents a robust distance measure, as demonstrated by accurate phylogenetic reconstructions. The genome conservation matrix for all presently sequenced organisms exhibits a remarkable ability to define evolutionary relationships across all taxonomic ranges. An assessment of taxonomic ranks with genome conservation shows that certain ranks are inadequately described and raises the possibility for a more precise and quantitative taxonomy in the future. All phylogenetic reconstructions are available at the genome phylogeny server: .     

系统发育基因组学研究进展

系统发育基因组学是利用全基因组数据构建系统发育树的新领域。全基因组数据能有效消除横向基因转移和类群间基因进化速率差异等因素对系统发育树的影响。根据所使用的全基因组数据的类型, 可以将系统发育基因组学方法分为以下5类：多基因联合建树方法, 基于基因含量的方法, 基于基因排列信息的方法, 基于序列短串含量特征信息的方法及基于代谢途径的方法。文章系统地总结了每一类方法的原理、速度、准确性、适用范围及在各个生物类群中的应用, 并对系统发育基因组学的前景及面临的挑战进行了概述。     

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.     

CTnscr94, a conjugative transposon found in enterobacteria.

Conjugational transposons are important for horizontal gene transfer in gram-positive and gram-negative bacteria, but have not been reported yet for enteric bacteria. Salmonella senftenberg 5494-57 has previously been shown to transfer by conjugation genes for a sucrose fermentation pathway which were located on a DNA element called scr-94. We report here that the corresponding scr genes for a phosphoenolpyruvate-dependent sucrose:phosphotransferase system and a sucrose metabolic pathway are located on a large (ca. 100 kb) conjugative transposon renamed CTnscr94. The self-transmissible element integrates at two specific attachment sites in a RecA-independent way into the chromosome of Escherichia coli K-12 strains. One site was identified within pheV, the structural gene for a tRNA(Phe). Sequencing of both ends of CTnscr94 revealed the presence of the 3' part of pheV on one end such that after integration of the element, a complete pheV gene is retained. CTnscr94 represents, to our knowledge, the first conjugational transposon found in enteric bacteria.     

Prediction of missing enzyme genes in a bacterial metabolic network. Reconstruction of the lysine-degradation pathway of Pseudomonas aeruginosa

The metabolic network is an important biological network which consists of enzymes and chemical compounds. However, a large number of metabolic pathways remains unknown, and most organism-specific metabolic pathways contain many missing enzymes. We present a novel method to identify the genes coding for missing enzymes using available genomic and chemical information from bacterial genomes. The proposed method consists of two steps: (a) estimation of the functional association between the genes with respect to chromosomal proximity and evolutionary association, using supervised network inference; and (b) selection of gene candidates for missing enzymes based on the original candidate score and the chemical reaction information encoded in the EC number. We applied the proposed methods to infer the metabolic network for the bacteria Pseudomonas aeruginosa from two genomic datasets: gene position and phylogenetic profiles. Next, we predicted several missing enzyme genes to reconstruct the lysine-degradation pathway in P. aeruginosa using EC number information. As a result, we identified PA0266 as a putative 5-aminovalerate aminotransferase (EC 2.6.1.48) and PA0265 as a putative glutarate semialdehyde dehydrogenase (EC 1.2.1.20). To verify our prediction, we conducted biochemical assays and examined the activity of the products of the predicted genes, PA0265 and PA0266, in a coupled reaction. We observed that the predicted gene products catalyzed the expected reactions; no activity was seen when both gene products were omitted from the reaction.     

Pleiotropy and GAL pathway degeneration in yeast

Traits that do not contribute to fitness are expected to be lost during the course of evolution, either as a result of selection or drift. The Leloir pathway of galactose metabolism (GAL) is an extensively studied metabolic pathway that degenerated on at least three independent occasions during the evolutionary diversification of yeasts, suggesting that the pathway is costly to maintain in environments that lack galactose. Here I test this hypothesis by competing GAL pathway deletion mutants of Saccharomyces cerevisiae against an isogenic strain with an intact GAL pathway under conditions where expression of the pathway is normally induced, repressed, or uninduced. These experiments do not support the hypothesis that pleiotropy drives GAL pathway degeneration, because mutations that knock out individual GAL genes do not tend to increase fitness in the absence of galactose. At a molecular level, this result can be explained by the fact that yeast uses inexpensive regulatory proteins to tightly regulate the expression of structural genes that are costly to express. I argue that these results have general relevance for our understanding of the fitness consequences of gene disruption in yeast.     

Diversity of DNA Sequences Among Restriction Endonucleases Producing Selenomonas ruminantium Isolates Detected by Enterobacterial Repetitive Intergenic Consensus Based Polymerase Chain Reaction (ERIC-PCR)

Enterobacterial repetitive intergenic consensus-based polymerase chain reaction (ERIC-PCR) was found useful for discrimination of rumen selenomonads. Simultaneous use of ERICIR and ERIC2 primers yielded strain-specific banding patterns. The patterns were compared using Dice similarity coefficients and a DNA relatedness dendrogram based on the unweighted pair group method using arithmetic averages (UPGMA) was constructed. Five clusters and four single strains were identified at a similarity level of 50%. Very weak grouping was observed for lactilytica and ruminantium subspecies ofSelenomonas ruminantium , indicating that lactate utilization has probably no taxonomic value. Restriction and modification phenotypes are weakly reflected in the dendrogram probably as the result of horizontal genetic transfer of genes encoding these phenotypic traits. While diverse in ERIC-PCR analysis, strains shown little variation in restriction fragment length polymorphism of amplified 16S-rRNA genes. All but one strain produced nearly identical profile indicating that majority of DNA diversity observed is due to epigenetic factors and not due to evolutionary divergence.