Adaptation to environmental temperature is a major determinant of molecular evolutionary rates in archaea

Methods to infer the ancestral conditions of life are commonly based on geological and paleontological analyses. Recently, several studies used genome sequences to gain information about past ecological conditions taking advantage of the property that the G+C and amino acid contents of bacterial and archaeal ribosomal DNA genes and proteins, respectively, are strongly influenced by the environmental temperature. The adaptation to optimal growth temperature (OGT) since the Last Universal Common Ancestor (LUCA) over the universal tree of life was examined, and it was concluded that LUCA was likely to have been a mesophilic organism and that a parallel adaptation to high temperature occurred independently along the two lineages leading to the ancestors of Bacteria on one side and of Archaea and Eukarya on the other side. Here, we focus on Archaea to gain a precise view of the adaptation to OGT over time in this domain. It has been often proposed on the basis of indirect evidence that the last archaeal common ancestor was a hyperthermophilic organism. Moreover, many results showed the influence of environmental temperature on the evolutionary dynamics of archaeal genomes: Thermophilic organisms generally display lower evolutionary rates than mesophiles. However, to our knowledge, no study tried to explain the differences of evolutionary rates for the entire archaeal domain and to investigate the evolution of substitution rates over time. A comprehensive archaeal phylogeny and a non homogeneous model of the molecular evolutionary process allowed us to estimate ancestral base and amino acid compositions and OGTs at each internal node of the archaeal phylogenetic tree. The last archaeal common ancestor is predicted to have been hyperthermophilic and adaptations to cooler environments can be observed for extant mesophilic species. Furthermore, mesophilic species present both long branches and high variation of nucleotide and amino acid compositions since the last archaeal common ancestor. The increase of substitution rates observed in mesophilic lineages along all their branches can be interpreted as an ongoing adaptation to colder temperatures and to new metabolisms. We conclude that environmental temperature is a major factor that governs evolutionary rates in Archaea.     

Widespread distribution of archaeal reverse gyrase in thermophilic bacteria suggests a complex history of vertical inheritance and lateral gene transfers

Reverse gyrase, an enzyme of uncertain funtion, is present in all hyperthermophilic archaea and bacteria. Previous phylogenetic studies have suggested that the gene for reverse gyrase has an archaeal origin and was transferred laterally (LGT) to the ancestors of the two bacterial hyperthermophilic phyla, Thermotogales and Aquificales. Here, we performed an in-depth analysis of the evolutionary history of reverse gyrase in light of genomic progress. We found genes coding for reverse gyrase in the genomes of several thermophilic bacteria that belong to phyla other than Aquificales and Thermotogales. Several of these bacteria are not, strictly speaking, hyperthermophiles because their reported optimal growth temperatures are below 80 degrees C. Furthermore, we detected a reverse gyrase gene in the sequence of the large plasmid of Thermus thermophilus strain HB8, suggesting a possible mechanism of transfer to the T. thermophilus strain HB8 involving plasmids and transposases. The archaeal part of the reverse gyrase tree is congruent with recent phylogenies of the archaeal domain based on ribosomal proteins or RNA polymerase subunits. Although poorly resolved, the complete reverse gyrase phylogeny suggests an ancient acquisition of the gene by bacteria via one or two LGT events, followed by its secondary distribution by LGT within bacteria. Finally, several genes of archaeal origin located in proximity to the reverse gyrase gene in bacterial genomes have bacterial homologues mostly in thermophiles or hyperthermophiles, raising the possibility that they were co-transferred with the reverse gyrase gene. Our new analysis of the reverse gyrase history strengthens the hypothesis that the acquisition of reverse gyrase may have been a crucial evolutionary step in the adaptation of bacteria to high-temperature environments. However, it also questions the role of this enzyme in thermophilic bacteria and the selective advantage its presence could provide.     

Accounting for evolutionary rate variation among sequence sites consistently changes universal phylogenies deduced from rRNA and protein-coding genes.

Phylogenetic analyses of gene and protein sequences have led to two major competing views of the universal phylogeny, the evolutionary tree relating the three kinds of living organisms, Bacteria, Archaea, and Eukarya. In the first scheme, called "the archaebacterial tree, " organisms of the same type are clustered together. In the second scenario, called "the eocyte tree," the archaeal phylum of Crenarchaeota is more closely related to eukaryotes than are other Archaea. A major property of the evolution of functional ribosomal and protein-encoding genes is that the rate of nucleotide and amino acid substitution varies across sequence sites. Here, using distance-based and maximum-likelihood methods, we show that universal phylogenies of ribosomal RNAs and RNA polymerases built by ignoring this variation are biased toward the archaebacterial tree because of attraction between long branches. In contrast, taking among-site rate variability into account gives support for the eocyte tree.     

The Evolutionary History of the Structure of 5S Ribosomal RNA

5S rRNA is the smallest nucleic acid component of the large ribosomal subunit, contributing to ribosomal assembly, stability, and function. Despite being a model for the study of RNA structure and RNA–protein interactions, the evolution of this universally conserved molecule remains unclear. Here, we explore the history of the three-domain structure of 5S rRNA using phylogenetic trees that are reconstructed directly from molecular structure. A total of 46 structural characters describing the geometry of 666 5S rRNAs were used to derive intrinsically rooted trees of molecules and molecular substructures. Trees of molecules revealed the tripartite nature of life. In these trees, superkingdom Archaea formed a paraphyletic basal group, while Bacteria and Eukarya were monophyletic and derived. Trees of molecular substructures supported an origin of the molecule in a segment that is homologous to helix I (α domain), its initial enhancement with helix III (β domain), and the early formation of the three-domain structure typical of modern 5S rRNA in Archaea. The delayed formation of the branched structure in Bacteria and Eukarya lends further support to the archaeal rooting of the tree of life. Remarkably, the evolution of molecular interactions between 5S rRNA and associated ribosomal proteins inferred from a census of domain structure in hundreds of genomes established a tight relationship between the age of 5S rRNA helices and the age of ribosomal proteins. Results suggest 5S rRNA originated relatively quickly but quite late in evolution, at a time when primordial metabolic enzymes and translation machinery were already in place. The molecule therefore represents a late evolutionary addition to the ribosomal ensemble that occurred prior to the early diversification of Archaea.     

Phylogenetic diversity and ecology of environmental Archaea

On the basis of culture studies, Archaea were thought to be synonymous with extreme environments. However, the large numbers of environmental rRNA gene sequences currently flooding into databases such as GenBank show that these organisms are present in almost all environments examined to date. Large sequence databases and new fast phylogenetic software allow more precise determination of the archaeal phylogenetic tree, but also indicate that our knowledge of archaeal diversity is incomplete. Although it is apparent that Archaea can be found in all environments, the chemistry of their ecological context is mostly unknown.     

Phylogenetic analyses of two "archaeal" genes in thermotoga maritima reveal multiple transfers between archaea and bacteria

The genome sequence of Thermotoga maritima revealed that 24% of its open reading frames (ORFs) showed the highest similarity scores to archaeal genes in BLAST analyses. Here we screened 16 strains from the genus Thermotoga and other related Thermotogales for the occurrence of two of these "archaeal" genes: the gene encoding the large subunit of glutamate synthase (gltB) and the myo-inositol 1P synthase gene (ino1). Both genes were restricted to the Thermotoga species within the Thermotogales. The distribution of the two genes, along with results from phylogenetic analyses, showed that they were acquired from Archaea during the divergence of the Thermotogales. Database searches revealed that three other bacteria-Dehalococcoides ethenogenes, Sinorhizobium meliloti, and Clostridium difficile-possess archaeal-type gltBs, and the phylogenetic analyses confirmed at least two lateral gene transfer (LGT) events between Bacteria and Archaea. These LGT events were also strongly supported by gene structure data, as the three domains in bacterial-type gltB are homologous to three independent ORFs in Archaea and Bacteria with archaeal-type gltBs. The ino1 gene has a scattered distribution among Bacteria, and apart from the Thermotoga strains it is found only in Aquifex aeolicus, D. ethenogenes, and some high-G+C Gram-positive bacteria. Phylogenetic analysis of the ino1 sequences revealed three highly supported prokaryotic clades, all containing a mixture of archaeal and bacterial sequences, and suggested that all bacterial ino1 genes had been recruited from archaeal donors. The Thermotoga strains and A. aeolicus acquired this gene independently from different archaeal species. Although transfer of genes from hyperthermophilic Archaea may have facilitated the evolution of bacterial hyperthermophily, between-domain transfers also affect mesophilic species. For hyperthermophiles, we hypothesize that LGT may be as much a consequence as the cause of adaptation to hyperthermophily.     

Phylogenetic Analyses of Archaeal Ribosomal DNA Sequences from Salt Pan Sediment of Mumbai,India

The archaeal diversity in salt pan sediment from Mumbai, India, was investigated by 16S rDNA-dependent molecular phylogeny. Small-subunit rRNA (16S rDNA) from salt pan sediment metagenome were amplified by polymerase chain reaction (PCR) using primers specific to the domain archaea. Thirty-two unique phylotypes were obtained by PCR-based RFLP of 16S rRNA genes using endonucleases Hae111 and Msp1, which were most suitable to score the genetic diversity. These phylotypes spanned a wide range within the domain Archaea including both Crenarchaeota and Euryarcheaota. However, none of the retrieved Crenarchaeota sequences could be grouped with previously cultured Crenarchaeota. Of all the Euryarcheaota sequences, three sequences were related to Haloarchaea, two were related to cultured Methanosarcina sp., and the remaining sequences were affiliated with uncultured Methanosarcina sp., Methanogenium sp., and Methanolobus sp. Most of the sequences determined were closely related to the sequences that had been previously obtained from metagenome of a variety of marine environments. The phylogenetic study of a site investigated for the first time revealed the presence of a highly diverse archaeal population and may represent novel sequences and organisms specially adapted to the salt pan sediment of Mumbai. These findings are of fundamental value for understanding the complexity of salt pan ecosystems.     

New Findings on Evolution of Metal Homeostasis Genes: Evidence from Comparative Genome Analysis of Bacteria and Archaea

In order to examine the natural history of metal homeostasis genes in prokaryotes, open reading frames with homology to characterized P_IB-type ATPases from the genomes of 188 bacteria and 22 archaea were investigated. Major findings were as follows. First, a high diversity in N-terminal metal binding motifs was observed. These motifs were distributed throughout bacterial and archaeal lineages, suggesting multiple loss and acquisition events. Second, the CopA locus separated into two distinct phylogenetic clusters, CopA1, which contained ATPases with documented Cu(I) influx activity, and CopA2, which contained both efflux and influx transporters and spanned the entire diversity of the bacterial domain, suggesting that CopA2 is the ancestral locus. Finally, phylogentic incongruences between 16S rRNA and P_IB-type ATPase gene trees identified at least 14 instances of lateral gene transfer (LGT) that had occurred among diverse microbes. Results from bootstrapped supported nodes indicated that (i) a majority of the transfers occurred among proteobacteria, most likely due to the phylogenetic relatedness of these organisms, and (ii) gram-positive bacteria with low moles percent G+C were often involved in instances of LGT. These results, together with our earlier work on the occurrence of LGT in subsurface bacteria (J. M. Coombs and T. Barkay, Appl. Environ. Microbiol. 70:1698-1707, 2004), indicate that LGT has had a minor role in the evolution of P_IB-type ATPases, unlike other genes that specify survival in metal-stressed environments. This study demonstrates how examination of a specific locus across microbial genomes can contribute to the understanding of phenotypes that are critical to the interactions of microbes with their environment.     

Comparative analysis of ribosomal proteins in complete genomes: an example of reductive evolution at the domain scale

A comprehensive investigation of ribosomal genes in complete genomes from 66 different species allows us to address the distribution of r-proteins between and within the three primary domains. Thirty-four r-protein families are represented in all domains but 33 families are specific to Archaea and Eucarya, providing evidence for specialisation at an early stage of evolution between the bacterial lineage and the lineage leading to Archaea and Eukaryotes. With only one specific r-protein, the archaeal ribosome appears to be a small-scale model of the eukaryotic one in terms of protein composition. However, the mechanism of evolution of the protein component of the ribosome appears dramatically different in Archaea. In Bacteria and Eucarya, a restricted number of ribosomal genes can be lost with a bias toward losses in intracellular pathogens. In Archaea, losses implicate 15% of the ribosomal genes revealing an unexpected plasticity of the translation apparatus and the pattern of gene losses indicates a progressive elimination of ribosomal genes in the course of archaeal evolution. This first documented case of reductive evolution at the domain scale provides a new framework for discussing the shape of the universal tree of life and the selective forces directing the evolution of prokaryotes.     

New findings on evolution of metal homeostasis genes: evidence from comparative genome analysis of bacteria and archaea

In order to examine the natural history of metal homeostasis genes in prokaryotes, open reading frames with homology to characterized P(IB)-type ATPases from the genomes of 188 bacteria and 22 archaea were investigated. Major findings were as follows. First, a high diversity in N-terminal metal binding motifs was observed. These motifs were distributed throughout bacterial and archaeal lineages, suggesting multiple loss and acquisition events. Second, the CopA locus separated into two distinct phylogenetic clusters, CopA1, which contained ATPases with documented Cu(I) influx activity, and CopA2, which contained both efflux and influx transporters and spanned the entire diversity of the bacterial domain, suggesting that CopA2 is the ancestral locus. Finally, phylogentic incongruences between 16S rRNA and P(IB)-type ATPase gene trees identified at least 14 instances of lateral gene transfer (LGT) that had occurred among diverse microbes. Results from bootstrapped supported nodes indicated that (i) a majority of the transfers occurred among proteobacteria, most likely due to the phylogenetic relatedness of these organisms, and (ii) gram-positive bacteria with low moles percent G+C were often involved in instances of LGT. These results, together with our earlier work on the occurrence of LGT in subsurface bacteria (J. M. Coombs and T. Barkay, Appl. Environ. Microbiol. 70:1698-1707, 2004), indicate that LGT has had a minor role in the evolution of P(IB)-type ATPases, unlike other genes that specify survival in metal-stressed environments. This study demonstrates how examination of a specific locus across microbial genomes can contribute to the understanding of phenotypes that are critical to the interactions of microbes with their environment.     

Vertical distribution and phylogenetic characterization of marine planktonic Archaea in the Santa Barbara Channel.

Newly described phylogenetic lineages within the domain Archaea have recently been found to be significant components of marine picoplankton assemblages. To better understand the ecology of these microorganisms, we investigated the relative abundance, distribution, and phylogenetic composition of Archaea in the Santa Barbara Channel. Significant amounts of archaeal rRNA and rDNA (genes coding for rRNA) were detected in all samples analyzed. The relative abundance of archaeal rRNA as measured by quantitative oligonucleotide hybridization experiments was low in surface waters but reached higher values (20 to 30% of prokaryotic rRNA) at depths below 100 m. Probes were developed for the two major groups of marine Archaea detected. rRNA originating from the euryarchaeal group (group II) was most abundant in surface waters, whereas rRNA from the crenarchaeal group (group I) dominated at depth. Clone libraries of PCR-amplified archaeal rRNA genes were constructed with samples from 0 and 200 m deep. Screening of libraries by hybridization with specific oligonucleotide probes, as well as subsequent sequencing of the cloned genes, indicated that virtually all archaeal rDNA clones recovered belonged to one of the two groups. The recovery of cloned rDNA sequence types in depth profiles exhibited the same trends as were observed in quantitative rRNA hybridization experiments. One representative of each of 18 distinct restriction fragment length polymorphism types was partially sequenced. Recovered sequences spanned most of the previously reported phylogenetic diversity detected in planktonic crenarchaeal and euryarchaeal groups. Several rDNA sequences appeared to be harbored in archaeal types which are widely distributed in marine coastal waters. In total, data suggest that marine planktonic crenarchaea and euryarchaea of temperate coastal habitats thrive in different zones of the water column. The relative rRNA abundance of the crenarchaeal group suggests that its members constitute a significant fraction of the prokaryotic biomass in subsurface coastal waters.     

Organisation of the S10, spc and alpha ribosomal protein gene clusters in prokaryotic genomes

Although it is well known that there is no long range colinearity in gene order in bacterial genomes, it is thought that there are several regions that are under strong structural constraints during evolution, in which gene order is extremely conserved. One such region is the str locus, containing the S10-spc-alpha operons. These operons contain genes coding for ribosomal proteins and for a number of housekeeping genes. We compared the organisation of these gene clusters in 111 sequenced prokaryotic genomes (99 bacterial and 12 archaeal genomes). We also compared the organisation to the phylogeny based on 16S ribosomal RNA gene sequences and the sequences of the ribosomal proteins L22, L16 and S14. Our data indicate that there is much variation in gene order and content in these gene clusters, both in bacterial as well as in archaeal genomes. Our data indicate that differential gene loss has occurred on multiple occasions during evolution. We also noted several discrepancies between phylogenetic trees based on 16S rRNA gene sequences and sequences of ribosomal proteins L16, L22 and S14, suggesting that horizontal gene transfer did play a significant role in the evolution of the S10-spc-alpha gene clusters.     

Phylogenetic analysis of Group I marine archaeal rRNA sequences emphasizes the hidden diversity within the primary group Archaea

Archaea form one of the three primary groups of extant life and are commonly associated with the extreme environments which many of their members inhabit. Currently, the Archaea are classified into two kingdoms, Crenarchaeota and Euryarchaeota, based on phylogenetic analysis of ribosomal RNA (rRNA) sequences. Molecular techniques allowing the retrieval and analysis of rRNA sequences from diverse environments are increasing our knowledge of archaeal diversity. This report describes the presence of marine Archaea in north-east Atlantic waters. Quantitative estimates indicated that the marine Archaea constitute 8 per cent of the total prokaryotic rRNA in Irish coastal waters. Phylogenetic analysis of the archaeal rRNA gene sequences revealed sufficient genetic diversity within Archaea to indicate that the current two-kingdom classification of Crenarchaeota and Euryarchaeota is restrictive.     

Electron acceptor-dependent identification of key anaerobic toluene degraders at a tar-oil-contaminated aquifer by Pyro-SIP

The diversity of archaeal communities growing in four hot springs (65-90 °C, pH 6.5) was assessed with 16S rRNA gene primers specific for the domain Archaea. Overall, mainly uncultured members of the Desulfurococcales, the Thermoproteales and the Korarchaeota, were identified. Based on this diversity, a set of chaperonin heat-shock protein (Hsp60) gene sequences from different archaeal species were aligned to design two degenerate primer sets for the amplification of the chaperonin gene: Ths and Kor (which can also detect the korarchaeotal chaperonin gene from one of the samples). A phylogenetic tree was constructed using the chaperonin sequences retrieved and other sequences from cultured representatives. The Alpha and Beta paralogs of the chaperonin gene were observed within the main clades and orthologs among them. Cultivated representatives from these clades were assigned to either paralog in the chaperonin tree. Uncultured representatives observed in the 16S rRNA gene analysis were found to be related to the Desulfurococcales. The topologies of the 16S rRNA gene and chaperonin phylogenetic trees were compared, and similar phylogenetic relationships were observed. Our results suggest that the chaperonin Hsp60 gene may be used as a phylogenetic marker for the clades found in this extreme environment.     

Archaeal diversity in waters from deep South African gold mines.

A culture-independent molecular analysis of archaeal communities in waters collected from deep South African gold mines was performed by performing a PCR-mediated terminal restriction fragment length polymorphism (T-RFLP) analysis of rRNA genes (rDNA) in conjunction with a sequencing analysis of archaeal rDNA clone libraries. The water samples used represented various environments, including deep fissure water, mine service water, and water from an overlying dolomite aquifer. T-RFLP analysis revealed that the ribotype distribution of archaea varied with the source of water. The archaeal communities in the deep gold mine environments exhibited great phylogenetic diversity; the majority of the members were most closely related to uncultivated species. Some archaeal rDNA clones obtained from mine service water and dolomite aquifer water samples were most closely related to environmental rDNA clones from surface soil (soil clones) and marine environments (marine group I [MGI]). Other clones exhibited intermediate phylogenetic affiliation between soil clones and MGI in the Crenarchaeota. Fissure water samples, derived from active or dormant geothermal environments, yielded archaeal sequences that exhibited novel phylogeny, including a novel lineage of Euryarchaeota. These results suggest that deep South African gold mines harbor novel archaeal communities distinct from those observed in other environments. Based on the phylogenetic analysis of archaeal strains and rDNA clones, including the newly discovered archaeal rDNA clones, the evolutionary relationship and the phylogenetic organization of the domain Archaea are reevaluated.     

A tree of life based on protein domain organizations

It is desirable to estimate a tree of life, a species tree including all available species in the 3 superkingdoms, Archaea, Bacteria, and Eukaryota, using not a limited number of genes but full-scale genome information. Here, we report a new method for constructing a tree of life based on protein domain organizations, that is, sequential order of domains in a protein, of all proteins detected in a genome of an organism. The new method is free from the identification of orthologous gene sets and therefore does not require the burdensome and error-prone computation. By pairwise comparisons of the repertoires of protein domain organizations of 17 archaeal, 136 bacterial, and 14 eukaryotic organisms, we computed evolutionary distances among them and constructed a tree of life. Our tree shows monophyly in Archaea, Bacteria, and Eukaryota and then monophyly in each of eukaryotic kingdoms and in most bacterial phyla. In addition, the branching pattern of the bacterial phyla in our tree is consistent with the widely accepted bacterial taxonomy and is very close to other genome-based trees. A couple of inconsistent aspects between the traditional trees and the genome-based trees including ours, however, would perhaps urge to revise the conventional view, particularly on the phylogenetic positions of hyperthermophiles.