Assessing evolutionary relationships among microbes from whole-genome analysis

Horizontally Transferred DNA in the Genome of the Fungus Pyricularia oryzae is Associated With Repressive Histone Modifications

Horizontal gene transfer (HGT) is a means of exchanging genetic material asexually. The process by which horizontally transferred genes are domesticated by the host genome is of great interest but is not well understood. In this study, we determined the telomere-to-telomere genome sequence of the wheat-infecting Pyricularia oryzae strain Br48. SNP analysis indicated that the Br48 strain is a hybrid of wheat- and Brachiaria-infecting strains by a sexual or parasexual cross. Comparative genomic analysis identified several megabase-scale "insertions" in the Br48 genome, some of which were possibly gained by HGT-related events from related species, such as P. pennisetigena or P. grisea. Notably, the mega-insertions often contained genes whose phylogeny is not congruent with the species phylogeny. Moreover, some of the genes have a close homolog even in distantly related organisms, such as basidiomycetes or prokaryotes, implying the involvement of multiple HGT events. Interestingly, the levels of the silent epigenetic marks H3K9me3 and H3K27me3 in a genomic region tended to be negatively correlated with the phylogenetic concordance of genes in the same region, suggesting that horizontally transferred DNA is preferentially targeted for epigenetic silencing. Indeed, the putative HGT-derived genes were activated when MoKmt6, the gene responsible for H3K27me3 modification, was deleted. Notably, these genes also tended to be up-regulated during infection, suggesting that they are now under host control and have contributed to establishing a fungal niche. In conclusion, this study suggests that epigenetic modifications have played an important role in the domestication of HGT-derived genes in the P. oryzae genome.

Evidence of repeated horizontal transfer of sterol C-5 desaturase encoding genes among dikarya fungi

Gene duplication and horizontal gene transfer (HGT) are two important but different forces for adaptive genome evolution. In eukaryotic organisms, gene duplication is considered to play a more important evolutionary role than HGT. However, certain fungal lineages have developed highly efficient mechanisms that avoid the occurrence of duplicated gene sequences within their genomes. While these mechanisms likely originated as a defense against harmful mobile genetic elements, they come with an evolutionary cost. A prominent example for a genome defense system is the RIP mechanism of the ascomycete fungus Neurospora crassa, which efficiently prevents sequence duplication within the genome and functional redundancy of the subsequent paralogs. Despite this tight control, the fungus possesses two functionally redundant sterol C-5 desaturase enzymes, ERG-10a and ERG-10b, that catalyze the same step during ergosterol biosynthesis. In this study, we addressed this conundrum by phylogenetic analysis of the two proteins and supporting topology tests. We obtained evidence that a primary HGT of a sterol C-5 desaturase gene from Tremellales (an order of Basidiomycota) into a representative of the Pezizomycotina (a subphylum of Ascomycota) is the origin of the ERG-10b sequence. The reconstructed phylogenies suggest that this HGT event was followed by multiple HGT events among other members of the Pezizomycotina, thereby generating a diverse group with members in the four classes Sordariomycetes, Xylonomycetes, Eurotiomycetes and Dothideomycetes, which all harbor the second sterol C-5 desaturase or maintained in some cases only the ERG-10b version of this enzyme. These results furnish an example for a gene present in numerous ascomycetous fungi but primarily acquired by an ancestral HGT event from another fungal phylum. Furthermore, these data indicate that HGT represents one mechanism to generate functional redundancy in organisms with a strict avoidance of gene duplications.

Molecular Evolution of Pseudomonas syringae Type III Secreted Effector Proteins

Diverse Gram-negative pathogens like Pseudomonas syringae employ type III secreted effector (T3SE) proteins as primary virulence factors that combat host immunity and promote disease. T3SEs can also be recognized by plant hosts and activate an effector triggered immune (ETI) response that shifts the interaction back toward plant immunity. Consequently, T3SEs are pivotal in determining the virulence potential of individual P. syringae strains, and ultimately help to restrict P. syringae pathogens to a subset of potential hosts that are unable to recognize their repertoires of T3SEs. While a number of effector families are known to be present in the P. syringae species complex, one of the most persistent challenges has been documenting the complex variation in T3SE contents across a diverse collection of strains. Using the entire pan-genome of 494 P. syringae strains isolated from more than 100 hosts, we conducted a global analysis of all known and putative T3SEs. We identified a total of 14,613 putative T3SEs, 4,636 of which were unique at the amino acid level, and show that T3SE repertoires of different P. syringae strains vary dramatically, even among strains isolated from the same hosts. We also find substantial diversification within many T3SE families, and in many cases find strong signatures of positive selection. Furthermore, we identify multiple gene gain and loss events for several families, demonstrating an important role of horizontal gene transfer (HGT) in the evolution of P. syringae T3SEs. These analyses provide insight into the evolutionary history of P. syringae T3SEs as they co-evolve with the host immune system, and dramatically expand the database of P. syringae T3SEs alleles.

Cyanobacterial phylogenetic analysis based on phylogenomics approaches render evolutionary diversification and adaptation: an overview of representative orders

Phylogenetic studies based on a definite set of marker genes usually reconstruct evolutionary relationships among the prokaryotic species. Based on specific target sequences, such studies represent variations and allow identification of similarities or dissimilarities in organisms. With the advent of completely sequenced genomes and accumulation of information on whole prokaryotic genomes, phylogenetic reconstructions should be considered more reliable if they are ideally based on entire genomes to resolve phylogenetic interest. We applied phylogenomics approaches taking into account completely sequenced cyanobacterial genomes to reconstruct underlying species that represented major taxonomic classes and belonged to distinctly different habitats (freshwater, marine, soils, and rocks). We did not rely on describing phylogeny of all representative class of cyanobacterial species on the basis of only ribosomal gene, 16S rDNA gene. In contrast, we analyzed combined molecular marker and phylogenomics approaches (genome alignment, gene content and gene order, composition vector and protein domain content) for accurately inferring phylogenetic relationship of species. We have shown that this approach reflects the impact of evolution on the organisms and considers connects with the ecological adaptation in cyanobacteria in different habitats. Analysis revealed that the members from marine habitat occupy different profile than those from freshwater. Impact of GC content and genomic repetitiveness over the diversification of cyanobacterial species and their possible role in adaptation was also reflected. Members occupying similar habitats cover more evolutionary distance together and also evolve various strategies for adaptation and survival either through genomic repetitiveness or preferences for genes of particular functions or modified GC content. Genomes undergo different changes for their adaptation in diverse habitats.

The MSDIN family in amanitin-producing mushrooms and evolution of the prolyl oligopeptidase genes

The biosynthetic pathway for amanitins and related cyclic peptides in deadly Amanita (Amanitaceae) mushrooms represents the first known ribosomal cyclic peptide pathway in the Fungi. Amanitins are found outside of the genus in distantly related agarics Galerina (Strophariaceae) and Lepiota (Agaricaceae). A long-standing question in the field persists: why is this pathway present in these phylogenetically disjunct agarics? Two deadly mushrooms, A. pallidorosea and A. subjunquillea, were deep sequenced, and sequences of biosynthetic genes encoding MSDINs (cyclic peptide precursor) and prolyl oligopeptidases (POPA and POPB) were obtained. The two Amanita species yielded 29 and 18 MSDINs, respectively. In addition, two MSDIN sequences were cloned from L. brunneoincarnata basidiomes. The toxin MSDIN genes encoding amatoxins or phallotoxins from the three genera were compared, and a phylogenetic tree constructed. Prolyl oligopeptidase B (POPB), a key enzyme in the biosynthetic pathway, was used in phylogenetic reconstruction to infer the evolutionary history of the genes. Phylogenies of POPB and POPA based on both coding and amino acid sequences showed very different results: while POPA genes clearly reflected the phylogeny of the host species, POPB did not; strikingly, it formed a well-supported monophyletic clade, despite that the species belong to different genera in disjunct families. POPA, a known house-keeping gene, was shown to be restricted in a branch containing only Amanita species and the phylogeny resembled that of those Amanita species. Phylogenetic analyses of MSDIN and POPB genes showed tight coordination and disjunct distribution. A POPB gene tree was compared with a corresponding species tree, and distances and substitution rates were compared. The result suggested POPB genes have significant smaller distances and rates than the house-keeping rpb2, discounting massive gene loss. Under this assumption, the incongruency between the gene tree and species tree was shown with strong support. Additionally, k-mer analyses consistently cluster Galerina and Amanita POPB genes, while Lepiota POPB is distinct. Our result suggests that horizontal gene transfer (HGT), at least between Amanita and Galerina, was involved in the acquisition of POPB genes, which may shed light on the evolution of the α-amanitin biosynthetic pathway.

Phenotypic and genotypic characterisation of an unique indigenous hypersaline unicellular cyanobacterium, Euhalothece sp.nov

A novel halotolerant species of cyanobacterium of the order Chroococcales was isolated from hypersaline estuary in Kwa-Zulu Natal, South Africa. A comprehensive polyphasic approach viz., cell morphology, pigment composition and complete genome sequence analysis was conducted to elucidate the taxonomic position of the isolated strain. The blue-green oval to rod-shaped cells were 14-18 μm in size, and contained a high amount of phycocyanin pigments. The strain was moderate thermotolerant/alkalitolerant halophile with the optimum conditions for growth at 35 °C, pH 8.5 and 120 g/l of NaCl. Based on 16S rRNA gene sequence phylogeny, the strain was related to members of the 'Euhalothece' subcluster (99%). The whole genome sequence was determined, and the annotated genes showed a 90% sequence similarity to the gas-vacuolate, spindle-shaped Dactylococcopsis salina PCC 8305. The size of the genome was determined to be 5,113,178 bp and contained 4332 protein-coding genes and 69 RNA genes with a G + C content of 46.7%. Genes encoding osmoregulation, oxidative stress, heat shock, persister cells, and UV-absorbing secondary metabolites, among others, were identified. Based on the phylogenetic analysis of the 16S rRNA gene sequences, physiological data, pigment compositions and genomic data, the strain is considered to represent a novel species of Euhalothece.

Drivers: A Biologically Contextualized, Cross-Inferential View of the Epidemiology of Neurodegenerative Disorders

Background:

Sutherland et al. (2011) suggested that, instead of risk factors for single neurodegenerative disorders (NDDs), there was a need to identify specific “drivers”, i.e., risk factors with impact on specific deposits, such as amyloid-β, tau, or α-synuclein, acting across entities.

Objectives and Methods:

Redefining drivers as “neither protein/gene- nor entity-specific features identifiable in the clinical and general epidemiology of conformational NDDs (CNDDs) as potential footprints of templating/spread/transfer mechanisms”, we conducted an analysis of the epidemiology of ten CNDDs, searching for patterns.

Results:

We identified seven potential drivers, each of which was shared by at least two CNDDs: 1) an age-at-exposure-related susceptibility to Creutzfeldt-Jakob disease (CJD) and several late-life CNDDs; 2) a relationship between age at onset, survival, and incidence; 3) shared genetic risk factors for CJD and late-life CNNDs; 4) partly shared personal (diagnostic, educational, behavioral, and social risk factors) predating clinical onset of late-life CNDDs; 5) two environmental risk factors, namely, surgery for sporadic CJD and amyotrophic lateral sclerosis, and Bordetella pertussis infection for Parkinson’s disease; 6) reticulo-endothelial system stressors or general drivers (andropause or premenopausal estrogen deficiency, APOEɛ4, and vascular risk factors) for late-life CNDDs such as dementia/Alzheimer’s disease, type-2 diabetes mellitus, and some sporadic cardiac and vascular degenerative diseases; and 7) a high, invariant incidence ratio of sporadic to genetic forms of mid- and late-life CNDDs, and type-2 diabetes mellitus.

Conclusion:

There might be a systematic epidemiologic pattern induced by specific proteins (PrP, TDP-43, SOD1, α-synuclein, amyloid-β, tau, Langerhans islet peptide, and transthyretin) or established combinations of these.

High-quality draft genome sequence of a new phytase-producing microorganism Pantoea sp. 3.5.1

Strain 3.5.1 was isolated from soils of the Republic of Tatarstan, Russia, on the basis of presence of a high phytate-degrading activity. Strains with such activities attract special interest because of its potential use as feed additives and natural manures. Strain 3.5.1 harbors a 99 % 16S rRNA nucleotide sequence similarity to different Pantoea species (P. vagans, P. ananatis, P. agglomerans, P. anthophila and Pantoea sp.) and exhibits unique biochemical properties that do not allow strain identification up to species. Moreover, the strain 3.5.1 shows a low ANI and MALDI-TOF Mass Spectrometry scores. Thus, it is likely that the strain 3.5.1 represents a new Pantoea species. Here, we present the genome sequence of Pantoea sp. strain 3.5.1. The 4,964,649 bp draft genome consists of 23 contigs with 4,556 protein-coding and 143 RNA genes. Genome sequencing and annotation revealed two phytase genes and putative regulatory genes controlling its activity.

Gene regulation by H-NS as a function of growth conditions depends on chromosomal position in Escherichia coli

Cellular adaptation to changing environmental conditions requires the coordinated regulation of expression of large sets of genes by global regulatory factors such as nucleoid associated proteins. Although in eukaryotic cells genomic position is known to play an important role in regulation of gene expression, it remains to be established whether in bacterial cells there is an influence of chromosomal position on the efficiency of these global regulators. Here we show for the first time that genome position can affect transcription activity of a promoter regulated by the histone-like nucleoid-structuring protein (H-NS), a global regulator of bacterial transcription and genome organization. We have used as a local reporter of H-NS activity the level of expression of a fluorescent reporter protein under control of an H-NS−regulated promoter (Phns) at different sites along the genome. Our results show that the activity of the Phns promoter depends on whether it is placed within the AT-rich regions of the genome that are known to be bound preferentially by H-NS. This modulation of gene expression moreover depends on the growth phase and the growth rate of the cells, reflecting the changes taking place in the relative abundance of different nucleoid proteins and the inherent heterogeneous organization of the nucleoid. Genomic position can thus play a significant role in the adaptation of the cells to environmental changes, providing a fitness advantage that can explain the selection of a gene’s position during evolution.

HGTector: an automated method facilitating genome-wide discovery of putative horizontal gene transfers

Background

First pass methods based on BLAST match are commonly used as an initial step to separate the different phylogenetic histories of genes in microbial genomes, and target putative horizontal gene transfer (HGT) events. This will continue to be necessary given the rapid growth of genomic data and the technical difficulties in conducting large-scale explicit phylogenetic analyses. However, these methods often produce misleading results due to their inability to resolve indirect phylogenetic links and their vulnerability to stochastic events.

Results

A new computational method of rapid, exhaustive and genome-wide detection of HGT was developed, featuring the systematic analysis of BLAST hit distribution patterns in the context of a priori defined hierarchical evolutionary categories. Genes that fall beyond a series of statistically determined thresholds are identified as not adhering to the typical vertical history of the organisms in question, but instead having a putative horizontal origin. Tests on simulated genomic data suggest that this approach effectively targets atypically distributed genes that are highly likely to be HGT-derived, and exhibits robust performance compared to conventional BLAST-based approaches. This method was further tested on real genomic datasets, including Rickettsia genomes, and was compared to previous studies. Results show consistency with currently employed categories of HGT prediction methods. In-depth analysis of both simulated and real genomic data suggests that the method is notably insensitive to stochastic events such as gene loss, rate variation and database error, which are common challenges to the current methodology. An automated pipeline was created to implement this approach and was made publicly available at: https://github.com/DittmarLab/HGTector. The program is versatile, easily deployed, has a low requirement for computational resources.

Conclusions

HGTector is an effective tool for initial or standalone large-scale discovery of candidate HGT-derived genes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-717) contains supplementary material, which is available to authorized users.

Whole genome phylogeny of Prochlorococcus marinus group of cyanobacteria: genome alignment and overlapping gene approach

Prochlorococcus is the smallest known oxygenic phototrophic marine cyanobacterium dominating the mid-latitude oceans. Physiologically and genetically distinct P. marinus isolates from many oceans in the world were assigned two different groups, a tightly clustered high-light (HL)-adapted and a divergent low-light (LL-) adapted clade. Phylogenetic analysis of this cyanobacterium on the basis of 16S rRNA and other conserved genes did not show consistency with its phenotypic behavior. We analyzed phylogeny of this genus on the basis of complete genome sequences through genome alignment, overlapping-gene content and gene-order approach. Phylogenetic tree of P. marinus obtained by comparing whole genome sequences in contrast to that based on 16S rRNA gene, corresponded well with the HL/LL ecotypic distinction of twelve strains and showed consistency with phenotypic classification of P. marinus. Evidence for the horizontal descent and acquisition of genes within and across the genus was observed. Many genes involved in metabolic functions were found to be conserved across these genomes and many were continuously gained by different strains as per their needs during the course of their evolution. Consistency in the physiological and genetic phylogeny based on whole genome sequence is established. These observations improve our understanding about the adaptation and diversification of these organisms under evolutionary pressure.

Systems biology approach to integrative comparative genomics

As more and more systems biology approaches are used to investigate the different types of biological macromolecules, increasing numbers of whole genomic studies are now available for a large array of organisms. Whether it is genomics, transcriptomics, proteomics, interactomics or metabolomics, the full complement of genomic information on all different levels can be juxtaposed between different organisms to reveal similarities or differences, and even to provide consensus models. At the intersection of comparative genomics and systems biology lies great possibility for discovery, analysis and prediction. This paper explores this nexus and the relationship from four general levels: DNA, RNA, protein and extragenomic. For each level, we provide an overview of the methods, discuss the potential challenges and survey the current research. Finally, we suggest some organizing principles and make proposals for new areas that will be important for future research.

Whole-proteome analysis of twelve species of alphaproteobacteria links four pathogens

Thousands of whole-genome and whole-proteome sequences have been made available through advances in sequencing technology, and sequences of millions more organisms will become available in the coming years. This wealth of genetic information will provide numerous opportunities to enhance our understanding of these organisms including a greater understanding of relationships among species. Researchers have used 16S rRNA and other gene sequences to study the evolutionary origins of bacteria, but these strategies do not provide insight into the sharing of genes among bacteria via horizontal transfer. In this work we use an open source software program called pClust to cluster proteins from the complete proteomes of twelve species of Alphaproteobacteria and generate a dendrogram from the resulting orthologous protein clusters. We compare the results with dendrograms constructed using the 16S rRNA gene and multiple sequence alignment of seven housekeeping genes. Analysis of the whole proteomes of these pathogens grouped Rickettsia typhi with three other animal pathogens whereas conventional sequence analysis failed to group these pathogens together. We conclude that whole-proteome analysis can give insight into relationships among species beyond their phylogeny, perhaps reflecting the effects of horizontal gene transfer and potentially providing insight into the functions of shared genes by means of shared phenotypes.

Stalking the fourth domain in metagenomic data: searching for, discovering, and interpreting novel, deep branches in marker gene phylogenetic trees

Background

Most of our knowledge about the ancient evolutionary history of organisms has been derived from data associated with specific known organisms (i.e., organisms that we can study directly such as plants, metazoans, and culturable microbes). Recently, however, a new source of data for such studies has arrived: DNA sequence data generated directly from environmental samples. Such metagenomic data has enormous potential in a variety of areas including, as we argue here, in studies of very early events in the evolution of gene families and of species.

Methodology/Principal Findings

We designed and implemented new methods for analyzing metagenomic data and used them to search the Global Ocean Sampling (GOS) Expedition data set for novel lineages in three gene families commonly used in phylogenetic studies of known and unknown organisms: small subunit rRNA and the recA and rpoB superfamilies. Though the methods available could not accurately identify very deeply branched ss-rRNAs (largely due to difficulties in making robust sequence alignments for novel rRNA fragments), our analysis revealed the existence of multiple novel branches in the recA and rpoB gene families. Analysis of available sequence data likely from the same genomes as these novel recA and rpoB homologs was then used to further characterize the possible organismal source of the novel sequences.

Conclusions/Significance

Of the novel recA and rpoB homologs identified in the metagenomic data, some likely come from uncharacterized viruses while others may represent ancient paralogs not yet seen in any cultured organism. A third possibility is that some come from novel cellular lineages that are only distantly related to any organisms for which sequence data is currently available. If there exist any major, but so-far-undiscovered, deeply branching lineages in the tree of life, we suggest that methods such as those described herein currently offer the best way to search for them.

Anchor-based whole genome phylogeny (ABWGP): a tool for inferring evolutionary relationship among closely related microorganisms [corrected]

Phenotypic behavior of a group of organisms can be studied using a range of molecular evolutionary tools that help to determine evolutionary relationships. Traditionally a gene or a set of gene sequences was used for generating phylogenetic trees. Incomplete evolutionary information in few selected genes causes problems in phylogenetic tree construction. Whole genomes are used as remedy. Now, the task is to identify the suitable parameters to extract the hidden information from whole genome sequences that truly represent evolutionary information. In this study we explored a random anchor (a stretch of 100 nucleotides) based approach (ABWGP) for finding distance between any two genomes, and used the distance estimates to compute evolutionary trees. A number of strains and species of Mycobacteria were used for this study. Anchor-derived parameters, such as cumulative normalized score, anchor order and indels were computed in a pair-wise manner, and the scores were used to compute distance/phylogenetic trees. The strength of branching was determined by bootstrap analysis. The terminal branches are clearly discernable using the distance estimates described here. In general, different measures gave similar trees except the trees based on indels. Overall the tree topology reflected the known biology of the organisms. This was also true for different strains of Escherichia coli. A new whole genome-based approach has been described here for studying evolutionary relationships among bacterial strains and species.

A phylogeny-driven genomic encyclopaedia of Bacteria and Archaea

Sequencing of bacterial and archaeal genomes has revolutionized our understanding of the many roles played by microorganisms. There are now nearly 1,000 completed bacterial and archaeal genomes available, most of which were chosen for sequencing on the basis of their physiology. As a result, the perspective provided by the currently available genomes is limited by a highly biased phylogenetic distribution. To explore the value added by choosing microbial genomes for sequencing on the basis of their evolutionary relationships, we have sequenced and analysed the genomes of 56 culturable species of Bacteria and Archaea selected to maximize phylogenetic coverage. Analysis of these genomes demonstrated pronounced benefits (compared to an equivalent set of genomes randomly selected from the existing database) in diverse areas including the reconstruction of phylogenetic history, the discovery of new protein families and biological properties, and the prediction of functions for known genes from other organisms. Our results strongly support the need for systematic 'phylogenomic' efforts to compile a phylogeny-driven 'Genomic Encyclopedia of Bacteria and Archaea' in order to derive maximum knowledge from existing microbial genome data as well as from genome sequences to come.

The case for biocentric microbiology

Microbiology is a relatively modern scientific discipline intended to objectively study microorganisms, including pathogens and nonpathogens. However, since its birth, this science has been negatively affected by anthropocentric convictions, including rational and irrational beliefs. Among these, for example, is the artificial separation between environmental and medical microbiology that weakens both disciplines. Anthropocentric microbiology also fails to properly answer questions concerning the evolution of microbial pathogenesis. Here, I argue that an exclusively biocentric microbiology is imperative for improving our understanding not only of the microbial world, but also of our own species, our guts, and the world around us.

Parsimony score of phylogenetic networks: hardness results and a linear-time heuristic

Phylogenies-the evolutionary histories of groups of organisms-play a major role in representing the interrelationships among biological entities. Many methods for reconstructing and studying such phylogenies have been proposed, almost all of which assume that the underlying history of a given set of species can be represented by a binary tree. Although many biological processes can be effectively modeled and summarized in this fashion, others cannot: recombination, hybrid speciation, and horizontal gene transfer result in networks of relationships rather than trees of relationships. In previous works, we formulated a maximum parsimony (MP) criterion for reconstructing and evaluating phylogenetic networks, and demonstrated its quality on biological as well as synthetic data sets. In this paper, we provide further theoretical results as well as a very fast heuristic algorithm for the MP criterion of phylogenetic networks. In particular, we provide a novel combinatorial definition of phylogenetic networks in terms of "forbidden cycles," and provide detailed hardness and hardness of approximation proofs for the "small" MP problem. We demonstrate the performance of our heuristic in terms of time and accuracy on both biological and synthetic data sets. Finally, we explain the difference between our model and a similar one formulated by Nguyen et al., and describe the implications of this difference on the hardness and approximation results.

Does gene translocation accelerate the evolution of laterally transferred genes?

Lateral gene transfer (LGT) and gene rearrangement are essential for shaping bacterial genomes during evolution. Separate attention has been focused on understanding the process of lateral gene transfer and the process of gene translocation. However, little is known about how gene translocation affects laterally transferred genes. Here we have examined gene translocations and lateral gene transfers in closely related genome pairs. The results reveal that translocated genes undergo elevated rates of evolution and gene translocation tends to take place preferentially in recently acquired genes. Translocated genes have a high probability to be truncated, suggesting that translocation followed by truncation/deletion might play an important role in the fast turnover of laterally transferred genes. Furthermore, more recently acquired genes have a higher proportion of genes on the leading strand, suggesting a strong strand bias of lateral gene transfer.

Lines of evidence for horizontal gene transfer of a phenazine producing operon into multiple bacterial species

Phenazines are secondary metabolites with broad-spectrum antibiotic activity against bacteria, fungi, and eukaryotes. In pseudomonad species, a conserved seven-gene phenazine operon (phzABCDEFG) is required for the conversion of chorismic acid to the broad-spectrum antibiotic phenazine-1-carboxylate. Previous analyses of genes involved in phenazine production from nonpseudomonad species uncovered a high degree of sequence similarity to pseudomonad homologues. The analyses undertaken in this study wished to eluciadate the evolutionary history of genes involved in the production of phenazines. Furthermore, I wanted to determine if the phenazine operon has been transferred through horizontal gene transfer. Analyses of GC content, codon usage patterns, frequency of 3:1 dinucleotides, sequence similarities, and phylogenetic reconstructions were undertaken to map the evolutionary history of phenazine genes from multiple bacterial species. Patchy phyletic distribution, high sequence similarities, and phylogenetic evidence infer that pseudomonad, Streptomyces cinnamonensis, Pantoea agglomerans, Burkholderia cepacia, Pectobacterium atrosepticum, Brevibacterium linens, and Mycobacterium abscessus species all contain a phenazine operon which has most likely been transferred among these species through horizontal gene transfer. The acquisition of an antibiotic-associated operon is significant, as it may increase the relative fitness of the recipient species.

Complete genome sequence of the aerobic CO-oxidizing thermophile Thermomicrobium roseum

In order to enrich the phylogenetic diversity represented in the available sequenced bacterial genomes and as part of an “Assembling the Tree of Life” project, we determined the genome sequence of Thermomicrobium roseum DSM 5159. T. roseum DSM 5159 is a red-pigmented, rod-shaped, Gram-negative extreme thermophile isolated from a hot spring that possesses both an atypical cell wall composition and an unusual cell membrane that is composed entirely of long-chain 1,2-diols. Its genome is composed of two circular DNA elements, one of 2,006,217 bp (referred to as the chromosome) and one of 919,596 bp (referred to as the megaplasmid). Strikingly, though few standard housekeeping genes are found on the megaplasmid, it does encode a complete system for chemotaxis including both chemosensory components and an entire flagellar apparatus. This is the first known example of a complete flagellar system being encoded on a plasmid and suggests a straightforward means for lateral transfer of flagellum-based motility. Phylogenomic analyses support the recent rRNA-based analyses that led to T. roseum being removed from the phylum Thermomicrobia and assigned to the phylum Chloroflexi. Because T. roseum is a deep-branching member of this phylum, analysis of its genome provides insights into the evolution of the Chloroflexi. In addition, even though this species is not photosynthetic, analysis of the genome provides some insight into the origins of photosynthesis in the Chloroflexi. Metabolic pathway reconstructions and experimental studies revealed new aspects of the biology of this species. For example, we present evidence that T. roseum oxidizes CO aerobically, making it the first thermophile known to do so. In addition, we propose that glycosylation of its carotenoids plays a crucial role in the adaptation of the cell membrane to this bacterium's thermophilic lifestyle. Analyses of published metagenomic sequences from two hot springs similar to the one from which this strain was isolated, show that close relatives of T. roseum DSM 5159 are present but have some key differences from the strain sequenced.

Evidence of recent interkingdom horizontal gene transfer between bacteria and Candida parapsilosis

BACKGROUND

To date very few incidences of interdomain gene transfer into fungi have been identified. Here, we used the emerging genome sequences of Candida albicans WO-1, Candida tropicalis, Candida parapsilosis, Clavispora lusitaniae, Pichia guilliermondii, and Lodderomyces elongisporus to identify recent interdomain HGT events. We refer to these as CTG species because they translate the CTG codon as serine rather than leucine, and share a recent common ancestor.

RESULTS

Phylogenetic and syntenic information infer that two C. parapsilosis genes originate from bacterial sources. One encodes a putative proline racemase (PR). Phylogenetic analysis also infers that there were independent transfers of bacterial PR enzymes into members of the Pezizomycotina, and protists. The second HGT gene in C. parapsilosis belongs to the phenazine F (PhzF) superfamily. Most CTG species also contain a fungal PhzF homolog. Our phylogeny suggests that the CTG homolog originated from an ancient HGT event, from a member of the proteobacteria. An analysis of synteny suggests that C. parapsilosis has lost the endogenous fungal form of PhzF, and subsequently reacquired it from a proteobacterial source. There is evidence that Schizosaccharomyces pombe and Basidiomycotina also obtained a PhzF homolog through HGT.

CONCLUSION

Our search revealed two instances of well-supported HGT from bacteria into the CTG clade, both specific to C. parapsilosis. Therefore, while recent interkingdom gene transfer has taken place in the CTG lineage, its occurrence is rare. However, our analysis will not detect ancient gene transfers, and we may have underestimated the global extent of HGT into CTG species.

Efficient parsimony-based methods for phylogenetic network reconstruction

MOTIVATION

Phylogenies--the evolutionary histories of groups of organisms-play a major role in representing relationships among biological entities. Although many biological processes can be effectively modeled as tree-like relationships, others, such as hybrid speciation and horizontal gene transfer (HGT), result in networks, rather than trees, of relationships. Hybrid speciation is a significant evolutionary mechanism in plants, fish and other groups of species. HGT plays a major role in bacterial genome diversification and is a significant mechanism by which bacteria develop resistance to antibiotics. Maximum parsimony is one of the most commonly used criteria for phylogenetic tree inference. Roughly speaking, inference based on this criterion seeks the tree that minimizes the amount of evolution. In 1990, Jotun Hein proposed using this criterion for inferring the evolution of sequences subject to recombination. Preliminary results on small synthetic datasets. Nakhleh et al. (2005) demonstrated the criterion's application to phylogenetic network reconstruction in general and HGT detection in particular. However, the naive algorithms used by the authors are inapplicable to large datasets due to their demanding computational requirements. Further, no rigorous theoretical analysis of computing the criterion was given, nor was it tested on biological data.

RESULTS

In the present work we prove that the problem of scoring the parsimony of a phylogenetic network is NP-hard and provide an improved fixed parameter tractable algorithm for it. Further, we devise efficient heuristics for parsimony-based reconstruction of phylogenetic networks. We test our methods on both synthetic and biological data (rbcL gene in bacteria) and obtain very promising results.

Evolution of sensory complexity recorded in a myxobacterial genome

Myxobacteria are single-celled, but social, eubacterial predators. Upon starvation they build multicellular fruiting bodies using a developmental program that progressively changes the pattern of cell movement and the repertoire of genes expressed. Development terminates with spore differentiation and is coordinated by both diffusible and cell-bound signals. The growth and development of Myxococcus xanthus is regulated by the integration of multiple signals from outside the cells with physiological signals from within. A collection of M. xanthus cells behaves, in many respects, like a multicellular organism. For these reasons M. xanthus offers unparalleled access to a regulatory network that controls development and that organizes cell movement on surfaces. The genome of M. xanthus is large (9.14 Mb), considerably larger than the other sequenced delta-proteobacteria. We suggest that gene duplication and divergence were major contributors to genomic expansion from its progenitor. More than 1,500 duplications specific to the myxobacterial lineage were identified, representing >15% of the total genes. Genes were not duplicated at random; rather, genes for cell-cell signaling, small molecule sensing, and integrative transcription control were amplified selectively. Families of genes encoding the production of secondary metabolites are overrepresented in the genome but may have been received by horizontal gene transfer and are likely to be important for predation.

Structural proteomics of minimal organisms: conservation of protein fold usage and evolutionary implications

BACKGROUND

Determining the complete repertoire of protein structures for all soluble, globular proteins in a single organism has been one of the major goals of several structural genomics projects in recent years.

RESULTS

We report that this goal has nearly been reached for several "minimal organisms"--parasites or symbionts with reduced genomes--for which over 95% of the soluble, globular proteins may now be assigned folds, overall 3-D backbone structures. We analyze the structures of these proteins as they relate to cellular functions, and compare conservation of fold usage between functional categories. We also compare patterns in the conservation of folds among minimal organisms and those observed between minimal organisms and other bacteria.

CONCLUSION

We find that proteins performing essential cellular functions closely related to transcription and translation exhibit a higher degree of conservation in fold usage than proteins in other functional categories. Folds related to transcription and translation functional categories were also overrepresented in minimal organisms compared to other bacteria.

Genome trees from conservation profiles

The concept of the genome tree depends on the potential evolutionary significance in the clustering of species according to similarities in the gene content of their genomes. In this respect, genome trees have often been identified with species trees. With the rapid expansion of genome sequence data it becomes of increasing importance to develop accurate methods for grasping global trends for the phylogenetic signals that mutually link the various genomes. We therefore derive here the methodological concept of genome trees based on protein conservation profiles in multiple species. The basic idea in this derivation is that the multi-component “presence-absence” protein conservation profiles permit tracking of common evolutionary histories of genes across multiple genomes. We show that a significant reduction in informational redundancy is achieved by considering only the subset of distinct conservation profiles. Beyond these basic ideas, we point out various pitfalls and limitations associated with the data handling, paving the way for further improvements. As an illustration for the methods, we analyze a genome tree based on the above principles, along with a series of other trees derived from the same data and based on pair-wise comparisons (ancestral duplication-conservation and shared orthologs). In all trees we observe a sharp discrimination between the three primary domains of life: Bacteria, Archaea, and Eukarya. The new genome tree, based on conservation profiles, displays a significant correspondence with classically recognized taxonomical groupings, along with a series of departures from such conventional clusterings.

Since Darwin's Origin of Species and Haeckel's Tree of Life, systematic biology has attempted to classify species into “family trees.” Genomics has provided a new framework permitting descriptions of sibling relations between species on the basis of their complete genetic blueprints. While trees based on single genes (rRNA), or limited numbers of genes have been useful, genome trees derived from complete genome comparisons should lead to more complete pictures of phylogenetic relations between various organisms. In order to reach such a global vision, procedures to establish sibling relationships should depend on an overall comparison that captures the evolutionary fates of proteins jointly in multiple genomes. This paper aims to establish a methodological basis to use genuine multidimensional procedures in the construction of genome trees. This approach completes the derivation of trees based on more classical techniques of pair-wise comparison between species. The authors survey classification schemes emerging from this approach, which either supports traditional views, such as the separation between the three phylogenetic domains Bacteria, Archaea, and Eukarya, or challenges them by suggesting, for example, intermingled clusterings of Proteobacteria with various other bacterial species.

Decoding the genomic tree of life

Genomes hold within them the record of the evolution of life on Earth. But genome fusions and horizontal gene transfer (HGT) seem to have obscured sufficiently the gene sequence record such that it is difficult to reconstruct the phylogenetic tree of life. HGT among prokaryotes is not random, however. Some genes (informational genes) are more difficult to transfer than others (operational genes). Furthermore, environmental, metabolic, and genetic differences among organisms restrict HGT, so that prokaryotes preferentially share genes with other prokaryotes having properties in common, including genome size, genome G+C composition, carbon utilization, oxygen utilization/sensitivity, and temperature optima, further complicating attempts to reconstruct the tree of life. A new method of phylogenetic reconstruction based on gene presence and absence, called conditioned reconstruction, has improved our prospects for reconstructing prokaryotic evolution. It is also able to detect past genome fusions, such as the fusion that appears to have created the first eukaryote. This genome fusion between a deep branching eubacterium, possibly an ancestor of the cyanobacterium and a proteobacterium, with an archaeal eocyte (crenarchaea), appears to be the result of an early symbiosis. Given new tools and new genes from relevant organisms, it should soon be possible to test current and future fusion theories for the origin of eukaryotes and to discover the general outlines of the prokaryotic tree of life.

Integron-associated Mobile Gene Cassettes Code for Folded Proteins: The Structure of Bal32a, a New Member of the Adaptable α+β Barrel Family

The wide-ranging physiology and large genetic variability observed for prokaryotes is largely attributed, not to the prokaryotic genome itself, but rather to mechanisms of lateral gene transfer. Cassette PCR has been used to sample the integron/gene cassette metagenome from different natural environments without laboratory cultivation of the host organism, and without prior knowledge of any target protein sequence. Since over 90% of cassette genes are unrelated to any sequence in the current databases, it is not clear whether these genes code for folded functional proteins. We have selected a sample of eight cassette-encoded genes with no known homologs; five have been isolated as soluble protein products and shown by biophysical techniques to be folded. In solution, at least three of these proteins organise as stable oligomeric assemblies. The tertiary structure of one of these, Bal32a derived from a contaminated soil site, has been solved by X-ray crystallography to 1.8 A resolution. From the three-dimensional structure, Bal32a is found to be a member of the highly adaptable alpha+beta barrel family of transport proteins and enzymes. In Bal32a, the barrel cavity is unusually deep and inaccessible to solvent. Polar side-chains in its interior are reminiscent of catalytic sites of limonene-1,2-epoxide hydrolase and nogalonic acid methyl ester cyclase. These studies demonstrate the viability of direct sampling of mobile DNA as a route for the discovery of novel proteins.

Detection and characterization of horizontal transfers in prokaryotes using genomic signature

Horizontal DNA transfer is an important factor of evolution and participates in biological diversity. Unfortunately, the location and length of horizontal transfers (HTs) are known for very few species. The usage of short oligonucleotides in a sequence (the so-called genomic signature) has been shown to be species-specific even in DNA fragments as short as 1 kb. The genomic signature is therefore proposed as a tool to detect HTs. Since DNA transfers originate from species with a signature different from those of the recipient species, the analysis of local variations of signature along recipient genome may allow for detecting exogenous DNA. The strategy consists in (i) scanning the genome with a sliding window, and calculating the corresponding local signature (ii) evaluating its deviation from the signature of the whole genome and (iii) looking for similar signatures in a database of genomic signatures. A total of 22 prokaryote genomes are analyzed in this way. It has been observed that atypical regions make up ∼6% of each genome on the average. Most of the claimed HTs as well as new ones are detected. The origin of putative DNA transfers is looked for among ∼12 000 species. Donor species are proposed and sometimes strongly suggested, considering similarity of signatures. Among the species studied, Bacillus subtilis, Haemophilus Influenzae and Escherichia coli are investigated by many authors and give the opportunity to perform a thorough comparison of most of the bioinformatics methods used to detect HTs.

Inter-genomic displacement via lateral gene transfer of bacterial trp operons in an overall context of vertical genealogy

Background

The growing conviction that lateral gene transfer plays a significant role in prokaryote genealogy opens up a need for comprehensive evaluations of gene-enzyme systems on a case-by-case basis. Genes of tryptophan biosynthesis are frequently organized as whole-pathway operons, an attribute that is expected to facilitate multi-gene transfer in a single step. We have asked whether events of lateral gene transfer are sufficient to have obscured our ability to track the vertical genealogy that underpins tryptophan biosynthesis.

Results

In 47 complete-genome Bacteria, the genes encoding the seven catalytic domains that participate in primary tryptophan biosynthesis were distinguished from any paralogs or xenologs engaged in other specialized functions. A reliable list of orthologs with carefully ascertained functional roles has thus been assembled and should be valuable as an annotation resource. The protein domains associated with primary tryptophan biosynthesis were then concatenated, yielding single amino-acid sequence strings that represent the entire tryptophan pathway. Lateral gene transfer of several whole-pathway trp operons was demonstrated by use of phylogenetic analysis. Lateral gene transfer of partial-pathway trp operons was also shown, with newly recruited genes functioning either in primary biosynthesis (rarely) or specialized metabolism (more frequently).

Conclusions

(i) Concatenated tryptophan protein trees are congruent with 16S rRNA subtrees provided that the genomes represented are of sufficiently close phylogenetic spacing. There are currently seven tryptophan congruency groups in the Bacteria. Recognition of a succession of others can be expected in the near future, but ultimately these should coalesce to a single grouping that parallels the 16S rRNA tree (except for cases of lateral gene transfer). (ii) The vertical trace of evolution for tryptophan biosynthesis can be deduced. The daunting complexities engendered by paralogy, xenology, and idiosyncrasies of nomenclature at this point in time have necessitated an expert-assisted manual effort to achieve a correct analysis. Once recognized and sorted out, paralogy and xenology can be viewed as features that enrich evolutionary histories.

Deriving the Genomic Tree of Life in the Presence of Horizontal Gene Transfer: Conditioned Reconstruction

The horizontal gene transfer (HGT) being inferred within prokaryotic genomes appears to be sufficiently massive that many scientists think it may have effectively obscured much of the history of life recorded in DNA. Here, we demonstrate that the tree of life can be reconstructed even in the presence of extensive HGT, provided the processes of genome evolution are properly modeled. We show that the dynamic deletions and insertions of genes that occur during genome evolution, including those introduced by HGT, may be modeled using techniques similar to those used to model nucleotide substitutions that occur during sequence evolution. In particular, we show that appropriately designed general Markov models are reasonable tools for reconstructing genome evolution. These studies indicate that, provided genomes contain sufficiently many genes and that the Markov assumptions are met, it is possible to reconstruct the tree of life. We also consider the fusion of genomes, a process not encountered in gene sequence evolution, and derive a method for the identification and reconstruction of genome fusion events. Genomic reconstructions of a well-defined classical four-genome problem, the root of the multicellular animals, show that the method, when used in conjunction with paralinear/logdet distances, performs remarkably well and is relatively unaffected by the recently discovered big genome artifact.

Sequence analysis of the mobile genome island pKLC102 of Pseudomonas aeruginosa C

The Pseudomonas aeruginosa plasmid pKLC102 coexists as a plasmid and a genome island in clone C strains. Whereas the related plasmid pKLK106 reversibly recombines with P. aeruginosa clone K chromosomes at one of the two tRNA(Lys) genes, pKLC102 is incorporated into the tRNA(Lys) gene only close to the pilA locus. Targeting of the other tRNA(Lys) copy in the chromosome is blocked by a 23,395-bp mosaic of truncated PAO open reading frames, transposons, and pKLC102 homologs. Annotation and phylogenetic analysis of the large 103,532-bp pKLC102 sequence revealed that pKLC102 is a hybrid of plasmid and phage origin. The plasmid lineage conferred oriV and genes for replication, partitioning, and conjugation, including a pil cluster encoding type IV thin sex pili and an 8,524-bp chvB glucan synthetase gene that is known to be a major determinant for host tropism and virulence. The phage lineage conferred integrase, att, and a syntenic set of conserved hypothetical genes also observed in the tRNA(Gly)-associated genome islands of P. aeruginosa clone C chromosomes. In subgroup C isolates from patients with cystic fibrosis, pKLC102 was irreversibly fixed into the chromosome by the insertion of the large 23,061-bp class I transposon TNCP23, which is a composite of plasmid, integron, and IS6100 elements. Intramolecular transposition of a copy of IS6100 led to chromosomal inversions and disruption of plasmid synteny. The case of pKLC102 in P. aeruginosa clone C documents the intraclonal evolution of a genome island from a mobile ancestor via a reversibly integrated state to irreversible incorporation and dissipation in the chromosome.

Horizontal gene transfer and phylogenetics

The initial analysis of complete genomes has suggested that horizontal gene transfer events are very frequent between microorganisms. This could potentially render the inference, and even the concept itself, of the organismal phylogeny impossible. However, a coherent phylogenetic pattern has recently emerged from an analysis of about a hundred genes, the so-called 'core', strongly suggesting that it is possible to infer the phylogeny of prokaryotes. Also, estimation of the frequency of horizontal gene transfers at the genome level in a phylogenetic context seems to indicate that it is rather low, although of significant biological impact. Nevertheless, it should be emphasized that the history of microorganisms cannot be properly represented by the phylogeny of the core, which represents only a tiny fraction of the genome. This history, even if horizontal gene transfers are rare, should be represented by a network surrounding the core phylogeny.

Evolution of the soluble diiron monooxygenases

Based on structural, biochemical, and genetic data, the soluble diiron monooxygenases can be divided into four groups: the soluble methane monooxygenases, the Amo alkene monooxygenase of Rhodococcus corallinus B-276, the phenol hydroxylases, and the four-component alkene/aromatic monooxygenases. The limited phylogenetic distribution of these enzymes among bacteria, together with available genetic evidence, indicates that they have been spread largely through horizontal gene transfer. Phylogenetic analyses reveal that the alpha- and beta-oxygenase subunits are paralogous proteins and were derived from an ancient gene duplication of a carboxylate-bridged diiron protein, with subsequent divergence yielding a catalytic alpha-oxygenase subunit and a structural beta-oxygenase subunit. The oxidoreductase and ferredoxin components of these enzymes are likely to have been acquired by horizontal transfer from ancestors common to unrelated diiron and Rieske center oxygenases and other enzymes. The cumulative results of phylogenetic reconstructions suggest that the alkene/aromatic monooxygenases diverged first from the last common ancestor for these enzymes, followed by the phenol hydroxylases, Amo alkene monooxygenase, and methane monooxygenases.

Extracting phylogenetic information from whole-genome sequencing projects: the lactic acid bacteria as a test case

The availability of an ever increasing number of complete genome sequences of diverse prokaryotic taxa has led to the introduction of novel approaches to infer phylogenetic relationships among bacteria. In the present study the sequences of the 16S rRNA gene and nine housekeeping genes were compared with the fraction of shared putative orthologous protein-encoding genes, conservation of gene order, dinucleotide relative abundance and codon usage among 11 genomes of species belonging to the lactic acid bacteria. In general there is a good correlation between the results obtained with various approaches, although it is clear that there is a stronger phylogenetic signal in some datasets than in others, and that different parameters have different taxonomic resolutions. It appears that trees based on different kinds of information derived from whole-genome sequencing projects do not provide much additional information about the phylogenetic relationships among bacterial taxa compared to more traditional alignment-based methods. Nevertheless, it is expected that the study of these novel forms of information will have its value in taxonomy, to determine which genes are shared, when genes or sets of genes were lost in evolutionary history, to detect the presence of horizontally transferred genes and/or confirm or enhance the phylogenetic signal derived from traditional methods. Although these conclusions are based on a relatively small dataset, they are largely in agreement with other studies and it is anticipated that similar trends will be observed when comparing other genomes.

On dynamics of overlapping genes in bacterial genomes

A genome-wide computational analysis of overlapping genes in various prokaryotic genomes was conducted to study the evolutionary rates of the formation and degradation of overlapping genes. Overlapping genes are adjacent genes that overlap partially or entirely. Overlapping gene pairs extracted from complete bacterial genome sequences were classified into three directional patterns, namely, 'convergent' (--><--), 'unidirectional' (-->-->), and 'divergent' (<---->). Divergent overlapping genes were rare compared to the other two patterns. We also compared the structures of overlapping genes that are not conserved in the homologous genes of different species. We found that overlapping genes are mainly generated from gene pairs with short intergenic regions. Overlapping genes caused by mutations at the 3'-end of genes occur more frequently than overlapping genes generated by mutations at the 5'-end of genes in closely related species. In contrast, gene overlaps due to mutations at the 5'-ends accumulate in distant species. These differences presumably correlate to the genomic divergence among species. The evolution of these overlapping gene structures is probably related to the evolutionary time scale, which suggests that the evolution of overlapping genes occurs at a universal mutation rate across species.

Traffic at the tmRNA gene

A partial screen for genetic elements integrated into completely sequenced bacterial genomes shows more significant bias in specificity for the tmRNA gene (ssrA) than for any type of tRNA gene. Horizontal gene transfer, a major avenue of bacterial evolution, was assessed by focusing on elements using this single attachment locus. Diverse elements use ssrA; among enterobacteria alone, at least four different integrase subfamilies have independently evolved specificity for ssrA, and almost every strain analyzed presents a unique set of integrated elements. Even elements using essentially the same integrase can be very diverse, as is a group with an ssrA-specific integrase of the P4 subfamily. This same integrase appears to promote damage routinely at attachment sites, which may be adaptive. Elements in arrays can recombine; one such event mediated by invertible DNA segments within neighboring elements likely explains the monophasic nature of Salmonella enterica serovar Typhi. One of a limited set of conserved sequences occurs at the attachment site of each enterobacterial element, apparently serving as a transcriptional terminator for ssrA. Elements were usually found integrated into tRNA-like sequence at the 3' end of ssrA, at subsites corresponding to those used in tRNA genes; an exception was found at the non-tRNA-like 3' end produced by ssrA gene permutation in cyanobacteria, suggesting that, during the evolution of new site specificity by integrases, tropism toward a conserved 3' end of an RNA gene may be as strong as toward a tRNA-like sequence. The proximity of ssrA and smpB, which act in concert, was also surveyed.

Evolutionary implications of microbial genome tetranucleotide frequency biases

We compared nucleotide usage pattern conservation for related prokaryotes by examining the representation of DNA tetranucleotide combinations in 27 representative microbial genomes. For each of the organisms studied, tetranucleotide usage departures from expectations (TUD) were shared between related organisms using both Markov chain analysis and a zero-order Markov method. Individual strains, multiple chromosomes, plasmids, and bacteriophages share TUDs within a species. TUDs varied between coding and noncoding DNA. Grouping prokaryotes based on TUD profiles resulted in relationships with important differences from those based on 16S rRNA phylogenies, which may reflect unequal rates of evolution of nucleotide usage patterns following divergence of particular organisms from a common ancestor. By both symmetrical tree distance and likelihood analysis, phylogenetic trees based on TUD profiles demonstrate a level of congruence with 16S rRNA trees similar to that of both RpoA and RecA trees. Congruence of these trees indicates that there exists phylogenetic signal in TUD patterns, most prominent in coding region DNA. Because relationships demonstrated in TUD-based analyses utilize whole genomes, they should be considered complementary to phylogenies based on single genetic elements, such as 16S rRNA.

Factors influencing the identification of transcription factor binding sites by cross-species comparison

As the number of sequenced genomes has grown, the questions of which species are most useful and how many genomes are sufficient for comparison have become increasingly important for comparative genomics studies. We have systematically addressed these questions with respect to phylogenetic footprinting of transcription factor (TF) binding sites in the gamma-proteobacteria, and have evaluated the statistical significance of our motif predictions. We used a study set of 166 Escherichia coli genes that have experimentally identified TF binding sites upstream of the gene, with orthologous data from nine additional gamma-proteobacteria for phylogenetic footprinting. Just three species were sufficient for approximately 74.0% of the motif predictions to correspond to the experimentally reported E. coli sites, and important characteristics to consider when choosing species were phylogenetic distance, genome size, and natural habitat. We also performed simulations using randomized data to determine the critical maximum a posteriori probability (MAP) values for statistical significance of our motif predictions (P = 0.05). Approximately 60% of motif predictions containing sites from just three species had average MAP values above these critical MAP values. The inclusion of a species very closely related to E. coli increased the number of statistically significant motif predictions, despite substantially increasing the critical MAP value.

Metabolic engineering of Lactococcus lactis: the impact of genomics and metabolic modelling

Lactic acid bacteria display a relatively simple and well described metabolism where the sugar source is converted mainly to lactic acid. Here we will shortly describe metabolic engineering strategies that led to the efficient re-routing of the lactococcal pyruvate metabolism to end-products other than lactic acid, including diacetyl and alanine. Moreover, we will review current metabolic engineering approaches that aim at increasing the flux through complex biosynthetic pathways, leading to exopolysaccharides and folic acid. Finally, the (future) impact of the developments in the area of genomics and corresponding high-throughput technologies will be discussed.

A phylogenomic approach to bacterial phylogeny: evidence of a core of genes sharing a common history

It has been claimed that complete genome sequences would clarify phylogenetic relationships between organisms, but up to now, no satisfying approach has been proposed to use efficiently these data. For instance, if the coding of presence or absence of genes in complete genomes gives interesting results, it does not take into account the phylogenetic information contained in sequences and ignores hidden paralogies by using a BLAST reciprocal best hit definition of orthology. In addition, concatenation of sequences of different genes as well as building of consensus trees only consider the few genes that are shared among all organisms. Here we present an attempt to use a supertree method to build the phylogenetic tree of 45 organisms, with special focus on bacterial phylogeny. This led us to perform a phylogenetic study of congruence of tree topologies, which allows the identification of a core of genes supporting similar species phylogeny. We then used this core of genes to infer a tree. This phylogeny presents several differences with the rRNA phylogeny, notably for the position of hyperthermophilic bacteria.

Pyrococcus genome comparison evidences chromosome shuffling-driven evolution

The genomes of three Pyrococcus species, P.abyssi, P.furiosus and P.horikoshii, were compared at the DNA level, taking advantage of our identification of their replication origins. Three types of rearrangements have been identified: (i) inversion and translation across the replication axis (origin/terminus), (ii) inversion and translocation restricted to a replichore (the half chromosome divided by the replication axis) and (iii) apparent mobility of long clusters of repeated sequences. Rearrangements restricted within a replichore were more common between P.furiosus and the two other Pyrococcus species than between P.horikoshii and P.abyssi. A strong correlation was found between 23 homologous insertion sequence elements, present only in P.furiosus, and recombined segment boundaries, suggesting that transposition events have been a major cause of genomic disruption in this species. Moreover, gene orientation bias was much more disrupted than strand composition biases in fragments that switched their orientation within a replichore upon recombination. This allowed us to conclude that one reversion and one translation occurred in P.abyssi after its divergence from P.horikoshii, and that a smaller segment has specifically recombined in P.furiosus. Whereas a majority of genes are transcribed in the same direction as DNA replication in P.horikoshii and P.abyssi, the colinearity of transcription and replication is only maintained for highly transcribed genes in P.furiosus. We discuss the implications of genomic rearrangements on gene orientation and composition biases, and their consequences on sequence evolution.

Evidence for lateral transfer of genes encoding ferredoxins, nitroreductases, NADH oxidase, and alcohol dehydrogenase 3 from anaerobic prokaryotes to Giardia lamblia and Entamoeba histolytica

Giardia lamblia and Entamoeba histolytica are amitochondriate, microaerophilic protists which use fermentation enzymes like those of bacteria to survive anaerobic conditions within the intestinal lumen. Genes encoding fermentation enzymes and related electron transport peptides (e.g., ferredoxins) in giardia organisms and amebae are hypothesized to be derived from either an ancient anaerobic eukaryote (amitochondriate fossil hypothesis), a mitochondrial endosymbiont (hydrogen hypothesis), or anaerobic bacteria (lateral transfer hypothesis). The goals here were to complete the molecular characterization of giardial and amebic fermentation enzymes and to determine the origins of the genes encoding them, when possible. A putative giardia [2Fe-2S]ferredoxin which had a hypothetical organelle-targeting sequence at its N terminus showed similarity to mitochondrial ferredoxins and the hydrogenosomal ferredoxin of Trichomonas vaginalis (another luminal protist). However, phylogenetic trees were star shaped, with weak bootstrap support, so we were unable to confirm or rule out the endosymbiotic origin of the giardia [2Fe-2S]ferredoxin gene. Putative giardial and amebic 6-kDa ferredoxins, ferredoxin-nitroreductase fusion proteins, and oxygen-insensitive nitroreductases each tentatively supported the lateral transfer hypothesis. Although there were not enough sequences to perform meaningful phylogenetic analyses, the unique common occurrence of these peptides and enzymes in giardia organisms, amebae, and the few anaerobic prokaryotes suggests the possibility of lateral transfer. In contrast, there was more robust phylogenetic evidence for the lateral transfer of G. lamblia genes encoding an NADH oxidase from a gram-positive coccus and a microbial group 3 alcohol dehydrogenase from thermoanaerobic prokaryotes. In further support of lateral transfer, the G. lamblia NADH oxidase and adh3 genes appeared to have an evolutionary history distinct from those of E. histolytica.

Multiple alignment of complete sequences (MACS) in the post-genomic era

Multiple alignment, since its introduction in the early seventies, has become a cornerstone of modern molecular biology. It has traditionally been used to deduce structure / function by homology, to detect conserved motifs and in phylogenetic studies. There has recently been some renewed interest in the development of multiple alignment techniques, with current opinion moving away from a single all-encompassing algorithm to iterative and / or co-operative strategies. The exploitation of multiple alignments in genome annotation projects represents a qualitative leap in the functional analysis process, opening the way to the study of the co-evolution of validated sets of proteins and to reliable phylogenomic analysis. However, the alignment of the highly complex proteins detected by today's advanced database search methods is a daunting task. In addition, with the explosion of the sequence databases and with the establishment of numerous specialized biological databases, multiple alignment programs must evolve if they are to successfully rise to the new challenges of the post-genomic era. The way forward is clearly an integrated system bringing together sequence data, knowledge-based systems and prediction methods with their inherent unreliability. The incorporation of such heterogeneous, often non-consistent, data will require major changes to the fundamental alignment algorithms used to date. Such an integrated multiple alignment system will provide an ideal workbench for the validation, propagation and presentation of this information in a format that is concise, clear and intuitive.