A korarchaeal genome reveals insights into the evolution of the Archaea

The origin of eukaryotes and their relationship with the Archaea: are we at a phylogenomic impasse?

The origin of eukaryotes and their evolutionary relationship with the Archaea is a major biological question and the subject of intense debate. In the context of the classical view of the universal tree of life, the Archaea and the Eukarya have a common ancestor, the nature of which remains undetermined. Alternative views propose instead that the Eukarya evolved directly from a bona fide archaeal lineage. Several recent large-scale phylogenomic studies using an array of approaches are divided in supporting either one or the other scenario, despite analysing largely overlapping data sets of universal genes. We examine the reasons for such a lack of consensus and consider how alternative approaches may enable progress in answering this fascinating and as-yet-unresolved question.

Evolution of diverse cell division and vesicle formation systems in Archaea

Recently a novel cell division system comprised of homologues of eukaryotic ESCRT-III (endosomal sorting complex required for transport III) proteins was discovered in the hyperthermophilic crenarchaeote Sulfolobus acidocaldarius. On the basis of this discovery, we undertook a comparative genomic analysis of the machineries for cell division and vesicle formation in Archaea. Archaea possess at least three distinct membrane remodelling systems: the FtsZ-based bacterial-type system, the ESCRT-III-based eukaryote-like system and a putative novel system that uses an archaeal actin-related protein. Many archaeal genomes encode assortments of components from different systems. Evolutionary reconstruction from these findings suggests that the last common ancestor of the extant Archaea possessed a complex membrane remodelling apparatus, different components of which were lost during subsequent evolution of archaeal lineages. By contrast, eukaryotes seem to have inherited all three ancestral systems.

Molecular signatures for the Crenarchaeota and the Thaumarchaeota

Crenarchaeotes found in mesophilic marine environments were recently placed into a new phylum of Archaea called the Thaumarchaeota. However, very few molecular characteristics of this new phylum are currently known which can be used to distinguish them from the Crenarchaeota. In addition, their relationships to deep-branching archaeal lineages are unclear. We report here detailed analyses of protein sequences from Crenarchaeota and Thaumarchaeota that have identified many conserved signature indels (CSIs) and signature proteins (SPs) (i.e., proteins for which all significant blast hits are from these groups) that are specific for these archaeal groups. Of the identified signatures 6 CSIs and 13 SPs are specific for the Crenarchaeota phylum; 6 CSIs and >250 SPs are uniquely found in various Thaumarchaeota (viz. Cenarchaeum symbiosum, Nitrosopumilus maritimus and a number of uncultured marine crenarchaeotes) and 3 CSIs and ~10 SPs are found in both Thaumarchaeota and Crenarchaeota species. Some of the molecular signatures are also present in Korarchaeum cryptofilum, which forms the independent phylum Korarchaeota. Although some of these molecular signatures suggest a distant shared ancestry between Thaumarchaeota and Crenarchaeota, our identification of large numbers of Thaumarchaeota-specific proteins and their deep branching between the Crenarchaeota and Euryarchaeota phyla in phylogenetic trees shows that they are distinct from both Crenarchaeota and Euryarchaeota in both genetic and phylogenetic terms. These observations support the placement of marine mesophilic archaea into the separate phylum Thaumarchaeota. Additionally, many CSIs and SPs have been found that are specific for different orders within Crenarchaeota (viz. Sulfolobales-3 CSIs and 169 SPs, Thermoproteales-5 CSIs and 25 SPs, Desulfurococcales-4 SPs, and Sulfolobales and Desulfurococcales-2 CSIs and 18 SPs). The signatures described here provide novel means for distinguishing the Crenarchaeota and the Thaumarchaeota and for the classification of related and novel species in different environments. Functional studies on these signature proteins could lead to discovery of novel biochemical properties that are unique to these groups of archaea.

Archaic chaos: intrinsically disordered proteins in Archaea

BACKGROUND

Many proteins or their regions known as intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) lack unique 3D structure in their native states under physiological conditions yet fulfill key biological functions. Earlier bioinformatics studies showed that IDPs and IDRs are highly abundant in different proteomes and carry out mostly regulatory functions related to molecular recognition and signal transduction. Archaea belong to an intriguing domain of life whose members, being microbes, are characterized by a unique mosaic-like combination of bacterial and eukaryotic properties and include inhabitants of some of the most extreme environments on the planet. With the expansion of the archaea genome data (more than fifty archaea species from five different phyla are known now), and with recent improvements in the accuracy of intrinsic disorder prediction, it is time to re-examine the abundance of IDPs and IDRs in the archaea domain.

RESULTS

The abundance of IDPs and IDRs in 53 archaea species is analyzed. The amino acid composition profiles of these species are generally quite different from each other. The disordered content is highly species-dependent. Thermoproteales proteomes have 14% of disordered residues, while in Halobacteria, this value increases to 34%. In proteomes of these two phyla, proteins containing long disordered regions account for 12% and 46%, whereas 4% and 26% their proteins are wholly disordered. These three measures of disorder content are linearly correlated with each other at the genome level. There is a weak correlation between the environmental factors (such as salinity, pH and temperature of the habitats) and the abundance of intrinsic disorder in Archaea, with various environmental factors possessing different disorder-promoting strengths. Harsh environmental conditions, especially those combining several hostile factors, clearly favor increased disorder content. Intrinsic disorder is highly abundant in functional Pfam domains of the archaea origin. The analysis based on the disordered content and phylogenetic tree indicated diverse evolution of intrinsic disorder among various classes and species of Archaea.

CONCLUSIONS

Archaea proteins are rich in intrinsic disorder. Some of these IDPs and IDRs likely evolve to help archaea to accommodate to their hostile habitats. Other archaean IDPs and IDRs possess crucial biological functions similar to those of the bacterial and eukaryotic IDPs/IDRs.

Enigmatic, ultrasmall, uncultivated Archaea

Metagenomics has provided access to genomes of as yet uncultivated microorganisms in natural environments, yet there are gaps in our knowledge-particularly for Archaea-that occur at relatively low abundance and in extreme environments. Ultrasmall cells (<500 nm in diameter) from lineages without cultivated representatives that branch near the crenarchaeal/euryarchaeal divide have been detected in a variety of acidic ecosystems. We reconstructed composite, near-complete approximately 1-Mb genomes for three lineages, referred to as ARMAN (archaeal Richmond Mine acidophilic nanoorganisms), from environmental samples and a biofilm filtrate. Genes of two lineages are among the smallest yet described, enabling a 10% higher coding density than found genomes of the same size, and there are noncontiguous genes. No biological function could be inferred for up to 45% of genes and no more than 63% of the predicted proteins could be assigned to a revised set of archaeal clusters of orthologous groups. Some core metabolic genes are more common in Crenarchaeota than Euryarchaeota, up to 21% of genes have the highest sequence identity to bacterial genes, and 12 belong to clusters of orthologous groups that were previously exclusive to bacteria. A small subset of 3D cryo-electron tomographic reconstructions clearly show penetration of the ARMAN cell wall and cytoplasmic membranes by protuberances extended from cells of the archaeal order Thermoplasmatales. Interspecies interactions, the presence of a unique internal tubular organelle [Comolli, et al. (2009) ISME J 3:159-167], and many genes previously only affiliated with Crenarchaea or Bacteria indicate extensive unique physiology in organisms that branched close to the time that Cren- and Euryarchaeotal lineages diverged.

The origin and early evolution of eukaryotes in the light of phylogenomics

Phylogenomics of eukaryote supergroups suggest a highly complex last common ancestor of eukaryotes and a key role of mitochondrial endosymbiosis in the origin of eukaryotes.

One bacterial cell, one complete genome

While the bulk of the finished microbial genomes sequenced to date are derived from cultured bacterial and archaeal representatives, the vast majority of microorganisms elude current culturing attempts, severely limiting the ability to recover complete or even partial genomes from these environmental species. Single cell genomics is a novel culture-independent approach, which enables access to the genetic material of an individual cell. No single cell genome has to our knowledge been closed and finished to date. Here we report the completed genome from an uncultured single cell of Candidatus Sulcia muelleri DMIN. Digital PCR on single symbiont cells isolated from the bacteriome of the green sharpshooter Draeculacephala minerva bacteriome allowed us to assess that this bacteria is polyploid with genome copies ranging from approximately 200–900 per cell, making it a most suitable target for single cell finishing efforts. For single cell shotgun sequencing, an individual Sulcia cell was isolated and whole genome amplified by multiple displacement amplification (MDA). Sanger-based finishing methods allowed us to close the genome. To verify the correctness of our single cell genome and exclude MDA-derived artifacts, we independently shotgun sequenced and assembled the Sulcia genome from pooled bacteriomes using a metagenomic approach, yielding a nearly identical genome. Four variations we detected appear to be genuine biological differences between the two samples. Comparison of the single cell genome with bacteriome metagenomic sequence data detected two single nucleotide polymorphisms (SNPs), indicating extremely low genetic diversity within a Sulcia population. This study demonstrates the power of single cell genomics to generate a complete, high quality, non-composite reference genome within an environmental sample, which can be used for population genetic analyzes.

Microbiology and geochemistry of Little Hot Creek, a hot spring environment in the Long Valley Caldera

A culture-independent community census was combined with chemical and thermodynamic analyses of three springs located within the Long Valley Caldera, Little Hot Creek (LHC) 1, 3, and 4. All three springs were approximately 80 degrees C, circumneutral, apparently anaerobic and had similar water chemistries. 16S rRNA gene libraries constructed from DNA isolated from spring sediment revealed moderately diverse but highly novel microbial communities. Over half of the phylotypes could not be grouped into known taxonomic classes. Bacterial libraries from LHC1 and LHC3 were predominantly species within the phyla Aquificae and Thermodesulfobacteria, while those from LHC4 were dominated by candidate phyla, including OP1 and OP9. Archaeal libraries from LHC3 contained large numbers of Archaeoglobales and Desulfurococcales, while LHC1 and LHC4 were dominated by Crenarchaeota unaffiliated with known orders. The heterogeneity in microbial populations could not easily be attributed to measurable differences in water chemistry, but may be determined by availability of trace amounts of oxygen to the spring sediments. Thermodynamic modeling predicted the most favorable reactions to be sulfur and nitrate respirations, yielding 40-70 kJ mol(-1) e(-) transferred; however, levels of oxygen at or below our detection limit could result in aerobic respirations yielding up to 100 kJ mol(-1) e(-) transferred. Important electron donors are predicted to be H(2), H(2)S, S(0), Fe(2+) and CH(4), all of which yield similar energies when coupled to a given electron acceptor. The results indicate that springs associated with the Long Valley Caldera contain microbial populations that show some similarities both to springs in Yellowstone and springs in the Great Basin.

Agmatine-conjugated cytidine in a tRNA anticodon is essential for AUA decoding in archaea

A modified base at the first (wobble) position of some tRNA anticodons is critical for deciphering the genetic code. In eukaryotes and eubacteria, AUA codons are decoded by tRNAsIle with modified bases pseudouridine (and/or inosine) and lysidine, respectively. The mechanism by which archaeal species translate AUA codons is unclear. We describe a polyamine-conjugated modified base, 2-agmatinylcytidine (agm(2)C or agmatidine), at the wobble position of archaeal tRNA(Ile) that decodes AUA codons specifically. We demonstrate that archaeal cells use agmatine to synthesize agm(2)C of tRNA(Ile). We also identified a new enzyme, tRNA(Ile)-agm(2)C synthetase (TiaS), that catalyzes agm(2)C formation in the presence of agmatine and ATP. Although agm(2)C is chemically similar to lysidine, TiaS constitutes a distinct class of enzyme from tRNA(Ile)-lysidine synthetase (TilS), suggesting that the decoding systems evolved convergently across domains.

Origin of the cell nucleus, mitosis and sex: roles of intracellular coevolution

BACKGROUND

The transition from prokaryotes to eukaryotes was the most radical change in cell organisation since life began, with the largest ever burst of gene duplication and novelty. According to the coevolutionary theory of eukaryote origins, the fundamental innovations were the concerted origins of the endomembrane system and cytoskeleton, subsequently recruited to form the cell nucleus and coevolving mitotic apparatus, with numerous genetic eukaryotic novelties inevitable consequences of this compartmentation and novel DNA segregation mechanism. Physical and mutational mechanisms of origin of the nucleus are seldom considered beyond the long-standing assumption that it involved wrapping pre-existing endomembranes around chromatin. Discussions on the origin of sex typically overlook its association with protozoan entry into dormant walled cysts and the likely simultaneous coevolutionary, not sequential, origin of mitosis and meiosis.

RESULTS

I elucidate nuclear and mitotic coevolution, explaining the origins of dicer and small centromeric RNAs for positionally controlling centromeric heterochromatin, and how 27 major features of the cell nucleus evolved in four logical stages, making both mechanisms and selective advantages explicit: two initial stages (origin of 30 nm chromatin fibres, enabling DNA compaction; and firmer attachment of endomembranes to heterochromatin) protected DNA and nascent RNA from shearing by novel molecular motors mediating vesicle transport, division, and cytoplasmic motility. Then octagonal nuclear pore complexes (NPCs) arguably evolved from COPII coated vesicle proteins trapped in clumps by Ran GTPase-mediated cisternal fusion that generated the fenestrated nuclear envelope, preventing lethal complete cisternal fusion, and allowing passive protein and RNA exchange. Finally, plugging NPC lumens by an FG-nucleoporin meshwork and adopting karyopherins for nucleocytoplasmic exchange conferred compartmentation advantages. These successive changes took place in naked growing cells, probably as indirect consequences of the origin of phagotrophy. The first eukaryote had 1-2 cilia and also walled resting cysts; I outline how encystation may have promoted the origin of meiotic sex. I also explain why many alternative ideas are inadequate.

CONCLUSION

Nuclear pore complexes are evolutionary chimaeras of endomembrane- and mitosis-related chromatin-associated proteins. The keys to understanding eukaryogenesis are a proper phylogenetic context and understanding organelle coevolution: how innovations in one cell component caused repercussions on others.

Agmatidine, a modified cytidine in the anticodon of archaeal tRNA(Ile), base pairs with adenosine but not with guanosine

Modification of the cytidine in the first anticodon position of the AUA decoding tRNA(Ile) (tRNA2(Ile)) of bacteria and archaea is essential for this tRNA to read the isoleucine codon AUA and to differentiate between AUA and the methionine codon AUG. To identify the modified cytidine in archaea, we have purified this tRNA species from Haloarcula marismortui, established its codon reading properties, used liquid chromatography-mass spectrometry (LC-MS) to map RNase A and T1 digestion products onto the tRNA, and used LC-MS/MS to sequence the oligonucleotides in RNase A digests. These analyses revealed that the modification of cytidine in the anticodon of tRNA2(Ile) adds 112 mass units to its molecular mass and makes the glycosidic bond unusually labile during mass spectral analyses. Accurate mass LC-MS and LC-MS/MS analysis of total nucleoside digests of the tRNA2(Ile) demonstrated the absence in the modified cytidine of the C2-oxo group and its replacement by agmatine (decarboxy-arginine) through a secondary amine linkage. We propose the name agmatidine, abbreviation C(+), for this modified cytidine. Agmatidine is also present in Methanococcus maripaludis tRNA2(Ile) and in Sulfolobus solfataricus total tRNA, indicating its probable occurrence in the AUA decoding tRNA(Ile) of euryarchaea and crenarchaea. The identification of agmatidine shows that bacteria and archaea have developed very similar strategies for reading the isoleucine codon AUA while discriminating against the methionine codon AUG.

Structure and function of the GINS complex, a key component of the eukaryotic replisome

High-fidelity chromosomal DNA replication is fundamental to all forms of cellular life and requires the complex interplay of a wide variety of essential and non-essential protein factors in a spatially and temporally co-ordinated manner. In eukaryotes, the GINS complex (from the Japanese go-ichi-ni-san meaning 5-1-2-3, after the four related subunits of the complex Sld5, Psf1, Psf2 and Psf3) was recently identified as a novel factor essential for both the initiation and elongation stages of the replication process. Biochemical analysis has placed GINS at the heart of the eukaryotic replication apparatus as a component of the CMG [Cdc45-MCM (minichromosome maintenance) helicase-GINS] complex that most likely serves as the replicative helicase, unwinding duplex DNA ahead of the moving replication fork. GINS homologues are found in the archaea and have been shown to interact directly with the MCM helicase and with primase, suggesting a central role for the complex in archaeal chromosome replication also. The present review summarizes current knowledge of the structure, function and evolution of the GINS complex in eukaryotes and archaea, discusses possible functions of the GINS complex and highlights recent results that point to possible regulation of GINS function in response to DNA damage.

The glycine-rich motif of Pyrococcus abyssi DNA polymerase D is critical for protein stability

A glycine-rich motif described as being involved in human polymerase delta proliferating cell nuclear antigen (PCNA) binding has also been identified in all euryarchaeal DNA polymerase D (Pol D) family members. We redefined the motif as the (G)-PYF box. In the present study, Pol D (G)-PYF box motif mutants from Pyrococcus abyssi were generated to investigate its role in functional interactions with the cognate PCNA. We demonstrated that this motif is not essential for interactions between PabPol D (P. abyssi Pol D) and PCNA, using surface plasmon resonance and primer extension studies. Interestingly, the (G)-PYF box is located in a hydrophobic region close to the active site. The (G)-PYF box mutants exhibited altered DNA binding properties. In addition, the thermal stability of all mutants was reduced compared to that of wild type, and this effect could be attributed to increased exposure of the hydrophobic region. These studies suggest that the (G)-PYF box motif mediates intersubunit interactions and that it may be crucial for the thermostability of PabPol D.

Binding to PCNA in Euryarchaeal DNA Replication requires two PIP motifs for DNA polymerase D and one PIP motif for DNA polymerase B

Replicative DNA polymerases possess a canonical C-terminal proliferating cell nuclear antigen (PCNA)-binding motif termed the PCNA-interacting protein (PIP) box. We investigated the role of the PIP box on the functional interactions of the two DNA polymerases, PabPol B (family B) and PabPol D (family D), from the hyperthermophilic euryarchaeon Pyrococcus abyssi, with its cognate PCNA. The PIP box was essential for interactions of PabPol B with PCNA, as shown by surface plasmon resonance and primer extension studies. In contrast, binding of PabPol D to PCNA was affected only partially by removing the PIP motif. We identified a second palindromic PIP box motif at the N-terminus of the large subunit of PabPol D that was required for the interactions of PabPol D with PCNA. Thus, two PIP motifs were needed for PabPol D for binding to PabPCNA. Moreover, the C-terminus of PabPCNA was essential for stimulation of PabPol D activity but not for stimulation of PabPol B activity. Neither DNA polymerase interacted with the PabPCNA interdomain connecting loop. Our data suggest that distinct processes are involved in PabPol D and PabPol B binding to PCNA, raising the possibility that Archaea require two mechanisms for recruiting replicative DNA polymerases at the replication fork.

When DNA replication and protein synthesis come together

In all organisms, DNA and protein are synthesized by dedicated, but unrelated, machineries that move along distinct templates with no apparent coordination. Therefore, connections between DNA replication and translation are a priori unexpected. However, recent findings support the existence of such connections throughout the three domains of life. In particular, we recently identified in archaeal genomes a conserved association between genes encoding DNA replication and ribosome-related proteins which all have eukaryotic homologs. We believe that this gene organization is biologically relevant and, moreover, that it suggests the existence of a mechanism coupling DNA replication and translation in Archaea and Eukarya.

Gain and loss of an intron in a protein-coding gene in Archaea: the case of an archaeal RNA pseudouridine synthase gene

Background

We previously found the first examples of splicing of archaeal pre-mRNAs for homologs of the eukaryotic CBF5 protein (also known as dyskerin in humans) in Aeropyrum pernix, Sulfolobus solfataricus, S. tokodaii, and S. acidocaldarirus, and also showed that crenarchaeal species in orders Desulfurococcales and Sulfolobales, except for Hyperthermus butylicus, Pyrodictium occultum, Pyrolobus fumarii, and Ignicoccus islandicus, contain the (putative) cbf5 intron. However, the exact timing of the intron insertion was not determined and verification of the putative secondary loss of the intron in some lineages was not performed.

Results

In the present study, we determined approximately two-thirds of the entire coding region of crenarchaeal Cbf5 sequences from 43 species. A phylogenetic analysis of our data and information from the available genome sequences suggested that the (putative) cbf5 intron existed in the common ancestor of the orders Desulfurococcales and Sulfolobales and that probably at least two independent lineages in the order Desulfurococcales lost the (putative) intron.

Conclusion

This finding is the first observation of a lineage-specific loss of a pre-mRNA intron in Archaea. As the insertion or deletion of introns in protein-coding genes in Archaea has not yet been seriously considered, our finding suggests the possible difficulty of accurately and completely predicting protein-coding genes in Archaea.

Flipping of alkylated DNA damage bridges base and nucleotide excision repair

Alkyltransferase-like proteins (ATLs) share functional motifs with the cancer chemotherapy target O(6)-alkylguanine-DNA alkyltransferase (AGT) and paradoxically protect cells from the biological effects of DNA alkylation damage, despite lacking the reactive cysteine and alkyltransferase activity of AGT. Here we determine Schizosaccharomyces pombe ATL structures without and with damaged DNA containing the endogenous lesion O(6)-methylguanine or cigarette-smoke-derived O(6)-4-(3-pyridyl)-4-oxobutylguanine. These results reveal non-enzymatic DNA nucleotide flipping plus increased DNA distortion and binding pocket size compared to AGT. Our analysis of lesion-binding site conservation identifies new ATLs in sea anemone and ancestral archaea, indicating that ATL interactions are ancestral to present-day repair pathways in all domains of life. Genetic connections to mammalian XPG (also known as ERCC5) and ERCC1 in S. pombe homologues Rad13 and Swi10 and biochemical interactions with Escherichia coli UvrA and UvrC combined with structural results reveal that ATLs sculpt alkylated DNA to create a genetic and structural intersection of base damage processing with nucleotide excision repair.

Identification of a crenarchaeal orthologue of Elf1: implications for chromatin and transcription in Archaea

The transcription machineries of Archaea and eukaryotes are similar in many aspects, but little is understood about archaeal chromatin and its role in transcription. Here, we describe the identification in hyperthermophilic Crenarchaeota and a Korarchaeon of an orthologue of the eukaryotic transcription elongation factor Elf1, which has been shown to function in chromatin structure maintenance of actively transcribed templates. Our discovery has implications for the relationship of chromatin and transcription in Archaea and the evolution of these processes in eukaryotes.

Phylogeny of prokaryotes: does it exist and why should we care?

Understanding microbial evolution is essential for gathering information on the most ancient events in the history of Life on our planet. Nevertheless, the idea that it is impossible to reconstruct the evolutionary history of prokaryotes because of horizontal gene transfer has become very popular. We review this important debate and how it can be solved.

Search for a 'Tree of Life' in the thicket of the phylogenetic forest

BACKGROUND

Comparative genomics has revealed extensive horizontal gene transfer among prokaryotes, a development that is often considered to undermine the 'tree of life' concept. However, the possibility remains that a statistical central trend still exists in the phylogenetic 'forest of life'.

RESULTS

A comprehensive comparative analysis of a 'forest' of 6,901 phylogenetic trees for prokaryotic genes revealed a consistent phylogenetic signal, particularly among 102 nearly universal trees, despite high levels of topological inconsistency, probably due to horizontal gene transfer. Horizontal transfers seemed to be distributed randomly and did not obscure the central trend. The nearly universal trees were topologically similar to numerous other trees. Thus, the nearly universal trees might reflect a significant central tendency, although they cannot represent the forest completely. However, topological consistency was seen mostly at shallow tree depths and abruptly dropped at the level of the radiation of archaeal and bacterial phyla, suggesting that early phases of evolution could be non-tree-like (Biological Big Bang). Simulations of evolution under compressed cladogenesis or Biological Big Bang yielded a better fit to the observed dependence between tree inconsistency and phylogenetic depth for the compressed cladogenesis model.

CONCLUSIONS

Horizontal gene transfer is pervasive among prokaryotes: very few gene trees are fully consistent, making the original tree of life concept obsolete. A central trend that most probably represents vertical inheritance is discernible throughout the evolution of archaea and bacteria, although compressed cladogenesis complicates unambiguous resolution of the relationships between the major archaeal and bacterial clades.

Evolution of DNA polymerases: an inactivated polymerase-exonuclease module in Pol epsilon and a chimeric origin of eukaryotic polymerases from two classes of archaeal ancestors

Background

Evolution of DNA polymerases, the key enzymes of DNA replication and repair, is central to any reconstruction of the history of cellular life. However, the details of the evolutionary relationships between DNA polymerases of archaea and eukaryotes remain unresolved.

Results

We performed a comparative analysis of archaeal, eukaryotic, and bacterial B-family DNA polymerases, which are the main replicative polymerases in archaea and eukaryotes, combined with an analysis of domain architectures. Surprisingly, we found that eukaryotic Polymerase ε consists of two tandem exonuclease-polymerase modules, the active N-terminal module and a C-terminal module in which both enzymatic domains are inactivated. The two modules are only distantly related to each other, an observation that suggests the possibility that Pol ε evolved as a result of insertion and subsequent inactivation of a distinct polymerase, possibly, of bacterial descent, upstream of the C-terminal Zn-fingers, rather than by tandem duplication. The presence of an inactivated exonuclease-polymerase module in Pol ε parallels a similar inactivation of both enzymatic domains in a distinct family of archaeal B-family polymerases. The results of phylogenetic analysis indicate that eukaryotic B-family polymerases, most likely, originate from two distantly related archaeal B-family polymerases, one form giving rise to Pol ε, and the other one to the common ancestor of Pol α, Pol δ, and Pol ζ. The C-terminal Zn-fingers that are present in all eukaryotic B-family polymerases, unexpectedly, are homologous to the Zn-finger of archaeal D-family DNA polymerases that are otherwise unrelated to the B family. The Zn-finger of Polε shows a markedly greater similarity to the counterpart in archaeal PolD than the Zn-fingers of other eukaryotic B-family polymerases.

Conclusion

Evolution of eukaryotic DNA polymerases seems to have involved previously unnoticed complex events. We hypothesize that the archaeal ancestor of eukaryotes encoded three DNA polymerases, namely, two distinct B-family polymerases and a D-family polymerase all of which contributed to the evolution of the eukaryotic replication machinery. The Zn-finger might have been acquired from PolD by the B-family form that gave rise to Pol ε prior to or in the course of eukaryogenesis, and subsequently, was captured by the ancestor of the other B-family eukaryotic polymerases. The inactivated polymerase-exonuclease module of Pol ε might have evolved by fusion with a distinct polymerase, rather than by duplication of the active module of Pol ε, and is likely to play an important role in the assembly of eukaryotic replication and repair complexes.

Reviewers

This article was reviewed by Patrick Forterre, Arcady Mushegian, and Chris Ponting. For the full reviews, please go to the Reviewers' Reports section.

The origins of phagocytosis and eukaryogenesis

Background

Phagocytosis, that is, engulfment of large particles by eukaryotic cells, is found in diverse organisms and is often thought to be central to the very origin of the eukaryotic cell, in particular, for the acquisition of bacterial endosymbionts including the ancestor of the mitochondrion.

Results

Comparisons of the sets of proteins implicated in phagocytosis in different eukaryotes reveal extreme diversity, with very few highly conserved components that typically do not possess readily identifiable prokaryotic homologs. Nevertheless, phylogenetic analysis of those proteins for which such homologs do exist yields clues to the possible origin of phagocytosis. The central finding is that a subset of archaea encode actins that are not only monophyletic with eukaryotic actins but also share unique structural features with actin-related proteins (Arp) 2 and 3. All phagocytic processes are strictly dependent on remodeling of the actin cytoskeleton and the formation of branched filaments for which Arp2/3 are responsible. The presence of common structural features in Arp2/3 and the archaeal actins suggests that the common ancestors of the archaeal and eukaryotic actins were capable of forming branched filaments, like modern Arp2/3. The Rho family GTPases that are ubiquitous regulators of phagocytosis in eukaryotes appear to be of bacterial origin, so assuming that the host of the mitochondrial endosymbiont was an archaeon, the genes for these GTPases come via horizontal gene transfer from the endosymbiont or in an earlier event.

Conclusion

The present findings suggest a hypothetical scenario of eukaryogenesis under which the archaeal ancestor of eukaryotes had no cell wall (like modern Thermoplasma) but had an actin-based cytoskeleton including branched actin filaments that allowed this organism to produce actin-supported membrane protrusions. These protrusions would facilitate accidental, occasional engulfment of bacteria, one of which eventually became the mitochondrion. The acquisition of the endosymbiont triggered eukaryogenesis, in particular, the emergence of the endomembrane system that eventually led to the evolution of modern-type phagocytosis, independently in several eukaryotic lineages.

Reviewers

This article was reviewed by Simonetta Gribaldo, Gaspar Jekely, and Pierre Pontarotti. For the full reviews, please go to the Reviewers' Reports section.

Integration of phenotypic metadata and protein similarity in Archaea using a spectral bipartitioning approach

In order to simplify and meaningfully categorize large sets of protein sequence data, it is commonplace to cluster proteins based on the similarity of those sequences. However, it quickly becomes clear that the sequence flexibility allowed a given protein varies significantly among different protein families. The degree to which sequences are conserved not only differs for each protein family, but also is affected by the phylogenetic divergence of the source organisms. Clustering techniques that use similarity thresholds for protein families do not always allow for these variations and thus cannot be confidently used for applications such as automated annotation and phylogenetic profiling. In this work, we applied a spectral bipartitioning technique to all proteins from 53 archaeal genomes. Comparisons between different taxonomic levels allowed us to study the effects of phylogenetic distances on cluster structure. Likewise, by associating functional annotations and phenotypic metadata with each protein, we could compare our protein similarity clusters with both protein function and associated phenotype. Our clusters can be analyzed graphically and interactively online.

Role of multiprotein bridging factor 1 in archaea: bridging the domains?

MBF1 (multiprotein bridging factor 1) is a highly conserved protein in archaea and eukaryotes. It was originally identified as a mediator of the eukaryotic transcription regulator BmFTZ-F1 (Bombyx mori regulator of fushi tarazu). MBF1 was demonstrated to enhance transcription by forming a bridge between distinct regulatory DNA-binding proteins and the TATA-box-binding protein. MBF1 consists of two parts: a C-terminal part that contains a highly conserved helix-turn-helix, and an N-terminal part that shows a clear divergence: in eukaryotes, it is a weakly conserved flexible domain, whereas, in archaea, it is a conserved zinc-ribbon domain. Although its function in archaea remains elusive, its function as a transcriptional co-activator has been deduced from thorough studies of several eukaryotic proteins, often indicating a role in stress response. In addition, MBF1 was found to influence translation fidelity in yeast. Genome context analysis of mbf1 in archaea revealed conserved clustering in the crenarchaeal branch together with genes generally involved in gene expression. It points to a role of MBF1 in transcription and/or translation. Experimental data are required to allow comparison of the archaeal MBF1 with its eukaryotic counterpart.

The making of tRNAs and more - RNase P and tRNase Z

Transfer-RNA (tRNA) molecules are essential players in protein biosynthesis. They are transcribed as precursors, which have to be extensively processed at both ends to become functional adaptors in protein synthesis. Two endonucleases that directly interact with the tRNA moiety, RNase P and tRNase Z, remove extraneous nucleotides on the molecule's 5'- and 3'-side, respectively. The ribonucleoprotein enzyme RNase P was identified almost 40 years ago and is considered a vestige from the "RNA world". Here, we present the state of affairs on prokaryotic RNase P, with a focus on recent findings on its role in RNA metabolism. tRNase Z was only identified 6 years ago, and we do not yet have a comprehensive understanding of its function. The current knowledge on prokaryotic tRNase Z in tRNA 3'-processing is reviewed here. A second, tRNase Z-independent pathway of tRNA 3'-end maturation involving 3'-exonucleases will also be discussed.

A genomic analysis of the archaeal system Ignicoccus hospitalis-Nanoarchaeum equitans

Sequencing of the complete genome of Ignicoccus hospitalis gives insight into its association with another species of Archaea, Nanoarchaeum equitans.

Background

The relationship between the hyperthermophiles Ignicoccus hospitalis and Nanoarchaeum equitans is the only known example of a specific association between two species of Archaea. Little is known about the mechanisms that enable this relationship.

Results

We sequenced the complete genome of I. hospitalis and found it to be the smallest among independent, free-living organisms. A comparative genomic reconstruction suggests that the I. hospitalis lineage has lost most of the genes associated with a heterotrophic metabolism that is characteristic of most of the Crenarchaeota. A streamlined genome is also suggested by a low frequency of paralogs and fragmentation of many operons. However, this process appears to be partially balanced by lateral gene transfer from archaeal and bacterial sources.

Conclusions

A combination of genomic and cellular features suggests highly efficient adaptation to the low energy yield of sulfur-hydrogen respiration and efficient inorganic carbon and nitrogen assimilation. Evidence of lateral gene exchange between N. equitans and I. hospitalis indicates that the relationship has impacted both genomes. This association is the simplest symbiotic system known to date and a unique model for studying mechanisms of interspecific relationships at the genomic and metabolic levels.

TBP domain symmetry in basal and activated archaeal transcription

The TATA box binding protein (TBP) is the platform for assembly of archaeal and eukaryotic transcription preinitiation complexes. Ancestral gene duplication and fusion events have produced the saddle-shaped TBP molecule, with its two direct-repeat subdomains and pseudo-two-fold symmetry. Collectively, eukaryotic TBPs have diverged from their present-day archaeal counterparts, which remain highly symmetrical. The similarity of the N- and C-halves of archaeal TBPs is especially pronounced in the Methanococcales and Thermoplasmatales, including complete conservation of their N- and C-terminal stirrups; along with helix H'1, the C-terminal stirrup of TBP forms the main interface with TFB/TFIIB. Here, we show that, in stark contrast to its eukaryotic counterparts, multiple substitutions in the C-terminal stirrup of Methanocaldococcus jannaschii (Mja) TBP do not completely abrogate basal transcription. Using DNA affinity cleavage, we show that, by assembling TFB through its conserved N-terminal stirrup, Mja TBP is in effect ambidextrous with regard to basal transcription. In contrast, substitutions in either its N- or the C-terminal stirrup abrogate activated transcription in response to the Lrp-family transcriptional activator Ptr2.

The deep archaeal roots of eukaryotes

The set of conserved eukaryotic protein-coding genes includes distinct subsets one of which appears to be most closely related to and, by inference, derived from archaea, whereas another one appears to be of bacterial, possibly, endosymbiotic origin. The “archaeal” genes of eukaryotes, primarily, encode components of information-processing systems, whereas the “bacterial” genes are predominantly operational. The precise nature of the archaeo–eukaryotic relationship remains uncertain, and it has been variously argued that eukaryotic informational genes evolved from the homologous genes of Euryarchaeota or Crenarchaeota (the major branches of extant archaea) or that the origin of eukaryotes lies outside the known diversity of archaea. We describe a comprehensive set of 355 eukaryotic genes of apparent archaeal origin identified through ortholog detection and phylogenetic analysis. Phylogenetic hypothesis testing using constrained trees, combined with a systematic search for shared derived characters in the form of homologous inserts in conserved proteins, indicate that, for the majority of these genes, the preferred tree topology is one with the eukaryotic branch placed outside the extant diversity of archaea although small subsets of genes show crenarchaeal and euryarchaeal affinities. Thus, the archaeal genes in eukaryotes appear to descend from a distinct, ancient, and otherwise uncharacterized archaeal lineage that acquired some euryarchaeal and crenarchaeal genes via early horizontal gene transfer.