Whole-genome analysis of photosynthetic prokaryotes

Insight in to the phylogeny of polyhydroxyalkanoate biosynthesis: Horizontal gene transfer

Polyhydroxyalkanoates (PHAs) are gaining more and more importance the world over due to their structural diversity and close analogy to plastics. Their biodegradability makes them extremely desirable substitutes for synthetic plastics. PHAs are produced in organisms under certain stress conditions. Here, we investigated 253 sequenced (completely and unfinished) genomes for the diversity and phylogenetics of the PHA biosynthesis. Discrepancies in the phylogenetic trees for phaA, phaB and phaC genes of the PHA biosynthesis have led to the suggestion that horizontal gene transfer (HGT) may be a major contributor for its evolution. Twenty four organisms belonging to diverse taxa were found to be involved in HGT. Among these, Bacillus cereus ATCC 14579 and Xanthomonas axonopodis pv. citri str. 306 seem to have acquired all the three genes through HGT events and have not been characterized so far as PHA producers. This study also revealed certain potential organisms such as Streptomyces coelicolor A3(2), Staphylococcus epidermidis ATCC 12228, Brucella suis 1330, Burkholderia sp., DSMZ 9242 and Leptospira interrogans serovar lai str. 56601, which can be transformed into novel PHA producers through recombinant DNA technology.

Abundant presence of the gamma-like Proteobacterial pufM gene in oxic seawater

Known anoxygenic photosynthetic bacteria (APB) affiliated to Gammaproteobacteria usually use anaerobic metabolism and are restricted to oxygen-free habitats. Here, we report abundant (average of 34.5%) presence of diverse APB related to gamma-like Proteobacteria in oxic oceanic surface water as indicated by the pufM gene, that encodes the M subunit of the light reaction centre complex. Thus, our sequences were most likely derived from aerobic anoxygenic phototrophs (AAnP). Two genetically distinct genotypes were revealed: one was from the oligotrophic North Pacific Ocean Gyre and the other, was from the trophic East China Sea and Bering Sea. The discovery of abundant presence of novel gamma-like Proteobacterial pufM gene in the oxic seawater extends the functional ecotypes of AAnP.

Searching for convergence in phylogenetic Markov chain Monte Carlo

Markov chain Monte Carlo (MCMC) is a methodology that is gaining widespread use in the phylogenetics community and is central to phylogenetic software packages such as MrBayes. An important issue for users of MCMC methods is how to select appropriate values for adjustable parameters such as the length of the Markov chain or chains, the sampling density, the proposal mechanism, and, if Metropolis-coupled MCMC is being used, the number of heated chains and their temperatures. Although some parameter settings have been examined in detail in the literature, others are frequently chosen with more regard to computational time or personal experience with other data sets. Such choices may lead to inadequate sampling of tree space or an inefficient use of computational resources. We performed a detailed study of convergence and mixing for 70 randomly selected, putatively orthologous protein sets with different sizes and taxonomic compositions. Replicated runs from multiple random starting points permit a more rigorous assessment of convergence, and we developed two novel statistics, delta and epsilon, for this purpose. Although likelihood values invariably stabilized quickly, adequate sampling of the posterior distribution of tree topologies took considerably longer. Our results suggest that multimodality is common for data sets with 30 or more taxa and that this results in slow convergence and mixing. However, we also found that the pragmatic approach of combining data from several short, replicated runs into a "metachain" to estimate bipartition posterior probabilities provided good approximations, and that such estimates were no worse in approximating a reference posterior distribution than those obtained using a single long run of the same length as the metachain. Precision appears to be best when heated Markov chains have low temperatures, whereas chains with high temperatures appear to sample trees with high posterior probabilities only rarely.

Occurrence of hydrogenases in cyanobacteria and anoxygenic photosynthetic bacteria: implications for the phylogenetic origin of cyanobacterial and algal hydrogenases

Hydrogenases are important enzymes in the energy metabolism of microorganisms. Therefore, they are widespread in prokaryotes. We analyzed the occurrence of hydrogenases in cyanobacteria and deduced a FeFe-hydrogenase in three different heliobacterial strains. This allowed the first phylogenetic analysis of the hydrogenases of all five major groups of photosynthetic bacteria (heliobacteria, green nonsulfur bacteria, green sulfur bacteria, photosynthetic proteobacteria, and cyanobacteria). In the case of both hydrogenases found in cyanobacteria (uptake and bidirectional), the green nonsulfur bacterium Chloroflexus aurantiacus was found to be the closest ancestor. Apart from a close relation between the archaebacterial and the green sulfur bacterial sulfhydrogenase, we could not find any evidence for horizontal gene transfer. Therefore, it would be most parsimonious if a Chloroflexus-like bacterium was the ancestor of Chloroflexus aurantiacus and cyanobacteria. After having transmitted both hydrogenase genes vertically to the different cyanobacterial species, either no, one, or both enzymes were lost, thus producing the current distribution. Our data and the available data from the literature on the occurrence of cyanobacterial hydrogenases show that the cyanobacterial uptake hydrogenase is strictly linked to the occurrence of the nitrogenase. Nevertheless, we did identify a nitrogen-fixing Synechococcus strain without an uptake hydrogenase. Since we could not find genes of a FeFe-hydrogenase in any of the tested cyanobacteria, although strains performing anoxygenic photosynthesis were also included in the analysis, a cyanobacterial origin of the contemporary FeFe-hydrogenase of algal plastids seems unlikely.

Prokaryotic photosynthesis and phototrophy illuminated

Genome sequencing projects are revealing new information about the distribution and evolution of photosynthesis and phototrophy. Although coverage of the five phyla containing photosynthetic prokaryotes (Chlorobi, Chloroflexi, Cyanobacteria, Proteobacteria and Firmicutes) is limited and uneven, genome sequences are (or soon will be) available for >100 strains from these phyla. Present knowledge of photosynthesis is almost exclusively based on data derived from cultivated species but metagenomic studies can reveal new organisms with novel combinations of photosynthetic and phototrophic components that have not yet been described. Metagenomics has already shown how the relatively simple phototrophy based upon rhodopsins has spread laterally throughout Archaea, Bacteria and eukaryotes. In this review, we present examples that reflect recent advances in phototroph biology as a result of insights from genome and metagenome sequencing.

The cyanobacterial genome core and the origin of photosynthesis

Comparative analysis of 15 complete cyanobacterial genome sequences, including "near minimal" genomes of five strains of Prochlorococcus spp., revealed 1,054 protein families [core cyanobacterial clusters of orthologous groups of proteins (core CyOGs)] encoded in at least 14 of them. The majority of the core CyOGs are involved in central cellular functions that are shared with other bacteria; 50 core CyOGs are specific for cyanobacteria, whereas 84 are exclusively shared by cyanobacteria and plants and/or other plastid-carrying eukaryotes, such as diatoms or apicomplexans. The latter group includes 35 families of uncharacterized proteins, which could also be involved in photosynthesis. Only a few components of cyanobacterial photosynthetic machinery are represented in the genomes of the anoxygenic phototrophic bacteria Chlorobium tepidum, Rhodopseudomonas palustris, Chloroflexus aurantiacus, or Heliobacillus mobilis. These observations, coupled with recent geological data on the properties of the ancient phototrophs, suggest that photosynthesis originated in the cyanobacterial lineage under the selective pressures of UV light and depletion of electron donors. We propose that the first phototrophs were anaerobic ancestors of cyanobacteria ("procyanobacteria") that conducted anoxygenic photosynthesis using a photosystem I-like reaction center, somewhat similar to the heterocysts of modern filamentous cyanobacteria. From procyanobacteria, photosynthesis spread to other phyla by way of lateral gene transfer.

Data mining for proteins characteristic of clades

A synapomorphy is a phylogenetic character that provides evidence of shared descent. Ideally a synapomorphy is ubiquitous within the clade of related organisms and nonexistent outside the clade, implying that it arose after divergence from other extant species and before the last common ancestor of the clade. With the recent proliferation of genetic sequence data, molecular synapomorphies have assumed great importance, yet there is no convenient means to search for them over entire genomes. We have developed a new program called Conserv, which can rapidly assemble orthologous sequences and rank them by various metrics, such as degree of conservation or divergence from out-group orthologs. We have used Conserv to conduct a largescale search for molecular synapomorphies for bacterial clades. The search discovered sequences unique to clades, such as Actinobacteria, Firmicutes and gamma-Proteobacteria, and shed light on several open questions, such as whether Symbiobacterium thermophilum belongs with Actinobacteria or Firmicutes. We conclude that Conserv can quickly marshall evidence relevant to evolutionary questions that would be much harder to assemble with other tools.

Rooting the tree of life by transition analyses

BACKGROUND

Despite great advances in clarifying the family tree of life, it is still not agreed where its root is or what properties the most ancient cells possessed--the most difficult problems in phylogeny. Protein paralogue trees can theoretically place the root, but are contradictory because of tree-reconstruction artefacts or poor resolution; ribosome-related and DNA-handling enzymes suggested one between neomura (eukaryotes plus archaebacteria) and eubacteria, whereas metabolic enzymes often place it within eubacteria but in contradictory places. Palaeontology shows that eubacteria are much more ancient than eukaryotes, and, together with phylogenetic evidence that archaebacteria are sisters not ancestral to eukaryotes, implies that the root is not within the neomura. Transition analysis, involving comparative/developmental and selective arguments, can polarize major transitions and thereby systematically exclude the root from major clades possessing derived characters and thus locate it; previously the 20 shared neomuran characters were thus argued to be derived, but whether the root was within eubacteria or between them and archaebacteria remained controversial.

RESULTS

I analyze 13 major transitions within eubacteria, showing how they can all be congruently polarized. I infer the first fully resolved prokaryote tree, with a basal stem comprising the new infrakingdom Glidobacteria (Chlorobacteria, Hadobacteria, Cyanobacteria), which is entirely non-flagellate and probably ancestrally had gliding motility, and two derived branches (Gracilicutes and Unibacteria/Eurybacteria) that diverged immediately following the origin of flagella. Proteasome evolution shows that the universal root is outside a clade comprising neomura and Actinomycetales (proteates), and thus lies within other eubacteria, contrary to a widespread assumption that it is between eubacteria and neomura. Cell wall and flagellar evolution independently locate the root outside Posibacteria (Actinobacteria and Endobacteria), and thus among negibacteria with two membranes. Posibacteria are derived from Eurybacteria and ancestral to neomura. RNA polymerase and other insertions strongly favour the monophyly of Gracilicutes (Proteobacteria, Planctobacteria, Sphingobacteria, Spirochaetes). Evolution of the negibacterial outer membrane places the root within Eobacteria (Hadobacteria and Chlorobacteria, both primitively without lipopolysaccharide): as all phyla possessing the outer membrane beta-barrel protein Omp85 are highly probably derived, the root lies between them and Chlorobacteria, the only negibacteria without Omp85, or possibly within Chlorobacteria.

CONCLUSION

Chlorobacteria are probably the oldest and Archaebacteria the youngest bacteria, with Posibacteria of intermediate age, requiring radical reassessment of dominant views of bacterial evolution. The last ancestor of all life was a eubacterium with acyl-ester membrane lipids, large genome, murein peptidoglycan walls, and fully developed eubacterial molecular biology and cell division. It was a non-flagellate negibacterium with two membranes, probably a photosynthetic green non-sulphur bacterium with relatively primitive secretory machinery, not a heterotrophic posibacterium with one membrane.

Quantitative exploration of the occurrence of lateral gene transfer by using nitrogen fixation genes as a case study

Lateral gene transfer (LGT) is now accepted as an important factor in the evolution of prokaryotes. Establishment of the occurrence of LGT is typically attempted by a variety of methods that includes the comparison of reconstructed phylogenetic trees, the search for unusual GC composition or codon usage within a genome, and identification of similarities between distant species as determined by best blast hits. We explore quantitative assessments of these strategies to study the prokaryotic trait of nitrogen fixation, the enzyme-catalyzed reduction of N(2) to ammonia. Phylogenies constructed on nitrogen fixation genes are not in agreement with the tree-of-life based on 16S rRNA but do not conclusively distinguish between gene loss and LGT hypotheses. Using a series of analyses on a set of complete genomes, our results distinguish two structurally distinct classes of MoFe nitrogenases whose distribution cuts across lines of vertical inheritance and makes us believe that a conclusive case for LGT has been made.

Genome-wide detection and analysis of homologous recombination among sequenced strains of Escherichia coli

BACKGROUND

Comparisons of complete bacterial genomes reveal evidence of lateral transfer of DNA across otherwise clonally diverging lineages. Some lateral transfer events result in acquisition of novel genomic segments and are easily detected through genome comparison. Other more subtle lateral transfers involve homologous recombination events that result in substitution of alleles within conserved genomic regions. This type of event is observed infrequently among distantly related organisms. It is reported to be more common within species, but the frequency has been difficult to quantify since the sequences under comparison tend to have relatively few polymorphic sites.

RESULTS

Here we report a genome-wide assessment of homologous recombination among a collection of six complete Escherichia coli and Shigella flexneri genome sequences. We construct a whole-genome multiple alignment and identify clusters of polymorphic sites that exhibit atypical patterns of nucleotide substitution using a random walk-based method. The analysis reveals one large segment (approximately 100 kb) and 186 smaller clusters of single base pair differences that suggest lateral exchange between lineages. These clusters include portions of 10% of the 3,100 genes conserved in six genomes. Statistical analysis of the functional roles of these genes reveals that several classes of genes are over-represented, including those involved in recombination, transport and motility.

CONCLUSION

We demonstrate that intraspecific recombination in E. coli is much more common than previously appreciated and may show a bias for certain types of genes. The described method provides high-specificity, conservative inference of past recombination events.

Comparative genomics of NAD biosynthesis in cyanobacteria

Biosynthesis of NAD(P) cofactors is of special importance for cyanobacteria due to their role in photosynthesis and respiration. Despite significant progress in understanding NAD(P) biosynthetic machinery in some model organisms, relatively little is known about its implementation in cyanobacteria. We addressed this problem by a combination of comparative genome analysis with verification experiments in the model system of Synechocystis sp. strain PCC 6803. A detailed reconstruction of the NAD(P) metabolic subsystem using the SEED genomic platform (http://theseed.uchicago.edu/FIG/index.cgi) helped us accurately annotate respective genes in the entire set of 13 cyanobacterial species with completely sequenced genomes available at the time. Comparative analysis of operational variants implemented in this divergent group allowed us to elucidate both conserved (de novo and universal pathways) and variable (recycling and salvage pathways) aspects of this subsystem. Focused genetic and biochemical experiments confirmed several conjectures about the key aspects of this subsystem. (i) The product of the slr1691 gene, a homolog of Escherichia coli gene nadE containing an additional nitrilase-like N-terminal domain, is a NAD synthetase capable of utilizing glutamine as an amide donor in vitro. (ii) The product of the sll1916 gene, a homolog of E. coli gene nadD, is a nicotinic acid mononucleotide-preferring adenylyltransferase. This gene is essential for survival and cannot be compensated for by an alternative nicotinamide mononucleotide (NMN)-preferring adenylyltransferase (slr0787 gene). (iii) The product of the slr0788 gene is a nicotinamide-preferring phosphoribosyltransferase involved in the first step of the two-step non-deamidating utilization of nicotinamide (NMN shunt). (iv) The physiological role of this pathway encoded by a conserved gene cluster, slr0787-slr0788, is likely in the recycling of endogenously generated nicotinamide, as supported by the inability of this organism to utilize exogenously provided niacin. Positional clustering and the co-occurrence profile of the respective genes across a diverse collection of cellular organisms provide evidence of horizontal transfer events in the evolutionary history of this pathway.

The Rieske Protein: A Case Study on the Pitfalls of Multiple Sequence Alignments and Phylogenetic Reconstruction

Previously published phylogenetic trees reconstructed on "Rieske protein" sequences frequently are at odds with each other, with those of other subunits of the parent enzymes and with small-subunit rRNA trees. These differences are shown to be at least partially if not completely due to problems in the reconstruction procedures. A major source of erroneous Rieske protein trees lies in the presence of a large, poorly conserved domain prone to accommodate very long insertions in well-defined structural hot spots substantially hampering multiple alignments. The remaining smaller domain, in contrast, is too conserved to allow distant phylogenies to be deduced with sufficient confidence. Three-dimensional structures of representatives from this protein family are now available from phylogenetically distant species and from diverse enzymes. Multiple alignments can thus be refined on the basis of these structures. We show that structurally guided alignments of Rieske proteins from Rieske-cytochrome b complexes and arsenite oxidases strongly reduce conflicts between resulting trees and those obtained on their companion enzyme subunits. Further problems encountered during this work, mainly consisting in database errors such as wrong annotations and frameshifts, are described. The obtained results are discussed against the background of hypotheses stipulating pervasive lateral gene transfer in prokaryotes.

Visualizing and assessing phylogenetic congruence of core gene sets: a case study of the gamma-proteobacteria

Here, we address a much-debated topic: is there or is there not an organismal tree of gamma-proteobacteria that can be unambiguously inferred from a core of shared genes? We apply several recently developed analytical methods to this problem, for the first time. Our heat map analyses of P values and of bootstrap bipartitions show the presence of conflicting phylogenetic signals among these core genes. Our synthesis reconstruction suggests that at least 10% of these genes have been laterally transferred during the divergence of the gamma-proteobacteria, and that for most of the rest, there is too little phylogenetic signal to permit firm conclusions about the mode of inheritance. Although there is clearly a central tendency in this data set (it is far from random), lateral gene transfers cannot be ruled out. Instead of an organismal tree, we propose that these core genes could be used to define a more subtle and partially reticulated pattern of relationships.

Aerobic anoxygenic photosynthesis genes and operons in uncultured bacteria in the Delaware River

Photosynthesis genes and operons of aerobic anoxygenic photosynthetic (AAP) bacteria have been examined in a variety of marine habitats, but genomic information about freshwater AAP bacteria is lacking. The goal of this study was to examine photosynthesis genes of AAP bacteria in the Delaware River. In a fosmid library, we found two clones bearing photosynthesis gene clusters with unique gene content and organization. Both clones contained 37 open reading frames, with most of those genes encoding known AAP bacterial proteins. The genes in one fosmid were most closely related to those of AAP bacteria in the Rhodobacter genus. The genes of the other clone were related to those of freshwater beta-proteobacteria. Both clones contained the acsF gene, which is required for aerobic bacteriochlorophyll synthesis, suggesting that these bacteria are not anaerobes. The beta-proteobacterial fosmid has the puf operon B-A-L-M-C and is the first example of an uncultured bacterium with this operon structure. The alpha-3-proteobacterial fosmid has a rare gene order (Q-B-A-L-M-X), previously observed only in the Rhodobacter genus. Phylogenetic analyses of photosynthesis genes revealed a possible freshwater cluster of AAP beta-proteobacteria. The data from both Delaware River clones suggest there are groups of freshwater or estuarine AAP bacteria distinct from those found in marine environments.

Photosynthesis: what color was its origin?

Recent studies using geological and molecular phylogenetic evidence suggest several alternative evolutionary scenarios for the origin of photosynthesis. The earliest photosynthetic group is variously thought to be heliobacteria, proteobacteria or a precursor of cyanobacteria, organisms whose photosynthetic pigments make them different colors.

Evidence of a large novel gene pool associated with prokaryotic genomic islands

Microbial genes that are “novel” (no detectable homologs in other species) have become of increasing interest as environmental sampling suggests that there are many more such novel genes in yet-to-be-cultured microorganisms. By analyzing known microbial genomic islands and prophages, we developed criteria for systematic identification of putative genomic islands (clusters of genes of probable horizontal origin in a prokaryotic genome) in 63 prokaryotic genomes, and then characterized the distribution of novel genes and other features. All but a few of the genomes examined contained significantly higher proportions of novel genes in their predicted genomic islands compared with the rest of their genome (Paired t test = 4.43E-14 to 1.27E-18, depending on method). Moreover, the reverse observation (i.e., higher proportions of novel genes outside of islands) never reached statistical significance in any organism examined. We show that this higher proportion of novel genes in predicted genomic islands is not due to less accurate gene prediction in genomic island regions, but likely reflects a genuine increase in novel genes in these regions for both bacteria and archaea. This represents the first comprehensive analysis of novel genes in prokaryotic genomic islands and provides clues regarding the origin of novel genes. Our collective results imply that there are different gene pools associated with recently horizontally transmitted genomic regions versus regions that are primarily vertically inherited. Moreover, there are more novel genes within the gene pool associated with genomic islands. Since genomic islands are frequently associated with a particular microbial adaptation, such as antibiotic resistance, pathogen virulence, or metal resistance, this suggests that microbes may have access to a larger “arsenal” of novel genes for adaptation than previously thought.

More than 250 microbial genomes have been sequenced to date. A significant proportion of the genes in these genomes have no apparent similarity to known genes and their functions are unknown (i.e., they appear to be novel). As the number of sequenced genomes increases, the number of these novel genes continues to increase. In this paper, the authors now show, through an analysis of a diverse range of prokaryotic genomes, that novel genes are more prevalent in regions called genomic islands. Genomic islands are clusters of genes in genomes that show evidence of horizontal origins. This study is notable since genomic islands disproportionately contain many genes of medical, agricultural, and environmental importance (e.g., animal and plant pathogen virulence factors, antibiotic resistance genes, phenolic degradation genes, etc.). The observation that high proportions of novel genes are also localized to genomic islands suggests that microbes may have access to a larger “arsenal” of novel genes for important adaptations than previously thought. These results also imply that there are different gene pools associated with recently horizontally transmitted genomic regions versus regions that are primarily vertically inherited. The authors suggest that further studies involving large-scale environmental genomic sampling are required to help characterize this understudied gene pool.

The structure of photosystem I and evolution of photosynthesis

Oxygenic photosynthesis is the principal producer of both oxygen and organic matter on earth. The primary step in this process--the conversion of sunlight into chemical energy--is driven by four multi-subunit membrane protein complexes named photosystem I, photosystem II, cytochrome b(6)f complex and F-ATPase. Photosystem I generates the most negative redox potential in nature and thus largely determines the global amount of enthalpy in living systems. The recent structural determination of PSI complexes from cyanobacteria and plants sheds light on the evolutionary forces that shaped oxygenic photosynthesis. The fortuitous formation of our solar system in a space plentiful of elements, our distance from the sun and the long time of uninterrupted evolution enabled the perfection of photosynthesis and the evolution of advanced organisms. The available structural information complements the knowledge gained from genomic and proteomic data to illustrate a more precise scenario for the evolution of life systems on earth.

An evolution-based analysis scheme to identify CO2/O2 specificity-determining factors for ribulose 1,5-bisphosphate carboxylase/oxygenase

Ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCo) catalyzes a rate-limiting step in photosynthetic carbon assimilation (reacting with CO2) and its competitive photo-respiratory carbon oxidation (reacting with O2). RuBisCo enzyme with an enhanced CO2/O2 specificity would boost the ability to make great progress in agricultural production and environmental management. RuBisCos in marine non-green algae, resulting from an earlier endo-symbiotic event, diverge greatly from those in green plants and cyanobacteria and, further, have the highest CO2/O2 specificity whereas RuBisCos in cyanobacteria have the lowest. We assumed that there exist different levels of CO2/O2 specificity-determining factors, corresponding to different evolutionary events and specificity levels. Based on this assumption, we devised a scheme to identify these substrate-determining factors. From this analysis, we are able to discover different categories of the CO2/O2 specificity-determining factors that show which residue substitutions account for (relatively) small specificity changes, as happened in green plants, or a tremendous enhancement, as observed in marine non-green algae. Therefore, the analysis can improve our understanding of molecular mechanisms in the substrate specificity development and prioritize candidate specificity-determining surface residues for site-directed mutagenesis.

Highways of gene sharing in prokaryotes

The extent to which lateral genetic transfer has shaped microbial genomes has major implications for the emergence of community structures. We have performed a rigorous phylogenetic analysis of >220,000 proteins from genomes of 144 prokaryotes to determine the contribution of gene sharing to current prokaryotic diversity, and to identify "highways" of sharing between lineages. The inferred relationships suggest a pattern of inheritance that is largely vertical, but with notable exceptions among closely related taxa, and among distantly related organisms that live in similar environments.

Horizontal gene transfer, genome innovation and evolution

To what extent is the tree of life the best representation of the evolutionary history of microorganisms? Recent work has shown that, among sets of prokaryotic genomes in which most homologous genes show extremely low sequence divergence, gene content can vary enormously, implying that those genes that are variably present or absent are frequently horizontally transferred. Traditionally, successful horizontal gene transfer was assumed to provide a selective advantage to either the host or the gene itself, but could horizontally transferred genes be neutral or nearly neutral? We suggest that for many prokaryotes, the boundaries between species are fuzzy, and therefore the principles of population genetics must be broadened so that they can be applied to higher taxonomic categories.

Evolutionary aspects of whole-genome biology

A decade of access to whole-genome sequences has been increasingly revealing about the informational network relating all living organisms. Although at one point there was concern that extensive horizontal gene transfer might hopelessly muddle phylogenies, it has not proved a severe hindrance. The melding of sequence and structural information is being used to great advantage, and the prospect exists that some of the earliest aspects of life on Earth can be reconstructed, including the invention of biosynthetic and metabolic pathways. Still, some fundamental phylogenetic problems remain, including determining the root--if there is one--of the historical relationship between Archaea, Bacteria and Eukarya.

The cobweb of life revealed by genome-scale estimates of horizontal gene transfer

With the availability of increasing amounts of genomic sequences, it is becoming clear that genomes experience horizontal transfer and incorporation of genetic information. However, to what extent such horizontal gene transfer (HGT) affects the core genealogical history of organisms remains controversial. Based on initial analyses of complete genomic sequences, HGT has been suggested to be so widespread that it might be the “essence of phylogeny” and might leave the treelike form of genealogy in doubt. On the other hand, possible biased estimation of HGT extent and the findings of coherent phylogenetic patterns indicate that phylogeny of life is well represented by tree graphs. Here, we reexamine this question by assessing the extent of HGT among core orthologous genes using a novel statistical method based on statistical comparisons of tree topology. We apply the method to 40 microbial genomes in the Clusters of Orthologous Groups database over a curated set of 297 orthologous gene clusters, and we detect significant HGT events in 33 out of 297 clusters over a wide range of functional categories. Estimates of positions of HGT events suggest a low mean genome-specific rate of HGT (2.0%) among the orthologous genes, which is in general agreement with other quantitative of HGT. We propose that HGT events, even when relatively common, still leave the treelike history of phylogenies intact, much like cobwebs hanging from tree branches.

A stastical approach applied to 297 orthologous gene clusters in 40 microbial genomes suggests a low rate of interspecies gene transfer. Species relationships can therefore be modeled with a tree structure.

The puf operon of the purple sulfur bacterium Amoebobacter purpureus: structure, transcription and phylogenetic analysis

The puf operon, encoding photosynthetic reaction center and light-harvesting genes, of the purple sulfur phototrophic bacterium Amoebobacter purpureus was cloned and sequenced. This revealed an unusual operon structure of the genes pufB1 A1 LMCB2 A2 B3 A3. The sequence represents the second complete puf operon available for Chromatiaceae. So far, additional sets of light-harvesting 1 (LH1) genes, pufB2 A2 and pufB3 A3 in the region downstream of pufC have only been described for Allochromatium vinosum. Along with reports of multiple LH1 polypeptides found in some Ectothiorhodospiraceae by direct protein sequencing, our results indicate that multiple LH1 genes may occur frequently in phototrophic gamma-proteobacteria. Phylogenetic analyses suggested a coevolution of the core puf genes pufB1 A1 LM. Separate analysis of the LH1 alpha and beta polypeptides revealed a high intraspecies relatedness for the secondary LH1beta polypeptides, possibly caused by functional constraints. In contrast, LH1alpha subunits of Amb. purpureus and Alc. vinosum are closely related (85% sequence identity) which could reflect horizontal gene transfer. RNA analyses suggested co-transcription of all puf genes in Amb. purpureus as a 5.5 kb primary transcript which appears to be more stable than the puf operon primary transcripts of purple non-sulfur bacteria. The 5' end of the transcript mapped to a putative promoter, which contains a -35 region located in an inverted repeat DNA sequence.

PentaPlot: a software tool for the illustration of genome mosaicism

Background

Dekapentagonal maps depict the phylogenetic relationships of five genomes in a visually appealing diagram and can be viewed as an alternative to a single evolutionary consensus tree. In particular, the generated maps focus attention on those gene families that significantly deviate from the consensus or plurality phylogeny. PentaPlot is a software tool that computes such dekapentagonal maps given an appropriate probability support matrix.

Results

The visualization with dekapentagonal maps critically depends on the optimal layout of unrooted tree topologies representing different evolutionary relationships among five organisms along the vertices of the dekapentagon. This is a difficult optimization problem given the large number of possible layouts. At its core our tool utilizes a genetic algorithm with demes and a local search strategy to search for the optimal layout. The hybrid genetic algorithm performs satisfactorily even in those cases where the chosen genomes are so divergent that little phylogenetic information has survived in the individual gene families.

Conclusion

PentaPlot is being made publicly available as an open source project at .

Displacement of ɛ-proteobacterial core genes by horizontally transferred homologous genes

The introduction of novel genes by horizontal gene transfer (HGT) is considered an alternative mechanism for genetic adaptation, leading to diversification and speciation. The goal of this study was to determine which genes that are present in all sequenced epsilon-proteobacterial genomes were acquired by HGT. In our approach we used BLAST analysis to reduce the number of genes that subsequently needed to be analysed using more in-depth phylogenetic methods, including neighbour-joining and maximum likelihood. Among the 991 core genes found in all five completed epsilon-proteobacterial genome sequences, we identified 30 genes that were probably acquired by HGT. It is proposed that these genes displaced an ancestral core gene with a similar function. Although it was not possible to identify putative donor taxa for all acquired genes, it was clear that genes were acquired from a wide range of Bacteria, including Spirochaetes, Firmicutes, Actinobacteria, mycoplasmas and several subdivisions of the Proteobacteria. We did not observe HGT from Archaea to the epsilon-Proteobacteria. The majority of acquired genes were operational genes involved in transport, metabolism, signal transduction and energy production and conversion.

Automatic selection of representative proteins for bacterial phylogeny

Background

Although there are now about 200 complete bacterial genomes in GenBank, deep bacterial phylogeny remains a difficult problem, due to confounding horizontal gene transfers and other phylogenetic "noise". Previous methods have relied primarily upon biological intuition or manual curation for choosing genomic sequences unlikely to be horizontally transferred, and have given inconsistent phylogenies with poor bootstrap confidence.

Results

We describe an algorithm that automatically picks "representative" protein families from entire genomes for use as phylogenetic characters. A representative protein family is one that, taken alone, gives an organismal distance matrix in good agreement with a distance matrix computed from all sufficiently conserved proteins. We then use maximum-likelihood methods to compute phylogenetic trees from a concatenation of representative sequences. We validate the use of representative proteins on a number of small phylogenetic questions with accepted answers. We then use our methodology to compute a robust and well-resolved phylogenetic tree for a diverse set of sequenced bacteria. The tree agrees closely with a recently published tree computed using manually curated proteins, and supports two proposed high-level clades: one containing Actinobacteria, Deinococcus, and Cyanobacteria ("Terrabacteria"), and another containing Planctomycetes and Chlamydiales.

Conclusion

Representative proteins provide an effective solution to the problem of selecting phylogenetic characters.

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.

Phylogenetic analyses of the core antenna domain: investigating the origin of photosystem I

Phototrophy, the conversion of light to biochemical energy, occurs throughout the Bacteria and plants, however, debate continues over how different phototrophic mechanisms and the bacteria that contain them are related. There are two types of phototrophic mechanisms in the Bacteria: reaction center type 1 (RC1) has core and core antenna domains that are parts of a single polypeptide, whereas reaction center type 2 (RC2) is composed of short core proteins without antenna domains. In cyanobacteria, RC2 is associated with separate core antenna proteins that are homologous to the core antenna domains of RC1. We reconstructed evolutionary relationships among phototrophic mechanisms based on a phylogeny of core antenna domains/proteins. Core antenna domains of 46 polypeptides were aligned, including the RC1 core proteins of heliobacteria, green sulfur bacteria, and photosystem I (PSI) of cyanobacteria and plastids, plus core antenna proteins of photosystem II (PSII) from cyanobacteria and plastids. Maximum likelihood, parsimony, and neighbor joining methods all supported a single phylogeny in which PSII core antenna proteins (PsbC, PsbB) arose within the cyanobacteria from duplications of the RC1-associated core antenna domains and accessory antenna proteins (IsiA, PcbA, PcbC) arose from duplications of PsbB. The data indicate an evolutionary history of RC1 in which an initially homodimeric reaction center was vertically transmitted to green sulfur bacteria, heliobacteria, and an ancestor of cyanobacteria. A heterodimeric RC1 (=PSI) then arose within the cyanobacterial lineage. In this scenario, the current diversity of core antenna domains/proteins is explained without a need to invoke horizontal transfer.

Evolutionary sequence analysis of complete eukaryote genomes

Background

Gene duplication and gene loss during the evolution of eukaryotes have hindered attempts to estimate phylogenies and divergence times of species. Although current methods that identify clusters of orthologous genes in complete genomes have helped to investigate gene function and gene content, they have not been optimized for evolutionary sequence analyses requiring strict orthology and complete gene matrices. Here we adopt a relatively simple and fast genome comparison approach designed to assemble orthologs for evolutionary analysis. Our approach identifies single-copy genes representing only species divergences (panorthologs) in order to minimize potential errors caused by gene duplication. We apply this approach to complete sets of proteins from published eukaryote genomes specifically for phylogeny and time estimation.

Results

Despite the conservative criterion used, 753 panorthologs (proteins) were identified for evolutionary analysis with four genomes, resulting in a single alignment of 287,000 amino acids. With this data set, we estimate that the divergence between deuterostomes and arthropods took place in the Precambrian, approximately 400 million years before the first appearance of animals in the fossil record. Additional analyses were performed with seven, 12, and 15 eukaryote genomes resulting in similar divergence time estimates and phylogenies.

Conclusion

Our results with available eukaryote genomes agree with previous results using conventional methods of sequence data assembly from genomes. They show that large sequence data sets can be generated relatively quickly and efficiently for evolutionary analyses of complete genomes.

Representative freshwater bacterioplankton isolated from Crater Lake, Oregon

High-throughput culturing (HTC) methods that rely on dilution to extinction in very-low-nutrient media were used to obtain bacterial isolates from Crater Lake, Oregon. 16S rRNA sequence determination and phylogenetic reconstruction were used to determine the potential ecological significance of isolated bacteria, both in Crater Lake and globally. Fifty-five Crater Lake isolates yielded 16 different 16S rRNA gene sequences. Thirty of 55 (55%) Crater Lake isolates had 16S rRNA gene sequences with 97% or greater similarity to sequences recovered previously from Crater Lake 16S rRNA gene clone libraries. Furthermore, 36 of 55 (65%) Crater Lake isolates were found to be members of widely distributed freshwater groups. These results confirm that HTC is a significant improvement over traditional isolation techniques that tend to enrich for microorganisms that do not predominate in their environment and rarely correlate with 16S rRNA gene clone library sequences. Although all isolates were obtained under dark, heterotrophic growth conditions, 2 of the 16 different groups showed evidence of photosynthetic capability as assessed by the presence of puf operon sequences, suggesting that photoheterotrophy may be a significant process in this oligotrophic, freshwater habitat.

Molecular basis of bacterial resistance to chloramphenicol and florfenicol

Chloramphenicol (Cm) and its fluorinated derivative florfenicol (Ff) represent highly potent inhibitors of bacterial protein biosynthesis. As a consequence of the use of Cm in human and veterinary medicine, bacterial pathogens of various species and genera have developed and/or acquired Cm resistance. Ff is solely used in veterinary medicine and has been introduced into clinical use in the mid-1990s. Of the Cm resistance genes known to date, only a small number also mediates resistance to Ff. In this review, we present an overview of the different mechanisms responsible for resistance to Cm and Ff with particular focus on the two different types of chloramphenicol acetyltransferases (CATs), specific exporters and multidrug transporters. Phylogenetic trees of the different CAT proteins and exporter proteins were constructed on the basis of a multisequence alignment. Moreover, information is provided on the mobile genetic elements carrying Cm or Cm/Ff resistance genes to provide a basis for the understanding of the distribution and the spread of Cm resistance--even in the absence of a selective pressure imposed by the use of Cm or Ff.

The complex architecture of oxygenic photosynthesis

Oxygenic photosynthesis is the principal producer of both oxygen and organic matter on earth. The primary step in this process - the conversion of sunlight into chemical energy - is driven by four, multisubunit, membrane-protein complexes that are known as photosystem I, photosystem II, cytochrome b(6)f and F-ATPase. Structural insights into these complexes are now providing a framework for the exploration not only of energy and electron transfer, but also of the evolutionary forces that shaped the photosynthetic apparatus.

Physcomitrella patens and Ceratodon purpureus, mosses as model organisms in photosynthesis studies

With the discovery of targeted gene replacement, moss biology has been rapidly advancing over the last 10 years. This study demonstrates the usefulness of moss as a model organism for plant photosynthesis research. The two mosses examined in this study, Physcomitrella patens and Ceratodon purpureus, are easily cultured through vegetative propagation. Growth tests were conducted to determine carbon sources suitable for maintaining heterotrophic growth while photosynthesis was blocked. Photosynthetic parameters examined in these plants indicated that the photosynthetic activity of Ceratodon and Physcomitrella is more similar to vascular plants than cyanobacteria or green algae. Ceratodon plants grown heterotrophically appeared etiolated in that the plants were taller and plastids did not differentiate thylakoid membranes. After returning to the light, the plants developed green, photosynthetically active chloroplasts. Furthermore, UV-induced mutagenesis was used to show that photosynthesis-deficient mutant Ceratodon plants could be obtained. After screening approximately 1000 plants, we obtained a number of mutants, which could be arranged into the following categories: high fluorescence, low fluorescence, fast and slow fluorescence quenching, and fast and slow greening. Our results indicate that in vivo biophysical analysis of photosynthetic activity in the mosses can be carried out which makes both mosses useful for photosynthesis studies, and Ceratodon best sustains perturbations in photosynthetic activity.

A genomic timescale of prokaryote evolution: insights into the origin of methanogenesis, phototrophy, and the colonization of land

Background

The timescale of prokaryote evolution has been difficult to reconstruct because of a limited fossil record and complexities associated with molecular clocks and deep divergences. However, the relatively large number of genome sequences currently available has provided a better opportunity to control for potential biases such as horizontal gene transfer and rate differences among lineages. We assembled a data set of sequences from 32 proteins (~7600 amino acids) common to 72 species and estimated phylogenetic relationships and divergence times with a local clock method.

Results

Our phylogenetic results support most of the currently recognized higher-level groupings of prokaryotes. Of particular interest is a well-supported group of three major lineages of eubacteria (Actinobacteria, Deinococcus, and Cyanobacteria) that we call Terrabacteria and associate with an early colonization of land. Divergence time estimates for the major groups of eubacteria are between 2.5–3.2 billion years ago (Ga) while those for archaebacteria are mostly between 3.1–4.1 Ga. The time estimates suggest a Hadean origin of life (prior to 4.1 Ga), an early origin of methanogenesis (3.8–4.1 Ga), an origin of anaerobic methanotrophy after 3.1 Ga, an origin of phototrophy prior to 3.2 Ga, an early colonization of land 2.8–3.1 Ga, and an origin of aerobic methanotrophy 2.5–2.8 Ga.

Conclusions

Our early time estimates for methanogenesis support the consideration of methane, in addition to carbon dioxide, as a greenhouse gas responsible for the early warming of the Earths' surface. Our divergence times for the origin of anaerobic methanotrophy are compatible with highly depleted carbon isotopic values found in rocks dated 2.8–2.6 Ga. An early origin of phototrophy is consistent with the earliest bacterial mats and structures identified as stromatolites, but a 2.6 Ga origin of cyanobacteria suggests that those Archean structures, if biologically produced, were made by anoxygenic photosynthesizers. The resistance to desiccation of Terrabacteria and their elaboration of photoprotective compounds suggests that the common ancestor of this group inhabited land. If true, then oxygenic photosynthesis may owe its origin to terrestrial adaptations.

'Conserved hypothetical' proteins: prioritization of targets for experimental study

Comparative genomics shows that a substantial fraction of the genes in sequenced genomes encodes 'conserved hypothetical' proteins, i.e. those that are found in organisms from several phylogenetic lineages but have not been functionally characterized. Here, we briefly discuss recent progress in functional characterization of prokaryotic 'conserved hypothetical' proteins and the possible criteria for prioritizing targets for experimental study. Based on these criteria, the chief one being wide phyletic spread, we offer two 'top 10' lists of highly attractive targets. The first list consists of proteins for which biochemical activity could be predicted with reasonable confidence but the biological function was predicted only in general terms, if at all ('known unknowns'). The second list includes proteins for which there is no prediction of biochemical activity, even if, for some, general biological clues exist ('unknown unknowns'). The experimental characterization of these and other 'conserved hypothetical' proteins is expected to reveal new, crucial aspects of microbial biology and could also lead to better functional prediction for medically relevant human homologs.

Seeing green bacteria in a new light: genomics-enabled studies of the photosynthetic apparatus in green sulfur bacteria and filamentous anoxygenic phototrophic bacteria

Based upon their photosynthetic nature and the presence of a unique light-harvesting antenna structure, the chlorosome, the photosynthetic green bacteria are defined as a distinctive group in the Bacteria. However, members of the two taxa that comprise this group, the green sulfur bacteria (Chlorobi) and the filamentous anoxygenic phototrophic bacteria ("Chloroflexales"), are otherwise quite different, both physiologically and phylogenetically. This review summarizes how genome sequence information facilitated studies of the biosynthesis and function of the photosynthetic apparatus and the oxidation of inorganic sulfur compounds in two model organisms that represent these taxa, Chlorobium tepidum and Chloroflexus aurantiacus. The genes involved in bacteriochlorophyll (BChl) c and carotenoid biosynthesis in these two organisms were identified by sequence homology with known BChl a and carotenoid biosynthesis enzymes, gene cluster analysis in Cfx. aurantiacus, and gene inactivation studies in Chl. tepidum. Based on these results, BChl a and BChl c biosynthesis is similar in the two organisms, whereas carotenoid biosynthesis differs significantly. In agreement with its facultative anaerobic nature, Cfx. aurantiacus in some cases apparently produces structurally different enzymes for heme and BChl biosynthesis, in which one enzyme functions under anoxic conditions and the other performs the same reaction under oxic conditions. The Chl. tepidum mutants produced with modified BChl c and carotenoid species also allow the functions of these pigments to be studied in vivo.

Homologs of plant PsbP and PsbQ proteins are necessary for regulation of photosystem ii activity in the cyanobacterium Synechocystis 6803

The mechanism of oxygen evolution by photosystem II (PSII) has remained highly conserved during the course of evolution from ancestral cyanobacteria to green plants. A cluster of manganese, calcium, and chloride ions, whose binding environment is optimized by PSII extrinsic proteins, catalyzes this water-splitting reaction. The accepted view is that in plants and green algae, the three extrinsic proteins are PsbO, PsbP, and PsbQ, whereas in cyanobacteria, they are PsbO, PsbV, and PsbU. Our previous proteomic analysis established the presence of a PsbQ homolog in the cyanobacterium Synechocystis 6803. The current study additionally demonstrates the presence of a PsbP homolog in cyanobacterial PSII. Both psbP and psbQ inactivation mutants exhibited reduced photoautotrophic growth as well as decreased water oxidation activity under CaCl(2)-depleted conditions. Moreover, purified PSII complexes from each mutant had significantly reduced activity. In cyanobacteria, one PsbQ is present per PSII complex, whereas PsbP is significantly substoichiometric. These findings indicate that both PsbP and PsbQ proteins are regulators that are necessary for the biogenesis of optimally active PSII in Synechocystis 6803. The new picture emerging from these data is that five extrinsic PSII proteins, PsbO, PsbP, PsbQ, PsbU, and PsbV, are present in cyanobacteria, two of which, PsbU and PsbV, have been lost during the evolution of green plants.

Speciation in Chlamydia: genomewide phylogenetic analyses identified a reliable set of acquired genes

Horizontal gene transfer (HGT), a process through which genomes acquire sequences from distantly related organisms, is believed to be a major source of genetic diversity in bacteria. A central question concerning the impact of HGT on bacterial genome evolution is the proportion of horizontally transferred sequences within genomes. This issue, however, remains unresolved because the various methods developed to detect potential HGT events identify different sets of genes. The present-day consensus is that phylogenetic analysis of individual genes is still the most objective and accurate approach for determining the occurrence and directionality of HGT. Here we present a genome-scale phylogenetic analysis of protein-encoding genes from five closely related Chlamydia, identifying a reliable set of sequences that have arisen via HGT since the divergence of the Chlamydia lineage. According to our knowledge, this is the first systematic phylogenetic inference-based attempt to establish a reliable set of acquired genes in a bacterial genome. Although Chlamydia are obligate intracellular parasites of higher eukaryotes, and thus suspected to be isolated from HGT more than the free-living species, our results show that their diversification has involved the introduction of foreign sequences into their genome. Furthermore, we also identified a complete set of genes that have undergone deletion, duplication, or rearrangement during this evolutionary period leading to the radiation of Chlamydia species. Our analysis may provide a deeper insight into how these medically important pathogens emerged and evolved from a common ancestor.

The Evolution of Modern Eukaryotic Phytoplankton

The community structure and ecological function of contemporary marine ecosystems are critically dependent on eukaryotic phytoplankton. Although numerically inferior to cyanobacteria, these organisms are responsible for the majority of the flux of organic matter to higher trophic levels and the ocean interior. Photosynthetic eukaryotes evolved more than 1.5 billion years ago in the Proterozoic oceans. However, it was not until the Mesozoic Era (251 to 65 million years ago) that the three principal phytoplankton clades that would come to dominate the modern seas rose to ecological prominence. In contrast to their pioneering predecessors, the dinoflagellates, coccolithophores, and diatoms all contain plastids derived from an ancestral red alga by secondary symbiosis. Here we examine the geological, geochemical, and biological processes that contributed to the rise of these three, distantly related, phytoplankton groups.

The evolutionary development of the protein complement of photosystem 2

During the transition from anoxygenic to oxygenic photosynthesis, the Type 2 reaction center underwent many changes, none so dramatic as the remarkable increase in complexity at the protein level, from only three or four subunits in the anoxygenic reaction center to possibly more than 25 in Photosystem 2 (PS2). The evolutionary source of most of these proteins is enigmatic, as they have no apparent homology to any other proteins in existing databases. However, some of the proteins in PS2 have apparent homologies to each other, suggesting ancient gene duplications have played an important role in the development of the complex. These homologies include the well-known examples of the D1 and D2 reaction center core proteins and the CP43 and CP47 core antenna proteins. In addition, PsbE and PsbF, the two subunits comprising cytochrome b-559, show homology to each other, suggesting that a homodimeric cytochrome preceded the heterodimeric one. Other potential homologies that appear to be statistically significant include PsbV with the N-terminal part of D1 and PsbT with PsbI. Most of the proteins that make up the photosynthetic apparatus bear no relation to any other proteins from any source. This suggests that a period of remarkable evolutionary innovation took place when the ability to make oxygen was invented. This was probably a response to the production of highly toxic oxygen and these new proteins served to protect and repair the photosynthetic apparatus from the harmful effects of oxygen.

Genome organization and localization of the pufLM genes of the photosynthesis reaction center in phylogenetically diverse marine Alphaproteobacteria

Genome organization, plasmid content and localization of the pufLM genes of the photosynthesis reaction center were studied by pulsed-field gel electrophoresis (PFGE) in marine phototrophic Alphaproteobacteria. Both anaerobic phototrophs (Rhodobacter veldkampii and Rhodobacter sphaeroides) and strictly aerobic anoxygenic phototrophs from the Roseobacter-Sulfitobacter-Silicibacter clade (Roseivivax halodurans, Roseobacter litoralis, Staleya guttiformis, Roseovarius tolerans, and five new strains isolated from dinoflagellate cultures) were investigated. The complete genome size was estimated for R. litoralis DSM6996(T) to be 4,704 kb, including three linear plasmids. All strains contained extrachromosomal elements of various conformations (linear or circular) and lengths (between 4.35 and 368 kb). In strain DFL-12, a member of a putative new genus isolated from a culture of the toxic dinoflagellate Prorocentrum lima, seven linear plasmids were found, together comprising 860 kb of genetic information. Hybridization with probes against the pufLM genes of the photosynthesis gene cluster after Southern transfer of the genomic DNAs showed these genes to be located on a linear plasmid of 91 kb in R. litoralis and on a linear plasmid of 120 kb in S. guttiformis, theoretically allowing their horizontal transfer. In all other strains, the pufLM genes were detected on the bacterial chromosome. The large number and significant size of the linear plasmids found especially in isolates from dinoflagellates might account for the metabolic versatility and presumed symbiotic association with eukaryotic hosts in these bacteria.

GeneOrder3.0: software for comparing the order of genes in pairs of small bacterial genomes

BACKGROUND

An increasing number of whole viral and bacterial genomes are being sequenced and deposited in public databases. In parallel to the mounting interest in whole genomes, the number of whole genome analyses software tools is also increasing. GeneOrder was originally developed to provide an analysis of genes between two genomes, allowing visualization of gene order and synteny comparisons of any small genomes. It was originally developed for comparing virus, mitochondrion and chloroplast genomes. This is now extended to small bacterial genomes of sizes less than 2 Mb.

RESULTS

GeneOrder3.0 has been developed and validated successfully on several small bacterial genomes (ca. 580 kb to 1.83 Mb) archived in the NCBI GenBank database. It is an updated web-based "on-the-fly" computational tool allowing gene order and synteny comparisons of any two small bacterial genomes. Analyses of several bacterial genomes show that a large amount of gene and genome re-arrangement occurs, as seen with earlier DNA software tools. This can be displayed at the protein level using GeneOrder3.0. Whole genome alignments of genes are presented in both a table and a dot plot. This allows the detection of evolutionary more distant relationships since protein sequences are more conserved than DNA sequences.

CONCLUSIONS

GeneOrder3.0 allows researchers to perform comparative analysis of gene order and synteny in genomes of sizes up to 2 Mb "on-the-fly."

AVAILABILITY

http://binf.gmu.edu/genometools.html and http://pasteur.atcc.org:8050/GeneOrder3.0.

THE RATE AND PATTERN OF CLADOGENESIS IN MICROBES

Theories of macroevolution rarely have been extended to include microbes; however, because microbes represent the most ancient and diverse assemblage of organismal diversity, such oversight limits our understanding of evolutionary history. Our analysis of phylogenetic trees for microbes suggests that macroevolution may differ between prokaryotes and both micro- and macroeukaryotes (mainly plants and animals). Phylogenetic trees inferred for prokaryotes and some microbial eukaryotes conformed to expectations assuming a constant rate of cladogenesis over time and among lineages: nevertheless, microbial eukaryote trees exhibited more variation in rates of cladogenesis than prokaryote trees. We hypothesize that the contrast of macroevolutionary dynamics between prokaryotes and many eukaryotes is due, at least in part, to differences in the prevalence of lateral gene transfer (LGT) between the two groups. Inheritance is predominantly, if not wholly, vertical within eukaryotes, a feature that allows for the emergence and maintenance of heritable variation among lineages. By contrast, frequent LGT in prokaryotes may ameliorate heritable variation in rate of cladogenesis resulting from the emergence of key innovations; thus, the inferred difference in macroevolution might reflect exclusivity of key innovations in eukaryotes and their promiscuous nature in prokaryotes.

A hybrid clustering approach to recognition of protein families in 114 microbial genomes

Background

Grouping proteins into sequence-based clusters is a fundamental step in many bioinformatic analyses (e.g., homology-based prediction of structure or function). Standard clustering methods such as single-linkage clustering capture a history of cluster topologies as a function of threshold, but in practice their usefulness is limited because unrelated sequences join clusters before biologically meaningful families are fully constituted, e.g. as the result of matches to so-called promiscuous domains. Use of the Markov Cluster algorithm avoids this non-specificity, but does not preserve topological or threshold information about protein families.

Results

We describe a hybrid approach to sequence-based clustering of proteins that combines the advantages of standard and Markov clustering. We have implemented this hybrid approach over a relational database environment, and describe its application to clustering a large subset of PDB, and to 328577 proteins from 114 fully sequenced microbial genomes. To demonstrate utility with difficult problems, we show that hybrid clustering allows us to constitute the paralogous family of ATP synthase F1 rotary motor subunits into a single, biologically interpretable hierarchical grouping that was not accessible using either single-linkage or Markov clustering alone. We describe validation of this method by hybrid clustering of PDB and mapping SCOP families and domains onto the resulting clusters.

Conclusion

Hybrid (Markov followed by single-linkage) clustering combines the advantages of the Markov Cluster algorithm (avoidance of non-specific clusters resulting from matches to promiscuous domains) and single-linkage clustering (preservation of topological information as a function of threshold). Within the individual Markov clusters, single-linkage clustering is a more-precise instrument, discerning sub-clusters of biological relevance. Our hybrid approach thus provides a computationally efficient approach to the automated recognition of protein families for phylogenomic analysis.

Visualization of the phylogenetic content of five genomes using dekapentagonal maps

Dekapentagonal maps depict phylogenetic information for orthologous genes present in five genomes, and provide a pre-screen for putatively horizontally transferred genes.

The methods presented here summarize phylogenetic relationships of genomes in visually appealing and informative figures. Dekapentagonal maps depict phylogenetic information for orthologous genes present in five genomes, and provide a pre-screen for putatively horizontally transferred genes. If the majority of individual gene phylogenies are unresolved, bipartition histograms provide a means of uncovering and analyzing the plurality consensus. Analyses of genomes representing five photosynthetic bacterial phyla and of the prokaryotic contributions to the eukaryotic cell illustrate the utility of the methods.

Characterization of the glucose-6-phosphate isomerase gene in Phytophthora infestans reveals the presence of multiple alleles

Glucose-6-phosphate isomerase (GPI) plays a key role in both glycolysis and gluconeogenesis. Isoforms of GPI are common, and therefore, its isozyme pattern is widely used to characterize isolates of Phytophthora infestans. Despite the importance of GPI in P. infestans studies, the gene encoding this enzyme has not yet been characterized. Furthermore, it has been suggested that P. infestans contains multiple copies of the gene but this hypothesis remains to be demonstrated. We have cloned and characterized GPI in various isolates of P. infestans as well as in several species of the genus Phytophthora. The gene contains 1671bp and encodes a protein with a predicted molecular weight of 60.8kDa. Multiple different alleles were identified and Southern analysis indicated certain P. infestans isolates carry several copies of the gene. Phylogenetic analysis revealed that P. infestans GPI is most closely related to sequences from Toxoplasma gondii, Arabidopsis thaliana, and Clarkia lewisii.