Phylogenomics of prokaryotic ribosomal proteins

Archaeal and bacterial ribosomes contain more than 50 proteins, including 34 that are universally conserved in the three domains of cellular life (bacteria, archaea, and eukaryotes). Despite the high sequence conservation, annotation of ribosomal (r-) protein genes is often difficult because of their short lengths and biased sequence composition. We developed an automated computational pipeline for identification of r-protein genes and applied it to 995 completely sequenced bacterial and 87 archaeal genomes available in the RefSeq database. The pipeline employs curated seed alignments of r-proteins to run position-specific scoring matrix (PSSM)-based BLAST searches against six-frame genome translations, mitigating possible gene annotation errors. As a result of this analysis, we performed a census of prokaryotic r-protein complements, enumerated missing and paralogous r-proteins, and analyzed the distributions of ribosomal protein genes among chromosomal partitions. Phyletic patterns of bacterial and archaeal r-protein genes were mapped to phylogenetic trees reconstructed from concatenated alignments of r-proteins to reveal the history of likely multiple independent gains and losses. These alignments, available for download, can be used as search profiles to improve genome annotation of r-proteins and for further comparative genomics studies.     

Global Patterns of Protein Domain Gain and Loss in Superkingdoms

Domains are modules within proteins that can fold and function independently and are evolutionarily conserved. Here we compared the usage and distribution of protein domain families in the free-living proteomes of Archaea, Bacteria and Eukarya and reconstructed species phylogenies while tracing the history of domain emergence and loss in proteomes. We show that both gains and losses of domains occurred frequently during proteome evolution. The rate of domain discovery increased approximately linearly in evolutionary time. Remarkably, gains generally outnumbered losses and the gain-to-loss ratios were much higher in akaryotes compared to eukaryotes. Functional annotations of domain families revealed that both Archaea and Bacteria gained and lost metabolic capabilities during the course of evolution while Eukarya acquired a number of diverse molecular functions including those involved in extracellular processes, immunological mechanisms, and cell regulation. Results also highlighted significant contemporary sharing of informational enzymes between Archaea and Eukarya and metabolic enzymes between Bacteria and Eukarya. Finally, the analysis provided useful insights into the evolution of species. The archaeal superkingdom appeared first in evolution by gradual loss of ancestral domains, bacterial lineages were the first to gain superkingdom-specific domains, and eukaryotes (likely) originated when an expanding proto-eukaryotic stem lineage gained organelles through endosymbiosis of already diversified bacterial lineages. The evolutionary dynamics of domain families in proteomes and the increasing number of domain gains is predicted to redefine the persistence strategies of organisms in superkingdoms, influence the make up of molecular functions, and enhance organismal complexity by the generation of new domain architectures. This dynamics highlights ongoing secondary evolutionary adaptations in akaryotic microbes, especially Archaea.     

Size comparisons among integral membrane transport protein homologues in bacteria, Archaea, and Eucarya

Integral membrane proteins from over 20 ubiquitous families of channels, secondary carriers, and primary active transporters were analyzed for average size differences between homologues from the three domains of life: Bacteria, Archaea, and Eucarya. The results showed that while eucaryotic homologues are consistently larger than their bacterial counterparts, archaeal homologues are significantly smaller. These size differences proved to be due primarily to variations in the sizes of hydrophilic domains localized to the N termini, the C termini, or specific loops between transmembrane alpha-helical spanners, depending on the family. Within the Eucarya domain, plant homologues proved to be substantially smaller than their animal and fungal counterparts. By contrast, extracytoplasmic receptors of ABC-type uptake systems in Archaea proved to be larger on average than those of their bacterial homologues, while cytoplasmic enzymes from different organisms exhibited little or no significant size differences. These observations presumably reflect evolutionary pressure and molecular mechanisms that must have been operative since these groups of organisms diverged from each other.     

The Root of the Universal Tree of Life Inferred from Anciently Duplicated Genes Encoding Components of the Protein-Targeting Machinery

The key protein of the signal recognition particle (termed SRP54 for Eucarya and Ffh for Bacteria) and the protein (termed SRα for Eucarya and Ftsy for bacteria) involved in the recognition and binding of the ribosome SRP nascent polypeptide complex are the products of an ancient gene duplication that appears to predate the divergence of all extant taxa. The paralogy of the genes encoding the two proteins (both of which are GTP triphosphatases) is argued by obvious sequence similarities between the N-terminal half of SRP54(Ffh) and the C-terminal half of SRα(Ftsy). This enables a universal phylogeny based on either protein to be rooted using the second protein as an outgroup. Phylogenetic trees inferred by various methods from an alignment (220 amino acid positions) of the shared SRP54(Ffh) and SRα(Ftsy) regions generate two reciprocally rooted universal trees corresponding to the two genes. The root of both trees is firmly positioned between Bacteria and Archaea/Eucarya, thus providing strong support for the notion (Iwabe et al. 1989; Gogarten et al. 1989) that the first bifurcation in the tree of life separated the lineage leading to Bacteria from a common ancestor to Archaea and Eucarya. None of the gene trees inferred from the two paralogues support a paraphyletic Archaea with the crenarchaeota as a sister group to Eucarya. Received: 19 March 1998 / Accepted: 5 June 1998     

Hot crenarchaeal viruses reveal deep evolutionary connections

The discovery of archaeal viruses provides insights into the fundamental biochemistry and evolution of the Archaea. Recent studies have identified a wide diversity of archaeal viruses within the hot springs of Yellowstone National Park and other high-temperature environments worldwide. These viruses are often morphologically unique and code for genes with little similarity to other known genes in the biosphere, a characteristic that has complicated efforts to trace their evolutionary history. Comparative genomics combined with structural analysis indicate that spindle-shaped virus lineages might be unique to the Archaea, whereas other icosahedral viruses might share a common lineage with viruses of Bacteria and Eukarya. These studies provide insights into the evolutionary history of viruses in all three domains of life.     

The origin and evolution of Archaea: a state of the art

Environmental surveys indicate that the Archaea are diverse and abundant not only in extreme environments, but also in soil, oceans and freshwater, where they may fulfil a key role in the biogeochemical cycles of the planet. Archaea display unique capacities, such as methanogenesis and survival at temperatures higher than 90 degrees C, that make them crucial for understanding the nature of the biota of early Earth. Molecular, genomics and phylogenetics data strengthen Woese's definition of Archaea as a third domain of life in addition to Bacteria and Eukarya. Phylogenomics analyses of the components of different molecular systems are highlighting a core of mainly vertically inherited genes in Archaea. This allows recovering a globally well-resolved picture of archaeal evolution, as opposed to what is observed for Bacteria and Eukarya. This may be due to the fact that no rapid divergence occurred at the emergence of present-day archaeal lineages. This phylogeny supports a hyperthermophilic and non-methanogenic ancestor to present-day archaeal lineages, and a profound divergence between two major phyla, the Crenarchaeota and the Euryarchaeota, that may not have an equivalent in the other two domains of life. Nanoarchaea may not represent a third and ancestral archaeal phylum, but a fast-evolving euryarchaeal lineage. Methanogenesis seems to have appeared only once and early in the evolution of Euryarchaeota. Filling up this picture of archaeal evolution by adding presently uncultivated species, and placing it back in geological time remain two essential goals for the future.     

Molecular evolutionary analyses of the Arabidopsis L7 ribosomal protein gene family

Cytoplasmic ribosomal protein (r-protein) genes in Arabidopsis thaliana are encoded by 80 multigene families that contain between two and seven members. Gene family members are typically similar at the protein sequence level, with the most divergent members of any gene family retaining 94% identity, on average. However, three Arabidopsis r-protein families - S15a, L7 and P2 - contain highly divergent family members. Here, we investigated the organization, structure, expression and molecular evolution of the L7 r-protein family. Phylogenetic analyses showed that L7 r-protein gene family members constitute two distinct phylogenetic groups. The first group including RPL7B, RPL7C and RPL7D has homologs in plants, animals and fungi. The second group represented by RPL7A is found in plants but has no orthologs from other fully-sequenced eukaryotic genomes. These two groups may have derived from a duplication event prior to the divergence of animals and plants. All four L7 r-protein genes are expressed and all exhibit a differential expression in inflorescence and flowers. RPL7A and RPL7B are less expressed than the other genes in all tissues analyzed. Molecular characterization of nucleic and protein sequences of L7 r-protein genes and analysis of their codon usage did not indicate any functional divergence. The probable evolution of an extra-ribosomal function of group 2 genes is discussed.     

艾丁嗜盐小盒菌B2菌株16S rRNA核苷酸序列

艾丁嗜盐小盒菌B2菌株(Haloarcula aidinensis, strain B2)16Sr RNA的核苷酸序列已以双脱氧核苷酸链终止法确定。该菌16Sr RNA显示出了典型的古生物类(Archaea)特性。虽然艾丁嗜盐小盒菌B2菌株在序列方面更接近细菌类(Bacteria)的16SrRNA,但它的序列也显示出与真核生物类(Eucarya)的某些特殊的相似性。在序列和结构方面,该菌与细菌类或真核生物类之间的相似程度要高于细菌类与真核生物类之间的相似程度。另外,该菌16SrRNA的序列与其它嗜盐菌序列相比较支持了以前的结论,即艾丁嗜盐小盒菌B2菌株应属于嗜盐小盒菌属(Haloarcula)的一新种。     

The structure and evolution of the ribosomal proteins encoded in the spc operon of the archaeon (Crenarchaeota) Sulfolobus acidocaldarius.

The genes for nine ribosomal proteins, L24, L5, S14, S8, L6, L18, S5, L30, and L15, have been isolated and sequenced from the spc operon in the archaeon (Crenarchaeota) Sulfolobus acidocaldarius, and the putative amino acid sequence of the proteins coded by these genes has been determined. In addition, three other genes in the spc operon, coding for ribosomal proteins S4E, L32E, and L19E (equivalent to rat ribosomal proteins S4, L32, and L19), were sequenced and the structure of the putative proteins was determined. The order of the ribosomal protein genes in the spc operon of the Crenarchaeota kingdom of Archaea is identical to that present in the Euryarchaeota kingdom of Archaea and also identical to that found in bacteria, except for the genes for r-proteins S4E, L32E, and L19E, which are absent in bacteria. Although AUG is the initiation codon in most of the spc genes, GUG (val) and UUG (leu) are also used as initiation codons in S. acidocaldarius. Over 70% of the codons in the Sulfolobus spc operon have A or U in the third position, reflecting the low GC content of Sulfolobus DNA. Phylogenetic analysis indicated that the archaeal r-proteins are a sister group of their eucaryotic counterparts but did not resolve the question of whether the Archaea is monophyletic, as suggested by the L6P, L15P, and L18P trees, or the question of whether the Crenarchaeota is separate from the Euryarchaeota and closer to the Eucarya, as suggested by the S8P, S5P, and L24P trees. In the case of the three Sulfolobus r-proteins that do not have a counterpart in the bacterial ribosome (S4E, L32E, and L19E), the archaeal r-proteins showed substantial identity to their eucaryotic equivalents, but in all cases the archaeal proteins formed a separate group from the eucaryotic proteins.     

Evolved RNA Secondary Structure and the Rooting of the Universal Tree of Life

The origin and diversification of RNA secondary structure were traced using cladistic methods. Structural components were coded as polarized and ordered multi-state characters, following a model of character state transformation outlined by considerations in statistical mechanics. Several classes of functional RNA were analyzed, including ribosomal RNA (rRNA). Considerable phylogenetic signal was present in their secondary structure. The intrinsically rooted phylogenies reconstructed from evolved RNA structure depicted those derived from nucleic acid sequence at all taxonomical levels, and grouped organisms in concordance with traditional classification, especially in the archaeal and eukaryal domains. Natural selection appears therefore to operate early in the information flow that originates in sequence and ends in an adapted phenotype. When examining the hierarchical classification of the living world, phylogenetic analysis of secondary structure of the small and large rRNA subunits reconstructed a universal tree of life that branched in three monophyletic groups corresponding to Eucarya, Archaea, and Bacteria, and was rooted in the eukaryotic branch. Ribosomal characters involved in the translational cycle could be easily traced and showed that transfer RNA (tRNA) binding domains in the large rRNA subunit evolved concurrently with the rest of the rRNA molecule. Results suggest it is equally parsimonious to consider that ancestral unicellular eukaryotes or prokaryotes gave rise to all extant life forms and provide a rare insight into the early evolution of nucleic acid and protein biosynthesis. Received: 13 September 2000 / Accepted: 27 August 2001     

Archaea and the prokaryote-to-eukaryote transition.

Since the late 1970s, determining the phylogenetic relationships among the contemporary domains of life, the Archaea (archaebacteria), Bacteria (eubacteria), and Eucarya (eukaryotes), has been central to the study of early cellular evolution. The two salient issues surrounding the universal tree of life are whether all three domains are monophyletic (i.e., all equivalent in taxanomic rank) and where the root of the universal tree lies. Evaluation of the status of the Archaea has become key to answering these questions. This review considers our cumulative knowledge about the Archaea in relationship to the Bacteria and Eucarya. Particular attention is paid to the recent use of molecular phylogenetic approaches to reconstructing the tree of life. In this regard, the phylogenetic analyses of more than 60 proteins are reviewed and presented in the context of their participation in major biochemical pathways. Although many gene trees are incongruent, the majority do suggest a sisterhood between Archaea and Eucarya. Altering this general pattern of gene evolution are two kinds of potential interdomain gene transferrals. One horizontal gene exchange might have involved the gram-positive Bacteria and the Archaea, while the other might have occurred between proteobacteria and eukaryotes and might have been mediated by endosymbiosis.     

A congruent phylogenomic signal places eukaryotes within the Archaea

Determining the relationships among the major groups of cellular life is important for understanding the evolution of biological diversity, but is difficult given the enormous time spans involved. In the textbook ‘three domains’ tree based on informational genes, eukaryotes and Archaea share a common ancestor to the exclusion of Bacteria. However, some phylogenetic analyses of the same data have placed eukaryotes within the Archaea, as the nearest relatives of different archaeal lineages. We compared the support for these competing hypotheses using sophisticated phylogenetic methods and an improved sampling of archaeal biodiversity. We also employed both new and existing tests of phylogenetic congruence to explore the level of uncertainty and conflict in the data. Our analyses suggested that much of the observed incongruence is weakly supported or associated with poorly fitting evolutionary models. All of our phylogenetic analyses, whether on small subunit and large subunit ribosomal RNA or concatenated protein-coding genes, recovered a monophyletic group containing eukaryotes and the TACK archaeal superphylum comprising the Thaumarchaeota, Aigarchaeota, Crenarchaeota and Korarchaeota. Hence, while our results provide no support for the iconic three-domain tree of life, they are consistent with an extended eocyte hypothesis whereby vital components of the eukaryotic nuclear lineage originated from within the archaeal radiation.     

Membrane lipids of mesophilic anaerobic bacteria thriving in peats have typical archaeal traits

The 16S ribosomal DNA based distinction between the bacterial and archaeal domains of life is strongly supported by the membrane lipid composition of the two domains; Bacteria generally contain dialkyl glycerol diester lipids, whereas Archaea produce isoprenoid dialkyl glycerol diether and membrane-spanning glycerol dialkyl glycerol tetraether (GDGT) lipids. Here we show that a new group of ecologically abundant membrane-spanning GDGT lipids, containing branched instead of isoprenoid carbon skeletons, are of a bacterial origin. This was revealed by examining the stereochemistry of the glycerol moieties of those branched tetraether membrane lipids, which was found to be the bacterial 1,2-di-O-alkyl-sn-glycerol stereoconfiguration and not the 2,3-di-O-alkyl-sn-glycerol stereoconfiguration as in archaeal membrane lipids. In addition, unequivocal evidence for the presence of cyclopentyl moieties in these bacterial membrane lipids was obtained by NMR. The biochemical traits of biosynthesis of tetraether membrane lipids and the formation of cyclopentyl moieties through internal cyclization, which were thought to be specific for the archaeal lineage of descent, thus also occur in the bacterial domain of life.     

Potential coexistence of both bacterial and eukaryotic small RNA biogenesis and functional related protein homologs in Archaea

RNA silencing plays crucial roles in both bacteria and eukaryotes, yet its machinery appears to differ in these two kingdoms. A couple of Argonaute protein homologs have been reported in some archaeal species in recent years. As Argonaute protein is the key component of eukaryotic RNA silencing pathways, such findings suggested the possibility of existence of eukaryotic RNA silencing like pathways in Archaea, which present the life forms between prokaryotes and eukaryotes. To further explore such hypothesis, we systematically screened 71 fully sequenced archaeal genomes, and identified some proteins containing homologous regions to the functional domains of eukaryotie RNA silencing pathway key proteins. The phylogenetic relationships of these proteins were analyzed. The conserved time-tional amino acids between archaeal and eukaryotic Piwi domains suggested their functional similarity. Our results provide new clues to the evolution of RNA silencing pathways.     

Bayesian phylogenetic analysis reveals two-domain topology of S-adenosylhomocysteine hydrolase protein sequences

S-Adenosylhomocysteine hydrolase (SahH) is involved in the degradation of the compound which inhibits methylation reactions. Using a Bayesian approach and other methods, we reconstructed a phylogenetic tree of amino acid sequences of this protein originating from all three major domains of living organisms. The SahH sequences formed two major branches: one composed mainly of Archaea and the other of eukaryotes and majority of bacteria, clearly contradicting the three-domain topology shown by small subunit rRNA gene. This topology suggests the occurrence of lateral transfer of this gene between the domains. Poor resolution of eukaryotes and bacteria excluded an ultimate conclusion in which out of the two domains this gene appeared first, however, the congruence of the secondary branches with SS rRNA and/or concatenated ribosomal protein datasets phylogenies suggested an "early" acquisition by some bacterial and eukaryotic phyla. Similarly, the branching pattern of Archaea reflected the phylogenies shown by SS rRNA and ribosomal proteins. SahH is widespread in Eucarya, albeit, due to reductive evolution, it is missing in the intracellular parasite Encephalitozoon cuniculi. On the other hand, the lack of affinity to the sequences from the alpha-Proteobacteria and cyanobacteria excludes a possibility of its acquisition in the course of mitochondrial or chloroplast endosymbioses. Unlike Archaea, most bacteria carry MTA/SAH nucleosidase, an enzyme involved also in metabolism of methylthioadenosine. However, the double function of MTA/SAH nucleosidase may be a barrier to ensure the efficient degradation of S-adenosylhomocysteine, specially when the intensity of methylation processes is high. This would explain the presence of S-adenosylhomocysteine hydrolase in the bacteria that have more complex metabolism. On the other hand, majority of obligate pathogenic bacteria due to simpler metabolism rely entirely on MTA/SAH nucleosidase. This could explain the observed phenetic pattern in which bacteria with larger (>6 Mb-million base pairs) genomes carry SAH hydrolase, whereas bacteria that have undergone reductive evolution usually carry MTA/SAH nucleosidase. This suggests that the presence or acquisition of S-adenosylhomocysteine hydrolase in bacteria may predispose towards higher metabolic, and in consequence, higher genomic complexity. The good examples are the phototrophic bacteria all of which carry this gene, however, the SahH phylogeny shows lack of congruence with SSU rRNA and photosyntethic genes, implying that the acquisition was independent and presumably preceded the acquisition of photosyntethic genes. The majority of cyanobacteria acquired this gene from Archaea, however, in some species the sahH gene was replaced by a copy from the beta- or gamma-Proteobacteria.     

Gene duplication and the evolution of group II chaperonins: implications for structure and function

Chaperonins are multisubunit protein-folding assemblies. They are composed of two distinct structural classes, which also have a characteristic phylogenetic distribution. Group I chaperonins (called GroEL/cpn60/hsp60) are present in Bacteria and eukaryotic organelles while group II chaperonins are found in Archaea (called the thermosome or TF55) and the cytoplasm of eukaryotes (called CCT or TriC). Gene duplication has been an important force in the evolution of group II chaperonins: Archaea possess one, two, or three homologous chaperonin subunit-encoding genes, and eight distinct CCT gene families (paralogs) have been described in eukaryotes. Phylogenetic analyses indicate that while the duplications in archaeal chaperonin genes have occurred numerous times independently in a lineage-specific fashion, the eight different CCT subunits found in eukaryotes are the products of duplications that occurred early and very likely only once in the evolution of the eukaryotic nuclear genome. Analyses of CCT sequences from diverse eukaryotic species reveal that each of the CCT subunits possesses a suite of invariant subunit-specific amino acid residues ("signatures"). When mapped onto the crystal structure of the archaeal chaperonin from Thermoplasma acidophilum, these signatures are located in the apical, intermediate, and equatorial domains. Regions that were found to be variable in length and/or amino acid sequence were localized primarily to the exterior of the molecule and, significantly, to the extreme tip of the apical domain (the "helical protrusion"). In light of recent biochemical and electron microscopic data describing specific CCT-substrate interactions, our results have implications for the evolution of subunit-specific functions in CCT.     

The Evolution of Histidine Biosynthesis in Archaea: Insights into the <Emphasis Type="Italic">his</Emphasis> Genes Structure and Organization in LUCA

The available sequences of genes encoding the enzymes associated with histidine biosynthesis suggest that this is an ancient metabolic pathway that was assembled prior to the diversification of Bacteria, Archaea, and Eucarya. Paralogous duplication, gene elongation, and fusion events of several different his genes have played a major role in shaping this biosynthetic route. We have analyzed the structure and organization of histidine biosynthetic genes from 55 complete archaeal genomes and combined it with phylogenetic inference in order to investigate the mechanisms responsible for the assembly of the his pathway and the origin of his operons. We show that a wide variety of different organizations of his genes exists in Archaea and that some his genes or entire his (sub-)operons have been likely transferred horizontally between Archaea and Bacteria. However, we show that, in most Archaea, his genes are monofunctional (except for hisD) and scattered throughout the genome, suggesting that his operons might have been assembled multiple times during evolution and that in some cases they are the result of recent evolutionary events. An evolutionary model for the structure and organization of his genes in LUCA is proposed.