An archaeal genomic signature

The discovery of archaea: from observed anomaly to consequential restructuring of the phylogenetic tree

Observational and experimental discoveries of new factual entities such as objects, systems, or processes, are major contributors to some advances in the life sciences. Yet, whereas discovery of theories was extensively deliberated by philosophers of science, very little philosophical attention was paid to the discovery of factual entities. This paper examines historical and philosophical aspects of the experimental discovery by Carl Woese of archaea, prokaryotes that comprise one of the three principal domains of the phylogenetic tree. Borrowing Kuhn's terminology, this discovery of a major biological entity was made during a 'normal science' project of building molecular taxonomy for prokaryotes. Unexpectedly, however, an observed anomaly instigated the discovery of archaea. Substantiation of the existence of the new archaeal entity and consequent reconstruction of the phylogenetic tree prompted replacement of a long-held model of a prokarya and eukarya bipartite tree of life by a new model of a tripartite tree comprising of bacteria, archaea, and eukarya. This paper explores the history and philosophical implications of the progression of Woese's project from normal science to anomaly-instigated model-changing discovery. It is also shown that the consequential discoveries of RNA splicing and of ribozymes were similarly prompted by unexpected irregularities during normal science activities. It is thus submitted that some discoveries of factual biological entities are triggered by unforeseen observational or experimental anomalies.

Phylogenomic Analysis of Metagenome-Assembled Genomes Deciphered Novel Acetogenic Nitrogen-Fixing Bathyarchaeota from Hot Spring Sediments

This study describes the phylogenomic analysis and metabolic insights of metagenome-assembled genomes (MAGs) retrieved from hot spring sediment samples. The metagenome-assembled sequences recovered three near-complete genomes belonging to the archaeal phylum. Analysis of genome-wide core genes and 16S rRNA-based phylogeny placed the ILS200 and ILS300 genomes within the uncultivated and largely understudied bathyarchaeal phylum, whereas ILS100 represented the phylum Thaumarchaeota. The average nucleotide identity (ANI) of the bin ILS100 was 76% with Nitrososphaeria_archaeon_isolate_SpSt-1069. However, the bins ILS200 and ILS300 showed ANI values of 75% and 70% with Candidatus_Bathyarchaeota_archaeon_isolate_DRTY-6_2_bin_115 and Candidatus_Bathyarchaeota_archaeon_BA1_ba1_01, respectively. The genomic potential of Bathyarchaeota bins ILS200 and ILS300 showed genes necessary for the Wood-Ljungdahl pathway, and the gene encoding the methyl coenzyme M reductase (mcr) complex essential for methanogenesis was absent. The metabolic potential of the assembled genomes included genes involved in nitrogen assimilation, including nitrogenase and the genes necessary for the urea cycle. The presence of these genes suggested the metabolic potential of Bathyarchaeota to fix nitrogen under extreme environments. In addition, the ILS200 and ILS300 genomes carried genes involved in the tricarboxylic acid (TCA) cycle, glycolysis, and degradation of organic carbons. Finally, we conclude that the reconstructed Bathyarchaeota bins are autotrophic acetogens and organo-heterotrophs.

IMPORTANCE We describe the Bathyarchaeota bins that are likely to be acetogens with a wide range of metabolic potential. These bins did not exhibit methanogenic machinery, suggesting methane production may not occur by all subgroup lineages of Bathyarchaeota. Phylogenetic analysis support that both ILS200 and ILS300 belonged to the Bathyarchaeota. The discovery of new bathyarchaeotal MAGs provides additional knowledge for understanding global carbon and nitrogen metabolism under extreme conditions.

What is an archaeon and are the Archaea really unique?

The recognition of the group Archaea as a major branch of the tree of life (ToL) prompted a new view of the evolution of biodiversity. The genomic representation of archaeal biodiversity has since significantly increased. In addition, advances in phylogenetic modeling of multi-locus datasets have resolved many recalcitrant branches of the ToL. Despite the technical advances and an expanded taxonomic representation, two important aspects of the origins and evolution of the Archaea remain controversial, even as we celebrate the 40th anniversary of the monumental discovery. These issues concern (i) the uniqueness (monophyly) of the Archaea, and (ii) the evolutionary relationships of the Archaea to the Bacteria and the Eukarya; both of these are relevant to the deep structure of the ToL. To explore the causes for this persistent ambiguity, I examine multiple datasets and different phylogenetic approaches that support contradicting conclusions. I find that the uncertainty is primarily due to a scarcity of information in standard datasets-universal core-genes datasets-to reliably resolve the conflicts. These conflicts can be resolved efficiently by comparing patterns of variation in the distribution of functional genomic signatures, which are less diffused unlike patterns of primary sequence variation. Relatively lower heterogeneity in distribution patterns minimizes uncertainties and supports statistically robust phylogenetic inferences, especially of the earliest divergences of life. This case study further highlights the limitations of primary sequence data in resolving difficult phylogenetic problems, and raises questions about evolutionary inferences drawn from the analyses of sequence alignments of a small set of core genes. In particular, the findings of this study corroborate the growing consensus that reversible substitution mutations may not be optimal phylogenetic markers for resolving early divergences in the ToL, nor for determining the polarity of evolutionary transitions across the ToL.

The Methanosarcina mazei MM2060 Gene Encodes a Bifunctional Kinase/Decarboxylase Enzyme Involved in Cobamide Biosynthesis

Cobamides (Cbas) are synthesized by many archaea, but some aspects of Cba biosynthesis in these microorganisms remain unclear. Here, we demonstrate that open reading frame MM2060 in the archaeum Methanosarcina mazei strain Gö1 encodes a bifunctional enzyme with l-threonine- O-3-phosphate (l-Thr-P) decarboxylase (EC 4.1.1.81) and l-Thr kinase activities (EC 2.7.1.177). In Salmonella enterica, where Cba biosynthesis has been extensively studied, the activities mentioned above are encoded by separate genes, namely, cobD and pduX, respectively. The activities associated with the MM2060 protein ( MmCobD) were validated in vitro and in vivo. In vitro, MmCobD used ATP and l-Thr as substrates and generated ADP, l-Thr-P, and ( R)-1-aminopropan-2-ol O-phosphate as products. Notably, MmCobD has a 111-amino acid C-terminal extension of unknown function, which contains a putative metal-binding motif. This C-terminal domain alone did not display activity either in vivo or in vitro. Although the C-terminal MmCobD domain was not required for l-Thr-P decarboxylase or l-Thr kinase activities in vivo, its absence negatively affected both activities. In vitro results suggested that this domain may have a regulatory or substrate-gating role. When purified under anoxic conditions, MmCobD displayed Michaelis-Menten kinetics and had a 1000-fold higher affinity for ATP and a catalytic efficiency 1300-fold higher than that of MmCobD purified under oxic conditions. To the best of our knowledge, MmCobD is the first example of a new class of l-Thr-P decarboxylases that also have l-Thr kinase activity. An archaeal protein with l-Thr kinase activity had not been identified prior to this work.

Comparative Analysis of Predicted Gene Expression among Crenarchaeal Genomes

Research into new methods for identifying highly expressed genes in anonymous genome sequences has been going on for more than 15 years. We presented here an alternative approach based on modified score of relative codon usage bias to identify highly expressed genes in crenarchaeal genomes. The proposed algorithm relies exclusively on sequence features for identifying the highly expressed genes. In this study, a comparative analysis of predicted highly expressed genes in five crenarchaeal genomes was performed using the score of Modified Relative Codon Bias Strength (MRCBS) as a numerical estimator of gene expression level. We found a systematic strong correlation between Codon Adaptation Index and MRCBS. Additionally, MRCBS correlated well with other expression measures. Our study indicates that MRCBS can consistently capture the highly expressed genes.

Anaerobic Metabolism in Haloferax Genus: Denitrification as Case of Study

A number of species of Haloferax genus (halophilic archaea) are able to grow microaerobically or even anaerobically using different alternative electron acceptors such as fumarate, nitrate, chlorate, dimethyl sulphoxide, sulphide and/or trimethylamine. This metabolic capability is also shown by other species of the Halobacteriaceae and Haloferacaceae families (Archaea domain) and it has been mainly tested by physiological studies where cell growth is observed under anaerobic conditions in the presence of the mentioned compounds. This work summarises the main reported features on anaerobic metabolism in the Haloferax, one of the better described haloarchaeal genus with significant potential uses in biotechnology and bioremediation. Special attention has been paid to denitrification, also called nitrate respiration. This pathway has been studied so far from Haloferax mediterranei and Haloferax denitrificans mainly from biochemical point of view (purification and characterisation of the enzymes catalysing the two first reactions). However, gene expression and gene regulation is far from known at the time of writing this chapter.

Environmental genes and genomes: understanding the differences and challenges in the approaches and software for their analyses

DNA-based taxonomic and functional profiling is widely used for the characterization of organismal communities across a rapidly increasing array of research areas that include the role of microbiomes in health and disease, biomonitoring, and estimation of both microbial and metazoan species richness. Two principal approaches are currently used to assign taxonomy to DNA sequences: DNA metabarcoding and metagenomics. When initially developed, each of these approaches mandated their own particular methods for data analysis; however, with the development of high-throughput sequencing (HTS) techniques they have begun to share many aspects in data set generation and processing. In this review we aim to define the current characteristics, goals and boundaries of each field, and describe the different software used for their analysis. We argue that an appreciation of the potential and limitations of each method can help underscore the improvements required by each field so as to better exploit the richness of current HTS-based data sets.

A backward view from 16S rRNA to archaea to the universal tree of life to progenotes: reminiscences of Carl Woese

I first became aware of Carl Woese in the mid-1970s when he and George Fox criticized a few of the 16S rRNA oligonucleotide sequences emerging from Strasbourg in the 10-12 y RNA sequencing project of the first 16S rRNA from Escherichia coli, some of which we were using for assembling RNA binding sites of ribosomal proteins. When I realized that they were attempting to sequence 16S rRNAs from a range of bacteria to classify them phylogenetically, I seriously questioned their sanity. Not because of the goal, which was admirable, but because of the sheer technical difficulty, and slowness, of sequencing large RNA molecules using the original Sanger RNA sequencing method, not to mention the health hazards of regularly preparing rRNA using 20-30 mCi [ (32)P]. My view changed radically, however, with their subsequent prediction of 5S rRNA secondary structures using a phylogenetic approach. Previously, the molecular biology community had been competing to generate the maximum numbers of base pairs in the model RNA molecule E. coli 5S RNA when Fox and Woese introduced the concept of compensatory base changes based on phylogeny for defining secondary structure and applied it to 5S RNA, they found evidence for only about 50% base pairing. This approach had previously been used for tRNA secondary structure predictions but its more general significance had never been acknowledged. Carl subsequently persuaded Harry Noller to apply the same method to predicting secondary structures of the large rRNAs.

[Gut microbiota in health and disease]

Gut microbiota is the community of live microorganisms residing in the digestive tract. There are many groups of researchers worldwide that are working at deciphering the collective genome of the human microbiota. Modern techniques for studying the microbiota have made us aware of an important number of nonculturable bacteria and of the relation between the microorganisms that live inside us and our homeostasis. The microbiota is essential for correct body growth, the development of immunity, and nutrition. Certain epidemics affecting humanity such as asthma and obesity may possibly be explained, at least partially, by alterations in the microbiota. Dysbiosis has been associated with a series of gastrointestinal disorders that include non-alcoholic fatty liver disease, celiac disease, and irritable bowel syndrome. The present article deals with the nomenclature, modern study techniques, and functions of gut microbiota, and its relation to health and disease.

Carl woese: from biophysics to evolutionary microbiology

This article is a tribute to Carl R. Woese, a biophysicist turned evolutionary microbiologist who passed away on December 30, 2012. We focus on his life, achievements, the discovery of Archaea and contributions to the development of molecular phylogeny. Further, the authors share their views and the lessons learnt from Woese's life with the microbiologists in India. We also emphasize the need for interdisciplinary collaboration and interaction for the progress and betterment of science.

Proteomic properties reveal phyloecological clusters of Archaea

In this study, we propose a novel way to describe the variety of environmental adaptations of Archaea. We have clustered 57 Archaea by using a non-redundant set of proteomic features, and verified that the clusters correspond to environmental adaptations to the archaeal habitats. The first cluster consists dominantly of hyperthermophiles and hyperthermoacidophilic aerobes. The second cluster joins together halophilic and extremely halophilic Archaea, while the third cluster contains mesophilic (mostly methanogenic) Archaea together with thermoacidophiles. The non-redundant subset of proteomic features was found to consist of five features: the ratio of charged residues to uncharged, average protein size, normalized frequency of beta-sheet, normalized frequency of extended structure and number of hydrogen bond donors. We propose this clustering to be termed phyloecological clustering. This approach could give additional insights into relationships among archaeal species that may be hidden by sole phylogenetic analysis.

Methanogenic archaea in subgingival sites: a review

Archaea are non-bacterial prokaryotes associated with oral microbiota in humans, but their roles in oral pathologies remain controversial. Several studies reported the molecular detection of methanogenic archaea from periodontitis, but the significance of this association has not been confirmed yet. An electronic search was therefore conducted in MEDLINE-Pubmed to identify all papers published in English connecting archaea and periodontal infections. Data analysis of the selected studies showed that five genera of methanogenic archaea have been detected in the subgingival microbiota, Methanobrevibacter oralis being the most frequently detected species in 41% of periodontitis patients and 55% of periodontal pockets compared to 6% of healthy subjects and 5% of periodontally-healthy sites (p < 10(-5) , Chi-squared test). Based on the five determination-criteria proposed by Socransky (association with disease, elimination of the organism, host response, animal pathogenicity and mechanisms of pathogenicity), M. oralis is a periodontal pathogen. The methanogenic archaea load correlating with periodontitis severity further supports the pathogenic role of methanogenic archaea in periodontitis. Therefore, detection and quantification of M. oralis in periodontal pockets could help the laboratory diagnosis and follow-up of periodontitis. Determining the origin, diversity and pathogenesis of archaea in periodontal infections warrants further investigations.

Selaginella Genome Analysis - Entering the "Homoplasy Heaven" of the MADS World

In flowering plants, arguably the most significant transcription factors regulating development are MADS-domain proteins, encoded by Type I and Type II MADS-box genes. Type II genes are divided into the MIKC(C) and MIKC* groups. In angiosperms, these types and groups play distinct roles in the development of female gametophytes, embryos, and seeds (Type I); vegetative and floral tissues in sporophytes (MIKC(C)); and male gametophytes (MIKC*), but their functions in other plants are largely unknown. The complete set of MADS-box genes has been described for several angiosperms and a moss, Physcomitrella patens. Our examination of the complete genome sequence of a lycophyte, Selaginella moellendorffii, revealed 19 putative MADS-box genes (13 Type I, 3 MIKC(C), and 3 MIKC*). Our results suggest that the most recent common ancestor of vascular plants possessed at least two Type I and two Type II genes. None of the S. moellendorffii MIKC(C) genes were identified as orthologs of any floral organ identity genes. This strongly corroborates the view that the clades of floral organ identity genes originated in a common ancestor of seed plants after the lineage that led to lycophytes had branched off, and that expansion of MIKC(C) genes in the lineage leading to seed plants facilitated the evolution of their unique reproductive organs. The number of MIKC* genes and the ratio of MIKC* to MIKC(C) genes is lower in S. moellendorffii and angiosperms than in P. patens, correlated with reduction of the gametophyte in vascular plants. Our data indicate that Type I genes duplicated and diversified independently within lycophytes and seed plants. Our observations on MADS-box gene evolution echo morphological evolution since the two lineages of vascular plants appear to have arrived independently at similar body plans. Our annotation of MADS-box genes in S. moellendorffii provides the basis for functional studies to reveal the roles of this crucial gene family in basal vascular plants.

The core and unique proteins of haloarchaea

Background

Since the first genome of a halophilic archaeon was sequenced in 2000, biologists have been advancing the understanding of genomic characteristics that allow for survival in the harsh natural environments of these organisms. An increase in protein acidity and GC-bias in the genome have been implicated as factors in tolerance to extreme salinity, desiccation, and high solar radiation. However, few previous attempts have been made to identify novel genes that would permit survival in such extreme conditions.

Results

With the recent release of several new complete haloarchaeal genome sequences, we have conducted a comprehensive comparative genomic analysis focusing on the identification of unique haloarchaeal conserved proteins that likely play key roles in environmental adaptation. Using bioinformatic methods, we have clustered 31,312 predicted proteins from nine haloarchaeal genomes into 4,455 haloarchaeal orthologous groups (HOGs). We assigned likely functions by association with established COG and KOG databases in NCBI. After identifying homologs in four additional haloarchaeal genomes, we determined that there were 784 core haloarchaeal protein clusters (cHOGs), of which 83 clusters were found primarily in haloarchaea. Further analysis found that 55 clusters were truly unique (tucHOGs) to haloarchaea and qualify as signature proteins while 28 were nearly unique (nucHOGs), the vast majority of which were coded for on the haloarchaeal chromosomes. Of the signature proteins, only one example with any predicted function, Ral, involved in desiccation/radiation tolerance in Halobacterium sp. NRC-1, was identified. Among the core clusters, 33% was predicted to function in metabolism, 25% in information transfer and storage, 10% in cell processes and signaling, and 22% belong to poorly characterized or general function groups.

Conclusion

Our studies have established conserved groups of nearly 800 protein clusters present in all haloarchaea, with a subset of 55 which are predicted to be accessory proteins that may be critical or essential for success in an extreme environment. These studies support core and signature genes and proteins as valuable concepts for understanding phylogenetic and phenotypic characteristics of coherent groups of organisms.

Major players on the microbial stage: why archaea are important

As microbiology undergoes a renaissance, fuelled in part by developments in new sequencing technologies, the massive diversity and abundance of microbes becomes yet more obvious. The Archaea have traditionally been perceived as a minor group of organisms forced to evolve into environmental niches not occupied by their more 'successful' and 'vigorous' counterparts, the bacteria. Here we outline some of the evidence gathered by an increasingly large and productive group of scientists that demonstrates not only that the Archaea contribute significantly to global nutrient cycling, but also that they compete successfully in 'mainstream' environments. Recent data suggest that the Archaea provide the major routes for ammonia oxidation in the environment. Archaea also have huge economic potential that to date has only been fully realized in the production of thermostable polymerases. Archaea have furnished us with key paradigms for understanding fundamentally conserved processes across all domains of life. In addition, they have provided numerous exemplars of novel biological mechanisms that provide us with a much broader view of the forms that life can take and the way in which micro-organisms can interact with other species. That this information has been garnered in a relatively short period of time, and appears to represent only a small proportion of what the Archaea have to offer, should provide further incentives to microbiologists to investigate the underlying biology of this fascinating domain.

Genome networks root the tree of life between prokaryotic domains

Eukaryotes arose from prokaryotes, hence the root in the tree of life resides among the prokaryotic domains. The position of the root is still debated, although pinpointing it would aid our understanding of the early evolution of life. Because prokaryote evolution was long viewed as a tree-like process of lineage bifurcations, efforts to identify the most ancient microbial lineage split have traditionally focused on positioning a root on a phylogenetic tree constructed from one or several genes. Such studies have delivered widely conflicting results on the position of the root, this being mainly due to methodological problems inherent to deep gene phylogeny and the workings of lateral gene transfer among prokaryotes over evolutionary time. Here, we report the position of the root determined with whole genome data using network-based procedures that take into account both gene presence or absence and the level of sequence similarity among all individual gene families that are shared across genomes. On the basis of 562,321 protein-coding gene families distributed across 191 genomes, we find that the deepest divide in the prokaryotic world is interdomain, that is, separating the archaebacteria from the eubacteria. This result resonates with some older views but conflicts with the results of most studies over the last decade that have addressed the issue. In particular, several studies have suggested that the molecular distinctness of archaebacteria is not evidence for their antiquity relative to eubacteria but instead stems from some kind of inherently elevated rate of archaebacterial sequence change. Here, we specifically test for such a rate elevation across all prokaryotic lineages through the analysis of all possible quartets among eight genes duplicated in all prokaryotes, hence the last common ancestor thereof. The results show that neither the archaebacteria as a group nor the eubacteria as a group harbor evidence for elevated evolutionary rates in the sampled genes, either in the recent evolutionary past or in their common ancestor. The interdomain prokaryotic position of the root is thus not attributable to lineage-specific rate variation.

Tight coupling of partial reactions in the acetyl-CoA decarbonylase/synthase (ACDS) multienzyme complex from Methanosarcina thermophila: acetyl C-C bond fragmentation at the a cluster promoted by protein conformational changes

Direct synthesis and cleavage of acetyl-CoA are carried out by the bifunctional CO dehydrogenase/acetyl-CoA synthase enzyme in anaerobic bacteria and by the acetyl-CoA decarbonylase/synthase (ACDS) multienzyme complex in Archaea. In both systems, a nickel- and Fe/S-containing active site metal center, the A cluster, catalyzes acetyl C-C bond formation/breakdown. Carbonyl group exchange of [1-(14)C]acetyl-CoA with unlabeled CO, a hallmark of CODH/ACS, is weakly active in ACDS, and exchange with CO(2) was up to 350 times faster, indicating tight coupling of CO release at the A cluster to CO oxidation to CO(2) at the C cluster in CO dehydrogenase. The basis for tight coupling was investigated by analysis of three recombinant A cluster proteins, ACDS beta subunit from Methanosarcina thermophila, acetyl-CoA synthase of Carboxydothermus hydrogenoformans (ACS(Ch)), and truncated ACS(Ch) lacking its 317-amino acid N-terminal domain. A comparison of acetyl-CoA synthesis kinetics, CO exchange, acetyltransferase, and A cluster Ni(+)-CO EPR characteristics demonstrated a direct role of the ACS N-terminal domain in promoting acetyl C-C bond fragmentation. Protein conformational changes, related to "open/closed" states previously identified crystallographically, were indicated to have direct effects on the coordination geometry and stability of the A cluster Ni(2+)-acetyl intermediate, controlling Ni(2+)-acetyl fragmentation and Ni(2+)(CO)(CH(3)) condensation. EPR spectral changes likely reflect variations in the Ni(+)-CO equatorial coordination environment in closed buried hydrophobic and open solvent-exposed states. The involvement of subunit-subunit interactions in ACDS, versus interdomain contacts in ACS, ensures that CO is not released from the ACDS beta subunit in the absence of appropriate interactions with the alpha(2)epsilon(2) CO dehydrogenase component. The resultant high efficiency CO transfer explains the low rate of CO exchange relative to CO(2).

An iron(II) dependent formamide hydrolase catalyzes the second step in the archaeal biosynthetic pathway to riboflavin and 7,8-didemethyl-8-hydroxy-5-deazariboflavin

The early steps in the biosynthesis of 7,8-didemethyl-8-hydroxy-5-deazariboflavin (Fo) and riboflavin in the archaea differ from the established eukaryotic and bacterial pathways. The archaeal pathway has been proposed to begin with an archaeal-specific GTP cyclohydrolase III that hydrolyzes the imidazole ring of GTP but does not remove the resulting formyl group from the formamide [Graham, D. E., Xu, H., and White, R. H. (2002) Biochemistry 41, 15074-15084 ]. This enzyme is different than the bacterial GTP cyclohydrolase II which catalyzes both reactions. Here we describe the identification and characterization of the formamide hydrolase that catalyzes the second step in the archaeal Fo and riboflavin biosynthetic pathway. The Methanocaldococcus jannaschii MJ0116 gene was cloned and heterologously expressed, and the resulting enzyme was shown to catalyze the formation of 2,5-diamino-6-ribosylamino-4(3H)-pyrimidinone 5'-phosphate (APy) and formate from 2-amino-5-formylamino-6-ribosylamino-4(3H)-pyrimidinone 5'-monophosphate (FAPy). The MJ0116-derived protein has been named ArfB to indicate that it catalyzes the second step in archaeal riboflavin and Fo biosynthesis. ArfB was found to require ferrous iron for activity although metal analysis by ICP indicated the presence of zinc as well as iron in the purified protein. The identification of this enzyme confirms the involvement of GTP cyclohydrolase III (ArfA) in archaeal riboflavin and Fo biosynthesis.

How the microbial world saved evolution from the scylla of molecular biology and the charybdis of the modern synthesis

In this commentary, we provide a personal overview of the conceptual history of microbiology and molecular biology over the course of the last hundred years, emphasizing the relationship of these fields to the problem of evolution. We argue that despite their apparent success, all three reached an impasse that arose from the influence of dogmatic or overly narrow perspectives. Finally, we describe how recent developments in microbiology are realizing Beijerinck's vision of a field that is fully integrated with molecular biology, microbial ecology, thereby challenging and extending current thinking in evolution.

Inversion-based genomic signatures

Background

Reconstructing complete ancestral genomes (at least in terms of their gene inventory and arrangement) is attracting much interest due to the rapidly increasing availability of whole genome sequences. While modest successes have been reported for mammalian and even vertebrate genomes, more divergent groups continue to pose a stiff challenge, mostly because current models of genomic evolution support too many choices.

Results

We describe a novel type of genomic signature based on rearrangements that characterizes evolutionary changes that must be common to all minimal rearrangement scenarios; by focusing on global patterns of rearrangements, such signatures bypass individual variations and sharply restrict the search space. We present the results of extensive simulation studies demonstrating that these signatures can be used to reconstruct accurate ancestral genomes and phylogenies even for widely divergent collections.

Conclusion

Focusing on genome triples rather than genomes pairs unleashes the full power of evolutionary analysis. Our genomic signature captures shared evolutionary events and thus can form the basis of a robust analysis and reconstruction of evolutionary history.

On the evolution of the tRNA-dependent amidotransferases, GatCAB and GatDE

Glutaminyl-tRNA synthetase and asparaginyl-tRNA synthetase evolved from glutamyl-tRNA synthetase and aspartyl-tRNA synthetase, respectively, after the split in the last universal communal ancestor (LUCA). Glutaminyl-tRNA(Gln) and asparaginyl-tRNA(Asn) were likely formed in LUCA by amidation of the mischarged species, glutamyl-tRNA(Gln) and aspartyl-tRNA(Asn), by tRNA-dependent amidotransferases, as is still the case in most bacteria and all known archaea. The amidotransferase GatCAB is found in both domains of life, while the heterodimeric amidotransferase GatDE is found only in Archaea. The GatB and GatE subunits belong to a unique protein family that includes Pet112 that is encoded in the nuclear genomes of numerous eukaryotes. GatE was thought to have evolved from GatB after the emergence of the modern lines of decent. Our phylogenetic analysis though places the split between GatE and GatB, prior to the phylogenetic divide between Bacteria and Archaea, and Pet112 to be of mitochondrial origin. In addition, GatD appears to have emerged prior to the bacterial-archaeal phylogenetic divide. Thus, while GatDE is an archaeal signature protein, it likely was present in LUCA together with GatCAB. Archaea retained both amidotransferases, while Bacteria emerged with only GatCAB. The presence of GatDE has favored a unique archaeal tRNA(Gln) that may be preventing the acquisition of glutaminyl-tRNA synthetase in Archaea. Archaeal GatCAB, on the other hand, has not favored a distinct tRNA(Asn), suggesting that tRNA(Asn) recognition is not a major barrier to the retention of asparaginyl-tRNA synthetase in many Archaea.

Enzyme-Driven Speciation: Crystallizing Archaea via Lipid Capture

As the origin(s) of life on Earth remains an open question, detailed characteristics about the "last universal ancestor" (LUA) continue to be obscured. Here we provide arguments that strengthen the bacterial-like nature of the LUA. Our view attempts to recreate the evolution of archaeal lipids, the major components of the distinctive membrane that encapsulates these ancient prokaryotes. We show that (S)- 3-O-geranylgeranylglyceryl phosphate synthase (GGGPS), a TIM-barrel protein that performs the committed step in archaeal lipid synthesis, likely evolved from the duplication and fusion of a (betaalpha)4 half-barrel ancestor. By comparison to the well-characterized HisA and HisF TIM-barrel proteins, we propose a time line for the invention of this diagnostic archaeal biosynthetic pathway. After excluding the possibility of horizontal gene transfer, we conclude that the evolutionary history of GGGPS mirrors the emergence of Archaea from the LUA. We illustrate aspects of this "lipid capture" model that support its likelihood in recreating key evolutionary events and, as our hypothesis is built on a single initiating event, we suggest that the appearance of GGGPS represents an example of enzyme-driven speciation.

Accurate phylogenetic classification of variable-length DNA fragments

Metagenome studies have retrieved vast amounts of sequence data from a variety of environments leading to new discoveries and insights into the uncultured microbial world. Except for very simple communities, the encountered diversity has made fragment assembly and the subsequent analysis a challenging problem. A taxonomic characterization of metagenomic fragments is required for a deeper understanding of shotgun-sequenced microbial communities, but success has mostly been limited to sequences containing phylogenetic marker genes. Here we present PhyloPythia, a composition-based classifier that combines higher-level generic clades from a set of 340 completed genomes with sample-derived population models. Extensive analyses on synthetic and real metagenome data sets showed that PhyloPythia allows the accurate classification of most sequence fragments across all considered taxonomic ranks, even for unknown organisms. The method requires no more than 100 kb of training sequence for the creation of accurate models of sample-specific populations and can assign fragments >or=1 kb with high specificity.

The cbiS gene of the archaeon Methanopyrus kandleri AV19 encodes a bifunctional enzyme with adenosylcobinamide amidohydrolase and alpha-ribazole-phosphate phosphatase activities

Here we report the initial biochemical characterization of the bifunctional alpha-ribazole-P (alpha-RP) phosphatase, adenosylcobinamide (AdoCbi) amidohydrolase CbiS enzyme from the hyperthermophilic methanogenic archaeon Methanopyrus kandleri AV19. The cbiS gene encodes a 39-kDa protein with two distinct segments, one of which is homologous to the AdoCbi amidohydrolase (CbiZ, EC 3.5.1.90) enzyme and the other of which is homologous to the recently discovered archaeal alpha-RP phosphatase (CobZ, EC 3.1.3.73) enzyme. CbiS function restored AdoCbi salvaging and alpha-RP phosphatase activity in strains of the bacterium Salmonella enterica where either step was blocked. The two halves of the cbiS genes retained their function in vivo when they were cloned separately. The CbiS enzyme was overproduced in Escherichia coli and was isolated to >95% homogeneity. High-performance liquid chromatography, UV-visible spectroscopy, and mass spectroscopy established alpha-ribazole and cobyric acid as the products of the phosphatase and amidohydrolase reactions, respectively. Reasons why the CbiZ and CobZ enzymes are fused in some archaea are discussed.

Tree disagreement: Measuring and testing incongruence in phylogenies

The branching patterns of phylogenetic trees often disagree even when they have been constructed using different portions of the same data. This phylogenetic discord (incongruence) can be explained by real differences in evolutionary process or history, but also may be due simply to random chance or sampling error. Techniques for measuring and testing the significance of phylogenetic incongruence are used widely in systematic biology, and are necessary when considering genome-scale datasets composed of multiple genes that may or may not have different histories. They are also applicable wherever tree algorithms are used for ordering and interpreting data (e.g., DNA microarrays). Here, I review the different incongruence tests and use them to test the phylogenetic discord of a potentially mobile genetic element (the widespread colonization Island) in the gamma-proteobacteria. I then consider how incongruence tests may be used as a starting point for phylogenetic analysis that accounts for horizontal transfer and duplication events as explanations for homoplasy.

The evolutionary history of Cys-tRNACys formation

The recent discovery of an alternate pathway for indirectly charging tRNA(Cys) has stimulated a re-examination of the evolutionary history of Cys-tRNA(Cys) formation. In the first step of the pathway, O-phosphoseryl-tRNA synthetase charges tRNA(Cys) with O-phosphoserine (Sep), a precursor of the cognate amino acid. In the following step, Sep-tRNA:Cys-tRNA synthase (SepCysS) converts Sep to Cys in a tRNA-dependent reaction. The existence of such a pathway raises several evolutionary questions, including whether the indirect pathway is a recent evolutionary invention, as might be implied from its localization to the Euryarchaea, or, as evidence presented here indicates, whether this pathway is more ancient, perhaps already in existence at the time of the last universal common ancestral state. A comparative phylogenetic approach is used, combining evolutionary information from protein sequences and structures, that takes both the signature of horizontal gene transfer and the recurrence of the full canonical phylogenetic pattern into account, to document the complete evolutionary history of cysteine coding and understand the nature of this process in the last universal common ancestral state. Resulting from the historical study of tRNA(Cys) aminoacylation and the integrative perspective of sequence, structure, and function are 3D models of O-phosphoseryl-tRNA synthetase and SepCysS, which provide experimentally testable predictions regarding the identity and function of key active-site residues in these proteins. The model of SepCysS is used to suggest a sulfhydrylation reaction mechanism, which is predicted to occur at the interface of a SepCysS dimer.

Evolutionary change - patterns and processes

The present review considered: (a) the factors that conditioned the early transition from non-life to life; (b) genome structure and complexity in prokaryotes, eukaryotes, and organelles; (c) comparative human chromosome genomics; and (d) the Brazilian contribution to some of these studies. Understanding the dialectical conflict between freedom and organization is fundamental to give meaning to the patterns and processes of organic evolution.

Constituents of SH1, a novel lipid-containing virus infecting the halophilic euryarchaeon Haloarcula hispanica

Recent studies have indicated that a number of bacterial and eukaryotic viruses that share a common architectural principle are related, leading to the proposal of an early common ancestor. A prediction of this model would be the discovery of similar viruses that infect archaeal hosts. Our main interest lies in icosahedral double-stranded DNA (dsDNA) viruses with an internal membrane, and we now extend our studies to include viruses infecting archaeal hosts. While the number of sequenced archaeal viruses is increasing, very little sequence similarity has been detected between bacterial and eukaryotic viruses. In this investigation we rigorously show that SH1, an icosahedral dsDNA virus infecting Haloarcula hispanica, possesses lipid structural components that are selectively acquired from the host pool. We also determined the sequence of the 31-kb SH1 genome and positively identified genes for 11 structural proteins, with putative identification of three additional proteins. The SH1 genome is unique and, except for a few open reading frames, shows no detectable similarity to other published sequences, but the overall structure of the SH1 virion and its linear genome with inverted terminal repeats is reminiscent of lipid-containing dsDNA bacteriophages like PRD1.

Archaeal genetics - the third way

For decades, archaea were misclassified as bacteria because of their prokaryotic morphology. Molecular phylogeny eventually revealed that archaea, like bacteria and eukaryotes, are a fundamentally distinct domain of life. Genome analyses have confirmed that archaea share many features with eukaryotes, particularly in information processing, and therefore can serve as streamlined models for understanding eukaryotic biology. Biochemists and structural biologists have embraced the study of archaea but geneticists have been more wary, despite the fact that genetic techniques for archaea are quite sophisticated. It is time for geneticists to start asking fundamental questions about our distant relatives.

A Methanocaldococcus jannaschii archaeal signature gene encodes for a 5-formaminoimidazole-4-carboxamide-1-beta-D-ribofuranosyl 5'-monophosphate synthetase. A new enzyme in purine biosynthesis

We have identified and characterized a new member of the ATP-grasp enzyme family that catalyzes the ATP- and formate-dependent formylation of 5-aminoimidazole-4-carboxamide-1-beta-D-ribofuranosyl 5'-monophosphate (AICAR) to 5-formaminoimidazole-4-carboxamide-1-beta-D-ribofuranosyl 5'-monophosphate (FAICAR) in the absence of folates. The enzyme, which we designate as PurP, is the product of the Methanocaldococcus jannaschii purP gene (MJ0136), which is a signature gene for Archaea. As is characteristic of reactions catalyzed by this family of enzymes, the other products of the reaction, ADP and P(i), were produced stoichiometrically with the amount of ATP, formate, and AICAR used. Formyl phosphate was found to substitute for ATP and formate in the reaction, yet the methylene analog, phosphonoacetaldehyde, was not an inhibitor or substrate for the reaction. The enzyme, along with PurO, which catalyzes the cyclization of FAICAR to inosine 5'-monophosphate, catalyzes the same overall transformation in purine biosynthesis as is accomplished by PurH in bacteria and eukaryotes. No homology exists between PurH and either PurO or PurP. 1H NMR and gas chromatography-mass spectrometry analysis of an M. jannaschii cell extract showed the presence of free formate that can be used by the enzyme for purine biosynthesis. This formate arises by the reduction of CO2 with hydrogen; this was demonstrated by incorporating 13C into the formate when M. jannaschii cell extracts were incubated with H13CO3- and hydrogen gas. The presence of this signature gene in all of the Archaea indicates the presence of a purine biosynthetic pathway proceeding in the absence of folate coenzymes.

The phosphoenolpyruvate carboxylase from Methanothermobacter thermautotrophicus has a novel structure

In Methanothermobacter thermautotrophicus, oxaloacetate synthesis is a major and essential CO(2)-fixation reaction. This methanogenic archaeon possesses two oxaloacetate-synthesizing enzymes, pyruvate carboxylase and phosphoenolpyruvate carboxylase. The phosphoenolpyruvate carboxylase from this organism was purified to homogeneity. The subunit size of this homotetrameric protein was 55 kDa, which is about half that of all known bacterial and eukaryotic phosphoenolpyruvate carboxylases (PPCs). The NH(2)-terminal sequence identified this enzyme as the product of MTH943, an open reading frame with no assigned function in the genome sequence. A BLAST search did not show an obvious sequence similarity between MTH943 and known PPCs, which are generally well conserved. This is the first report of a new type of phosphoenolpyruvate carboxylase that we call PpcA ("A" for "archaeal"). Homologs to PpcA were present in most archaeal genomic sequences, but only in three bacterial (Clostridium perfringens, Oenococcus oeni, and Leuconostoc mesenteroides) and no eukaryotic genomes. PpcA was the only recognizable oxaloacetate-producing enzyme in Methanopyrus kandleri, a hydrothermal vent organism. Each PpcA-containing organism lacked a PPC homolog. The activity of M. thermautotrophicus PpcA was not influenced by acetyl coenzyme A and was about 50 times less sensitive to aspartate than the Escherichia coli PPC. The catalytic core (including His(138), Arg(587), and Gly(883)) of the E. coli PPC was partly conserved in PpcA, but three of four aspartate-binding residues (Lys(773), Arg(832), and Asn(881)) were not. PPCs probably evolved from PpcA through a process that added allosteric sites to the enzyme. The reverse is also equally possible.

A new biology for a new century

Biology today is at a crossroads. The molecular paradigm, which so successfully guided the discipline throughout most of the 20th century, is no longer a reliable guide. Its vision of biology now realized, the molecular paradigm has run its course. Biology, therefore, has a choice to make, between the comfortable path of continuing to follow molecular biology's lead or the more invigorating one of seeking a new and inspiring vision of the living world, one that addresses the major problems in biology that 20th century biology, molecular biology, could not handle and, so, avoided. The former course, though highly productive, is certain to turn biology into an engineering discipline. The latter holds the promise of making biology an even more fundamental science, one that, along with physics, probes and defines the nature of reality. This is a choice between a biology that solely does society's bidding and a biology that is society's teacher.

Development of additional selectable markers for the halophilic archaeon Haloferax volcanii based on the leuB and trpA genes

Since most archaea are extremophilic and difficult to cultivate, our current knowledge of their biology is confined largely to comparative genomics and biochemistry. Haloferax volcanii offers great promise as a model organism for archaeal genetics, but until now there has been a lack of a wide variety of selectable markers for this organism. We describe here isolation of H. volcanii leuB and trpA genes encoding 3-isopropylmalate dehydrogenase and tryptophan synthase, respectively, and development of these genes as a positive selection system. DeltaleuB and DeltatrpA mutants were constructed in a variety of genetic backgrounds and were shown to be auxotrophic for leucine and tryptophan, respectively. We constructed both integrative and replicative plasmids carrying the leuB or trpA gene under control of a constitutive promoter. The use of these selectable markers in deletion of the lhr gene of H. volcanii is described.

The genome of Nanoarchaeum equitans: insights into early archaeal evolution and derived parasitism

The hyperthermophile Nanoarchaeum equitans is an obligate symbiont growing in coculture with the crenarchaeon Ignicoccus. Ribosomal protein and rRNA-based phylogenies place its branching point early in the archaeal lineage, representing the new archaeal kingdom Nanoarchaeota. The N. equitans genome (490,885 base pairs) encodes the machinery for information processing and repair, but lacks genes for lipid, cofactor, amino acid, or nucleotide biosyntheses. It is the smallest microbial genome sequenced to date, and also one of the most compact, with 95% of the DNA predicted to encode proteins or stable RNAs. Its limited biosynthetic and catabolic capacity indicates that N. equitans' symbiotic relationship to Ignicoccus is parasitic, making it the only known archaeal parasite. Unlike the small genomes of bacterial parasites that are undergoing reductive evolution, N. equitans has few pseudogenes or extensive regions of noncoding DNA. This organism represents a basal archaeal lineage and has a highly reduced genome.

Experimental determination and system level analysis of essential genes in Escherichia coli MG1655

Defining the gene products that play an essential role in an organism's functional repertoire is vital to understanding the system level organization of living cells. We used a genetic footprinting technique for a genome-wide assessment of genes required for robust aerobic growth of Escherichia coli in rich media. We identified 620 genes as essential and 3,126 genes as dispensable for growth under these conditions. Functional context analysis of these data allows individual functional assignments to be refined. Evolutionary context analysis demonstrates a significant tendency of essential E. coli genes to be preserved throughout the bacterial kingdom. Projection of these data over metabolic subsystems reveals topologic modules with essential and evolutionarily preserved enzymes with reduced capacity for error tolerance.

Horizontal gene transfer: a critical view

It has been suggested that horizontal gene transfer (HGT) is the "essence of phylogeny." In contrast, much data suggest that this is an exaggeration resulting in part from a reliance on inadequate methods to identify HGT events. In addition, the assumption that HGT is a ubiquitous influence throughout evolution is questionable. Instead, rampant global HGT is likely to have been relevant only to primitive genomes. In modern organisms we suggest that both the range and frequencies of HGT are constrained most often by selective barriers. As a consequence those HGT events that do occur most often have little influence on genome phylogeny. Although HGT does occur with important evolutionary consequences, classical Darwinian lineages seem to be the dominant mode of evolution for modern organisms.

How big is the iceberg of which organellar genes in nuclear genomes are but the tip?

As more and more complete bacterial and archaeal genome sequences become available, the role of lateral gene transfer (LGT) in shaping them becomes more and more clear. Over the long term, it may be the dominant force, affecting most genes in most prokaryotes. We review the history of LGT, suggesting reasons why its prevalence and impact were so long dismissed. We discuss various methods purporting to measure the extent of LGT, and evidence for and against the notion that there is a core of never-exchanged genes shared by all genomes, from which we can deduce the "true" organismal tree. We also consider evidence for, and implications of, LGT between prokaryotes and phagocytic eukaryotes.

Prokaryotic evolution in light of gene transfer

Accumulating prokaryotic gene and genome sequences reveal that the exchange of genetic information through both homology-dependent recombination and horizontal (lateral) gene transfer (HGT) is far more important, in quantity and quality, than hitherto imagined. The traditional view, that prokaryotic evolution can be understood primarily in terms of clonal divergence and periodic selection, must be augmented to embrace gene exchange as a creative force, itself responsible for much of the pattern of similarities and differences we see between prokaryotic microbes. Rather than replacing periodic selection on genetic diversity, gene loss, and other chromosomal alterations as important players in adaptive evolution, gene exchange acts in concert with these processes to provide a rich explanatory paradigm-some of whose implications we explore here. In particular, we discuss (1) the role of recombination and HGT in giving phenotypic "coherence" to prokaryotic taxa at all levels of inclusiveness, (2) the implications of these processes for the reconstruction and meaning of "phylogeny," and (3) new views of prokaryotic adaptation and diversification based on gene acquisition and exchange.

On the evolution of cells

A theory for the evolution of cellular organization is presented. The model is based on the (data supported) conjecture that the dynamic of horizontal gene transfer (HGT) is primarily determined by the organization of the recipient cell. Aboriginal cell designs are taken to be simple and loosely organized enough that all cellular componentry can be altered and/or displaced through HGT, making HGT the principal driving force in early cellular evolution. Primitive cells did not carry a stable organismal genealogical trace. Primitive cellular evolution is basically communal. The high level of novelty required to evolve cell designs is a product of communal invention, of the universal HGT field, not intralineage variation. It is the community as a whole, the ecosystem, which evolves. The individual cell designs that evolved in this way are nevertheless fundamentally distinct, because the initial conditions in each case are somewhat different. As a cell design becomes more complex and interconnected a critical point is reached where a more integrated cellular organization emerges, and vertically generated novelty can and does assume greater importance. This critical point is called the "Darwinian Threshold" for the reasons given.

Archaeal phylogeny based on ribosomal proteins

Until recently, phylogenetic analyses of Archaea have mainly been based on ribosomal RNA (rRNA) sequence comparisons, leading to the distinction of the two major archaeal phyla: the Euryarchaeota and the Crenarchaeota. Here, thanks to the recent sequencing of several archaeal genomes, we have constructed a phylogeny based on the fusion of the sequences of the 53 ribosomal proteins present in most of the archaeal species. This phylogeny was remarkably congruent with the rRNA phylogeny, suggesting that both reflected the actual phylogeny of the domain Archaea even if some nodes remained unresolved. In both cases, the branches leading to hyperthermophilic species were short, suggesting that the evolutionary rate of their genes has been slowed down by structural constraints related to environmental adaptation. In addition, to estimate the impact of lateral gene transfer (LGT) on our tree reconstruction, we used a new method that revealed that 8 genes out of the 53 ribosomal proteins used in our study were likely affected by LGT. This strongly suggested that a core of 45 nontransferred ribosomal protein genes existed in Archaea that can be tentatively used to infer the phylogeny of this domain. Interestingly, the tree obtained using only the eight ribosomal proteins likely affected by LGT was not very different from the consensus tree, indicating that LGT mainly brought random phylogenetic noise. The major difference involves organisms living in similar environments, suggesting that LGTs are mainly directed by the physical proximity of the organisms rather than by their phylogenetic proximity.

Evolutionary analysis by whole-genome comparisons

A total of 37 complete genome sequences of bacteria, archaea, and eukaryotes were compared. The percentage of orthologous genes of each species contained within any of the other 36 genomes was established. In addition, the mean identity of the orthologs was calculated. Several conclusions result: (i) a greater absolute number of orthologs of a given species is found in larger species than in smaller ones; (ii) a greater percentage of the orthologous genes of smaller genomes is contained in other species than is the case for larger genomes, which corresponds to a larger proportion of essential genes; (iii) before species can be specifically related to one another in terms of gene content, it is first necessary to correct for the size of the genome; (iv) eukaryotes have a significantly smaller percentage of bacterial orthologs after correction for genome size, which is consistent with their placement in a separate domain; (v) the archaebacteria are specifically related to one another but are not significantly different in gene content from the bacteria as a whole; (vi) determination of the mean identity of all orthologs (involving hundreds of gene comparisons per genome pair) reduces the impact of errors in misidentification of orthologs and to misalignments, and thus it is far more reliable than single gene comparisons; (vii) however, there is a maximum amount of change in protein sequences of 37% mean identity, which limits the use of percentage sequence identity to the lower taxa, a result which should also be true for single gene comparisons of both proteins and rRNA; (viii) most of the species that appear to be specifically related based upon gene content also appear to be specifically related based upon the mean identity of orthologs; (ix) the genes of a majority of species considered in this study have diverged too much to allow the construction of all-encompassing evolutionary trees. However, we have shown that eight species of gram-negative bacteria, six species of gram-positive bacteria, and eight species of archaebacteria are specifically related in terms of gene content, mean identity of orthologs, or both.

Bootstrap, Bayesian probability and maximum likelihood mapping: exploring new tools for comparative genome analyses

BACKGROUND

Horizontal gene transfer (HGT) played an important role in shaping microbial genomes. In addition to genes under sporadic selection, HGT also affects housekeeping genes and those involved in information processing, even ribosomal RNA encoding genes. Here we describe tools that provide an assessment and graphic illustration of the mosaic nature of microbial genomes.

RESULTS

We adapted the Maximum Likelihood (ML) mapping to the analyses of all detected quartets of orthologous genes found in four genomes. We have automated the assembly and analyses of these quartets of orthologs given the selection of four genomes. We compared the ML-mapping approach to more rigorous Bayesian probability and Bootstrap mapping techniques. The latter two approaches appear to be more conservative than the ML-mapping approach, but qualitatively all three approaches give equivalent results. All three tools were tested on mitochondrial genomes, which presumably were inherited as a single linkage group.

CONCLUSIONS

In some instances of interphylum relationships we find nearly equal numbers of quartets strongly supporting the three possible topologies. In contrast, our analyses of genome quartets containing the cyanobacterium Synechocystis sp. indicate that a large part of the cyanobacterial genome is related to that of low GC Gram positives. Other groups that had been suggested as sister groups to the cyanobacteria contain many fewer genes that group with the Synechocystis orthologs. Interdomain comparisons of genome quartets containing the archaeon Halobacterium sp. revealed that Halobacterium sp. shares more genes with Bacteria that live in the same environment than with Bacteria that are more closely related based on rRNA phylogeny. Many of these genes encode proteins involved in substrate transport and metabolism and in information storage and processing. The performed analyses demonstrate that relationships among prokaryotes cannot be accurately depicted by or inferred from the tree-like evolution of a core of rarely transferred genes; rather prokaryotic genomes are mosaics in which different parts have different evolutionary histories. Probability mapping is a valuable tool to explore the mosaic nature of genomes.

The origin of the eukaryotic cell: a genomic investigation

We have collected a set of 347 proteins that are found in eukaryotic cells but have no significant homology to proteins in Archaea and Bacteria. We call these proteins eukaryotic signature proteins (ESPs). The dominant hypothesis for the formation of the eukaryotic cell is that it is a fusion of an archaeon with a bacterium. If this hypothesis is accepted then the three cellular domains, Eukarya, Archaea, and Bacteria, would collapse into two cellular domains. We have used the existence of this set of ESPs to test this hypothesis. The evidence of the ESPs implicates a third cell (chronocyte) in the formation of the eukaryotic cell. The chronocyte had a cytoskeleton that enabled it to engulf prokaryotic cells and a complex internal membrane system where lipids and proteins were synthesized. It also had a complex internal signaling system involving calcium ions, calmodulin, inositol phosphates, ubiquitin, cyclin, and GTP-binding proteins. The nucleus was formed when a number of archaea and bacteria were engulfed by a chronocyte. This formation of the nucleus would restore the three cellular domains as the Chronocyte was not a cell that belonged to the Archaea or to the Bacteria.