Ancient origin of the divergent forms of leucyl-tRNA synthetases in the Halobacteriales
© Andam et al.; licensee BioMed Central Ltd. 2012
Received: 4 October 2011
Accepted: 27 April 2012
Published: 13 June 2012
Horizontal gene transfer (HGT) has greatly impacted the genealogical history of many lineages, particularly for prokaryotes, with genes frequently moving in and out of a line of descent. Many genes that were acquired by a lineage in the past likely originated from ancestral relatives that have since gone extinct. During the course of evolution, HGT has played an essential role in the origin and dissemination of genetic and metabolic novelty.
Three divergent forms of leucyl-tRNA synthetase (LeuRS) exist in the archaeal order Halobacteriales, commonly known as haloarchaea. Few haloarchaeal genomes have the typical archaeal form of this enzyme and phylogenetic analysis indicates it clusters within the Euryarchaeota as expected. The majority of sequenced halobacterial genomes possess a bacterial form of LeuRS. Phylogenetic reconstruction puts this larger group of haloarchaea at the base of the bacterial domain. The most parsimonious explanation is that an ancient transfer of LeuRS took place from an organism related to the ancestor of the bacterial domain to the haloarchaea. The bacterial form of LeuRS further underwent gene duplications and/or gene transfers within the haloarchaea, with some genomes possessing two distinct types of bacterial LeuRS. The cognate tRNALeu also reveals two distinct clusters for the haloarchaea; however, these tRNALeu clusters do not coincide with the groupings found in the LeuRS tree, revealing that LeuRS evolved independently of its cognate tRNA.
The study of leucyl-tRNA synthetase in haloarchaea illustrates the importance of gene transfer originating in lineages that went extinct since the transfer occurred. The haloarchaeal LeuRS and tRNALeu did not co-evolve.
The archaeal order Halobacteriales, commonly known as haloarchaea, consists of extremely halophilic, aerobic or facultative anaerobic organisms currently classified into 29 genera (http://www.the-icsp.org/taxa/halobacterlist.htm) These organisms are the dominant taxa in hypersaline ecosystems, such as salterns, salt and soda lakes and coastal areas, in which NaCl concentrations can reach 150–350 g/L . Members of the Halobacteriales are known to undergo frequent HGT and recombination [2–4]. The recently identified methylaspartate cycle for acetyl-CoA assimilation in haloarchaea consists of enzymes acquired through HGT. The pre-existing genes acquired from different bacterial donors were originally involved in various metabolic processes . Analyses of the bacteriorhodopsin and halorhodopsin proteins in the haloarchaea also suggest that HGT has played a role in their evolution .
The Halobacteriales are usually considered to have evolved from a group of halophilic methanogens. Phylogenies based on rRNA, concatenated proteins, and proteins involved in transcription and translation show the Halobacteriales as a sister group to the Methanosarcinales [7, 8]. However, whole-genome-based phylogenetic reconstructions often placed them at the base of the archaeal domain [9, 10], which might reflect the high number of genes in the haloarchaea that are of bacterial origin. Alternatively, they could be from extinct archaeal lineages that left a “fossil” in the molecular record.
Very few studies have provided evidence for ancient transfers from now-extinct lineages that existed prior to or during the time of last universal common ancestor (LUCA) of all life, or of each three domains. An example is the case of the rare pyrrolysyl-tRNA synthetase (PylRS) that charges the tRNAPyl with the non-canonical amino acid pyrrolysine (Pyl) . This rare enzyme has a restricted distribution, to date found only in members of the archaeal order Methanosarcinales, the firmicute Desulfitobacterium hafniense and a Deltaproteobacterium endosymbiont . In relation to the other aaRS, PylRS is placed as a deep-branching lineage within the aaRS subclass IIb, emerging prior to the most recent common ancestor of the bacterial and archaeal/eukaryal domains . The phylogenetic distribution of this enzyme suggests that these extant taxa acquired PylRS through several HGT episodes from an ancient, most likely extinct, lineage . A rare form of seryl-tRNA synthetase (SerRS), to date only found in some methanogens, based on phylogenetic reconstruction was suggested to have been acquired through HGT from a deep branching lineage . The patchy distribution of another uncommon Class II aaRS, O-phosphoseryl-tRNA synthetase (SepRS), is also indicative of ancient HGT. SepRS charges tRNACys with O-phosphoserine (Sep), a precursor of cysteine (Cys), to form Sep-tRNACys and is then converted to Cys-tRNACys. Phylogenetic analyses show that SepRS predates the duplication event that gave rise to the two phenylalanyl-tRNA synthetases (PheRS) subunits and also arose before the divergence of the organismal LUCA .
A challenge in the analyses of genetic contributions of ancient lineages to existing genomes is the absence of information about donor lineages because majority of them are already extinct. Genes that arose prior to the time of LUCA are expected to exhibit high divergence from their homologs. This would reflect an extremely long coalescence time to a most recent common molecular ancestor, occurring well before the organismal common ancestor . In the case of the PylRS , the rare form of SerRS , and SepRS , we can infer from their phylogenetic histories that they likely were already present prior to or during the time of the organismal LUCA.
In this study, we show the existence of two forms of LeuRS in the Halobacteriales that arose through ancient HGT. The bacterial form of LeuRS in the haloarchaea was likely acquired from a relative of the ancestor of the bacterial domain and further underwent gene duplication, transfer and divergence within the haloarchaea. We also discuss the impact of ancient HGT events in generating genetic diversity in present-day lineages.
Results and discussion
Two major clades of Halobacteriales in the LeuRS phylogeny
aaRS are ancient enzymes that catalyze the attachment of tRNA with its cognate amino acid during the translation process. This function is essential in maintaining the fidelity of the genetic code and all 20 aminoacyl-tRNA species are essential for all living organisms. Although aaRSs are part of the conserved "information processing and storage" gene set, aaRS are frequently transferred across species boundaries and even between domains [18–20], most likely due to the limited interactions with other biomolecules .
A single protein can contain parts that differ in phylogeny and substitution rates. We used GARD (Genetic Algorithm for Recombination Detection ) to investigate if different parts of the LeuRSs in haloarchaea have different histories. Using MUSCLE  and SATé  alignments, GARD determined breakpoints corresponding to position 780 (MUSCLE) and 628 (SATé) in the Halogeometricum LeuRS sequence, respectively. Further inspection of the multiple sequence alignment revealed that most of the phylogenetic information distinguishing the archaeal and bacterial type LeuRSs is contained in the larger amino terminal part of the alignment. This part contains the domain that catalyzes the esterification between leucine and tRNA, and contains many positions universally conserved between the domains. The carboxy terminal part of the alignment encodes the tRNA recognition domain. While GARD found a significant difference between the tree topologies determined for the two parts of the multiple sequence alignment, in both phylogenies reconstructed separately for the two parts the SATé alignment, the LeuRS-B sequences group at the base of the bacterial homologs, whereas LeuRS-A group with the euryarchaeal homologs (see Additional file 2: Figure S2). The role of the two parts of LeuRS in interacting with tRNALeu are illustrated in Additional file 3: Figure S3. Using the breakpoint from the GARD analysis of the MUSCLE alignment resulted in a carboxy terminal portion that was too short for reliable phylogenetic reconstruction. It is noteworthy that in the maximum likelihood phylogeny for this short fragment all haloarchaea grouped together, albeit with a bootstrap support value of only 47%. As most of the haloarchaeal fragments failed a chi-square test for compositional homogeneity, this finding may reflect a shared compositional bias in the haloarchaeal sequences, although the possibility that the carboxyterminal part of LeuRS might have a different evolutionary history from the rest of the enzyme cannot be excluded.
To explore the possibility that placement of the haloarchaeal LeuRS-B reflects an artifact created through long branch attraction, we calculated the pairwise distances between representatives of the bacterial LeuRS (Salinibacter ruber and Halanaerobium prevalens), archaeal LeuRS (Haloferax volcanii, Halogeometricum borinquense, Methanocorpusculum labreanum, Pyrococcus furiosus), haloarchaeal LeuRS-B (the two LeuRS-B copies in Halomicrobium mukohataei and Haloterrigena turkmenica) and the outgroup (Isoleucyl-tRNA synthetase from Methanopyrus kandleri and Thermotoga maritima). Mean pairwise distances from the outgroup do not show significant differences (0.5364 ± 0.0511 for the archaeal LeuRS, 0.3915 ± 0.0268 for the bacterial LeuRS, and 0.4038 ± 0.0791 for the haloarchaeal LeuRS-B). Analysis of compositional homogeneity using chi-square test as implemented in the program TREE-PUZZLE  indicated that the LeuRS-B sequences do not have atypical composition (P > 0.05). We do not find evidence that the placement of haloarchaeal LeuRS-B at the base of the bacterial homologs is due to an artifact created by these sequences being more divergent or having a different composition, and we find no indication of a close association of Halobacteriales LeuRS-B sequences with any specific bacterial or archaeal group. Nevertheless, artifacts created in the alignment certainly have the potential to increase apparent support values, thus a placement of the LeuRS-B sequences within the cluster of bacterial homologs cannot be excluded.
We performed more detailed phylogenetic analyses of the two haloarchaeal clusters and their closest relatives to determine the phylogenetic relationships among the members of each group (Figure 2). We analyzed 14 haloarchaeal genomes that were available in the NCBI completed microbial genome database. Out of these, only three genomes carry the LeuRS-A form – Haloferax volcanii, Halogeometricum borinquense and Haladaptatus paucihalophilus. Their sequences show close affinities to members of the Methanomicrobiales and Methanobacteriales (Figure 2a). The bacterial version LeuRS-B exhibits a more complicated picture (Figure 2b). Two highly-supported clusters can be observed, which we refer to as B’ and B”. In five of the genomes included in this study (Natrialba magadii, Haloterrigena turkmenica, Halomicrobium mukohataei, Haloarcula marismortui and Halorhabdus utahensis), both B’ and B” are present. Two possible scenarios can explain the observed distribution of LeuRS-B. The observation that B' and B'' group together at the base of the bacteria indicates their divergence occurred either in the donating lineage, or following the transfer. The two distinct scenarios are (a) the B form was already present in the haloarchaeal ancestor; versus (b) the B form was later acquired, but spread to different haloarchaeal groups through biased gene transfer .
A second line of support for HGT of the two B forms comes from parametric bootstrapping analysis as implemented in LGT3State . In this test, the null model requires that HGT is absent in the evolution of the LeuRS-B genes and that gene loss events can explain the distribution patterns. This model implies that the most recent Halobacteriales ancestor carried both types of LeuRS-B. The second model assumes gene losses and gains of the alternative forms can occur, that is, a genome carrying LeuRS-B’ can gain the LeuRS-B”, and vice versa, resulting in a genome with both types present, from which one type may eventually be lost. Using the LGT3State program , we generated 1000 bootstrap distributions under the gene loss only model. Thus, we have 1000 datasets reflecting the outcomes under the null model, which are compared to the real data. The distribution of the 1000 likelihood values gives us a measure of what to expect under the null hypothesis. The log-likelihood values obtained for the bootstrapped samples evaluated under the HGT model ranged from −43.2 to −49.6, and are much lower than the log-likelihood values when assuming the HGT model for the original tree (−6.35). Hence, we can reject the gene loss only model with a significance level of P<0.001.
Interestingly, we also observed that no genome possesses only the B” form (Figure 2b), i.e., B” is always found to co-exist with the B’. For the genomes that carry the two B copies, maintenance of the two functionally identical enzymes likely confers a selective advantage to the host. In bacteria, differential sensitivity of multiple copies of aaRS with redundant functions may benefit the organism against naturally occurring antibiotics . The antibiotic capabilities of Archaea have only recently been investigated. Peptide antibiotics produced by some members of the Archaea, referred to as archaeocins, have been identified from haloarchaea and Sulfolobus and were reported to exhibit cross-kingdom toxicity . A recent study showed that methanogenic archaea exhibit differences in susceptibility to various antibiotics, such as ampicillin, streptomycin, gentamicin, rifampicin, ofloxacin, tetracycline . It is also possible that there is a difference in the functional efficiency of the two LeuRS-B forms, with B” being less efficient in aminoacylating some of its cognate tRNAs. This may be similar to the intragenomic heterogeneity in the ribosomal operons of Haloarcula marismortui, which exhibit differences in gene expression under different environmental conditions . Alternatively, the functioning enzyme may consist of a B'B'' heterodimer, allowing more degrees of freedom to accommodate destabilizing mutations , as observed in Aquifex aeolicus[32, 33]; the transition from a homo- to a heterodimer initially might not have been adaptive, but the resulting heterodimer nevertheless may be under strong purifying selection . However, the latter scenario is unlikely as the genes encoding the B' and B" forms are located in different parts of the genomes (Figure 3).
Haladaptatus paucihalophilus possesses both the A and the B’ form of LeuRS (Figure 2). Both copies are located adjacent to each other and are divergently transcribed. Two of its flanking genes (coding for a thermosome subunit and alanine dehydrogenase) are also found in the genomic neighborhood of leuS in the other two haloarchaea that possess only the A form (Haloferax volcanii and Halogeometricum borinquense; Figure 3). This is compatible with the scenario that Haladaptatus originally had the A form and has subsequently acquired the B’ form through HGT from another haloarchaeon.
The archaeal and bacterial forms of LeuRS are significantly distinct from each other (Additional file 4: Table S1). The identities between the A and B forms range from 21-26%, reflecting the very deep divergence that gave rise to these two forms. In contrast, the two LeuRS-B forms exhibit 46–53% identity between the two B-types suggesting a more recent divergence event.
Scattered distribution of the different LeuRS in the Halobacteriales
Previous studies have reported the challenge of using the 16S rRNA phylogeny to determine the evolutionary relationships of the Halobacteriales . Two factors have been implicated: the presence of multiple divergent copies of this gene in a single genome in many haloarchaeal species and that recombination of the rRNA gene occurs frequently between species . Paralogous copies of rRNA operons in these organisms have been reported to show more than 5% divergence , and identical sequences have been found in strains that are otherwise clearly differentiated, making it difficult to establish accurate Halobacteriales relationships.
Mapping the presence and absence of the three LeuRS in the MLSA tree shows that all species belonging to clades I and III possess both B forms of the bacterial LeuRS. Given that LeuRS genes were frequently transferred within the haloarchaea, we do not interpret the co-occurrence of the B' and B" forms as shared derived character for clade I and III. For the archaeal version (LeuRS-A), we observed a dispersed distribution, mostly in branches that appear to have diverged more recently. If we consider the MLSA tree as a suitable representation of the species phylogeny of this group, and only take into account the distribution of LeuRS types within this group, then the initial assumption would be that the ancestor of the Halobacteriales possessed the bacterial form of LeuRS. However, another more likely scenario is that the presence of the archaeal version of the enzyme (LeuRS-A) is the ancestral state in the Halobacteriales. The clustering of the haloarchaeal LeuRS-A cluster within the euryarchaeal homologs, specifically with those from methanogens, would indicate shared ancestry [21, 22], and the archaeal LeuRS would be vertically inherited by the Halobacteriales. The single divergence event that gave rise to the B' and B'' forms likely took place early in the evolution of the Halobacteriales, followed by the spread or retention of both forms of LeuRS-B within the order.
Assuming that the Halobacteriales ancestor originally possessed the archaeal form acquired through vertical inheritance from the common ancestor of all Archaea, it later on gained the bacterial LeuRS through horizontal transfer from a deep branching bacterial lineage, possibly still unsampled or now extinct. The finding that the haloarchaeal LeuRS-B diverged before the homologs found in bacteria suggests that either the lineage donating LeuRS-B to the haloarchaea or the haloarchaea themselves coexisted with the bacterial most recent common ancestor. More than one lineage could have carried the bacterial version of LeuRS before it was transferred to the haloarchaea; however, provided that the deep branching of the haloarchaeal LeuRS form B is not an artifact, all the scenarios imply that the bacterial version now residing in the haloarchaea, coexisted with the ancestor of the bacterial domain. Following transfer to the haloarchaea, the bacterial form spread among the majority of the Halobacteriales through vertical inheritance and HGT biased toward close relatives [14, 41], with some species possessing one form while in others, both forms of the bacterial LeuRS are retained.
We then compared the LeuRS-A (Figure 2a) and LeuRS-B (Figure 2b) with the MLSA tree (Figure 4) to see if there are any conflicting topologies between them. For LeuRS-A, we observed similarity regarding the placement of the three species. Haloferax and Halogeometricum group together, and Haladaptatus is found at the base (Figure 2a). The topology of the LeuRS-B” tree was also similar to the MLSA tree, except for the placement of Halorhabdus (Figure 2b). This, however, is not highly supported and therefore we cannot draw any conclusion from it. In LeuRS-B”, the groupings of Natrialba and Haloterrigena, and of Haloarcula and Halomicrobium are similar to what we found in the MLSA tree. In comparing the LeuRS-B’ and the MLSA tree, we also observed the same clustering of the above mentioned two pairs of haloarchaea. An important conflict, however, is the phylogenetic position of Halomicrobium; the MLSA tree places it in clade III, while in the LeuRS-B’ tree, its position is highly supported at the base of the clade II (Figure 2a). Within clade III of the MLSA tree, Haloarcula and Halomicrobium have a closer relationship than with Halorhabdus. Hence, the LeuRS-B’ topology indicates a transfer from clade II to Halomicrobium. Another possible conflict is that of Natronomonas, which clusters with the clade II species in the LeuRS tree.
Topologies of the MLSA tree and each of the LeuRS trees indicate that (1) the Halobacteriales came to possess the archaeal form through common ancestry with the rest of the Archaea that was eventually lost in a majority of the Halobacteriales, and (2) the bacterial LeuRS types were vertically and horizontally inherited within the group. We can be certain that at least one HGT event took place – the transfer from a deep branching, currently unsampled bacterial lineage diverging most likely before the bacterial common ancestor to the Halobacteriales.
Archaeal tRNALeu phylogeny shows two groups of haloarchaea
Superficially similar to the LeuRS tree, two main groups of Halobacteriales are found in the tRNALeu tree (Figure 5). However, the distribution of the haloarchaea into the two groups differs significantly from that found in their corresponding synthetase tree. In the LeuRS tree, the smaller group of Halobacteriales consists of Haloferax, Halogeometricum and Haladaptatus, and the majority is found in a bigger cluster distinct from it (Figure 2). In contrast, the three genera mentioned above do not group together in the tRNALeu tree (Figure 5). One cluster consists of Haloferax and Haladaptatus, together with Haloarcula, Halobacterium, Halorhabdus, Natronomonas, Haloquadratum, Natrialba and Halorubrum. A second cluster is comprised of Halogeometricum, Haloterrigena, Halomicrobium and Halalkalicoccus.
The discovery of the conflicting groupings of haloarchaea in the LeuRS and the tRNALeu phylogenies begs the question of the evolution regarding LeuRS-tRNALeu metabolic interaction in these organisms. Our results suggest that the evolutionary route that the haloarchaeal tRNALeu took was independent of the evolution of the aaRS that aminoacylates it. This implies that the LeuRS and tRNALeu can be horizontally acquired independently, and one does not seem to strongly restrict the evolution of the other. tRNAs are often involved in HGT, with many found in close proximity to mobile elements and genomic islands . The lack of co-evolution we find for tRNALeu and LeuRS is in contrast to the finding that human but not E. coli TyrRS could complement yeast whose TyrRS gene had been disrupted . However, this reported "species specificity" was found to be due to a small peptide element in TyrRS, whose modification allowed the switching of species-specific aminoacylation across taxonomic domains .
The horizontal acquisition of aaRS of the same specificity might reflect a stochastic event in the evolution of these ancient enzymes. Numerous HGT events have been reported in many aaRS of different amino acid specificity, and these involved transfers at different taxonomic levels [18–20]. If these enzymes have been undergoing horizontal transfers in many extant lineages without affecting the evolution of their cognate tRNA, we cannot exclude the possibility that these transfers occurred without any impact to their aminoacylation capacities. Hence, the frequent transfers and current distribution of aaRS may instead reflect neutral stochastic transfers  and replacements. On the other hand, different aaRS forms in some instances were shown to provide differential sensitivity to naturally occurring antibiotics (see discussion in ). The possibility of selection through antibiotic resistance is seen in duplicate forms of same-specificity aaRS in Bacteria [47–49], and was suggested as a possible driving force behind the replacement of aaRS homeoalleles . However, this hypothesis still requires further investigation.
Conclusions: the impact of ancient HGTs on present-day lineages
Shared ancestry alone cannot explain the incredible variety in the genetic material that lineages possess. With HGT, organismal evolution becomes a patchwork of genes from varied sources. As a lineage evolves through time, it receives and loses genes and gene fragments, not only from its close relatives through biased gene transfer [14, 41] but also from the mobilome , from distantly related taxa, and from organisms that existed alongside a particular extant lineage and that have now become extinct. Ancient lineages are an important source of genetic diversity in extant taxa. Through vertical inheritance alone, this molecular diversity would likely take millions of years to achieve. When transfers occur from deep branching lineages and the recipient passes the horizontally acquired genetic material to its descendants, the ancient genes are maintained in the genomes or pan-genomes of existing organisms even if the original donors went extinct since the transfer occurred.
Highly divergent genes that have patchy distributions in extant lineages provide strong evidence for ancient HGTs, as observed not only in the haloarchaeal LeuRS but also reported in PylRS , SepRS , and the rare forms of Ser  and ThrRS . In these cases, the genes still exhibit some degree of similarity with their homologs, allowing reliable phylogenetic reconstruction. More challenging are genes that have no recognizable homolog in other existing lineages as is the case for the thousands of intriguing orphan genes (or ORFans ) and gene families in extant genomes, whose evolutionary histories remain unaccounted for because they are present in only a small group of closely related organisms . Assuming that most lineages that ever existed are now extinct [54, 55], it is remarkable that evolution preserves some of these genes as molecular “fossils” . However, the relative contributions from fast evolving phages and other components of the mobilome , unsampled or extinct lineages, and gene creation from previously non-coding DNA  remain to be ascertained.
The horizontal transmission of bacterial LeuRS to the Halobacteriales provides evidence for prokaryotic lineages that existed in the distant past and for their position in the Tree/Net of Life. The extremely long branch that in most molecular phylogenies separates LUCA from the base of the bacterial domain may have been populated by lineages that existed in the past. The results of the LeuRS analyses may be interpreted as evidence for the fourth domain of life, as reported in . At present, we can only deduce few and tentative characteristics of these ancient, deep-branching lineages; however, future work may identify other donations made by these ancient lineages, possibly leading to a better characterization of these long extinct cousins of modern bacteria.
Protein sequences of LeuRS from the three domains were retrieved by BLASTP searches of the non-redundant protein database and the BLAST microbial genome database from the National Center for Biotechnology Information (NCBI) website . For the global phylogenetic analysis, 325 LeuRS sequences were used. Sequences were aligned using the MUSCLE algorithm  with default parameters. Maximum likelihood phylogenetic reconstruction of the LeuRS sequences was performed using PhyML v3.0  with 100 bootstrap replicates, WAG  substitution model, estimated portions of invariable sites, four substitution-rate categories, estimated Γ distribution parameter, estimated amino acid frequencies, and NJ starting tree. Maximum likelihood distances were calculated using the programs PUZZLEBOOT v1.03  and TREE-PUZZLE  using the WAG  + Γ + I model to calculate pairwise maximum likelihood distances and NEIGHBOR  for tree reconstruction. Posterior probability values were generated using MrBayes v. 3.1.2 , with a fixed WAG  amino acid substitution model using four rate categories approximating a Γ distribution, four chains and a random starting tree. We used a specified number of generations for each aaRS analyses (145,000 for the haloarchaeal LeuRS type A and 150,000 for the haloarchaeal LeuRS type B) sampling every 100th generation. The first 25% of the sampled generations were removed from the analysis as burn-in. Inspection of the convergence parameter and log likelihood values reveals that the exploration of the tree space has reached a plateau.
For comparison, we used a multilocus sequence analysis (MLSA) approach that comprises five housekeeping genes that encode for V-type ATP synthase subunit B (AtpB), elongation factor 2 (EF-2), DNA repair and recombination protein (RadA), RNA polymerase subunit B’ (RpoB’) and preprotein translocase subunit (SecY) . The sequence data for each protein were partitioned using MrBayes v.3.1.2  and phylogenetic estimates were calculated from the different data partitions. Tree reconstruction and calculation of posterior probability values were generated using MrBayes v.3.1.2 .
DNA sequences encoding the archaeal tRNALeu were obtained by BLASTN searches. tRNA tree reconstruction and bootstrapping were performed using PhyML v3.0  with estimated portions of invariable sites, four substitution-rate categories, estimated ts/tv ratio, estimated Γ distribution parameter, estimated amino acid frequencies, BioNJ starting tree, 100 bootstrap replicates and GTR  nucleotide substitution model. Posterior probability values for the tRNALeu tree were generated using MrBayes v. 3.1.2 , with a fixed GTR nucleotide substitution model using four rate categories approximating a Γ distribution, four chains, a random starting tree, 50,000,000 generations sampling every 100th generation. The first 25% of the sampled generations were removed from the analysis as burn-in. Branch lengths and topologies of all phylograms were calculated with PhyML v3.0 . Inspection of the convergence parameter and log likelihood values reveals that the exploration of the tree space has reached a plateau. The substitution models used for each approach were determined using ProtTest  and jModelTest .
For GARD analyses, a smaller dataset was selected, containing 13 bacterial, 9 haloarchaeal, and 14 sequences from other archaea (Additional file 2: Figure S2). Analyses were performed using GARD as implemented on the datamonkey . The selection of the most appropriate substitution model (BLOSSUM62) was performed using the model selection program provided on the server. Using an alignment with MUSCLE as starting point, the sequences were realigned using SATé 2.03  selecting MUSCLE for merger and ProGammaIBLOSSUM62 as substitution model. Protein structure files were downloaded from the RCSB Protein Data Bank  and visualized using the Swiss PDB viewer . Positions in the Thermus thermophilus and Pyrococcus horikoshii structures corresponding to the identified breakpoint were identified using the multiple sequence alignment.
Genomic synteny among several members of the Halobacteriales and other Archaea was analyzed to identify the genes surrounding the leuS gene. This was done by aligning the genomes using the Integrated Microbial Genomes software tool provided by the U.S. Department of Energy Joint Genome Institute (http://img.jgi.doe.gov/cgi-bin/w/main.ci).
The authors would like to thank Pascal Lapierre and the Biotechnology Bioservices Center of the University of Connecticut for technical support, and Gregory P. Fournier (MIT) for stimulating discussions and suggestions. This work was supported by the National Science Foundation Grant (DEB 0830024) to JPG. RTP was supported by the National Science Foundation (DEB 0919290 and 0830024) and the U.S.Israel Binational Science Foundation (award number 2007043). The structure depicted on the pdf-cover is from .
- Oren A: The order Halobacteriales. Prokaryotes. Edited by: Dworkin M, Falkow S, Rosenberg E, Schleifer K, Stackebrandt E. 2006, Springer Science+Business Media, New York, 113-164. 2006View ArticleGoogle Scholar
- Cuadros-Orellana S, Martin-Cuadrado AB, Legault B, D'Auria G, Zhaxybayeva O, Papke RT, Rodriguez-Valera F: Genomic plasticity in prokaryotes: the case of the square haloarchaeon. ISME J. 2007, 1 (3): 235-245.PubMedView ArticleGoogle Scholar
- Papke RT, Koenig JE, Rodriguez-Valera F, Doolittle WF: Frequent recombination in a saltern population of Halorubrum. Science. 2004, 306 (5703): 1928-1929.PubMedGoogle Scholar
- Rhodes M, Spear J, Oren A, House C: Differences in lateral gene transfer in hypersaline versus thermal environments. BMC Evol Biol. 2011, 11 (1): 199.PubMedPubMed CentralView ArticleGoogle Scholar
- Khomyakova M, Bukmez O, Thomas LK, Erb TJ, Berg IA: A methylaspartate cycle in haloarchaea. Science. 2011, 331 (6015): 334-337.PubMedView ArticleGoogle Scholar
- Sharma AK, Walsh DA, Bapteste E, Rodriguez-Valera F, Ford Doolittle W, Papke RT: Evolution of rhodopsin ion pumps in haloarchaea. BMC Evol Biol. 2007, 7: 79.PubMedPubMed CentralView ArticleGoogle Scholar
- Brochier-Armanet C, Forterre P, Gribaldo S: Phylogeny and evolution of the Archaea: one hundred genomes later. Curr Opin Microbiol. 2011, 14 (3): 274-281.PubMedView ArticleGoogle Scholar
- Woese CR: Bacterial evolution. Microbiol Rev. 1987, 51 (2): 221-271.PubMedPubMed CentralGoogle Scholar
- Wolf YI, Rogozin IB, Grishin NV, Koonin EV: Genome trees and the tree of life. Trends Genet. 2002, 18 (9): 472-479.PubMedView ArticleGoogle Scholar
- Korbel JO, Snel B, Huynen MA, Bork P: SHOT: a web server for the construction of genome phylogenies. Trends Genet. 2002, 18 (3): 158-162.PubMedView ArticleGoogle Scholar
- Srinivasan G, James CM, Krzycki JA: Pyrrolysine encoded by UAG in Archaea: charging of a UAG-decoding specialized tRNA. Science. 2002, 296 (5572): 1459-1462.PubMedView ArticleGoogle Scholar
- Fournier GP, Huang J, Gogarten JP: Horizontal gene transfer from extinct and extant lineages: biological innovation and the coral of life. Philos Transac R Soc Lond B Biol Sci. 2009, 364: 2229-2239.View ArticleGoogle Scholar
- Kavran JM, Gundllapalli S, O'Donoghue P, Englert M, Soll D, Steitz TA: Structure of pyrrolysyl-tRNA synthetase, an archaeal enzyme for genetic code innovation. Proc Natl Acad Sci USA. 2007, 104 (27): 11268-11273.PubMedPubMed CentralView ArticleGoogle Scholar
- Andam CP, Gogarten JP: Biased gene transfer in microbial evolution. Nat Rev Microbiol. 2011, 9 (7): 543-555.PubMedView ArticleGoogle Scholar
- Sauerwald A, Zhu W, Major TA, Roy H, Palioura S, Jahn D, Whitman WB, Yates JR, Ibba M, Soll D: RNA-dependent cysteine biosynthesis in archaea. Science. 2005, 307 (5717): 1969-1972.PubMedView ArticleGoogle Scholar
- O'Donoghue P, Sethi A, Woese CR, Luthey-Schulten ZA: The evolutionary history of Cys-tRNACys formation. Proc Natl Acad Sci USA. 2005, 102 (52): 19003-19008.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhaxybayeva O, Gogarten JP: Cladogenesis, coalescence and the evolution of the three domains of life. Trends Genet. 2004, 20 (4): 182-187.PubMedView ArticleGoogle Scholar
- Woese CR, Olsen GJ, Ibba M, Soll D: Aminoacyl-tRNA synthetases, the genetic code, and the evolutionary process. Microbiol Mol Biol Rev. 2000, 64: 202-236.PubMedPubMed CentralView ArticleGoogle Scholar
- Wolf YI, Aravind L, Grishin NV, Koonin EV: Evolution of aminoacyl-tRNA synthetases–analysis of unique domain architectures and phylogenetic trees reveals a complex history of horizontal gene transfer events. Genome Res. 1999, 9 (8): 689-710.PubMedGoogle Scholar
- Luque I, Riera-Alberola ML, Andujar A, de Ochoa Alda JA: Intraphylum diversity and complex evolution of cyanobacterial aminoacyl-tRNA synthetases. Mol Biol Evol. 2008, 25 (11): 2369-2389.PubMedView ArticleGoogle Scholar
- Brochier C, Forterre P, Gribaldo S: An emerging phylogenetic core of Archaea: phylogenies of transcription and translation machineries converge following addition of new genome sequences. BMC Evol Biol. 2005, 5 (1): 36.PubMedPubMed CentralView ArticleGoogle Scholar
- Brochier-Armanet C, Boussau B, Gribaldo S, Forterre P: Mesophilic crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota. Nat Rev Microbiol. 2008, 6: 245-252.PubMedView ArticleGoogle Scholar
- Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SDW: Automated phylogenetic detection of recombination using a genetic algorithm. Mol Biol Evol. 2006, 23 (10): 1891-1901.PubMedView ArticleGoogle Scholar
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucl Acids Res. 2004, 32 (5): 1792-1797.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu K, Warnow TJ, Holder MT, Nelesen SM, Yu J, Stamatakis AP, Linder CR: SATe-II: very fast and accurate simultaneous estimation of multiple sequence alignments and phylogenetic trees. Syst Biol. 2012, 61 (1): 90-106.PubMedView ArticleGoogle Scholar
- Stern A, Mayrose I, Penn O, Shaul S, Gophna U, Pupko T: An evolutionary analysis of lateral gene transfer in thymidylate synthase enzymes. Syst Biol. 2010, 59 (2): 212-225.PubMedPubMed CentralView ArticleGoogle Scholar
- Kondrashov FA, Kondrashov AS: Role of selection in fixation of gene duplications. J Theor Biol. 2006, 239 (2): 141-151.PubMedView ArticleGoogle Scholar
- O'Connor E, Shand R: Halocins and sulfolobicins: the emerging story of archaeal protein and peptide antibiotics. J Ind Microbiol Biotech. 2001, 28: 23-31.View ArticleGoogle Scholar
- Dridi B, Fardeau ML, Ollivier B, Raoult D, Drancourt M: The antimicrobial resistance pattern of cultured human methanogens reflects the unique phylogenetic position of archaea. J Antimicrob Chemother. 2011, 66 (9): 2038-2044.PubMedView ArticleGoogle Scholar
- Lopez-Lopez A, Benlloch S, Bonfa M, Rodriguez-Valera F, Mira A: Intragenomic 16S rDNA divergence in Haloarcula marismortui is an adaptation to different temperatures. J Mol Evol. 2007, 65 (6): 687-696.PubMedView ArticleGoogle Scholar
- Gray MW, Lukes J, Archibald JM, Keeling PJ, Doolittle WF: Irremediable Complexity?. Science. 2010, 330 (6006): 920-921.PubMedView ArticleGoogle Scholar
- Xu MG, Chen JF, Martin F, Zhao MW, Eriani G, Wang ED: Leucyl-tRNA synthetase consisting of two subunits from hyperthermophilic bacteria Aquifex aeolicus. J Biol Chem. 2002, 277 (44): 41590-41596.PubMedView ArticleGoogle Scholar
- Olieric N, Bey G, Nierengarten H, Wang ED, Moras D, Eriani G, Cavarelli J: Expression, purification, and characterization of a new heterotetramer structure of leucyl-tRNA synthetase from Aquifex aeolicus in Escherichia coli. Protein Expr Purif. 2006, 47 (1): 1-9.PubMedView ArticleGoogle Scholar
- Stoltzfus A: On the possibility of constructive neutral evolution. J Mol Evol. 1999, 49 (2): 169-181.PubMedView ArticleGoogle Scholar
- Boucher Y, Douady CJ, Sharma AK, Kamekura M, Doolittle WF: Intragenomic heterogeneity and intergenomic recombination among haloarchaeal rRNA genes. J Bacteriol. 2004, 186 (12): 3980-3990.PubMedPubMed CentralView ArticleGoogle Scholar
- Papke RT, Zhaxybayeva O, Feil EJ, Sommerfeld K, Muise D, Doolittle WF: Searching for species in haloarchaea. Proc Natl Acad Sci USA. 2007, 104 (35): 14092-14097.PubMedPubMed CentralView ArticleGoogle Scholar
- Walsh D, Bapteste E, Kamekura M, Doolittle W: Evolution of the RNA polymerase B' subunit gene (rpoB') in Halobacteriales: a complementary molecular marker to the SSU rRNA gene. Mol Biol Evol. 2004, 21: 2340-2351.PubMedView ArticleGoogle Scholar
- Minegishi H, Kamekura M, Itoh T, Echigo A, Usami R, Hasimoto T: Further refinement of the phylogeny of the Halobacteriaceae based on the full-length RNA polymerase subunit B' (rpoB')gene. Int J Syst Evol Microbiol. 2010, 60: 2398-2408.PubMedView ArticleGoogle Scholar
- Papke R, White E, Reddy P, Weigel G, Kamekura M, Minegishi H, Usami R, Ventosa A: A multilocus sequence analysis (MLSA) approach to Halobacteriales phylogeny and taxonomy. Int J Syst Evol Microbiol. 2011, 61 (12): 2984-2995.PubMedView ArticleGoogle Scholar
- Anderson I, Scheuner C, Goker M, Mavromatis K, Hooper SD, Porat I, Klenk HP, Ivanova N, Kyrpides N: Novel insights into the diversity of catabolic metabolism from ten haloarchaeal genomes. PLoS One. 2011, 6 (5): e20237.PubMedPubMed CentralView ArticleGoogle Scholar
- Andam CP, Williams D, Gogarten JP: Biased gene transfer mimics patterns created through shared ancestry. Proc Natl Acad Sci USA. 2010, 107: 10679-10684.PubMedPubMed CentralView ArticleGoogle Scholar
- Widmann J, Harris J, Lozupone C, Wolfson A, Knight R: Stable tRNA-based phylogenies using only 76 nucleotides. RNA. 2010, 16: 1469-1477.PubMedPubMed CentralView ArticleGoogle Scholar
- Ou H, Chen L, Lonnen J, Chaudhuri R, Thani A, Smith R, Garton N, Hinton J, Pallen M, Barer M, et al: A novel strategy for the identification of genomic islands by comparative analysis of the contents and contexts of tRNA sites in closely related bacteria. Nucleic Acids Res. 2006, 34: e3.PubMedPubMed CentralView ArticleGoogle Scholar
- Wakasugi K, Quinn CL, Tao N, Schimmel P: Genetic code in evolution: switching species-specific aminoacylation with a peptide transplant. EMBO J. 1998, 17 (1): 297-305.PubMedPubMed CentralView ArticleGoogle Scholar
- Gogarten JP, Townsend JP: Horizontal gene transfer, genome innovation and evolution. Nat Rev Microbiol. 2005, 3 (9): 679-687.PubMedView ArticleGoogle Scholar
- Andam CP, Fournier GP, Gogarten JP: Multilevel populations and the evolution of antibiotic resistance through horizontal gene transfer. FEMS Microbiol Rev. 2011, 35 (5): 756-767.PubMedView ArticleGoogle Scholar
- Gentry DR, Ingraham KA, Stanhope MJ, Rittenhouse S, Jarvest RL, O'Hanlon PJ, Brown JR, Holmes DJ: Variable sensitivity to bacterial methionyl-tRNA synthetase inhibitors reveals subpopulations of Streptococcus pneumoniae with two distinct methionyl-tRNA synthetase genes. Antimicrob Agents Chemother. 2003, 47 (6): 1784-1789.PubMedPubMed CentralView ArticleGoogle Scholar
- Brown JR, Gentry D, Becker JA, Ingraham K, Holmes DJ, Stanhope MJ: Horizontal transfer of drug-resistant aminoacyl-transfer-RNA synthetases of anthrax and Gram-positive pathogens. EMBO Rep. 2003, 4 (7): 692-698.PubMedPubMed CentralView ArticleGoogle Scholar
- Vecchione JJ, Sello JK: A novel tryptophanyl-tRNA synthetase gene confers high-level resistance to indolmycin. Antimicrob Agents Chemother. 2009, 53 (9): 3972-3980.PubMedPubMed CentralView ArticleGoogle Scholar
- Frost LS, Leplae R, Summers AO, Toussaint A: Mobile genetic elements: the agents of open source evolution. Nat Rev Microbiol. 2005, 3 (9): 722-732.PubMedView ArticleGoogle Scholar
- Andam CP, Gogarten JP: Biased gene transfer and its implications for the concept of lineage. Biol Direct. 2011, 6 (1): 47.PubMedPubMed CentralView ArticleGoogle Scholar
- Siew N, Fischer D: Structural biology sheds light on the puzzle of genomic ORFans. J Mol Biol. 2004, 342: 369-373.PubMedView ArticleGoogle Scholar
- Lapierre P, Gogarten JP: Estimating the size of the bacterial pan-genome. Trends Genet. 2009, 25 (3): 107-110.PubMedView ArticleGoogle Scholar
- Gogarten JP, Fournier GP, Zhaxybayeva O: Gene transfer and the reconstruction of life’s early history from genomic data. Space Science Reviews. 2008, 135 (1): 115-131.View ArticleGoogle Scholar
- Fournier GP, Huang J, Gogarten JP: Horizontal gene transfer from extinct and extant lineages: biological innovation and the coral of life. Philos Trans R Soc Lond B Biol Sci. 2009, 364 (1527): 2229-2239.PubMedPubMed CentralView ArticleGoogle Scholar
- Daubin V, Ochman H: Bacterial genomes as new gene homes: the genealogy of ORFans in E. coli. Genome Res. 2004, 14 (6): 1036-1042.PubMedPubMed CentralView ArticleGoogle Scholar
- Tautz D, Domazet-Loso T: The evolutionary origin of orphan genes. Nat Rev Genet. 2011, 12 (10): 692-702.PubMedView ArticleGoogle Scholar
- Boyer M, Madoui MA, Gimenez G, La Scola B, Raoult D: Phylogenetic and phyletic studies of informational genes in genomes highlight existence of a 4 domain of life including giant viruses. PLoS One. 2010, 5 (12): e15530.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView ArticleGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704.PubMedView ArticleGoogle Scholar
- Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol. 2001, 18 (5): 691-699.PubMedView ArticleGoogle Scholar
- Holder ME, Roger A: PUZZLEBOOT. 1999, Marine Biological Laboratory, Woods HoleGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504.PubMedView ArticleGoogle Scholar
- Felsenstein J: PHYLIP Manual v1.03. 1993, University Herbarium, University of California, BerkeleyGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574.PubMedView ArticleGoogle Scholar
- Lanave C, Preparata G, Saccone C, Serio G: A new method for calculating evolutionary substitution rates. J Mol Evol. 1984, 20 (1): 86-93.PubMedView ArticleGoogle Scholar
- Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105.PubMedView ArticleGoogle Scholar
- Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256.PubMedView ArticleGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28 (1): 235-242.PubMedPubMed CentralView ArticleGoogle Scholar
- Guex N, Peitsch MC: SWISS-MODEL and the Swiss-PdbViewer: an environment for comparative protein modeling. Electrophoresis. 1997, 18 (15): 2714-2723.PubMedView ArticleGoogle Scholar
- Fukunaga R, Yokoyama S: The crystal structure of Leucyl-tRNA synthetase and tRNA(leucine) complex PDB ID: 1wz2. http://dx.doi.org/10.2210/pdb1wz2/pdbIn RCSB Protein Data Bank.
- Tukalo M, Yaremchuk A, Fukunaga R, Yokoyama S, Cusack S: The crystal structure of leucyl-tRNA synthetase complexed with tRNALeu in the post-transfer-editing conformation. Nat Struct Mol Biol. 2005, 12 (10): 923-930.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.