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.     

Phylogenetic and phyletic studies of informational genes in genomes highlight existence of a 4 domain of life including giant viruses

The discovery of Mimivirus, with its very large genome content, made it possible to identify genes common to the three domains of life (Eukarya, Bacteria and Archaea) and to generate controversial phylogenomic trees congruent with that of ribosomal genes, branching Mimivirus at its root. Here we used sequences from metagenomic databases, Marseillevirus and three new viruses extending the Mimiviridae family to generate the phylogenetic trees of eight proteins involved in different steps of DNA processing. Compared to the three ribosomal defined domains, we report a single common origin for Nucleocytoplasmic Large DNA Viruses (NCLDV), DNA processing genes rooted between Archaea and Eukarya, with a topology congruent with that of the ribosomal tree. As for translation, we found in our new viruses, together with Mimivirus, five proteins rooted deeply in the eukaryotic clade. In addition, comparison of informational genes repertoire based on phyletic pattern analysis supports existence of a clade containing NCLDVs clearly distinct from that of Eukarya, Bacteria and Archaea. We hypothesize that the core genome of NCLDV is as ancient as the three currently accepted domains of life.     

Comparison of microbial populations in model and natural rumens using 16S ribosomal RNA-targeted probes

A model rumen system, dual-flow continuous culture fermenters, was evaluated by two comparative criteria in two experiments using ribosomal (r)RNA-targeted DNA probes to compare key microbial groups in samples. The initial experiment measured temporal changes in population structure during adaptation of ruminal microbial populations in fermenters over 240 h. The fermenter inoculum contained 34.9% Bacteria, 60.1% Eukarya and 6.8% Archaea measured as a fraction of total small subunit (SSU) rRNA quantified using a universal probe. The cellulolytic bacterial genus Fibrobacter comprised 9.5% of total SSU rRNA in the inoculum. After 240 h of fermenter operation, the average abundance was 80.9% Bacteria, 6.1% Eukarya, 5.1% Archaea and Fibrobacter genus accounted for 6.6% of the total SSU rRNA. Divergence between ruminal and fermenter population structure was evaluated in the second experiment and samples were classified as ruminal, inoculum or fermenter (96, 120, 144 and 168 h of fermenter operation). Fermenter samples had higher relative abundances of Bacteria (84.5%) and Archaea (2.1%) and lower relative abundances of Eukarya (1.8%) than ruminal samples (average 48.0% Bacteria, 1.3% Archaea and 61.5% Eukarya). The relative abundance of Fibrobacter was similar in all samples, averaging 2.5%. The ruminal and fermenter samples had similar proportions of F. succinogenes and F. succinogenes subgroup 3 (as a percentage of Fibrobacter SSU rRNA). Fibrobacter succinogenes subgroup 1 and F. intestinalis proportions of Fibrobacter were lower in fermenter samples (8.2% and 0.7% respectively) than in ruminal samples (28.4% and 2.2% respectively). Fermenters were able to maintain a core prokaryotic community structure similar to the native microbial community in the rumen. Although protozoa populations were lost, maintenance of Fibrobacter and archaeal populations indicated that the model system supported a functional community structure similar to the rumen. This model rumen system may serve as a suitable tool for studying aspects of ruminal microbial ecology and may resolve some of the relationships between microbial community structure and function by providing control of experimental conditions.     

The Origin and Evolution of tRNA Inferred from Phylogenetic Analysis of Structure

The evolutionary history of the two structural and functional domains of tRNA is controversial but harbors the secrets of early translation and the genetic code. To explore the origin and evolution of tRNA, we reconstructed phylogenetic trees directly from molecular structure. Forty-two structural characters describing the geometry of 571 tRNAs and three statistical parameters describing thermodynamic and mechanical features of molecules quantitatively were used to derive phylogenetic trees of molecules and molecular substructures. Trees of molecules failed to group tRNA according to amino acid specificity and did not reveal the tripartite nature of life, probably due to loss of phylogenetic signal or because tRNA diversification predated organismal diversification. Trees of substructures derived from both structural and statistical characters support the origin of tRNA in the acceptor arm and the hypothesis that the top half domain composed of acceptor and pseudouridine (TΨC) arms is more ancient than the bottom half domain composed of dihydrouridine (DHU) and anticodon arms. This constitutes the cornerstone of the genomic tag hypothesis that postulates tRNAs were ancient telomeres in the RNA world. The trees of substructures suggest a model for the evolution of the major functional and structural components of tRNA. In this model, short RNA hairpins with stems homologous to the acceptor arm of present day tRNAs were extended with regions homologous to TΨC and anticodon arms. The DHU arm was then incorporated into the resulting three-stemmed structure to form a proto-cloverleaf structure. The variable region was the last structural addition to the molecular repertoire of evolving tRNA substructures. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The Origin of Modern 5S rRNA: A Case of Relating Models of Structural History to Phylogenetic Data

Evolutionary models of molecular structures must incorporate molecular information at different levels of structural complexity and must be phrased within a phylogenetic perspective. In this regard, phylogenetic trees of substructures that are reconstructed from molecular features that contribute to order and thermodynamic stability show that a gradual model of evolution of 5S rRNA structure is more parsimonious than models that invoke large segmental duplications of the molecule. The search for trees of substructures that are most parsimonious, by their very nature, defines an objective strategy to select models of molecular change that best fit structural data. When combined with additional data, such as the age of protein domains that interact with RNA substructures, these trees can be used to falsify unlikely hypotheses.     

The archaeal origins of the eukaryotic translational system

Among the 78 eukaryotic ribosomal proteins, eleven are specific to Eukarya, 33 are common only to Archaea and Eukarya and 34 are homologous (at least in part) to those of both Bacteria and Archaea. Several other translational proteins are common only to Eukarya and Archaea (e.g., IF2a, SRP19, etc.), whereas others are shared by the three phyla (e.g., EFTu/EF1A and SRP54). Although this and other analyses strongly support an archaeal origin for a substantial fraction of the eukaryotic translational machinery, especially the ribosomal proteins, there have been numerous unique and ubiquitous additions to the eukaryotic translational system besides the 11 unique eukaryotic ribosomal proteins. These include peptide additions to most of the 67 archaeal homolog proteins, rRNA insertions, the 5.8S RNA and the Alu extension to the SRP RNA. Our comparative analysis of these and other eukaryotic features among the three different cellular phylodomains supports the idea that an archaeal translational system was most likely incorporated by means of endosymbiosis into a host cell that was neither bacterial nor archaeal in any modern sense. Phylogenetic analyses provide support for the timing of this acquisition coinciding with an ancient bottleneck in prokaryotic diversity.     

Detailed analysis of RNA-protein interactions within the bacterial ribosomal protein L5/5S rRNA complex

The crystal structure of ribosomal protein L5 from Thermus thermophilus complexed with a 34-nt fragment comprising helix III and loop C of Escherichia coli 5S rRNA has been determined at 2.5 A resolution. The protein specifically interacts with the bulged nucleotides at the top of loop C of 5S rRNA. The rRNA and protein contact surfaces are strongly stabilized by intramolecular interactions. Charged and polar atoms forming the network of conserved intermolecular hydrogen bonds are located in two narrow planar parallel layers belonging to the protein and rRNA, respectively. The regions, including these atoms conserved in Bacteria and Archaea, can be considered an RNA-protein recognition module. Comparison of the T. thermophilus L5 structure in the RNA-bound form with the isolated Bacillus stearothermophilus L5 structure shows that the RNA-recognition module on the protein surface does not undergo significant changes upon RNA binding. In the crystal of the complex, the protein interacts with another RNA molecule in the asymmetric unit through the beta-sheet concave surface. This protein/RNA interface simulates the interaction of L5 with 23S rRNA observed in the Haloarcula marismortui 50S ribosomal subunit.     

Evolution according to large ribosomal subunit RNA

Evolutionary trees were constructed, by distance methods, from an alignment of 225 complete large subunit (LSU) rRNA sequences, representing Eucarya, Archaea, Bacteria, plastids, and mitochondria. A comparison was made with trees based on sets of small subunit (SSU) rRNA sequences. Trees constructed on the set of 172 species and organelles for which the sequences of both molecules are known had a very similar topology, at least with respect to the divergence order of large taxa such as the eukaryotic kingdoms and the bacterial divisions. However, since there are more than ten times as many SSU as LSU rRNA sequences, it is possible to select many SSU rRNA sequence sets of equivalent size but different species composition. The topologies of these trees showed considerable differences according to the particular species set selected.The effect of the dataset and of different distance correction methods on tree topology was tested for both LSU and SSU rRNA by repetitive random sampling of a single species from each large taxon. The impact of the species set on the topology of the resulting consensus trees is much lower using LSU than using SSU rRNA. This might imply that LSU rRNA is a better molecule for studying wide-range relationships. The mitochondria behave clearly as a monophyletic group, clustering with the Proteobacteria. Gram-positive bacteria appear as two distinct groups, which are found clustered together in very few cases. Archaea behave as if monophyletic in most cases, but with a low confidence.Abbreviations LSU rRNA large subunit ribosomal RNA - SSU rRNA small subunit ribosomal RNA - JC Jukes and Cantor - JN Jin and Nei Correspondence to: R. De Wachter     

Genomics and early cellular evolution. The origin of the DNA world.

The sequencing of several genomes from each of the three domains of life (Archaea, Bacteria and Eukarya) has provided a huge amount of data that can be used to gain insight about early cellular evolution. Some features of the universal tree of life based on rRNA polygenies have been confirmed, such as the division of the cellular living world into three domains. The monophyly of each domain is supported by comparative genomics. However, the hyperthermophilic nature of the 'last universal common ancestor' (LUCA) is not confirmed. Comparative genomics has revealed that gene transfers have been (and still are) very frequent in genome evolution. Nevertheless, a core of informational genes appears more resistant to transfer, testifying for a close relationship between archaeal and eukaryal informational processes. This observation can be explained either by a common unique history between Archaea and Eukarya or by an atypical evolution of these systems in Bacteria. At the moment, comparative genomics still does not allow to choose between a simple LUCA, possibly with an RNA genome, or a complex LUCA, with a DNA genome and informational mechanisms similar to those of Archaea and Eukarya. Further comparative studies on informational mechanisms in the three domains should help to resolve this critical question. The role of viruses in the origin and evolution of DNA genomes also appears an area worth of active investigations. I suggest here that DNA and DNA replication mechanisms appeared first in the virus world before being transferred into cellular organisms.     

Identification and role of functionally important motifs in the 970 loop of Escherichia coli 16S ribosomal RNA

The 970 loop (helix 31) of Escherichia coli 16S ribosomal RNA contains two modified nucleotides, m²G966 and m⁵C967. Positions A964, A969, and C970 are conserved among the Bacteria, Archaea, and Eukarya. The nucleotides present at positions 965, 966, 967, 968, and 971, however, are only conserved and unique within each domain. All organisms contain a modified nucleoside at position 966, but the type of the modification is domain specific. Biochemical and structure studies have placed this loop near the P site and have shown it to be involved in the decoding process and in binding the antibiotic tetracycline. To identify the functional components of this ribosomal RNA hairpin, the eight nucleotides of the 970 loop of helix 31 were subjected to saturation mutagenesis and 107 unique functional mutants were isolated and analyzed. Nonrandom nucleotide distributions were observed at each mutated position among the functional isolates. Nucleotide identity at positions 966 and 969 significantly affects ribosome function. Ribosomes with single mutations of m²G966 or m⁵C967 produce more protein in vivo than do wild-type ribosomes. Overexpression of initiation factor 3 specifically restored wild-type levels of protein synthesis to the 966 and 967 mutants, suggesting that modification of these residues is important for initiation factor 3 binding and for the proper initiation of protein synthesis.     

Structural motifs of the bacterial ribosomal proteins S20, S18 and S16 that contact rRNA present in the eukaryotic ribosomal proteins S25, S26 and S27A,respectively

The majority of constitutive proteins in the bacterial 30S ribosomal subunit have orthologues in Eukarya and Archaea. The eukaryotic counterparts for the remainder (S6, S16, S18 and S20) have not been identified. We assumed that amino acid residues in the ribosomal proteins that contact rRNA are to be constrained in evolution and that the most highly conserved of them are those residues that are involved in forming the secondary protein structure. We aligned the sequences of the bacterial ribosomal proteins from the S20p, S18p and S16p families, which make multiple contacts with rRNA in the Thermus thermophilus 30S ribosomal subunit (in contrast to the S6p family), with the sequences of the unassigned eukaryotic small ribosomal subunit protein families. This made it possible to reveal that the conserved structural motifs of S20p, S18p and S16p that contact rRNA in the bacterial ribosome are present in the ribosomal proteins S25e, S26e and S27Ae, respectively. We suggest that ribosomal protein families S20p, S18p and S16p are homologous to the families S25e, S26e and S27Ae, respectively.     

Global phylogeny determined by the combination of protein domains in proteomes

The majority of proteins consist of multiple domains that are either repeated or combined in defined order. In this study, we survey the combination of protein domains defined at fold and fold superfamily levels in 185 genomes belonging to organisms that have been fully sequenced and introduce a method that reconstructs rooted phylogenomic trees from the content and arrangement of domains in proteins at a genomic level. We find that the majority of domain combinations were unique to Archaea, Bacteria, or Eukarya, suggesting most combinations originated after life had diversified. Domain repeat and domain repeat within multidomain proteins increased notably in eukaryotes, mainly at the expense of single-domain and domain-pair proteins. This increase was mostly confined to Metazoa. We also find an unbalanced sharing of domain combinations which suggests that Eukarya is more closely related to Bacteria than to Archaea, an observation that challenges the widely assumed eukaryote-archaebacterial sisterhood relationship. The occurrence and abundance of the molecular repertoire (interactome) of domain combinations was used to generate phylogenomic trees. These global interactome-based phylogenies described organismal histories satisfactorily, revealing the tripartite nature of life, and supporting controversial evolutionary patterns, such as the Coelomata hypothesis, the grouping of plants and animals, and the Gram-positive origin of bacteria. Results suggest strongly that the process of domain combination is not random but curved by evolution, rejecting the null hypothesis of domain modules combining in the absence of natural selection or an optimality criterion.     

具有真细菌和真核生物融合特征的古生菌转录系统

古生菌是一类区别于真细菌和真核生物的第三域生命形式 ,转录是生物体遗传信息传递系统中的一个中心环节。近年来研究结果表明 ,古生菌的转录系统具有真细菌和真核生物的融合特征 :古生菌的基本转录装置包括RNA聚合酶、基本转录因子、启动子元件等与真核生物相似 ;而古生菌的转录调控机制却更加类似于真细菌 ,在古生菌中发现并鉴定了许多类似于真细菌的转录调控蛋白。另外古生菌还具有某些独特的转录调控方式     

Eukaryotes first: how could that be?

In the half century since the formulation of the prokaryote : eukaryote dichotomy, many authors have proposed that the former evolved from something resembling the latter, in defiance of common (and possibly common sense) views. In such ‘eukaryotes first’ (EF) scenarios, the last universal common ancestor is imagined to have possessed significantly many of the complex characteristics of contemporary eukaryotes, as relics of an earlier ‘progenotic’ period or RNA world. Bacteria and Archaea thus must have lost these complex features secondarily, through ‘streamlining’. If the canonical three-domain tree in which Archaea and Eukarya are sisters is accepted, EF entails that Bacteria and Archaea are convergently prokaryotic. We ask what this means and how it might be tested.     

Protein S20 Binds Two 16S rRNA Sites as Assembly Is Initiated

Ribosomal protein S20 is a primary binding protein that bridges the 5′ domain and the 3′ minor domain of the 16S ribosomal RNA (rRNA) in the 30S ribosomal subunit. Using time-dependent dimethyl sulfate modification, we have determined that as it is bound to 16S rRNA, protein S20 causes rapid protection of bases A246, A274, A279, and A282 in the stem region of helix 11 in the 5′ domain and moderately fast modifications of helix 44 bases A1433 and A1434 in the 3′ minor domain. At a later time, enhancements occur with bases A181and A190 in helix 9, bases A325 and A327 in helix 13, and base C264 at the distal end of helix 11 in the 5′ domain of 16S rRNA. The modifications that occur in the stem region of helix 11 are distant from the binding site of protein S20, as determined from the crystal structure. Simultaneous addition of protein S17 with S20 to the complex significantly alters the modifications caused by protein S20 in the stem region of helix 11 but does not alter the remaining modifications. Our results indicate that protein S20 is binding to at least two alternate 16S rRNA sites during the early assembly process.     

An HflX-type GTPase from Sulfolobus solfataricus binds to the 50S ribosomal subunit in all nucleotide-bound states

HflX GTPases are found in all three domains of life, the Bacteria, Archaea, and Eukarya. HflX from Escherichia coli has been shown to bind to the 50S ribosomal subunit in a nucleotide-dependent manner, and this interaction strongly stimulates its GTPase activity. We recently determined the structure of an HflX ortholog from the archaeon Sulfolobus solfataricus (SsoHflX). It revealed the presence of a novel HflX domain that might function in RNA binding and is linked to a canonical G domain. This domain arrangement is common to all archaeal, bacterial, and eukaryotic HflX GTPases. This paper shows that the archaeal SsoHflX, like its bacterial orthologs, binds to the 50S ribosomal subunit. This interaction does not depend on the presence of guanine nucleotides. The HflX domain is sufficient for ribosome interaction. Binding appears to be restricted to free 50S ribosomal subunits and does not occur with 70S ribosomes engaged in translation. The fingerprint (1)H-(15)N heteronuclear correlation nuclear magnetic resonance (NMR) spectrum of SsoHflX reveals a large number of well-resolved resonances that are broadened upon binding to the 50S ribosomal subunit. The GTPase activity of SsoHflX is stimulated by crude fractions of 50S ribosomal subunits, but this effect is lost with further high-salt purification of the 50S ribosomal subunits, suggesting that the stimulation depends on an extrinsic factor bound to the 50S ribosomal subunit. Our results reveal common properties but also marked differences between archaeal and bacterial HflX proteins.     

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.     

Universal and domain-specific sequences in 23S–28S ribosomal RNA identified by computational phylogenetics

Comparative analysis of ribosomal RNA (rRNA) sequences has elucidated phylogenetic relationships. However, this powerful approach has not been fully exploited to address ribosome function. Here we identify stretches of evolutionarily conserved sequences, which correspond with regions of high functional importance. For this, we developed a structurally aligned database, FLORA (full-length organismal rRNA alignment) to identify highly conserved nucleotide elements (CNEs) in 23S–28S rRNA from each phylogenetic domain (Eukarya, Bacteria, and Archaea). Universal CNEs (uCNEs) are conserved in sequence and structural position in all three domains. Those in regions known to be essential for translation validate our approach. Importantly, some uCNEs reside in areas of unknown function, thus identifying novel sequences of likely great importance. In contrast to uCNEs, domain-specific CNEs (dsCNEs) are conserved in just one phylogenetic domain. This is the first report of conserved sequence elements in rRNA that are domain-specific; they are largely a eukaryotic phenomenon. The locations of the eukaryotic dsCNEs within the structure of the ribosome suggest they may function in nascent polypeptide transit through the ribosome tunnel and in tRNA exit from the ribosome. Our findings provide insights and a resource for ribosome function studies.