Ancient horizontal gene transfer

Examining Ancient Inter-domain Horizontal Gene Transfer

Details of the genomic changes that occurred in the ancestors of Eukarya, Archaea and Bacteria are elusive. Ancient interdomain horizontal gene transfer (IDHGT) amongst the ancestors of these three domains has been difficult to detect and analyze because of the extreme degree of divergence of genes in these three domains and because most evidence for such events are poorly supported. In addition, many researchers have suggested that the prevalence of IDHGT events early in the evolution of life would most likely obscure the patterns of divergence of major groups of organisms let alone allow the tracking of horizontal transfer at this level. In order to approach this problem, we mined the E. coli genome for genes with distinct paralogs. Using the 1,268 E. coli K-12 genes with 40% or higher similarity level to a paralog elsewhere in the E. coli genome we detected 95 genes found exclusively in Bacteria and Archaea and 86 genes found in Bacteria and Eukarya. These genes form the basis for our analysis of IDHGT. We also applied a newly developed statistical test (the node height test), to examine the robustness of these inferences and to corroborate the phylogenetically identified cases of ancient IDHGT. Our results suggest that ancient inter domain HGT is restricted to special cases, mostly involving symbiosis in eukaryotes and specific adaptations in prokaryotes. Only three genes in the Bacteria + Eukarya class (Deoxyxylulose-5-phosphate synthase (DXPS), fructose 1,6-phosphate aldolase class II protein and glucosamine-6-phosphate deaminase) and three genes–in the Bacteria + Archaea class (ABC-type FE3+ -siderophore transport system, ferrous iron transport protein B, and dipeptide transport protein) showed evidence of ancient IDHGT. However, we conclude that robust estimates of IDHGT will be very difficult to obtain due to the methodological limitations and the extreme sequence saturation of the genes suspected of being involved in IDHGT.

The mosaic genome of Anaeromyxobacter dehalogenans strain 2CP-C suggests an aerobic common ancestor to the delta-proteobacteria

Anaeromyxobacter dehalogenans strain 2CP-C is a versaphilic delta-Proteobacterium distributed throughout many diverse soil and sediment environments. 16S rRNA gene phylogenetic analysis groups A. dehalogenans together with the myxobacteria, which have distinguishing characteristics including strictly aerobic metabolism, sporulation, fruiting body formation, and surface motility. Analysis of the 5.01 Mb strain 2CP-C genome substantiated that this organism is a myxobacterium but shares genotypic traits with the anaerobic majority of the delta-Proteobacteria (i.e., the Desulfuromonadales). Reflective of its respiratory versatility, strain 2CP-C possesses 68 genes coding for putative c-type cytochromes, including one gene with 40 heme binding motifs. Consistent with its relatedness to the myxobacteria, surface motility was observed in strain 2CP-C and multiple types of motility genes are present, including 28 genes for gliding, adventurous (A-) motility and 17 genes for type IV pilus-based motility (i.e., social (S-) motility) that all have homologs in Myxococcus xanthus. Although A. dehalogenans shares many metabolic traits with the anaerobic majority of the delta-Proteobacteria, strain 2CP-C grows under microaerophilic conditions and possesses detoxification systems for reactive oxygen species. Accordingly, two gene clusters coding for NADH dehydrogenase subunits and two cytochrome oxidase gene clusters in strain 2CP-C are similar to those in M. xanthus. Remarkably, strain 2CP-C possesses a third NADH dehydrogenase gene cluster and a cytochrome cbb₃ oxidase gene cluster, apparently acquired through ancient horizontal gene transfer from a strictly anaerobic green sulfur bacterium. The mosaic nature of the A. dehalogenans strain 2CP-C genome suggests that the metabolically versatile, anaerobic members of the delta-Proteobacteria may have descended from aerobic ancestors with complex lifestyles.

Whole genome-based assessment of the taxonomic position of the arthropod pathogenic bacterium Rickettsiella grylli

Rickettsiella grylli is an intracellular bacterial pathogen of aquatic and terrestrial arthropods. Previous determination of its 16S rRNA-encoding sequence has led to the taxonomic classification of the genus Rickettsiella in the class Gammaproteobacteria, order Legionellales, family Coxiellaceae, i.e. in close vicinity to vertebrate pathogenic bacteria of the genera Coxiella and Legionella. Here we use the additional information available from the recently published first whole genome sequence from this genus to evaluate critically the taxonomic classification of R. grylli beyond the 16S rRNA gene level. Using phylogenetic reconstruction, together with significance testing on a data basis defined by a core set of 211 previously identified families of protein-encoding genes, together with a reanalysis of 16S rRNA gene data, the present study firmly corroborates the assignment of this species to both the class Gammaproteobacteria and the order Legionellales. However, the results obtained from concatenated and single protein, single protein-encoding gene, and 16S rRNA gene data demonstrate a similar phylogenetic distance of R. grylli to both the Coxiellaceae and the Legionellaceae and are, therefore, inconsistent with its current family-level classification. Consequently, a respective reorganization of the order Legionellales is proposed.

Comparative genomics of human-like Schistosoma japonicum genes indicates a putative mechanism for host-parasite relationship

Schistosoma japonicum causes schistosomiasis in humans and livestock in the Asia-Pacific region. We assembled more than 43,700 S. japonicum expressed sequence tags and conducted comparative genomic analyses between S. japonicum and its human host. Some schistosome genes showed exceptionally high similarity in nucleotide sequence to their human homologues, of which five exhibited anomalous phylogeny and human codon usage bias. The most plausible explanation for their presence is horizontal gene transfer from host to parasite. Functional evidence suggests that S. japonicum might exploit host endocrine and immune signals for cell development and maturation via these host-like genes.

Evidence of multiple horizontal transfers of the long terminal repeat retrotransposon RIRE1 within the genus Oryza

Horizontal gene transfer, defined as the transmission of genetic material between reproductively isolated species, has been considered for a long time to be a rare phenomenon. Most well-documented cases of horizontal gene transfer have been described in prokaryotes or in animals and they often involve transposable elements. The most abundant class of transposable elements in plant genomes are the long terminal repeat (LTR) retrotransposons. Because of their propensity to increase their copy number while active, LTR retrotransposons can have a significant impact on genomics changes during evolution. In a previous study, we showed that in the wild rice species Oryza australiensis, 60% of the genome is composed of only three families of LTR retrotransposons named RIRE1, Wallabi and Kangourou. In the present study, using both in silico and experimental approaches, we show that one of these three families, RIRE1, has been transferred horizontally between O. australiensis and seven other reproductively isolated Oryza species. This constitutes a new case of horizontal transfer in plants.

Identification of multiple independent horizontal gene transfers into poxviruses using a comparative genomics approach

Background

Poxviruses are important pathogens of humans, livestock and wild animals. These large dsDNA viruses have a set of core orthologs whose gene order is extremely well conserved throughout poxvirus genera. They also contain many genes with sequence and functional similarity to host genes which were probably acquired by horizontal gene transfer.

Although phylogenetic trees can indicate the occurrence of horizontal gene transfer and even uncover multiple events, their use may be hampered by uncertainties in both the topology and the rooting of the tree. We propose to use synteny conservation around the horizontally transferred gene (HTgene) to distinguish between single and multiple events.

Results

Here we devise a method that incorporates comparative genomic information into the investigation of horizontal gene transfer, and we apply this method to poxvirus genomes. We examined the synteny conservation around twenty four pox genes that we identified, or which were reported in the literature, as candidate HTgenes. We found support for multiple independent transfers into poxviruses for five HTgenes. Three of these genes are known to be important for the survival of the virus in or out of the host cell and one of them increases susceptibility to some antiviral drugs.

Conclusion

In related genomes conserved synteny information can provide convincing evidence for multiple independent horizontal gene transfer events even in the absence of a robust phylogenetic tree for the HTgene.

Topological variation in single-gene phylogenetic trees

A large-scale phylogenetic study of the human lineage dramatically points up the problems of using single genes to build phylogenetic trees.

A recent large-scale phylogenomic study has shown the great degree of topological variation that can be found among eukaryotic phylogenetic trees constructed from single genes, highlighting the problems that can be associated with gene sampling in phylogenetic studies.

Two closely related pathways of nicotine catabolism in Arthrobacter nicotinovorans and Nocardioides sp. strain JS614

A virtually identical nicotine catabolic pathway including the heterotrimeric molybdenum enzyme nicotine and 6-hydroxy-pseudo-oxynicotine dehydrogenase, 6-hydroxy-L: -nicotine oxidase, 2,6-dihydroxy-pseudo-oxynicotine hydrolase, and 2,6-dihydroxypyridine hydroxylase have been identified in A. nicotinovorans and Nocardioides sp. JS614. Enzymes catalyzing the same reactions and similar protein antigens were detected in the extracts of the two microorganisms. Nicotine blue and methylamine, two end products of nicotine catabolism were detected in the growth medium of both bacterial species. Nicotine catabolic genes are clustered on pAO1 in A. nicotinovorans, but located chromosomally in Nocardioides sp. JS614.

Wolbachia genome integrated in an insect chromosome: evolution and fate of laterally transferred endosymbiont genes

Recent accumulation of microbial genome data has demonstrated that lateral gene transfers constitute an important and universal evolutionary process in prokaryotes, while those in multicellular eukaryotes are still regarded as unusual, except for endosymbiotic gene transfers from mitochondria and plastids. Here we thoroughly investigated the bacterial genes derived from a Wolbachia endosymbiont on the nuclear genome of the beetle Callosobruchus chinensis. Exhaustive PCR detection and Southern blot analysis suggested that approximately 30% of Wolbachia genes, in terms of the gene repertoire of wMel, are present on the insect nuclear genome. Fluorescent in situ hybridization located the transferred genes on the proximal region of the basal short arm of the X chromosome. Molecular evolutionary and other lines of evidence indicated that the transferred genes are probably derived from a single lateral transfer event. The transferred genes were, for the length examined, structurally disrupted, freed from functional constraints, and transcriptionally inactive. Hence, most, if not all, of the transferred genes have been pseudogenized. Notwithstanding this, the transferred genes were ubiquitously detected from Japanese and Taiwanese populations of C. chinensis, while the number of the transferred genes detected differed between the populations. The transferred genes were not detected from congenic beetle species, indicating that the transfer event occurred after speciation of C. chinensis, which was estimated to be one or several million years ago. These features of the laterally transferred endosymbiont genes are compared with the evolutionary patterns of mitochondrial and plastid genome fragments acquired by nuclear genomes through recent endosymbiotic gene transfers.

Distinct functions for the two PsbP-like proteins PPL1 and PPL2 in the chloroplast thylakoid lumen of Arabidopsis

PsbP, an extrinsic subunit of photosystem II (PSII), is a nuclear-encoded protein that optimizes the water-splitting reaction in vivo. In addition to PsbP, higher plants have two nuclear-encoded genes for PsbP homologs (PsbP-like proteins [PPLs]) that show significant sequence similarity to a cyanobacterial PsbP homolog (cyanoP); however, the function of PPLs in higher plants has not yet been elucidated. In this study, we characterized Arabidopsis (Arabidopsis thaliana) mutants lacking either of two PPLs, PPL1 and PPL2. Phylogenetic analysis suggests that PPL1 would be an ortholog of cyanoP, and PPL2 and PsbP may have a paralogous relationship with PPL1. Analysis on mRNA expression profiles showed that PPL1 expressed under stress conditions and PPL2 coexpressed with the subunits of chloroplast NAD(P)H dehydrogenase (NDH) complex. Consistent with these suggestions, PSII activity in a ppl1 mutant was more sensitive to high-intensity light than wild type, and the recovery of photoinhibited PSII activity was delayed in ppl1 plants. Therefore, PPL1 is required for efficient repair of photodamaged PSII. Furthermore, the stoichiometric level and activity of the chloroplast NDH complex in thylakoids were severely decreased in a ppl2 mutant, demonstrating that PPL2 is a novel thylakoid lumenal factor required for accumulation of the chloroplast NDH complex. These results suggest that during endosymbiosis and subsequent gene transfer to the host nucleus, cyanoP from ancient cyanobacteria evolved into PPL1, PPL2, and PsbP, and each of them has a distinct role in photosynthetic electron transfer in Arabidopsis.

Global extent of horizontal gene transfer

Horizontal gene transfer (HGT) is thought to play an important role in the evolution of species and innovation of genomes. There have been many convincing evidences for HGT for specific genes or gene families, but there has been no estimate of the global extent of HGT. Here, we present a method of identifying HGT events within a given protein family and estimate the global extent of HGT in all curated protein domain families ( approximately 8,000) listed in the Pfam database. The results suggest four conclusions: (i) for all protein domain families in Pfam, the fixation of genes horizontally transferred is not a rampant phenomenon between organisms with substantial phylogenetic separations (1.1-9.7% of Pfam families surveyed at three taxonomic ranges studied show indication of HGT); (ii) however, at the level of domains, >50% of Archaea have one or more protein domains acquired by HGT, and nearly 30-50% of Bacteria did the same when examined at three taxonomic ranges. But, the equivalent value for Eukarya is <10%; (iii) HGT will have very little impact in the construction of organism phylogeny, when the construction methods use whole genomes, large numbers of common genes, or SSU rRNAs; and (iv) there appears to be no strong preference of HGT for protein families of particular cellular or molecular functions.

Plasmids and rickettsial evolution: insight from Rickettsia felis

BACKGROUND

The genome sequence of Rickettsia felis revealed a number of rickettsial genetic anomalies that likely contribute not only to a large genome size relative to other rickettsiae, but also to phenotypic oddities that have confounded the categorization of R. felis as either typhus group (TG) or spotted fever group (SFG) rickettsiae. Most intriguing was the first report from rickettsiae of a conjugative plasmid (pRF) that contains 68 putative open reading frames, several of which are predicted to encode proteins with high similarity to conjugative machinery in other plasmid-containing bacteria.

METHODOLOGY/PRINCIPAL FINDINGS

Using phylogeny estimation, we determined the mode of inheritance of pRF genes relative to conserved rickettsial chromosomal genes. Phylogenies of chromosomal genes were in agreement with other published rickettsial trees. However, phylogenies including pRF genes yielded different topologies and suggest a close relationship between pRF and ancestral group (AG) rickettsiae, including the recently completed genome of R. bellii str. RML369-C. This relatedness is further supported by the distribution of pRF genes across other rickettsiae, as 10 pRF genes (or inactive derivatives) also occur in AG (but not SFG) rickettsiae, with five of these genes characteristic of typical plasmids. Detailed characterization of pRF genes resulted in two novel findings: the identification of oriV and replication termination regions, and the likelihood that a second proposed plasmid, pRFdelta, is an artifact of the original genome assembly.

CONCLUSION/SIGNIFICANCE

Altogether, we propose a new rickettsial classification scheme with the addition of a fourth lineage, transitional group (TRG) rickettsiae, that is unique from TG and SFG rickettsiae and harbors genes from possible exchanges with AG rickettsiae via conjugation. We offer insight into the evolution of a plastic plasmid system in rickettsiae, including the role plasmids may have played in the acquirement of virulence traits in pathogenic strains, and the likely origin of plasmids within the rickettsial tree.

In planta horizontal transfer of a major pathogenicity effector gene

Xanthomonas citri pv. citri is a clonal group of strains that causes citrus canker disease and appears to have originated in Asia. A phylogenetically distinct clonal group that causes identical disease symptoms on susceptible citrus, X. citri pv. aurantifolii, arose more recently in South America. Genomes of X. citri pv. aurantifolii strains carry two DNA fragments that hybridize to pthA, an X. citri pv. citri gene which encodes a major type III pathogenicity effector protein that is absolutely required to cause citrus canker. Marker interruption mutagenesis and complementation revealed that X. citri pv. aurantifolii strain B69 carried one functional pthA homolog, designated pthB, that was required to cause cankers on citrus. Gene pthB was found among 38 open reading frames on a 37,106-bp plasmid, designated pXcB, which was sequenced and annotated. No additional pathogenicity effectors were found on pXcB, but 11 out of 38 open reading frames appeared to encode a type IV transfer system. pXcB transferred horizontally in planta, without added selection, from B69 to a nonpathogenic X. citri pv. citri (pthA::Tn5) mutant strain, fully restoring canker. In planta transfer efficiencies were very high (>0.1%/recipient) and equivalent to those observed for agar medium with antibiotic selection, indicating that pthB conferred a strong selective advantage to the recipient strain. A single pathogenicity effector that can confer a distinct selective advantage in planta may both facilitate plasmid survival following horizontal gene transfer and account for the origination of phylogenetically distinct groups of strains causing identical disease symptoms.

The mycorrhiza helper bacteria revisited

In natural conditions, mycorrhizal fungi are surrounded by complex microbial communities, which modulate the mycorrhizal symbiosis. Here, the focus is on the so-called mycorrhiza helper bacteria (MHB). This concept is revisited, and the distinction is made between the helper bacteria, which assist mycorrhiza formation, and those that interact positively with the functioning of the symbiosis. After considering some examples of MHB from the literature, the ecological and evolutionary implications of the relationships of MHB with mycorrhizal fungi are discussed. The question of the specificity of the MHB effect is addressed, and an assessment is made of progress in understanding the mechanisms of the MHB effect, which has been made possible through the development of genomics. Finally, clear evidence is presented suggesting that some MHB promote the functioning of the mycorrhizal symbiosis. This is illustrated for three critical functions of practical significance: nutrient mobilization from soil minerals, fixation of atmospheric nitrogen, and protection of plants against root pathogens. The review concludes with discussion of future research priorities regarding the potentially very fruitful concept of MHB.

Origin and evolution of the mitochondrial aminoacyl-tRNA synthetases

Many theories favor a fusion of 2 prokaryotic genomes for the origin of the Eukaryotes, but there are disagreements on the origin, timing, and cellular structures of the cells involved. Equally controversial is the source of the nuclear genes for mitochondrial proteins, although the alpha-proteobacterial contribution to the mitochondrial genome is well established. Phylogenetic inferences show that the nuclearly encoded mitochondrial aminoacyl-tRNA synthetases (aaRSs) occupy a position in the tree that is not close to any of the currently sequenced alpha-proteobacterial genomes, despite cohesive and remarkably well-resolved alpha-proteobacterial clades in 12 of the 20 trees. Two or more alpha-proteobacterial clusters were observed in 8 cases, indicative of differential loss of paralogous genes or horizontal gene transfer. Replacement and retargeting events within the nuclear genomes of the Eukaryotes was indicated in 10 trees, 4 of which also show split alpha-proteobacterial groups. A majority of the mitochondrial aaRSs originate from within the bacterial domain, but none specifically from the alpha-Proteobacteria. For some aaRS, the endosymbiotic origin may have been erased by ongoing gene replacements on the bacterial as well as the eukaryotic side. For others that accurately resolve the alpha-proteobacterial divergence patterns, the lack of affiliation with mitochondria is more surprising. We hypothesize that the ancestral eukaryotic gene pool hosted primordial "bacterial-like" genes, to which a limited set of alpha-proteobacterial genes, mostly coding for components of the respiratory chain complexes, were added and selectively maintained.

Origin of the bacterial SET domain genes: vertical or horizontal?

The presence of Supressor of variegation-Enhanser of zeste-Trithorax (SET) domain genes in bacteria is a current paradigm for lateral genetic exchange between eukaryotes and prokaryotes. Because a major function of SET domain proteins is the chemical modification of chromatin and bacteria do not have chromatin, there is no apparent functional requirement for the existence of bacterial SET domain genes. Consequently, their finding in only a small fraction of pathogenic and symbiotic bacteria was taken as evidence that bacteria have obtained the SET domain genes from their hosts. Furthermore, it was proposed that the products of the genes would, most likely, be involved in bacteria-host interactions. The broadened scope of sequenced bacterial genomes to include also free-living and environmental species provided a larger sample to analyze the bacterial SET domain genes. By phylogenetic analysis, examination of individual chromosomal regions for signs of insertion, and evaluating the chromosomal versus SET domain genes' GC contents, we provide evidence that SET domain genes have existed in the bacterial domain of life independently of eukaryotes. The bacterial genes have undergone an evolution of their own unconnected to the evolution of the eukaryotic SET domain genes. Initial finding of SET domain genes in predominantly pathogenic and symbiotic bacteria resulted, most probably, from a biased sample. However, a lateral transfer of SET domain genes may have occurred between some bacteria and a family of Archaea. A model for the evolution and distribution of SET domain genes in bacteria is proposed.

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.

The tree of one percent

As lateral gene transfer among prokaryotes and endosymbiotic gene transfer (from organelles) among eukaryotes are fundamentally not tree-like in nature, biologists need to depart from the notion that all genomes are related by a single bifurcating tree.

Two significant evolutionary processes are fundamentally not tree-like in nature - lateral gene transfer among prokaryotes and endosymbiotic gene transfer (from organelles) among eukaryotes. To incorporate such processes into the bigger picture of early evolution, biologists need to depart from the preconceived notion that all genomes are related by a single bifurcating tree.

Horizontal gene transfer of a plastid gene in the non-photosynthetic flowering plants Orobanche and Phelipanche (Orobanchaceae)

Plastid sequences are among the most widely used in phylogenetic and phylogeographic studies in flowering plants, where they are usually assumed to evolve like non-recombining, uniparentally transmitted, single-copy genes. Among others, this assumption can be violated by intracellular gene transfer (IGT) within cells or by the exchange of genes across mating barriers (horizontal gene transfer, HGT). We report on HGT of a plastid region including rps2, trnL-F, and rbcL in a group of non-photosynthetic flowering plants. Species of the parasitic broomrape genus Phelipanche harbor two copies of rps2, a plastid ribosomal gene, one corresponding to the phylogenetic position of the respective species, the other being horizontally acquired from the related broomrape genus Orobanche. While the vertically transmitted copies probably reside within the plastid genome, the localization of the horizontally acquired copies is not known. With both donor and recipient being parasitic plants, a possible pathway for the exchange of genetic material is via a commonly attacked host.

Comparative genomics of the class 4 histone deacetylase family indicates a complex evolutionary history

Background

Histone deacetylases are enzymes that modify core histones and play key roles in transcriptional regulation, chromatin assembly, DNA repair, and recombination in eukaryotes. Three types of related histone deacetylases (classes 1, 2, and 4) are widely found in eukaryotes, and structurally related proteins have also been found in some prokaryotes. Here we focus on the evolutionary history of the class 4 histone deacetylase family.

Results

Through sequence similarity searches against sequenced genomes and expressed sequence tag data, we identified members of the class 4 histone deacetylase family in 45 eukaryotic and 37 eubacterial species representative of very distant evolutionary lineages. Multiple phylogenetic analyses indicate that the phylogeny of these proteins is, in many respects, at odds with the phylogeny of the species in which they are found. In addition, the eukaryotic members of the class 4 histone deacetylase family clearly display an anomalous phyletic distribution.

Conclusion

The unexpected phylogenetic relationships within the class 4 histone deacetylase family and the anomalous phyletic distribution of these proteins within eukaryotes might be explained by two mechanisms: ancient gene duplication followed by differential gene losses and/or horizontal gene transfer. We discuss both possibilities in this report, and suggest that the evolutionary history of the class 4 histone deacetylase family may have been shaped by horizontal gene transfers.

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.

The slow road to the eukaryotic genome

The eukaryotic genome is a mosaic of eubacterial and archaeal genes in addition to those unique to itself. The mosaic may have arisen as the result of two prokaryotes merging their genomes, or from genes acquired from an endosymbiont of eubacterial origin. A third possibility is that the eukaryotic genome arose from successive events of lateral gene transfer over long periods of time. This theory does not exclude the endosymbiont, but questions whether it is necessary to explain the peculiar set of eukaryotic genes. We use phylogenetic studies and reconstructions of ancestral first appearances of genes on the prokaryotic phylogeny to assess evidence for the lateral gene transfer scenario. We find that phylogenies advanced to support fusion can also arise from a succession of lateral gene transfer events. Our reconstructions of ancestral first appearances of genes reveal that the various genes that make up the eukaryotic mosaic arose at different times and in diverse lineages on the prokaryotic tree, and were not available in a single lineage. Successive events of lateral gene transfer can explain the unusual mosaic structure of the eukaryotic genome, with its content linked to the immediate adaptive value of the genes its acquired. Progress in understanding eukaryotes may come from identifying ancestral features such as the eukaryotic splicesome that could explain why this lineage invaded, or created, the eukaryotic niche.

Evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase family of the enolase superfamily

Understanding how proteins evolve to provide both exquisite specificity and proficient activity is a fundamental problem in biology that has implications for protein function prediction and protein engineering. To study this problem, we analyzed the evolution of structure and function in the o-succinylbenzoate synthase/N-acylamino acid racemase (OSBS/NAAAR) family, part of the mechanistically diverse enolase superfamily. Although all characterized members of the family catalyze the OSBS reaction, this family is extraordinarily divergent, with some members sharing <15% identity. In addition, a member of this family, Amycolatopsis OSBS/NAAAR, is promiscuous, catalyzing both dehydration and racemization. Although the OSBS/NAAAR family appears to have a single evolutionary origin, no sequence or structural motifs unique to this family could be identified; all residues conserved in the family are also found in enolase superfamily members that have different functions. Based on their species distribution, several uncharacterized proteins similar to Amycolatopsis OSBS/NAAAR appear to have been transmitted by lateral gene transfer. Like Amycolatopsis OSBS/NAAAR, these might have additional or alternative functions to OSBS because many are from organisms lacking the pathway in which OSBS is an intermediate. In addition to functional differences, the OSBS/NAAAR family exhibits surprising structural variations, including large differences in orientation between the two domains. These results offer several insights into protein evolution. First, orthologous proteins can exhibit significant structural variation, and specificity can be maintained with little conservation of ligand-contacting residues. Second, the discovery of a set of proteins similar to Amycolatopsis OSBS/NAAAR supports the hypothesis that new protein functions evolve through promiscuous intermediates. Finally, a combination of evolutionary, structural, and sequence analyses identified characteristics that might prime proteins, such as Amycolatopsis OSBS/NAAAR, for the evolution of new activities.

Presence/absence polymorphism for alternative pathogenicity islands in Pseudomonas viridiflava, a pathogen of Arabidopsis

The contribution of arms race dynamics to plant-pathogen coevolution has been called into question by the presence of balanced polymorphisms in resistance genes of Arabidopsis thaliana, but less is known about the pathogen side of the interaction. Here we investigate structural polymorphism in pathogenicity islands (PAIs) in Pseudomonas viridiflava, a prevalent bacterial pathogen of A. thaliana. PAIs encode the type III secretion system along with its effectors and are essential for pathogen recognition in plants. P. viridiflava harbors two structurally distinct and highly diverged PAI paralogs (T- and S-PAI) that are integrated in different chromosome locations in the P. viridiflava genome. Both PAIs are segregating as presence/absence polymorphisms such that only one PAI ([T-PAI, nablaS-PAI] and [nablaT-PAI, S-PAI]) is present in any individual cell. A worldwide population survey identified no isolate with neither or both PAI. T-PAI and S-PAI genotypes exhibit virulence differences and a host-specificity tradeoff. Orthologs of each PAI can be found in conserved syntenic locations in other Pseudomonas species, indicating vertical phylogenetic transmission in this genus. Molecular evolutionary analysis of PAI sequences also argues against "recent" horizontal transfer. Spikes in nucleotide divergence in flanking regions of PAI and nabla-PAI alleles suggest that the dual PAI polymorphism has been maintained in this species under some form of balancing selection. Virulence differences and host specificities are hypothesized to be responsible for the maintenance of the dual PAI system in this bacterial pathogen.

A Phylogenetic and Structural Analysis of Truncated Hemoglobins

Truncated hemoglobins (trHbs) are heme proteins found in bacteria, plants, and unicellular eukaryotes. They are distantly related to vertebrate hemoglobins and are typically shorter than these by 20-40 residues. The multiple amino acid deletions, insertions, and replacements result in distinctive alterations of the canonical globin fold and a wide range of chemical properties. An early phylogenetic analysis categorized trHbs into three groups, I (trHbN), II (trHbO), and III (trHbP). Here, we revisit this analysis with 111 trHbs. We find that trHbs are orthologous within each group and paralogous across the groups. Group I globins form the most disparate set and separate into two divergent subgroups. Group II is comparatively homogeneous, whereas Group III displays the highest level of overall conservation. In Group I and Group II globins, for which some ligand binding and structural data are available, an improved description of probable protein-ligand interactions is achieved. Other conservation trends are either confirmed (essential glycines in loops), refined (lining of ligand access tunnel), or newly identified (helix start signal). The Group III globins, so far uncharacterized, exhibit recognizable heme cavity residues while lacking some of the residues thought to be important to the trHb fold. An analysis of the phylogenetic trees of each group provides a plausible scenario for the emergence of trHbs, by which the Group II trHb gene was the original gene, and the Group I trHb and Group III trHb genes were obtained via duplication and transfer events.

Lateral gene transfer of a dermonecrotic toxin between spiders and bacteria

MOTIVATION

Spiders in the genus Loxosceles, including the notoriously toxic brown recluse, cause severe necrotic skin lesions owing to the presence of a venom enzyme called sphingomyelinase D (SMaseD). This enzyme activity is unknown elsewhere in the animal kingdom but is shared with strains of pathogenic Corynebacteria that cause various illnesses in farm animals. The presence of the same toxic activity only in distantly related organisms poses an interesting and medically important question in molecular evolution.

RESULTS

We use superpositions of recently determined structures and sequence comparisons to infer that both bacterial and spider SMaseDs originated from a common, broadly conserved domain family, the glycerophosphoryl diester phosphodiesterases. We also identify a unique sequence/structure motif present in both SMaseDs but not in the ancestral family, supporting SMaseD origin through a single divergence event in either bacteria or spiders, followed by lateral gene transfer from one lineage to the other.

The last common ancestor: what's in a name?

Twenty completely sequenced cellular genomes from the three major domains were analyzed using twice one-way BLAST searches in order to define the set of the most conserved protein-encoding sequences to characterize the gene complement of the last common ancestor of extant life. The resulting set is dominated by different putative ATPases, and by molecules involved in gene expression and RNA metabolism. DEAD-type RNA helicase and enolase genes, which are known to be part of the RNA degradosome, are as conserved as many transcription and translation genes. This suggests the early evolution of a control mechanism for gene expression at the RNA level, providing additional support to the hypothesis that during early cellular evolution RNA molecules played a more prominent role. Conserved sequences related to biosynthetic pathways include those encoding putative phosphoribosyl pyrophosphate synthase and thioredoxin, which participate in nucleotide metabolism. Although the information contained in the available databases corresponds only to a minor portion of biological diversity, the sequences reported here are likely to be part of an essential and highly conserved pool of proteins domains common to all organisms.

Clustered genes related to sulfate respiration in uncultured prokaryotes support the theory of their concomitant horizontal transfer

The dissimilatory reduction of sulfate is an ancient metabolic process central to today's biogeochemical cycling of sulfur and carbon in marine sediments. Until now its polyphyletic distribution was most parsimoniously explained by multiple horizontal transfers of single genes rather than by a not-yet-identified "metabolic island." Here we provide evidence that the horizontal transfer of a gene cluster may indeed be responsible for the patchy distribution of sulfate-reducing prokaryotes (SRP) in the phylogenetic tree. We isolated three DNA fragments (32 to 41 kb) from uncultured, closely related SRP from DNA directly extracted from two distinct marine sediments. Fosmid ws39f7, and partially also fosmids ws7f8 and hr42c9, harbored a core set of essential genes for the dissimilatory reduction of sulfate, including enzymes for the reduction of sulfur intermediates and synthesis of the prosthetic group of the dissimilatory sulfite reductase. Genome comparisons suggest that encoded membrane proteins universally present among SRP are critical for electron transfer to cytoplasmic enzymes. In addition, novel, conserved hypothetical proteins that are likely involved in dissimilatory sulfate reduction were identified. Based on comparative genomics and previously published experimental evidence, a more comprehensive model of dissimilatory sulfate reduction is presented. The observed clustering of genes involved in dissimilatory sulfate reduction has not been previously found. These findings strongly support the hypothesis that genes responsible for dissimilatory sulfate reduction were concomitantly transferred in a single event among prokaryotes. The acquisition of an optimized gene set would enormously facilitate a successful implementation of a novel pathway.

A case of promiscuity: Agrobacterium's endless hunt for new partners

Agrobacterium tumefaciens is a phytopathogenic bacterium that induces the 'crown gall' disease in plants by transfer and integration of a segment of its tumor-inducing (Ti) plasmid DNA into the genome of numerous plant species that represent most of the higher plant families. Recently, it has been shown that, under laboratory conditions, the host range of Agrobacterium can be extended to non-plant eukaryotic organisms. These include yeast, filamentous fungi, cultivated mushrooms and human cultured cells. In this article, we present Agrobacterium-mediated transformation of non-plant organisms as a source of new protocols for genetic transformation, as a unique tool for genomic studies (insertional mutagenesis or targeted DNA integration) and as a useful model system to study bacterium-host cell interactions. Moreover, better knowledge of the DNA-transfer mechanisms from bacteria to eukaryotic organisms can also help in understanding horizontal gene transfer--a driving force throughout biological evolution.

Perspective on mutagenesis and repair: the standard model and alternate modes of mutagenesis

The basic ideas of replication, mutagenesis, and repair have outlined a picture of how point mutations occur that has provided a valuable framework for theory and experiment, much as the Standard Model of particle physics has done for our concept of fundamental particles. However, alternative modes of mutagenesis are being defined that are changing our perspective of the "Standard Model" of mutagenesis, requiring an expanded model. The genome is now envisioned as being in dynamic equilibrium between a multitude of forces for mutational change and forces that counteract such change. By maintaining a delicate balance between these forces, cells avoid unwanted or excessive mutations. Yet, cells allow mutagenesis to occur under certain conditions. We can define an emerging paradigm. Namely, mechanisms exist that can direct point mutations to specific designated genes or regions of genes. In some cases, this is achieved by specific enzymes, and in other cases high mutability is programmed into the sequence of certain genes to help generate diversity. In yet additional cases, general mutability is increased under stress, and selective forces allow the recovery of favorable mutants.

Genome phylogenies indicate a meaningful alpha-proteobacterial phylogeny and support a grouping of the mitochondria with the Rickettsiales

Placement of the mitochondrial branch on the tree of life has been problematic. Sparse sampling, the uncertainty of how lateral gene transfer might overwrite phylogenetic signals, and the uncertainty of phylogenetic inference have all contributed to the issue. Here we address this issue using a supertree approach and completed genomic sequences. We first determine that a sensible alpha-proteobacterial phylogenetic tree exists and that it can confidently be inferred using orthologous genes. We show that congruence across these orthologous gene trees is significantly better than might be expected by random chance. There is some evidence of horizontal gene transfer within the alpha-proteobacteria, but it appears to be restricted to a minority of genes ( approximately 23%) most of whom ( approximately 74%) can be categorized as operational. This means that placement of the mitochondrion should not be excessively hampered by interspecies gene transfer. We then show that there is a consistently strong signal for placement of the mitochondrion on this tree and that this placement is relatively insensitive to methodological approach or data set. A concatenated alignment was created consisting of 15 mitochondrion-encoded proteins that are unlikely to have undergone any lateral gene transfer in the timeline under consideration. This alignment infers that the sister group of the mitochondria, for the taxa that have been sampled, is the order Rickettsiales.

Assessing Protein Co-evolution in the Context of the Tree of Life Assists in the Prediction of the Interactome

The identification of the whole set of protein interactions taking place in an organism is one of the main tasks in genomics, proteomics and systems biology. One of the computational techniques used by many investigators for studying and predicting protein interactions is the comparison of evolutionary histories (phylogenetic trees), under the hypothesis that interacting proteins would be subject to a similar evolutionary pressure resulting in a similar topology of the corresponding trees. Here, we present a new approach to predict protein interactions from phylogenetic trees, which incorporates information on the overall evolutionary histories of the species (i.e. the canonical "tree of life") in order to correct by the expected background similarity due to the underlying speciation events. We test the new approach in the largest set of annotated interacting proteins for Escherichia coli. This assessment of co-evolution in the context of the tree of life leads to a highly significant improvement (P(N) by sign test approximately 10E-6) in predicting interaction partners with respect to the previous technique, which does not incorporate information on the overall speciation tree. For half of the proteins we found a real interactor among the 6.4% top scores, compared with the 16.5% by the previous method. We applied the new method to the whole E.coli proteome and propose functions for some hypothetical proteins based on their predicted interactors. The new approach allows us also to detect non-canonical evolutionary events, in particular horizontal gene transfers. We also show that taking into account these non-canonical evolutionary events when assessing the similarity between evolutionary trees improves the performance of the method predicting interactions.

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.

Taxonomic utility of a phylogenetic analysis of phosphoglycerate kinase proteins of Archaea, Bacteria, and Eukaryota: Insights by Bayesian analyses

We studied 131 protein sequences of the essentially ubiquitous glycolytic enzyme 3-phosphoglycerate kinase (3-PGK) by Bayesian analyses in three Domains: 15 Archaea, 83 Bacteria, and 33 Eukaryota. The posterior distribution of phylogenetic trees developed were based on a uniform prior, the WAG model of protein evolution, Metropolis-Hastings sampling in a Markov chain Monte Carlo analysis, and a package of diagnostics to critically evaluate the validity of the analyses. The 15 Archaea separated with high posterior probability. The archaean Phyla Euryarchaeota and the apparently Euryarchaeota derived Crenarchaeota were monophyletic. The 33 Eukaryota separated into two main groups: the non-chlorophyllous forms with coherent sub-groupings of Euglenozoa, Alveolata, Fungi, and Metazoa and all the chlorophyllous species studied: the Plantae (Viridaeplantae), chlorophyllous Stramenopiles, and the chlorophyllous Bacteria. This association supports other opinions concerning the related lineage of cyanobacteria and the Plantae. The 3-PGK sequences from 83 Bacteria in almost every instance associated by their recognized taxal group: alpha-, beta-, gamma-, epsilon-proteobacteria, Chlamydia, Actinobacteridae, and Firmicutes. Firmicutes sequences were subdivided into three apparently monophyletic groups: the anaerobic Clostridia, the spore-forming Bacillales and a group containing the Mollicutes, Lactobacillales and non-spore-forming Bacillales. The 3-PGK-gene tree assemblage was notable both for its pervasive clustering in three Domains according to recognized taxonomic groupings of Class, Order, Family, and Genus. The 3-PGK enzyme or 3-PGK-like activity may have played a central role in the metabolism of the Universal Ancestor.

Stochastic models for horizontal gene transfer: taking a random walk through tree space

Horizontal gene transfer (HGT) plays a critical role in evolution across all domains of life with important biological and medical implications. I propose a simple class of stochastic models to examine HGT using multiple orthologous gene alignments. The models function in a hierarchical phylogenetic framework. The top level of the hierarchy is based on a random walk process in "tree space" that allows for the development of a joint probabilistic distribution over multiple gene trees and an unknown, but estimable species tree. I consider two general forms of random walks. The first form is derived from the subtree prune and regraft (SPR) operator that mirrors the observed effects that HGT has on inferred trees. The second form is based on walks over complete graphs and offers numerically tractable solutions for an increasing number of taxa. The bottom level of the hierarchy utilizes standard phylogenetic models to reconstruct gene trees given multiple gene alignments conditional on the random walk process. I develop a well-mixing Markov chain Monte Carlo algorithm to fit the models in a Bayesian framework. I demonstrate the flexibility of these stochastic models to test competing ideas about HGT by examining the complexity hypothesis. Using 144 orthologous gene alignments from six prokaryotes previously collected and analyzed, Bayesian model selection finds support for (1) the SPR model over the alternative form, (2) the 16S rRNA reconstruction as the most likely species tree, and (3) increased HGT of operational genes compared to informational genes.

Paradigm change in evolutionary microbiology

Thomas Kuhn had little to say about scientific change in biological science, and biologists are ambivalent about how applicable his framework is for their disciplines. We apply Kuhn's account of paradigm change to evolutionary microbiology, where key Darwinian tenets are being challenged by two decades of findings from molecular phylogenetics. The chief culprit is lateral gene transfer, which undermines the role of vertical descent and the representation of evolutionary history as a tree of life. To assess Kuhn's relevance to this controversy, we add a social analysis of the scientists involved to the historical and philosophical debates. We conclude that while Kuhn's account may capture aspects of the pattern (or outcome) of an episode of scientific change, he has little to say about how the process of generating new understandings is occurring in evolutionary microbiology. Once Kuhn's application is limited to that of an initial investigative probe into how scientific problem-solving occurs, his disciplinary scope becomes broader.

Gene transfer to plants by diverse species of bacteria

Agrobacterium is widely considered to be the only bacterial genus capable of transferring genes to plants. When suitably modified, Agrobacterium has become the most effective vector for gene transfer in plant biotechnology. However, the complexity of the patent landscape has created both real and perceived obstacles to the effective use of this technology for agricultural improvements by many public and private organizations worldwide. Here we show that several species of bacteria outside the Agrobacterium genus can be modified to mediate gene transfer to a number of diverse plants. These plant-associated symbiotic bacteria were made competent for gene transfer by acquisition of both a disarmed Ti plasmid and a suitable binary vector. This alternative to Agrobacterium-mediated technology for crop improvement, in addition to affording a versatile 'open source' platform for plant biotechnology, may lead to new uses of natural bacteria-plant interactions to achieve plant transformation.

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.

Discovery of a free-living chlorophyll d-producing cyanobacterium with a hybrid proteobacterial/cyanobacterial small-subunit rRNA gene

Chlorophyll d-producing cyanobacteria are a recently described group of phototrophic bacteria that is a major focus of photosynthesis research, previously known only from marine environments in symbiosis with eukaryotes. We have discovered a free-living member of this group from a eutrophic hypersaline lake. Phylogenetic analyses indicated these strains are closely related to each other but not to prochlorophyte cyanobacteria that also use an alternative form of chlorophyll as the major light-harvesting pigment. We have also demonstrated that these bacteria acquired a fragment of the small-subunit rRNA gene encoding a conserved hairpin in the bacterial ribosome from a proteobacterial donor at least 10 million years before the present. Thus, our most widely used phylogenetic marker can be a mosaic of sequence fragments with widely divergent evolutionary histories.

Pleomorphic mammalian tumor-derived bacteria self-organize as multicellular mammalian eukaryotic-like organisms: morphogenetic properties in vitro, possible origins, and possible roles in mammalian `tumor ecologies'

Highly pleomorphic bacteria have regularly been isolated from mammalian tumors and leukemic bloods. Here, it is shown that highly pleomorphic, cell-wall deficient bacteria derived from a mammalian tumor self-organize in vitro into mammalian tissue-like morphogenetic patterns consisting of multicellular tissue-like sheets and capillary-like networks. It is proposed that these pleomorphic mammalian tumor-derived (MTD) bacteria, during morphogenesis, express mammalian tissue morphogenesis-related genes that were acquired through eukaryote-to-prokaryote DNA transfer. Similar pleomorphic MTD bacteria might play important roles as symbiotic multicellular mammalian eukaryotic-like organisms in mammalian 'tumor ecologies' that include malignant and nonmalignant mammalian eukaryotic cells. From a mammalian tumor ecology perspective, eradication of tumors in some mammalian hosts may depend upon the elimination of pleomorphic MTD bacteria self-organized as symbiotic multicellular mammalian eukaryotic-like organisms. Further investigations of the extraordinary mammalian eukaryotic-like multicellularity of these bacteria may yield fundamental insights into the evolution of multicellularity and multicellular development and may challenge basic assumptions regarding cellular evolution.

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.

Phylogenetic inference in Rafflesiales: the influence of rate heterogeneity and horizontal gene transfer

Background

The phylogenetic relationships among the holoparasites of Rafflesiales have remained enigmatic for over a century. Recent molecular phylogenetic studies using the mitochondrial matR gene placed Rafflesia, Rhizanthes and Sapria (Rafflesiaceae s. str.) in the angiosperm order Malpighiales and Mitrastema (Mitrastemonaceae) in Ericales. These phylogenetic studies did not, however, sample two additional groups traditionally classified within Rafflesiales (Apodantheaceae and Cytinaceae). Here we provide molecular phylogenetic evidence using DNA sequence data from mitochondrial and nuclear genes for representatives of all genera in Rafflesiales.

Results

Our analyses indicate that the phylogenetic affinities of the large-flowered clade and Mitrastema, ascertained using mitochondrial matR, are congruent with results from nuclear SSU rDNA when these data are analyzed using maximum likelihood and Bayesian methods. The relationship of Cytinaceae to Malvales was recovered in all analyses. Relationships between Apodanthaceae and photosynthetic angiosperms varied depending upon the data partition: Malvales (3-gene), Cucurbitales (matR) or Fabales (atp1). The latter incongruencies suggest that horizontal gene transfer (HGT) may be affecting the mitochondrial gene topologies. The lack of association between Mitrastema and Ericales using atp1 is suggestive of HGT, but greater sampling within eudicots is needed to test this hypothesis further.

Conclusions

Rafflesiales are not monophyletic but composed of three or four independent lineages (families): Rafflesiaceae, Mitrastemonaceae, Apodanthaceae and Cytinaceae. Long-branch attraction appears to be misleading parsimony analyses of nuclear small-subunit rDNA data, but model-based methods (maximum likelihood and Bayesian analyses) recover a topology that is congruent with the mitochondrial matR gene tree, thus providing compelling evidence for organismal relationships. Horizontal gene transfer appears to be influencing only some taxa and some mitochondrial genes, thus indicating that the process is acting at the single gene (not whole genome) level.

Genome-wide molecular clock and horizontal gene transfer in bacterial evolution

We describe a simple theoretical framework for identifying orthologous sets of genes that deviate from a clock-like model of evolution. The approach used is based on comparing the evolutionary distances within a set of orthologs to a standard intergenomic distance, which was defined as the median of the distribution of the distances between all one-to-one orthologs. Under the clock-like model, the points on a plot of intergenic distances versus intergenomic distances are expected to fit a straight line. A statistical technique to identify significant deviations from the clock-like behavior is described. For several hundred analyzed orthologous sets representing three well-defined bacterial lineages, the alpha-Proteobacteria, the gamma-Proteobacteria, and the Bacillus-Clostridium group, the clock-like null hypothesis could not be rejected for approximately 70% of the sets, whereas the rest showed substantial anomalies. Subsequent detailed phylogenetic analysis of the genes with the strongest deviations indicated that over one-half of these genes probably underwent a distinct form of horizontal gene transfer, xenologous gene displacement, in which a gene is displaced by an ortholog from a different lineage. The remaining deviations from the clock-like model could be explained by lineage-specific acceleration of evolution. The results indicate that although xenologous gene displacement is a major force in bacterial evolution, a significant majority of orthologous gene sets in three major bacterial lineages evolved in accordance with the clock-like model. The approach described here allows rapid detection of deviations from this mode of evolution on the genome scale.

Statistical Analysis of the Spatial Distribution of Operons in the Transcriptional Regulation Network of Escherichia coli

We have performed a statistical analysis of the spatial distribution of operons along the DNA in the transcriptional regulation network of Escherichia coli. The analysis reveals that pairs of operons that regulate each other and those that are co-regulated tend to lie much closer to one another than would be expected for a random network. Moreover, these pairs of operons tend to be transcribed in diverging directions. This spatial arrangement of operons allows the upstream regulatory domains to overlap and interfere with each other and our analysis also demonstrates the statistical significance of this motif of overlapping operons. Overlapping operons afford additional regulatory control, such as the correlated or anticorrelated expression of operons. We show by a mean-field analysis of a feed-forward loop that overlapping operons can drastically enhance the performance of gene regulatory networks. Our results suggest that regulatory control can provide a selective pressure that drives operons together in the course of evolution.

Multivariate entropy distance method for prokaryotic gene identification

A new simple method is found for efficient and accurate identification of coding sequences in prokaryotic genome. The method employs a Shannon description of artificial language for DNA sequences. It consists in translating a DNA sequence into a pseudo-amino acid sequence with 20 fundamental words according to the universal genetic code. With an entropy-density profile (EDP), the method maps a sequence of finite length to a vector and then analyzes its position in the 20-dimensional phase space depending on its nature. It is found that the ratio of the relative distance to an averaged coding and non-coding EDP over a small number (up to one) of open reading frames (ORFs) can serve as a good coding potential. An iterative algorithm is designed for finding a set of "root" sequences using this coding potential. A multivariate entropy distance (MED) algorithm is then proposed for the identification of prokaryotic genes; it has a feature to combine the use of a coding potential and an EDP-based sequence similarity analysis. The current version of MED is unsupervised, parameter-free and simple to implement. It is demonstrated to be able to detect 95-99% genes with 10-30% of additional genes when tested against the RefSeq database of NCBI and to detect 97.5-99.8% of confirmed genes with known functions. It is also shown to be able to find a set of (functionally known) genes that are missed by other well-known gene finding algorithms. All measurements show that the MED algorithm reaches a similar performance level as the algorithms like GeneMark and Glimmer for prokaryotic gene prediction.