Phylogenetic analyses of cyanobacterial genomes: quantification of horizontal gene transfer events

Evolutionary dynamics of light-independent protochlorophyllide oxidoreductase genes in the secondary plastids of cryptophyte algae

Plastid genes encoding light-independent protochlorophyllide oxidoreductase (LIPOR) subunits were isolated from cryptophyte algae, the first example of such genes in plastids of secondary endosymbiotic origin. The presence of functional and nonfunctional copies of LIPOR genes in cryptophytes suggests that light-independent chlorophyll biosynthesis is a nonessential pathway in these organisms.

Structural analysis of a non-ribosomal halogenated cyclic peptide and its putative operon from Microcystis: implications for evolution of cyanopeptolins

The structure of the major peptide produced by Microcystis cf. wesenbergii NIVA-CYA 172/5, the halogenated heptapeptide cyanopeptolin-984, was determined using LC/MS/MS. A gene cluster encoding a peptide synthetase putatively producing a cyanopeptolin was cloned from the same strain and sequenced. The cluster consists of four genes encoding peptide synthetases and one gene encoding a halogenase. Two additional ORFs transcribed in the opposite direction were found in the 5' flanking sequence; one of these encodes an ABC transporter. The overall organization of the cyanopeptolin synthetase operon (mcn) resembles a previously analysed anabaenopeptilide synthetase operon (apd) from Anabaena strain 90. Phylogenetic analyses of the individual domains from Mcn, Apd and other cyanobacterial peptide synthetases showed clustering of the adenylation domains according to function irrespective of operon origin - indicating strong functional constraints across peptide synthetases. In contrast, the condensation and thiolation domains to a large extent grouped according to operon affiliation or position in the respective operons. Phylogenetic analyses of condensation domains indicated that N-terminal domains and domains that condense L-amino acids and D-amino acids, respectively, form three separate groups. Although recombination events are likely to be involved in the evolution of mcn, no clear evidence of genetic recombination between the two cyanopeptolin gene clusters was found. Within the genus Microcystis, microcystin and cyanopeptolin synthetases have an evolutionary history of genomic coexistence. However, the data indicated that the two classes of peptide synthetase gene clusters have evolved independently.

Comparison of cyanopeptolin genes in Planktothrix, Microcystis, and Anabaena strains: evidence for independent evolution within each genus

The major cyclic peptide cyanopeptolin 1138, produced by Planktothrix strain NIVA CYA 116, was characterized and shown to be structurally very close to the earlier-characterized oscillapeptin E. A cyanopeptolin gene cluster likely to encode the corresponding peptide synthetase was sequenced from the same strain. The 30-kb oci gene cluster contains two novel domains previously not detected in nonribosomal peptide synthetase gene clusters (a putative glyceric acid-activating domain and a sulfotransferase domain), in addition to seven nonribosomal peptide synthetase modules. Unlike in two previously described cyanopeptolin gene clusters from Anabaena and Microcystis, a halogenase gene is not present. The three cyanopeptolin gene clusters show similar gene and domain arrangements, while the binding pocket signatures deduced from the adenylation domain sequences and the additional tailoring domains vary. This suggests loss and gain of tailoring domains within each genus, after the diversification of the three clades, as major events leading to the present diversity. The ABC transporter genes associated with the cyanopeptolin gene clusters form a monophyletic clade and accordingly are likely to have evolved as part of the functional unit. Phylogenetic analyses of adenylation and condensation domains, including domains from cyanopeptolins and microcystins, show a closer similarity between the Planktothrix and Microcystis cyanopeptolin domains than between these and the Anabaena domain. No clear evidence of recombination between cyanopeptolins and microcystins could be detected. There were no strong indications of horizontal gene transfer of cyanopeptolin gene sequences across the three genera, supporting independent evolution within each genus.

Extensive horizontal transfer of core genome genes between two Lactobacillus species found in the gastrointestinal tract

BACKGROUND

While genes that are conserved between related bacterial species are usually thought to have evolved along with the species, phylogenetic trees reconstructed for individual genes may contradict this picture and indicate horizontal gene transfer. Individual trees are often not resolved with high confidence, however, and in that case alternative trees are generally not considered as contradicting the species tree, although not confirming it either. Here we conduct an in-depth analysis of 401 protein phylogenetic trees inferred with varying levels of confidence for three lactobacilli from the acidophilus complex. At present the relationship between these bacteria, isolated from environments as diverse as the gastrointestinal tract (Lactobacillus acidophilus and Lactobacillus johnsonii) and yogurt (Lactobacillus delbrueckii ssp. bulgaricus), is ambiguous due to contradictory phenotypical and 16S rRNA based classifications.

RESULTS

Among the 401 phylogenetic trees, those that could be reconstructed with high confidence support the 16S-rRNA tree or one alternative topology in an astonishing 3:2 ratio, while the third possible topology is practically absent. Lowering the confidence threshold for trees to be taken into consideration does not significantly affect this ratio, and therefore suggests that gene transfer may have affected as much as 40% of the core genome genes. Gene function bias suggests that the 16S rRNA phylogeny of the acidophilus complex, which indicates that L. acidophilus and L. delbrueckii ssp. bulgaricus are the closest related of these three species, is correct. A novel approach of comparison of interspecies protein divergence data employed in this study allowed to determine that gene transfer most likely took place between the lineages of the two species found in the gastrointestinal tract.

CONCLUSION

This case-study reports an unprecedented level of phylogenetic incongruence, presumably resulting from extensive horizontal gene transfer. The data give a first indication of the large extent of gene transfer that may take place in the gastrointestinal tract and its accumulated effect. For future studies, our results should encourage a careful weighing of data on phylogenetic tree topology, confidence and distribution to conclude on the absence or presence and extent of horizontal gene transfer.

The ancestor of the Paulinella chromatophore obtained a carboxysomal operon by horizontal gene transfer from a Nitrococcus-like gamma-proteobacterium

BACKGROUND

Paulinella chromatophora is a freshwater filose amoeba with photosynthetic endosymbionts (chromatophores) of cyanobacterial origin that are closely related to free-living Prochlorococcus and Synechococcus species (PS-clade). Members of the PS-clade of cyanobacteria contain a proteobacterial form 1A RubisCO (ribulose-1,5-bisphosphate carboxylase/oxygenase) that was acquired by horizontal gene transfer (HGT) of a carboxysomal operon. In rDNA-phylogenies, the Paulinella chromatophore diverged basal to the PS-clade, raising the question whether the HGT occurred before or after the split of the chromatophore ancestor.

RESULTS

Phylogenetic analyses of the almost complete rDNA operon with an improved taxon sampling containing most known cyanobacterial lineages recovered the Paulinella chromatophore as sister to the complete PS-clade. The sequence of the complete carboxysomal operon of Paulinella was determined. Analysis of RubisCO large subunit (rbcL) sequences revealed that Paulinella shares the proteobacterial form 1A RubisCO with the PS-clade. The gamma-proteobacterium Nitrococcus mobilis was identified as sister of the Paulinella chromatophore and the PS-clade in the RubisCO phylogeny. Gene content and order in the carboxysomal operon correlates well with the RubisCO phylogeny demonstrating that the complete carboxysomal operon was acquired by the common ancestor of the Paulinella chromatophore and the PS-clade through HGT. The carboxysomal operon shows a significantly elevated AT content in Paulinella, which in the rbcL gene is confined to third codon positions. Combined phylogenies using rbcL and the rDNA-operon resulted in a nearly fully resolved tree of the PS-clade.

CONCLUSION

The HGT of the carboxysomal operon predated the divergence of the chromatophore ancestor from the PS-clade. Following HGT and divergence of the chromatophore ancestor, diversification of the PS-clade into at least three subclades occurred. The gamma-proteobacterium Nitrococcus mobilis represents the closest known relative to the donor of the carboxysomal operon. The isolated position of the Paulinella chromatophore in molecular phylogenies as well as its elevated AT content suggests that the Paulinella chromatophore has already undergone typical steps in the reductive evolution of an endosymbiont.

The power of phylogenetic approaches to detect horizontally transferred genes

BACKGROUND

Horizontal gene transfer plays an important role in evolution because it sometimes allows recipient lineages to adapt to new ecological niches. High genes transfer frequencies were inferred for prokaryotic and early eukaryotic evolution. Does horizontal gene transfer also impact phylogenetic reconstruction of the evolutionary history of genomes and organisms? The answer to this question depends at least in part on the actual gene transfer frequencies and on the ability to weed out transferred genes from further analyses. Are the detected transfers mainly false positives, or are they the tip of an iceberg of many transfer events most of which go undetected by current methods?

RESULTS

Phylogenetic detection methods appear to be the method of choice to infer gene transfers, especially for ancient transfers and those followed by orthologous replacement. Here we explore how well some of these methods perform using in silico transfers between the terminal branches of a gamma proteobacterial, genome based phylogeny. For the experiments performed here on average the AU test at a 5% significance level detects 90.3% of the transfers and 91% of the exchanges as significant. Using the Robinson-Foulds distance only 57.7% of the exchanges and 60% of the donations were identified as significant. Analyses using bipartition spectra appeared most successful in our test case. The power of detection was on average 97% using a 70% cut-off and 94.2% with 90% cut-off for identifying conflicting bipartitions, while the rate of false positives was below 4.2% and 2.1% for the two cut-offs, respectively. For all methods the detection rates improved when more intervening branches separated donor and recipient.

CONCLUSION

Rates of detected transfers should not be mistaken for the actual transfer rates; most analyses of gene transfers remain anecdotal. The method and significance level to identify potential gene transfer events represent a trade-off between the frequency of erroneous identification (false positives) and the power to detect actual transfer events.

Pattern pluralism and the Tree of Life hypothesis

Darwin claimed that a unique inclusively hierarchical pattern of relationships between all organisms based on their similarities and differences [the Tree of Life (TOL)] was a fact of nature, for which evolution, and in particular a branching process of descent with modification, was the explanation. However, there is no independent evidence that the natural order is an inclusive hierarchy, and incorporation of prokaryotes into the TOL is especially problematic. The only data sets from which we might construct a universal hierarchy including prokaryotes, the sequences of genes, often disagree and can seldom be proven to agree. Hierarchical structure can always be imposed on or extracted from such data sets by algorithms designed to do so, but at its base the universal TOL rests on an unproven assumption about pattern that, given what we know about process, is unlikely to be broadly true. This is not to say that similarities and differences between organisms are not to be accounted for by evolutionary mechanisms, but descent with modification is only one of these mechanisms, and a single tree-like pattern is not the necessary (or expected) result of their collective operation. Pattern pluralism (the recognition that different evolutionary models and representations of relationships will be appropriate, and true, for different taxa or at different scales or for different purposes) is an attractive alternative to the quixotic pursuit of a single true TOL.

A simulation test bed for hypotheses of genome evolution

MOTIVATION

Microbial genomes undergo evolutionary processes such as gene family expansion and contraction, variable rates and patterns of sequence substitution and lateral genetic transfer. Simulation tools are essential for both the generation of data under different evolutionary models and the validation of analytical methods on such data. However, meaningful investigation of phenomena such as lateral genetic transfer requires the simultaneous consideration of many underlying evolutionary processes.

RESULTS

We have developed EvolSimulator, a software package that combines non-stationary sequence and gene family evolution together with models of lateral genetic transfer, within a customizable birth-death model of speciation and extinction. Here, we examine simulated data sets generated with EvolSimulator using existing statistical techniques from the evolutionary literature, showing in detail each component of the simulation strategy.

AVAILABILITY

Source code, manual and other information are freely available at www.bioinformatics.org.au/evolsim.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Refuting phylogenetic relationships

BACKGROUND

Phylogenetic methods are philosophically grounded, and so can be philosophically biased in ways that limit explanatory power. This constitutes an important methodologic dimension not often taken into account. Here we address this dimension in the context of concatenation approaches to phylogeny.

RESULTS

We discuss some of the limits of a methodology restricted to verificationism, the philosophy on which gene concatenation practices generally rely. As an alternative, we describe a software which identifies and focuses on impossible or refuted relationships, through a simple analysis of bootstrap bipartitions, followed by multivariate statistical analyses. We show how refuting phylogenetic relationships could in principle facilitate systematics. We also apply our method to the study of two complex phylogenies: the phylogeny of the archaea and the phylogeny of the core of genes shared by all life forms. While many groups are rejected, our results left open a possible proximity of N. equitans and the Methanopyrales, of the Archaea and the Cyanobacteria, and as well the possible grouping of the Methanobacteriales/Methanoccocales and Thermosplasmatales, of the Spirochaetes and the Actinobacteria and of the Proteobacteria and firmicutes.

CONCLUSION

It is sometimes easier (and preferable) to decide which species do not group together than which ones do. When possible topologies are limited, identifying local relationships that are rejected may be a useful alternative to classical concatenation approaches aiming to find a globally resolved tree on the basis of weak phylogenetic markers.

REVIEWERS

This article was reviewed by Mark Ragan, Eugene V Koonin and J Peter Gogarten.