Comparative genomics of Archaea: how much have we learned in six years, and what's next?

Microbial communities of halite deposits and other hypersaline environments

Large regions of Earth's surface are underlain by salt deposits that evaporated from ancient oceans and are populated by extreme halophilic microbes. While the microbiology of ancient evaporites has been well studied, the ecology of halite deposits and more recently formed NaCl "salticle" stalactite structures (speleothems) in a Triassic halite mine are less well characterized. The microbiome of Kilroot Salt Mine was profiled using conventional and enhanced culturing techniques. From this, 89 halophilic archaeal isolates from six known genera, and 55 halophilic or halotolerant bacterial isolates from 18 genera were obtained. Culture-independent metagenomic approaches also revealed that culturing techniques were inadvertently biased toward specific taxa, and the need for optimized isolation procedures are required to enhance cultivation diversity. Speleothems formed from saturated brines are unique structures that have the potential to entomb haloarchaea cells for thousands of years within fluid inclusions. The presence of such fluid inclusions, alongside the high abundance of genes related to glycerol metabolism, biofilm formation, and persister cell formation is highly suggestive of an environmental niche that could promote longevity and survivability. Finally, previous studies reporting the discovery of novel biocatalysts from the Kilroot mine microbiome, suggests that this environment may be an untapped source of chemical diversity with high biodiscovery potential.

Archaea: An Agro-Ecological Perspective

Microorganisms inhabiting bulk soil and rhizosphere play an important role in soil biogeochemical cycles leading to enhanced plant growth and productivity. In this context, the role of bacteria is well established, however, not much reports are available about the role archaea plays in this regard. Literature suggests that archaea also play a greater role in nutrient cycling of carbon, nitrogen, sulfur, and other minerals, possess various plant growth promoting attributes, and can impart tolerance to various abiotic stresses (especially osmotic and oxidative) in areas of high salinity, low and high temperatures and hydrogen ion concentrations. Thermoacidophilic archaea have been found to potentially involve in bioleaching of mineral ores and bioremediation of chemical pollutants and aromatic compounds. Looking at immense potential of archaea in promoting plant growth, alleviating abiotic stresses, and remediating contaminated sites, detailed studies are required to establish their role in different ecological processes, and their interactions in rhizosphere with plant and other microflora (bacteria and fungi) in different ecosystems. In this review, a brief discussion on archaea from the agro-ecological point of view is presented.

LUCA to LECA, the Lucacene: A model for the gigayear delay from the first prokaryote to eukaryogenesis

It is puzzling why life on Earth consisted of prokaryotes for up to 2.5 ± 0.5 billion years (Gy) before the appearance of the first eukaryotes. This period, from LUCA (Last Universal Common Ancestor) to LECA (Last Eucaryotic Common Ancestor), we have named the Lucacene, to suggest all prokaryotic descendants of LUCA before the appearance of LECA. Here we present a simple model based on horizontal gene transfer (HGT). It is the process of HGT from Bacteria to Archaea and its reverse that we wish to simulate and estimate its duration until eukaryogenesis. Rough quantitation of its parameters shows that the model may explain the long duration of the Lucacene.

Comparative Genomics Screen Identifies Microbe-Associated Molecular Patterns from 'Candidatus Liberibacter' spp. That Elicit Immune Responses in Plants

Citrus huanglongbing (HLB), caused by phloem-limited 'Candidatus Liberibacter' bacteria, is a destructive disease threatening the worldwide citrus industry. The mechanisms of pathogenesis are poorly understood and no efficient strategy is available to control HLB. Here, we used a comparative genomics screen to identify candidate microbe-associated molecular patterns (MAMPs) from 'Ca. Liberibacter' spp. We identified the core genome from multiple 'Ca. Liberibacter' pathogens, and searched for core genes with signatures of positive selection. We hypothesized that genes encoding putative MAMPs would evolve to reduce recognition by the plant immune system, while retaining their essential functions. To efficiently screen candidate MAMP peptides, we established a high-throughput microtiter plate-based screening assay, particularly for citrus, that measured reactive oxygen species (ROS) production, which is a common immune response in plants. We found that two peptides could elicit ROS production in Arabidopsis and Nicotiana benthamiana. One of these peptides elicited ROS production and defense gene expression in HLB-tolerant citrus genotypes, and induced MAMP-triggered immunity against the bacterial pathogen Pseudomonas syringae. Our findings identify MAMPs that boost immunity in citrus and could help prevent or reduce HLB infection.

Phyletic Distribution and Lineage-Specific Domain Architectures of Archaeal Two-Component Signal Transduction Systems

The two-component signal transduction (TCS) machinery is a key mechanism of sensing environmental changes in the prokaryotic world. TCS systems have been characterized thoroughly in bacteria but to a much lesser extent in archaea. Here, we provide an updated census of more than 2,000 histidine kinases and response regulators encoded in 218 complete archaeal genomes, as well as unfinished genomes available from metagenomic data. We describe the domain architectures of the archaeal TCS components, including several novel output domains, and discuss the evolution of the archaeal TCS machinery. The distribution of TCS systems in archaea is strongly biased, with high levels of abundance in haloarchaea and thaumarchaea but none detected in the sequenced genomes from the phyla Crenarchaeota, Nanoarchaeota, and Korarchaeota The archaeal sensor histidine kinases are generally similar to their well-studied bacterial counterparts but are often located in the cytoplasm and carry multiple PAS and/or GAF domains. In contrast, archaeal response regulators differ dramatically from the bacterial ones. Most archaeal genomes do not encode any of the major classes of bacterial response regulators, such as the DNA-binding transcriptional regulators of the OmpR/PhoB, NarL/FixJ, NtrC, AgrA/LytR, and ActR/PrrA families and the response regulators with GGDEF and/or EAL output domains. Instead, archaea encode multiple copies of response regulators containing either the stand-alone receiver (REC) domain or combinations of REC with PAS and/or GAF domains. Therefore, the prevailing mechanism of archaeal TCS signaling appears to be via a variety of protein-protein interactions, rather than direct transcriptional regulation.IMPORTANCE Although the Archaea represent a separate domain of life, their signaling systems have been assumed to be closely similar to the bacterial ones. A study of the domain architectures of the archaeal two-component signal transduction (TCS) machinery revealed an overall similarity of archaeal and bacterial sensory modules but substantial differences in the signal output modules. The prevailing mechanism of archaeal TCS signaling appears to involve various protein-protein interactions rather than direct transcription regulation. The complete list of histidine kinases and response regulators encoded in the analyzed archaeal genomes is available online at http://www.ncbi.nlm.nih.gov/Complete_Genomes/TCSarchaea.html.

New insights into marine group III Euryarchaeota, from dark to light

Marine Euryarchaeota remain among the least understood major components of marine microbial communities. Marine group II Euryarchaeota (MG-II) are more abundant in surface waters (4-20% of the total prokaryotic community), whereas marine group III Euryarchaeota (MG-III) are generally considered low-abundance members of deep mesopelagic and bathypelagic communities. Using genome assembly from direct metagenome reads and metagenomic fosmid clones, we have identified six novel MG-III genome sequence bins from the photic zone (Epi1-6) and two novel bins from deep-sea samples (Bathy1-2). Genome completeness in those genome bins varies from 44% to 85%. Photic-zone MG-III bins corresponded to novel groups with no similarity, and significantly lower GC content, when compared with previously described deep-MG-III genome bins. As found in many other epipelagic microorganisms, photic-zone MG-III bins contained numerous photolyase and rhodopsin genes, as well as genes for peptide and lipid uptake and degradation, suggesting a photoheterotrophic lifestyle. Phylogenetic analysis of these photolyases and rhodopsins as well as their genomic context suggests that these genes are of bacterial origin, supporting the hypothesis of an MG-III ancestor that lived in the dark ocean. Epipelagic MG-III occur sporadically and in relatively small proportions in marine plankton, representing only up to 0.6% of the total microbial community reads in metagenomes. None of the reconstructed epipelagic MG-III genomes were present in metagenomes from aphotic zone depths or from high latitude regions. Most low-GC bins were highly enriched at the deep chlorophyll maximum zones, with the exception of Epi1, which appeared evenly distributed throughout the photic zone worldwide.

Comparative genomic analysis of evolutionarily conserved but functionally uncharacterized membrane proteins in archaea: Prediction of novel components of secretion, membrane remodeling and glycosylation systems

A systematic comparative genomic analysis of all archaeal membrane proteins that have been projected to the last archaeal common ancestor gene set led to the identification of several novel components of predicted secretion, membrane remodeling, and protein glycosylation systems. Among other findings, most crenarchaea have been shown to encode highly diverged orthologs of the membrane insertase YidC, which is nearly universal in bacteria, eukaryotes, and euryarchaea. We also identified a vast family of archaeal proteins, including the C-terminal domain of N-glycosylation protein AglD, as membrane flippases homologous to the flippase domain of bacterial multipeptide resistance factor MprF, a bifunctional lysylphosphatidylglycerol synthase and flippase. Additionally, several proteins were predicted to function as membrane transporters. The results of this work, combined with our previous analyses, reveal an unexpected diversity of putative archaeal membrane-associated functional systems that remain to be functionally characterized. A more general conclusion from this work is that the currently available collection of archaeal (and bacterial) genomes could be sufficient to identify (almost) all widespread functional modules and develop experimentally testable predictions of their functions.

A new class of marine Euryarchaeota group II from the Mediterranean deep chlorophyll maximum

We have analyzed metagenomic fosmid clones from the deep chlorophyll maximum (DCM), which, by genomic parameters, correspond to the 16S ribosomal RNA (rRNA)-defined marine Euryarchaeota group IIB (MGIIB). The fosmid collections associated with this group add up to 4 Mb and correspond to at least two species within this group. From the proposed essential genes contained in the collections, we infer that large sections of the conserved regions of the genomes of these microbes have been recovered. The genomes indicate a photoheterotrophic lifestyle, similar to that of the available genome of MGIIA (assembled from an estuarine metagenome in Puget Sound, Washington Pacific coast), with a proton-pumping rhodopsin of the same kind. Several genomic features support an aerobic metabolism with diversified substrate degradation capabilities that include xenobiotics and agar. On the other hand, these MGIIB representatives are non-motile and possess similar genome size to the MGIIA-assembled genome, but with a lower GC content. The large phylogenomic gap with other known archaea indicates that this is a new class of marine Euryarchaeota for which we suggest the name Thalassoarchaea. The analysis of recruitment from available metagenomes indicates that the representatives of group IIB described here are largely found at the DCM (ca. 50 m deep), in which they are abundant (up to 0.5% of the reads), and at the surface mostly during the winter mixing, which explains formerly described 16S rRNA distribution patterns. Their uneven representation in environmental samples that are close in space and time might indicate sporadic blooms.

Comparative genomics of the marine bacterial genus Glaciecola reveals the high degree of genomic diversity and genomic characteristic for cold adaptation

To what extent the genomes of different species belonging to one genus can be diverse and the relationship between genomic differentiation and environmental factor remain unclear for oceanic bacteria. With many new bacterial genera and species being isolated from marine environments, this question warrants attention. In this study, we sequenced all the type strains of the published species of Glaciecola, a recently defined cold-adapted genus with species from diverse marine locations, to study the genomic diversity and cold-adaptation strategy in this genus.The genome size diverged widely from 3.08 to 5.96 Mb, which can be explained by massive gene gain and loss events. Horizontal gene transfer and new gene emergence contributed substantially to the genome size expansion. The genus Glaciecola had an open pan-genome. Comparative genomic research indicated that species of the genus Glaciecola had high diversity in genome size, gene content and genetic relatedness. This may be prevalent in marine bacterial genera considering the dynamic and complex environments of the ocean. Species of Glaciecola had some common genomic features related to cold adaptation, which enable them to thrive and play a role in biogeochemical cycle in the cold marine environments.

Carbohydrate metabolism in Archaea: current insights into unusual enzymes and pathways and their regulation

The metabolism of Archaea, the third domain of life, resembles in its complexity those of Bacteria and lower Eukarya. However, this metabolic complexity in Archaea is accompanied by the absence of many "classical" pathways, particularly in central carbohydrate metabolism. Instead, Archaea are characterized by the presence of unique, modified variants of classical pathways such as the Embden-Meyerhof-Parnas (EMP) pathway and the Entner-Doudoroff (ED) pathway. The pentose phosphate pathway is only partly present (if at all), and pentose degradation also significantly differs from that known for bacterial model organisms. These modifications are accompanied by the invention of "new," unusual enzymes which cause fundamental consequences for the underlying regulatory principles, and classical allosteric regulation sites well established in Bacteria and Eukarya are lost. The aim of this review is to present the current understanding of central carbohydrate metabolic pathways and their regulation in Archaea. In order to give an overview of their complexity, pathway modifications are discussed with respect to unusual archaeal biocatalysts, their structural and mechanistic characteristics, and their regulatory properties in comparison to their classic counterparts from Bacteria and Eukarya. Furthermore, an overview focusing on hexose metabolic, i.e., glycolytic as well as gluconeogenic, pathways identified in archaeal model organisms is given. Their energy gain is discussed, and new insights into different levels of regulation that have been observed so far, including the transcript and protein levels (e.g., gene regulation, known transcription regulators, and posttranslational modification via reversible protein phosphorylation), are presented.

Updated clusters of orthologous genes for Archaea: a complex ancestor of the Archaea and the byways of horizontal gene transfer

Background

Collections of Clusters of Orthologous Genes (COGs) provide indispensable tools for comparative genomic analysis, evolutionary reconstruction and functional annotation of new genomes. Initially, COGs were made for all complete genomes of cellular life forms that were available at the time. However, with the accumulation of thousands of complete genomes, construction of a comprehensive COG set has become extremely computationally demanding and prone to error propagation, necessitating the switch to taxon-specific COG collections. Previously, we reported the collection of COGs for 41 genomes of Archaea (arCOGs). Here we present a major update of the arCOGs and describe evolutionary reconstructions to reveal general trends in the evolution of Archaea.

Results

The updated version of the arCOG database incorporates 91% of the pangenome of 120 archaea (251,032 protein-coding genes altogether) into 10,335 arCOGs. Using this new set of arCOGs, we performed maximum likelihood reconstruction of the genome content of archaeal ancestral forms and gene gain and loss events in archaeal evolution. This reconstruction shows that the last Common Ancestor of the extant Archaea was an organism of greater complexity than most of the extant archaea, probably with over 2,500 protein-coding genes. The subsequent evolution of almost all archaeal lineages was apparently dominated by gene loss resulting in genome streamlining. Overall, in the evolution of Archaea as well as a representative set of bacteria that was similarly analyzed for comparison, gene losses are estimated to outnumber gene gains at least 4 to 1. Analysis of specific patterns of gene gain in Archaea shows that, although some groups, in particular Halobacteria, acquire substantially more genes than others, on the whole, gene exchange between major groups of Archaea appears to be largely random, with no major ‘highways’ of horizontal gene transfer.

Conclusions

The updated collection of arCOGs is expected to become a key resource for comparative genomics, evolutionary reconstruction and functional annotation of new archaeal genomes. Given that, in spite of the major increase in the number of genomes, the conserved core of archaeal genes appears to be stabilizing, the major evolutionary trends revealed here have a chance to stand the test of time.

Reviewers

This article was reviewed by (for complete reviews see the Reviewers’ Reports section): Dr. PLG, Prof. PF, Dr. PL (nominated by Prof. JPG).

A pursuit of lineage-specific and niche-specific proteome features in the world of archaea

Background

Archaea evoke interest among researchers for two enigmatic characteristics –a combination of bacterial and eukaryotic components in their molecular architectures and an enormous diversity in their life-style and metabolic capabilities. Despite considerable research efforts, lineage- specific/niche-specific molecular features of the whole archaeal world are yet to be fully unveiled. The study offers the first large-scale in silico proteome analysis of all archaeal species of known genome sequences with a special emphasis on methanogenic and sulphur-metabolising archaea.

Results

Overall amino acid usage in archaea is dominated by GC-bias. But the environmental factors like oxygen requirement or thermal adaptation seem to play important roles in selection of residues with no GC-bias at the codon level. All methanogens, irrespective of their thermal/salt adaptation, show higher usage of Cys and have relatively acidic proteomes, while the proteomes of sulphur-metabolisers have higher aromaticity and more positive charges. Despite of exhibiting thermophilic life-style, korarchaeota possesses an acidic proteome. Among the distinct trends prevailing in COGs (Cluster of Orthologous Groups of proteins) distribution profiles, crenarchaeal organisms display higher intra-order variations in COGs repertoire, especially in the metabolic ones, as compared to euryarchaea. All methanogens are characterised by a presence of 22 exclusive COGs.

Conclusions

Divergences in amino acid usage, aromaticity/charge profiles and COG repertoire among methanogens and sulphur-metabolisers, aerobic and anaerobic archaea or korarchaeota and nanoarchaeota, as elucidated in the present study, point towards the presence of distinct molecular strategies for niche specialization in the archaeal world.

Nmag_2608, an extracellular ubiquitin-like domain-containing protein from the haloalkaliphilic archaeon Natrialba magadii

Ubiquitin-like proteins (Ubls) and ubiquitin-like domain-containing proteins (Ulds) found in both eukaryotes and prokaryotes display an ubiquitin fold. We previously characterized a 124-amino acid polypeptide (P400) from the haloalkaliphilic archaeon Natrialba magadii having structural homology with ubiquitin family proteins. The reported N. magadii's genome allowed the identification of the Nmag_2608 gene for the protein containing P400, which belongs to specific orthologs of halophilic organisms. It was found that Nmag_2608 has an N-terminal signal peptide with a lipobox motif characteristic of bacterial lipoproteins. Also, it presents partial identity with the ubiquitin-like domain-containing proteins, soluble ligand binding β-grasp proteins. Western blots and heterologous expression tests in E. coli evidenced that Nmag_2608 is processed and secreted outside the cell, where it could perform its function. The analysis of Nmag_2608 expression in N. magadii's cells suggests a co-transcription with the adjoining Nmag_2609 gene encoding a protein of the cyclase family. Also, the transcript level decreased in cells grown in low salinity and starved. To conclude, this work reports for the first time an extracellular archaeal protein with an ubiquitin-like domain.

The archaeal 'TACK' superphylum and the origin of eukaryotes

Although most hypotheses to explain the emergence of the eukaryotic lineage are conflicting, some consensus exists concerning the requirement of a genomic fusion between archaeal and bacterial components. Recent phylogenomic studies have provided support for eocyte-like scenarios in which the alleged 'archaeal parent' of the eukaryotic cell emerged from the Crenarchaeota/Thaumarchaeota. Here, we provide evidence for a scenario in which this archaeal parent emerged from within the 'TACK' superphylum that comprises the Thaumarchaeota, Crenarchaeota and Korarchaeota, as well as the recently proposed phylum 'Aigarchaeota'. In support of this view, functional and comparative genomics studies have unearthed an increasing number of features that are uniquely shared by the TACK superphylum and eukaryotes, including proteins involved in cytokinesis, membrane remodeling, cell shape determination and protein recycling.

Crenarchaeal CdvA forms double-helical filaments containing DNA and interacts with ESCRT-III-like CdvB

BACKGROUND

The phylum Crenarchaeota lacks the FtsZ cell division hallmark of bacteria and employs instead Cdv proteins. While CdvB and CdvC are homologues of the eukaryotic ESCRT-III and Vps4 proteins, implicated in membrane fission processes during multivesicular body biogenesis, cytokinesis and budding of some enveloped viruses, little is known about the structure and function of CdvA. Here, we report the biochemical and biophysical characterization of the three Cdv proteins from the hyperthermophilic archaeon Metallospherae sedula.

METHODOLOGY/PRINCIPAL FINDINGS

Using sucrose density gradient ultracentrifugation and negative staining electron microscopy, we evidenced for the first time that CdvA forms polymers in association with DNA, similar to known bacterial DNA partitioning proteins. We also observed that, in contrast to full-lengh CdvB that was purified as a monodisperse protein, the C-terminally deleted CdvB construct forms filamentous polymers, a phenomenon previously observed with eukaryotic ESCRT-III proteins. Based on size exclusion chromatography data combined with detection by multi-angle laser light scattering analysis, we demonstrated that CdvC assembles, in a nucleotide-independent way, as homopolymers resembling dodecamers and endowed with ATPase activity in vitro. The interactions between these putative cell division partners were further explored. Thus, besides confirming the previous observations that CdvB interacts with both CdvA and CdvC, our data demonstrate that CdvA/CdvB and CdvC/CdvB interactions are not mutually exclusive.

CONCLUSIONS/SIGNIFICANCE

Our data reinforce the concept that Cdv proteins are closely related to the eukaryotic ESCRT-III counterparts and suggest that the organization of the ESCRT-III machinery at the Crenarchaeal cell division septum is organized by CdvA an ancient cytoskeleton protein that might help to coordinate genome segregation.

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.

Hot Transcriptomics

DNA microarray technology allows for a quick and easy comparison of complete transcriptomes, resulting in improved molecular insight in fluctuations of gene expression. After emergence of the microarray technology about a decade ago, the technique has now matured and has become routine in many molecular biology laboratories. Numerous studies have been performed that have provided global transcription patterns of many organisms under a wide range of conditions. Initially, implementation of this high-throughput technology has lead to high expectations for ground breaking discoveries. Here an evaluation is performed of the insight that transcriptome analysis has brought about in the field of hyperthermophilic archaea. The examples that will be discussed have been selected on the basis of their impact, in terms of either biological insight or technological progress.

Modeling of Escherichia coli Endonuclease V structure in complex with DNA

Endonuclease V (EndoV) is a metal-dependent DNA repair enzyme involved in removal of deaminated bases (e.g., deoxyuridine, deoxyinosine, and deoxyxanthosine), with pairing specificities different from the original bases. Homologs of EndoV are present in all major phyla from bacteria to humans and their function is quite well analyzed. EndoV has been combined with DNA ligase to develop an enzymatic method for mutation scanning and has been engineered to obtain variants with different substrate specificities that serve as improved tools in mutation recognition and cancer mutation scanning. However, little is known about the structure and mechanism of substrate DNA binding by EndoV. Here, we present the results of a bioinformatic analysis and a structural model of EndoV from Escherichia coli in complex with DNA. The structure was obtained by a combination of fold-recognition, comparative modeling, de novo modeling and docking methods. The modeled structure provides a convenient tool to study protein sequence-structure-function relationships in EndoV and to engineer its further variants.

Genomic context analysis in Archaea suggests previously unrecognized links between DNA replication and translation

Specific functional interactions of proteins involved in DNA replication and/or DNA repair or transcription might occur in Archaea, suggesting a previously unrecognized regulatory network coupling DNA replication and translation, which might also exist in Eukarya.

Background

Comparative analysis of genomes is valuable to explore evolution of genomes, deduce gene functions, or predict functional linking between proteins. Here, we have systematically analyzed the genomic environment of all known DNA replication genes in 27 archaeal genomes to infer new connections for DNA replication proteins from conserved genomic associations.

Results

Two distinct sets of DNA replication genes frequently co-localize in archaeal genomes: the first includes the genes for PCNA, the small subunit of the DNA primase (PriS), and Gins15; the second comprises the genes for MCM and Gins23. Other genomic associations of genes encoding proteins involved in informational processes that may be functionally relevant at the cellular level have also been noted; in particular, the association between the genes for PCNA, transcription factor S, and NudF. Surprisingly, a conserved cluster of genes coding for proteins involved in translation or ribosome biogenesis (S27E, L44E, aIF-2 alpha, Nop10) is almost systematically contiguous to the group of genes coding for PCNA, PriS, and Gins15. The functional relevance of this cluster encoding proteins conserved in Archaea and Eukarya is strongly supported by statistical analysis. Interestingly, the gene encoding the S27E protein, also known as metallopanstimulin 1 (MPS-1) in human, is overexpressed in multiple cancer cell lines.

Conclusion

Our genome context analysis suggests specific functional interactions for proteins involved in DNA replication between each other or with proteins involved in DNA repair or transcription. Furthermore, it suggests a previously unrecognized regulatory network coupling DNA replication and translation in Archaea that may also exist in Eukarya.

Mesophilic Crenarchaeota: proposal for a third archaeal phylum, the Thaumarchaeota

The archaeal domain is currently divided into two major phyla, the Euryarchaeota and Crenarchaeota. During the past few years, diverse groups of uncultivated mesophilic archaea have been discovered and affiliated with the Crenarchaeota. It was recently recognized that these archaea have a major role in geochemical cycles. Based on the first genome sequence of a crenarchaeote, Cenarchaeum symbiosum, we show that these mesophilic archaea are different from hyperthermophilic Crenarchaeota and branch deeper than was previously assumed. Our results indicate that C. symbiosum and its relatives are not Crenarchaeota, but should be considered as a third archaeal phylum, which we propose to name Thaumarchaeota (from the Greek 'thaumas', meaning wonder).

The Very Early Stages of Biological Evolution and the Nature of the Last Common Ancestor of the Three Major Cell Domains

Quantitative estimates of the gene complement of the last common ancestor of all extant organisms, that is, the cenancestor, may be hindered by ancient horizontal gene transfer events and polyphyletic gene losses, as well as by biases in genome databases and methodological artifacts. Nevertheless, most reports agree that the last common ancestor resembled extant prokaryotes. A significant number of the highly conserved genes are sequences involved in the synthesis, degradation, and binding of RNA, including transcription and translation. Although the gene complement of the cenancestor includes sequences that may have originated in different epochs, the extraordinary conservation of RNA-related sequences supports the hypothesis that the last common ancestor was an evolutionary outcome of the so-called RNA/protein world. The available evidence suggests that the cenancestor was not a hyperthermophile, but it is currently not possible to assess its ecological niche or its mode of energy acquisition and carbon sources.

Clusters of orthologous genes for 41 archaeal genomes and implications for evolutionary genomics of archaea

Background

An evolutionary classification of genes from sequenced genomes that distinguishes between orthologs and paralogs is indispensable for genome annotation and evolutionary reconstruction. Shortly after multiple genome sequences of bacteria, archaea, and unicellular eukaryotes became available, an attempt on such a classification was implemented in Clusters of Orthologous Groups of proteins (COGs). Rapid accumulation of genome sequences creates opportunities for refining COGs but also represents a challenge because of error amplification. One of the practical strategies involves construction of refined COGs for phylogenetically compact subsets of genomes.

Results

New Archaeal Clusters of Orthologous Genes (arCOGs) were constructed for 41 archaeal genomes (13 Crenarchaeota, 27 Euryarchaeota and one Nanoarchaeon) using an improved procedure that employs a similarity tree between smaller, group-specific clusters, semi-automatically partitions orthology domains in multidomain proteins, and uses profile searches for identification of remote orthologs. The annotation of arCOGs is a consensus between three assignments based on the COGs, the CDD database, and the annotations of homologs in the NR database. The 7538 arCOGs, on average, cover ~88% of the genes in a genome compared to a ~76% coverage in COGs. The finer granularity of ortholog identification in the arCOGs is apparent from the fact that 4538 arCOGs correspond to 2362 COGs; ~40% of the arCOGs are new. The archaeal gene core (protein-coding genes found in all 41 genome) consists of 166 arCOGs. The arCOGs were used to reconstruct gene loss and gene gain events during archaeal evolution and gene sets of ancestral forms. The Last Archaeal Common Ancestor (LACA) is conservatively estimated to possess 996 genes compared to 1245 and 1335 genes for the last common ancestors of Crenarchaeota and Euryarchaeota, respectively. It is inferred that LACA was a chemoautotrophic hyperthermophile that, in addition to the core archaeal functions, encoded more idiosyncratic systems, e.g., the CASS systems of antivirus defense and some toxin-antitoxin systems.

Conclusion

The arCOGs provide a convenient, flexible framework for functional annotation of archaeal genomes, comparative genomics and evolutionary reconstructions. Genomic reconstructions suggest that the last common ancestor of archaea might have been (nearly) as advanced as the modern archaeal hyperthermophiles. ArCOGs and related information are available at: .

Reviewers

This article was reviewed by Peer Bork, Patrick Forterre, and Purificacion Lopez-Garcia.

Protein disulfide oxidoreductases and the evolution of thermophily: was the last common ancestor a heat-loving microbe?

Protein disulfide oxidoreductases (PDOs) are redox enzymes that catalyze dithiol-disulfide exchange reactions. Their sequences and structure reveal the presence of two thioredoxin fold units, each of which is endowed with a catalytic site CXXC motif. PDOs are the outcome of an ancient gene duplication event. They have been described in a number of thermophilic and hyperthermophilic species, where they play a critical role in the structural stabilization of intracellular proteins. PDOs are homologous to both the N-terminal domain of the bacterial alkyl hydroperoxide reductase (AhpF) and to the eukaryotic protein disulfide isomerase (PDI). Phylogenetic analysis of PDOs suggests that they first evolved in the crenarchaeota, spreading from them into the Bacteria via the euryarchaeota. These results imply that the last common ancestor (LCA) of all extant living beings lacked a PDO and argue, albeit weakly, against a thermophilic LCA.

Comparative analysis of predicted gene expression among deep-sea genomes

Deep-sea species live in an environment that is specifically characterized by extreme temperature and hydrostatic pressure. In this work, predicted highly expressed (PHX) genes are comparatively analyzed for six deep-sea microbes, which allows us to pinpoint the common highly expressed genes shared by them. The relationships between gene expression level and some basic properties such as genomic G + C content, optimal growth temperature (OGT), and environmental hydrostatic pressure of the six deep-sea species are also investigated. We find that the percentage of PHX genes out of a whole genome positively correlates to OGT for the deep-sea genomes, whereas such positive correlation seems not to exist between environmental hydrostatic pressure and percentage of PHX genes. Moreover, there exists a negative correlation between genomic G + C content and diversity of gene expression level for the deep-sea genomes, which is in sharp contrast to land-living microbes. We report the top 20 PHX genes for the six deep-sea genomes and find no common highly expressed genes shared by them except for ribosomal proteins, transcription factors, and translation factors. Our present work proffers a paradigm for studying the relationship between environmental factors and microbes' predicted gene expression level.

Crystal structure of the human GINS complex

The GINS complex mediates the assembly of the MCM2-7 (minichromosome maintenance) complex with proteins in a replisome progression complex. The eukaryotic GINS complex is composed of Sld5, Psf1, Psf2, and Psf3, which must be assembled for cell proliferation. We determined the crystal structure of the human GINS complex: GINS forms an elliptical shape with a small central channel. The structures of Sld5 and Psf2 resemble those of Psf1 and Psf3, respectively. In addition, the N-terminal and C-terminal domains of Sld5/Psf1 are permuted in Psf2/Psf3, which suggests that the four proteins have evolved from a common ancestor. Using a structure-based mutational analysis, we identified the functionally critical surface regions of the GINS complex.

Genome and proteome of long-chain alkane degrading Geobacillus thermodenitrificans NG80-2 isolated from a deep-subsurface oil reservoir

The complete genome sequence of Geobacillus thermodenitrificans NG80-2, a thermophilic bacillus isolated from a deep oil reservoir in Northern China, consists of a 3,550,319-bp chromosome and a 57,693-bp plasmid. The genome reveals that NG80-2 is well equipped for adaptation into a wide variety of environmental niches, including oil reservoirs, by possessing genes for utilization of a broad range of energy sources, genes encoding various transporters for efficient nutrient uptake and detoxification, and genes for a flexible respiration system including an aerobic branch comprising five terminal oxidases and an anaerobic branch comprising a complete denitrification pathway for quick response to dissolved oxygen fluctuation. The identification of a nitrous oxide reductase gene has not been previously described in Gram-positive bacteria. The proteome further reveals the presence of a long-chain alkane degradation pathway; and the function of the key enzyme in the pathway, the long-chain alkane monooxygenase LadA, is confirmed by in vivo and in vitro experiments. The thermophilic soluble monomeric LadA is an ideal candidate for treatment of environmental oil pollutions and biosynthesis of complex molecules.

Identification of genes encoding tRNA modification enzymes by comparative genomics

As the molecular adapters between codons and amino acids, transfer-RNAs are pivotal molecules of the genetic code. The coding properties of a tRNA molecule do not reside only in its primary sequence. Posttranscriptional nucleoside modifications, particularly in the anticodon loop, can modify cognate codon recognition, affect aminoacylation properties, or stabilize the codon-anticodon wobble base pairing to prevent ribosomal frameshifting. Despite a wealth of biophysical and structural knowledge of the tRNA modifications themselves, their pathways of biosynthesis had been until recently only partially characterized. This discrepancy was mainly due to the lack of obvious phenotypes for tRNA modification-deficient strains and to the difficulty of the biochemical assays used to detect tRNA modifications. However, the availability of hundreds of whole-genome sequences has allowed the identification of many of these missing tRNA-modification genes. This chapter reviews the methods that were used to identify these genes with a special emphasis on the comparative genomic approaches. Methods that link gene and function but do not rely on sequence homology will be detailed, with examples taken from the tRNA modification field.

High precision multi-genome scale reannotation of enzyme function by EFICAz

BACKGROUND

The functional annotation of most genes in newly sequenced genomes is inferred from similarity to previously characterized sequences, an annotation strategy that often leads to erroneous assignments. We have performed a reannotation of 245 genomes using an updated version of EFICAz, a highly precise method for enzyme function prediction.

RESULTS

Based on our three-field EC number predictions, we have obtained lower-bound estimates for the average enzyme content in Archaea (29%), Bacteria (30%) and Eukarya (18%). Most annotations added in KEGG from 2005 to 2006 agree with EFICAz predictions made in 2005. The coverage of EFICAz predictions is significantly higher than that of KEGG, especially for eukaryotes. Thousands of our novel predictions correspond to hypothetical proteins. We have identified a subset of 64 hypothetical proteins with low sequence identity to EFICAz training enzymes, whose biochemical functions have been recently characterized and find that in 96% (84%) of the cases we correctly identified their three-field (four-field) EC numbers. For two of the 64 hypothetical proteins: PA1167 from Pseudomonas aeruginosa, an alginate lyase (EC 4.2.2.3) and Rv1700 of Mycobacterium tuberculosis H37Rv, an ADP-ribose diphosphatase (EC 3.6.1.13), we have detected annotation lag of more than two years in databases. Two examples are presented where EFICAz predictions act as hypothesis generators for understanding the functional roles of hypothetical proteins: FLJ11151, a human protein overexpressed in cancer that EFICAz identifies as an endopolyphosphatase (EC 3.6.1.10), and MW0119, a protein of Staphylococcus aureus strain MW2 that we propose as candidate virulence factor based on its EFICAz predicted activity, sphingomyelin phosphodiesterase (EC 3.1.4.12).

CONCLUSION

Our results suggest that we have generated enzyme function annotations of high precision and recall. These predictions can be mined and correlated with other information sources to generate biologically significant hypotheses and can be useful for comparative genome analysis and automated metabolic pathway reconstruction.

Origin of the bacterial SET domain genes: vertical or horizontal?

The presence of Supressor of variegation-Enhanser of zeste-Trithorax (SET) domain genes in bacteria is a current paradigm for lateral genetic exchange between eukaryotes and prokaryotes. Because a major function of SET domain proteins is the chemical modification of chromatin and bacteria do not have chromatin, there is no apparent functional requirement for the existence of bacterial SET domain genes. Consequently, their finding in only a small fraction of pathogenic and symbiotic bacteria was taken as evidence that bacteria have obtained the SET domain genes from their hosts. Furthermore, it was proposed that the products of the genes would, most likely, be involved in bacteria-host interactions. The broadened scope of sequenced bacterial genomes to include also free-living and environmental species provided a larger sample to analyze the bacterial SET domain genes. By phylogenetic analysis, examination of individual chromosomal regions for signs of insertion, and evaluating the chromosomal versus SET domain genes' GC contents, we provide evidence that SET domain genes have existed in the bacterial domain of life independently of eukaryotes. The bacterial genes have undergone an evolution of their own unconnected to the evolution of the eukaryotic SET domain genes. Initial finding of SET domain genes in predominantly pathogenic and symbiotic bacteria resulted, most probably, from a biased sample. However, a lateral transfer of SET domain genes may have occurred between some bacteria and a family of Archaea. A model for the evolution and distribution of SET domain genes in bacteria is proposed.

Characterization of an ATP-dependent DNA ligase from the acidophilic archaeon "Ferroplasma acidarmanus" Fer1

Analysis of the genome of "Ferroplasma acidarmanus" Fer1, an archaeon that is an extreme acidophile, identified an open reading frame encoding a putative ATP-dependent DNA ligase, which we termed FaLig. The deduced amino acid sequence of FaLig contains 595 amino acids, with a predicted molecular mass of 67.8 kDa. "F. acidarmanus" Fer1 is classified as a Euryarchaeote, but phylogenetic analysis using amino acid sequences showed that FaLig is more similar to DNA ligases from Crenarchaeota, suggesting that lateral transfer of these genes has occurred among archaea. The gene sequence encoding FaLig was cloned into a bacterial expression vector harbouring an upstream His-tag to aid purification. Conditions for expression and purification from Escherichia coli were identified and recombinant FaLig was confirmed to be an ATP-dependent DNA ligase. Optimal conditions for nick-joining by the protein were pH 6-7, 0.5 mM ATP, in the presence of either Mg(2+) or Mn(2+). Using a range of nicked, double-stranded nucleic acids, ligation was detected with the same substrates as previously determined for other DNA ligases. Although FaLig is the DNA ligase from one of the most extreme acidophilic organism yet studied, this characterization suggests that its biochemical mechanism is analogous to that of enzymes from other cellular systems.

LC-MS/MS based proteomic analysis and functional inference of hypothetical proteins in Desulfovibrio vulgaris

High efficiency capillary liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to examine the proteins extracted from Desulfovibrio vulgaris cells across six treatment conditions. While our previous study provided a proteomic overview of the cellular metabolism based on proteins with known functions [W. Zhang, M.A. Gritsenko, R.J. Moore, D.E. Culley, L. Nie, K. Petritis, E.F. Strittmatter, D.G. Camp II, R.D. Smith, F.J. Brockman, A proteomic view of the metabolism in Desulfovibrio vulgaris determined by liquid chromatography coupled with tandem mass spectrometry, Proteomics 6 (2006) 4286-4299], this study describes the global detection and functional inference for hypothetical D. vulgaris proteins. Using criteria that a given peptide of a protein is identified from at least two out of three independent LC-MS/MS measurements and that for any protein at least two different peptides are identified among the three measurements, 129 open reading frames (ORFs) originally annotated as hypothetical proteins were found to encode expressed proteins. Functional inference for the conserved hypothetical proteins was performed by a combination of several non-homology based methods: genomic context analysis, phylogenomic profiling, and analysis of a combination of experimental information, including peptide detection in cells grown under specific culture conditions and cellular location of the proteins. Using this approach we were able to assign possible functions to 20 conserved hypothetical proteins. This study demonstrated that a combination of proteomics and bioinformatics methodologies can provide verification of the expression of hypothetical proteins and improve genome annotation.

Distribution, structure and diversity of "bacterial" genes encoding two-component proteins in the Euryarchaeota

The publicly available annotated archaeal genome sequences (23 complete and three partial annotations, October 2005) were searched for the presence of potential two-component open reading frames (ORFs) using gene category lists and BLASTP. A total of 489 potential two-component genes were identified from the gene category lists and BLASTP. Two-component genes were found in 14 of the 21 Euryarchaeal sequences (October 2005) and in neither the Crenarchaeota nor the Nanoarchaeota. A total of 20 predicted protein domains were identified in the putative two-component ORFs that, in addition to the histidine kinase and receiver domains, also includes sensor and signalling domains. The detailed structure of these putative proteins is shown, as is the distribution of each class of two-component genes in each species. Potential members of orthologous groups have been identified, as have any potential operons containing two or more two-component genes. The number of two-component genes in those Euryarchaeal species which have them seems to be linked more to lifestyle and habitat than to genome complexity, with most examples being found in Methanospirillum hungatei, Haloarcula marismortui, Methanococcoides burtonii and the mesophilic Methanosarcinales group. The large numbers of two-component genes in these species may reflect a greater requirement for internal regulation. Phylogenetic analysis of orthologous groups of five different protein classes, three probably involved in regulating taxis, suggests that most of these ORFs have been inherited vertically from an ancestral Euryarchaeal species and point to a limited number of key horizontal gene transfer events.

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.

Control of ribosome turnover during growth of the haloalkaliphilic archaeon Natronococcus occultus

The metabolism of ribosomes during growth of the haloalkaliphilic archaeon Natronococcus occultus was examined. The ribosome content was higher during exponential growth and diminished to 35% of the maximum in the stationary stage. The incorporation of H3-orotic acid and C14-uracil into rRNA was higher during exponential growth. After that, it decreased to 39% of the maximum in the stationary stage. The labeling of non-ribosomal RNA took place almost exclusively in the exponential stage. From loss of radioactivity, the half-life of rRNA was 11.43, 14.85, 5.28 and 7.14 h during the initial, exponential, late exponential and stationary growth stages, respectively. These results suggested that increased synthesis combined with diminished degradation were responsible for the high ribosome content displayed by Ncc. occultus during exponential growth. In contrast, diminished synthesis together with increased degradation provoked its posterior loss.

A comparative categorization of protein function encoded in bacterial or archeal genomic islands

Genomes of prokaryotes harbor genomic islands (GIs), which are frequently acquired via horizontal gene transfer (HGT). Here I present an analysis of GIs with respect to gene-encoded functions. GIs were identified by statistical analysis of codon usage and clustering. Genes classified as putatively alien (pA) or putatively native (pN) were categorized according to the COG database. Among pA and pN genes, the distribution of COG functions and classes were studied for different groupings of prokaryotes. Groups were formed according to taxonomical relation or habitats. In all groups, genes related to class L (replication, recombination, and repair) were statistically significantly overrepresented in GIs. GIs of bacteria and archaea showed a distinct pattern of preferences. In archeal GIs, genes belonging to COG class M (cell wall/membrane/envelope biogenesis) or Q (secondary metabolites biosynthesis, transport, and catabolism) were more frequent. In bacterial GIs, genes of classes U (intracellular trafficking, secretion, and vesicular transport), N (cell motility), and V (defense mechanisms) were predominant. Underrepresentation was strongest for genes belonging to class J (translation, ribosomal structure, and biogenesis). Among single COG functions overrepresented in GIs were transferases and transporters. In both superkingdoms, HGT enhances genomic content by meeting demands that are independent of the studied habitats. These findings are in agreement with the complexity theory, which predicts the preferential import of operational genes. However, only specific subsets of operational genes were enriched in GIs. Modification of the cell envelope, cell motility, secretion, and protection of cellular DNA are major issues in HGT.

Congruence of evidence for a Methanopyrus-proximal root of life based on transfer RNA and aminoacyl-tRNA synthetase genes

Among 60 organisms, the intraspecies genetic distances between tRNAs cognate for different amino acids, between the initiator and elongator tRNAs for Met, and between potentially paralogous pairs of aminoacyl-tRNA synthetases are found to be at a minimum within the Methanopyrus kandleri genome. These results indicate an exact congruence between the evidence from tRNA and aminoacyl-tRNA synthetase genes locating the root of life closest to this organism.

Comparative genomics of Thermus thermophilus and Deinococcus radiodurans: divergent routes of adaptation to thermophily and radiation resistance

Background

Thermus thermophilus and Deinococcus radiodurans belong to a distinct bacterial clade but have remarkably different phenotypes. T. thermophilus is a thermophile, which is relatively sensitive to ionizing radiation and desiccation, whereas D. radiodurans is a mesophile, which is highly radiation- and desiccation-resistant. Here we present an in-depth comparison of the genomes of these two related but differently adapted bacteria.

Results

By reconstructing the evolution of Thermus and Deinococcus after the divergence from their common ancestor, we demonstrate a high level of post-divergence gene flux in both lineages. Various aspects of the adaptation to high temperature in Thermus can be attributed to horizontal gene transfer from archaea and thermophilic bacteria; many of the horizontally transferred genes are located on the single megaplasmid of Thermus. In addition, the Thermus lineage has lost a set of genes that are still present in Deinococcus and many other mesophilic bacteria but are not common among thermophiles. By contrast, Deinococcus seems to have acquired numerous genes related to stress response systems from various bacteria. A comparison of the distribution of orthologous genes among the four partitions of the Deinococcus genome and the two partitions of the Thermus genome reveals homology between the Thermus megaplasmid (pTT27) and Deinococcus megaplasmid (DR177).

Conclusion

After the radiation from their common ancestor, the Thermus and Deinococcus lineages have taken divergent paths toward their distinct lifestyles. In addition to extensive gene loss, Thermus seems to have acquired numerous genes from thermophiles, which likely was the decisive contribution to its thermophilic adaptation. By contrast, Deinococcus lost few genes but seems to have acquired many bacterial genes that apparently enhanced its ability to survive different kinds of environmental stresses. Notwithstanding the accumulation of horizontally transferred genes, we also show that the single megaplasmid of Thermus and the DR177 megaplasmid of Deinococcus are homologous and probably were inherited from the common ancestor of these bacteria.

Yeast (ADPribosyl)ation: revisiting a controversial question

The controversy about the occurrence of an (ADPribosyl)ating activity in yeast is still standing up. Here we discuss this topic on the basis of results obtained with classic experiments proposed over years as basis to characterize an (ADPribosyl)ation system in any organism. Independent results obtained in two different laboratories were in line with each other and went towards the occurrence of an active (ADPribosyl)ating system in Saccharomyces cerevisiae. In fact data collected from nuclear preparations of cultured cells matched those from baker's yeast and lyophilized yeast cells. Yeast (ADPribosyl)ating enzyme is a protein of 80-90 kDa, as determined by electrophoresis on polyacrylamide gel in sodium dodecyl sulphate, followed by immunoblotting with antibodies against anti-poly(ADPribose) polymerase catalytic site. It synthesizes products, that, after digestion with phosphodiesterase, co-migrates mainly with phosphoribosyl adenosine monophosphate after thin layer chromatography on silica gel plate.

Evolutionary and functional genomics of the Archaea

In the past two years, archaeal genomics has achieved several breakthroughs. On the evolutionary front the most exciting development was the sequencing and analysis of the genome of Nanoarchaeum equitans, a tiny parasitic organism that has only approximately 540 genes. The genome of Nanoarchaeum shows signs of extreme rearrangement including the virtual absence of conserved operons and the presence of several split genes. Nanoarchaeum is distantly related to other archaea, and it has been proposed to represent a deep archaeal branch that is distinct from Euryarchaeota and Crenarchaeota. This would imply that many features of its gene repertoire and genome organization might be ancestral. However, additional genome analysis has provided a more conservative suggestion - that Nanoarchaeum is a highly derived euryarchaeon. Also there have been substantial developments in functional genomics, including the discovery of the elusive aminoacyl-tRNA synthetase that is involved in both the biosynthesis of cysteine and its incorporation into proteins in methanogens, and the first experimental validation of the predicted archaeal exosome.

Characterization of a thermophilic ATP-dependent DNA ligase from the euryarchaeon Pyrococcus horikoshii

Archaea encode a DNA ligase composed of a C-terminal catalytic domain typical of ATP-dependent ligases plus an N-terminal domain similar to that found in eukaryotic cellular and poxvirus DNA ligases. All archaeal DNA ligases characterized to date have ATP-dependent adenylyltransferase and nick-joining activities. However, recent reports of dual-specificity ATP/NAD+ ligases in two Thermococcus species and Pyrococcus abyssi and an ATP/ADP ligase in Aeropyrum pernix raise the prospect that certain archaeal enzymes might exemplify an undifferentiated ancestral stage in the evolution of ligase substrate specificity. Here we analyze the biochemical properties of Pyrococcus horikoshii DNA ligase. P. horikoshii ligase catalyzes auto-adenylylation and nick sealing in the presence of a divalent cation and ATP; it is unable to utilize NAD+ or ADP to promote ligation in lieu of ATP. P. horikoshii ligase is thermophilic in vitro, with optimal adenylyltransferase activity at 90 degrees C and nick-joining activity at 70 to 90 degrees C. P. horikoshii ligase resembles the ligases of Methanobacterium thermautotrophicum and Sulfolobus shibatae in its strict specificity for ATP.

Discovering novel biology by in silico archaeology

Archaea are prokaryotes that evolved in parallel with bacteria. Since the discovery of the distinct status of the Archaea, extensive physiological and biochemical research has been conducted to elucidate the molecular basis of their remarkable lifestyle and their unique biology. Here, we discuss how in-depth comparative genomics has been used to improve the annotation of archaeal genomes. Combined with experimental verification, bioinformatic analysis contributes to the ongoing discovery of novel metabolic conversions and control mechanisms, and as such to a better understanding of the intriguing biology of the Archaea.

Ancestral paralogs and pseudoparalogs and their role in the emergence of the eukaryotic cell

Gene duplication is a crucial mechanism of evolutionary innovation. A substantial fraction of eukaryotic genomes consists of paralogous gene families. We assess the extent of ancestral paralogy, which dates back to the last common ancestor of all eukaryotes, and examine the origins of the ancestral paralogs and their potential roles in the emergence of the eukaryotic cell complexity. A parsimonious reconstruction of ancestral gene repertoires shows that 4137 orthologous gene sets in the last eukaryotic common ancestor (LECA) map back to 2150 orthologous sets in the hypothetical first eukaryotic common ancestor (FECA) [paralogy quotient (PQ) of 1.92]. Analogous reconstructions show significantly lower levels of paralogy in prokaryotes, 1.19 for archaea and 1.25 for bacteria. The only functional class of eukaryotic proteins with a significant excess of paralogous clusters over the mean includes molecular chaperones and proteins with related functions. Almost all genes in this category underwent multiple duplications during early eukaryotic evolution. In structural terms, the most prominent sets of paralogs are superstructure-forming proteins with repetitive domains, such as WD-40 and TPR. In addition to the true ancestral paralogs which evolved via duplication at the onset of eukaryotic evolution, numerous pseudoparalogs were detected, i.e. homologous genes that apparently were acquired by early eukaryotes via different routes, including horizontal gene transfer (HGT) from diverse bacteria. The results of this study demonstrate a major increase in the level of gene paralogy as a hallmark of the early evolution of eukaryotes.

Evolution of a microbial nitrilase gene family: a comparative and environmental genomics study

BACKGROUND

Completed genomes and environmental genomic sequences are bringing a significant contribution to understanding the evolution of gene families, microbial metabolism and community eco-physiology. Here, we used comparative genomics and phylogenetic analyses in conjunction with enzymatic data to probe the evolution and functions of a microbial nitrilase gene family. Nitrilases are relatively rare in bacterial genomes, their biological function being unclear.

RESULTS

We examined the genetic neighborhood of the different subfamily genes and discovered conserved gene clusters or operons associated with specific nitrilase clades. The inferred evolutionary transitions that separate nitrilases which belong to different gene clusters correlated with changes in their enzymatic properties. We present evidence that Darwinian adaptation acted during one of those transitions and identified sites in the enzyme that may have been under positive selection.

CONCLUSION

Changes in the observed biochemical properties of the nitrilases associated with the different gene clusters are consistent with a hypothesis that those enzymes have been recruited to a novel metabolic pathway following gene duplication and neofunctionalization. These results demonstrate the benefits of combining environmental genomic sampling and completed genomes data with evolutionary and biochemical analyses in the study of gene families. They also open new directions for studying the functions of nitrilases and the genes they are associated with.